|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
%1! q5 B/ u9 r2 d/ }8 V# X4 c$ \
syms x y b b1 n m b2 p c g r;. x1 W2 F( T* \3 i3 C9 _
v11=y*(-b1-m+r)+(1-y)*(-b1-b2+n);
6 y, V/ R Y) V3 Dv12=y*r-(1-y)*b2; x, P7 D+ A. v6 f& q" [+ h+ T$ b) ~
v1=simplify(x*v11-(1-x)*v12)
Y( O9 Q$ S6 g( r( X5 G: y& xfx=simplify(x*(v11-v1))
" f3 z2 ~0 Y ^2 ?* P: z%2
* O" ~" n1 h+ o6 F I. i/ Nv21=x*(-c+p+m+g)+(1-x)*(-c+p+g);
7 G* A% k& n) O3 k/ q( n! h, Av22=x*(p-n)+(1-x)*p;
+ k8 `: x0 ^6 {0 a8 ov2=simplify(y*v21-(1-y)*v22)
) t5 B( U+ I( z8 h9 B; g$ ~4 C% I0 qfy=simplify(y*(v21-v2))8 C4 S! m& S; ]- q1 E
%构建雅可比矩阵
0 N4 d C/ s8 _ I" Udisp=(['雅可比矩阵']);
( d) W1 G9 k) y; Q2 i8 EA = [diff(fx,x) diff(fx,y);% Y5 q. X. B$ e& Q7 p# H2 Y
diff(fy,x) diff(fy,y)]
9 l* f/ C3 x) @: ?3 k%2求均衡点 6 }$ m) I) k8 B9 N; P
equ=[fx==0,fy==0]; $ b9 k# B2 Z# Q
answ=solve(equ,[x,y]);
) V. _/ `/ A2 R5 P" t2 fdisp('均衡点:');A1=[answ.x,answ.y]: o1 _4 B$ o* s; h
%均衡点个数
" ~$ V' i1 ^. M9 c5 L* A! o6 o& ^disp(['均衡点个数:']);length(answ.x) |
$ d/ U3 w2 Y8 T/ Q' x+ S |
|
/1 
发表于 2022-11-11 11:00
收藏
淘帖
支持！
反对！