|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
%1. S4 H) q4 s) C/ @! A9 K) Z
syms x y b b1 n m b2 p c g r;- x; O6 w" [0 K: r) Y7 E
v11=y*(-b1-m+r)+(1-y)*(-b1-b2+n);9 f/ Q! z) T) Q1 n
v12=y*r-(1-y)*b2;
1 J: O- Z! S8 d7 a! Mv1=simplify(x*v11-(1-x)*v12)
% ~1 {, j) h2 R* ^( qfx=simplify(x*(v11-v1))4 J$ g5 h. x+ C* ^) j2 n' f1 m) \
%2, @3 t) n. g( f
v21=x*(-c+p+m+g)+(1-x)*(-c+p+g);0 n$ b8 ^2 E5 |/ e& h- O
v22=x*(p-n)+(1-x)*p;
; D4 d! N' N3 Y/ X/ e2 v& L1 hv2=simplify(y*v21-(1-y)*v22)
% z% z) [+ |, i$ `fy=simplify(y*(v21-v2))
H6 U5 R2 I0 S& X" U$ ~%构建雅可比矩阵
3 a6 p7 s) _9 w1 |, k5 Mdisp=(['雅可比矩阵']);/ S- d1 {, D+ S4 i" p ]1 s5 K
A = [diff(fx,x) diff(fx,y);, W8 y0 L2 [4 u
diff(fy,x) diff(fy,y)]
2 l$ K* j& J1 q%2求均衡点 ; y* }( r0 {2 j, w
equ=[fx==0,fy==0]; 4 T& F* ^2 U, U% L) W% h' E# |
answ=solve(equ,[x,y]); 8 R" `/ s* i" k8 E8 F
disp('均衡点:');A1=[answ.x,answ.y]* {/ B a' A- V& [+ U4 w
%均衡点个数
; S9 a$ ~2 T# n; y& U/ W+ ~disp(['均衡点个数:']);length(answ.x) | % C9 C. ?, U& e' B
|
|