EDA365电子论坛网

标题: 求二阶矩阵的均衡点，跑不出来 [打印本页]
作者: justlikethisis    时间: 2022-11-11 11:00
标题: 求二阶矩阵的均衡点，跑不出来
%1& U2 ^& W+ t  p) K1 @; L1 Z
syms x y b b1 n m b2 p c g r;" \# x  ?: b# ?- K% d, S" t
v11=y*(-b1-m+r)+(1-y)*(-b1-b2+n);
% c2 b) E2 M* s( S( O7 zv12=y*r-(1-y)*b2;
, n- x  b1 Q' Y- l! ov1=simplify(x*v11-(1-x)*v12)7 f; A  }, h% \/ b6 z( B4 u
fx=simplify(x*(v11-v1))+ b  L) d; R2 y7 ^
%2
' _) m0 Y: O. B% h, K2 W1 f9 I6 Nv21=x*(-c+p+m+g)+(1-x)*(-c+p+g);
' e7 l$ W1 h, g8 K* ^$ e* jv22=x*(p-n)+(1-x)*p;+ ?: g( v3 O6 p( B, R4 [, _
v2=simplify(y*v21-(1-y)*v22)% w! \: K0 o7 k  A! m
fy=simplify(y*(v21-v2))+ Y  f6 A& C6 s* n8 Y
%构建雅可比矩阵5 ?+ Z% k* J7 H
disp=(['雅可比矩阵']);
3 h) V+ _4 Q5 ^6 m8 I& JA = [diff(fx,x) diff(fx,y);
. r( }2 l# r# R0 ?    diff(fy,x) diff(fy,y)]2 e; B9 q2 b& ~9 Z+ `2 S
%2求均衡点                                                                                % P: I2 Z# T6 n3 v0 x) ]0 S+ k
equ=[fx==0,fy==0];                                                               / N. ]" m  U! y) c% Y2 M% q
answ=solve(equ,[x,y]);                                                               : W& a# ]+ s6 J9 C3 `9 `5 ]
disp('均衡点:');A1=[answ.x,answ.y]
( D3 ]8 E( i7 x%均衡点个数                                                                                 
1 F# }' w/ v( f0 g5 J* e! T8 Tdisp(['均衡点个数:']);length(answ.x)
  t; ], n- C2 ~9 J4 v
作者: fantasyqqq    时间: 2022-11-11 13:09
这里 disp=(['雅可比矩阵']);，把等号删了



欢迎光临 EDA365电子论坛网 (https://bbs.eda365.com/) Powered by Discuz! X3.2