找回密码
 注册
关于网站域名变更的通知
查看: 569|回复: 1
打印 上一主题 下一主题

求二阶矩阵的均衡点,跑不出来

[复制链接]

该用户从未签到

跳转到指定楼层
1#
发表于 2022-11-11 11:00 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式

EDA365欢迎您登录!

您需要 登录 才可以下载或查看,没有帐号?注册

x
%1
8 w. V* y* s3 U1 R) S! d1 rsyms x y b b1 n m b2 p c g r;
  |1 B6 W! }) _2 C4 X: s. ]6 m& qv11=y*(-b1-m+r)+(1-y)*(-b1-b2+n);
" G/ n' k) e0 D9 Y0 Yv12=y*r-(1-y)*b2;
% {3 r4 P: T. `# V) }v1=simplify(x*v11-(1-x)*v12)
4 S/ X: t1 z: E* v& g* nfx=simplify(x*(v11-v1))
, A' Z+ T# y) _" s0 q%2; F# I" ]5 O) U8 i
v21=x*(-c+p+m+g)+(1-x)*(-c+p+g);. i, ]* w) R, G* L9 N4 E
v22=x*(p-n)+(1-x)*p;
& I' ^) X( k4 x7 v4 kv2=simplify(y*v21-(1-y)*v22)
  O2 g6 m1 {' W+ A) {8 e. Lfy=simplify(y*(v21-v2))- ?( m5 H2 i& \* ?/ J9 _
%构建雅可比矩阵& c7 G# c$ A9 g3 t8 u! F
disp=(['雅可比矩阵']);9 t4 I3 ?' z# Q9 X
A = [diff(fx,x) diff(fx,y);
/ p2 y9 d: U4 I' I1 ^' U+ X- |2 O    diff(fy,x) diff(fy,y)], k" M3 X1 b' C' B" @3 r$ I
%2求均衡点                                                                                7 P; q& K( X3 e) {: D2 h* r7 W) Z
equ=[fx==0,fy==0];                                                               & o+ a' T- M# r4 S
answ=solve(equ,[x,y]);                                                               
* g, Q0 F2 _2 g) O$ ndisp('均衡点:');A1=[answ.x,answ.y]7 [( B+ ^4 ^5 y: v0 r$ @
%均衡点个数                                                                                 
# R: B# k  }; {. W* \disp(['均衡点个数:']);length(answ.x)

- k$ T- T* x0 w( {
  • TA的每日心情
    开心
    2022-11-22 15:53
  • 签到天数: 2 天

    [LV.1]初来乍到

    2#
    发表于 2022-11-11 13:09 | 只看该作者
    这里 disp=(['雅可比矩阵']);,把等号删了
    您需要登录后才可以回帖 登录 | 注册

    本版积分规则

    关闭

    推荐内容上一条 /1 下一条

    EDA365公众号

    关于我们|手机版|EDA365电子论坛网 ( 粤ICP备18020198号-1 )

    GMT+8, 2025-10-28 00:10 , Processed in 0.125000 second(s), 23 queries , Gzip On.

    深圳市墨知创新科技有限公司

    地址:深圳市南山区科技生态园2栋A座805 电话:19926409050

    快速回复 返回顶部 返回列表