|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
%1
, e! B. L" {! V# V1 Q asyms x y b b1 n m b2 p c g r;5 F% b* @! X# Z# l \' ~
v11=y*(-b1-m+r)+(1-y)*(-b1-b2+n);
- v% j+ @+ Y' [' r) T* i2 `v12=y*r-(1-y)*b2;! L, R0 M; u: h+ E; R( t
v1=simplify(x*v11-(1-x)*v12)1 h, n5 J3 ?5 O1 s- b
fx=simplify(x*(v11-v1))' Y+ l' ^" z; x# g/ ^' F% a. [1 \5 v
%2
$ v8 R+ s, F! a0 L$ @- Lv21=x*(-c+p+m+g)+(1-x)*(-c+p+g);
, r% l6 ~% }5 n9 A" j1 iv22=x*(p-n)+(1-x)*p;3 F+ q8 ~9 ^6 t' D" Z Y: ^( n
v2=simplify(y*v21-(1-y)*v22)8 y6 h& C7 ^0 `: G9 ~2 c( w; F' k
fy=simplify(y*(v21-v2))8 P& f. Y. c( |% T! F) G1 ~( n
%构建雅可比矩阵$ o2 R7 R0 R4 O
disp=(['雅可比矩阵']);1 z7 n4 w6 ?9 R( j8 g3 _1 A
A = [diff(fx,x) diff(fx,y);
V% I" X9 K+ T' X diff(fy,x) diff(fy,y)]
8 s9 y$ m2 q# Q%2求均衡点
2 v) u' F% O. N. |" mequ=[fx==0,fy==0]; ! D8 @, g% e2 K& R) ~* B/ q
answ=solve(equ,[x,y]); 7 w p% e& d7 z1 H
disp('均衡点:');A1=[answ.x,answ.y]
' A } T+ _. r%均衡点个数 $ j/ M; K7 _" I& k
disp(['均衡点个数:']);length(answ.x) |
: g1 O7 u# L, z3 `- | |
|
/1 
发表于 2022-11-11 11:00
收藏
淘帖
支持！
反对！