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

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

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x
%1
( Q: q5 E% F5 p8 C: A, Asyms x y b b1 n m b2 p c g r;. }% C2 d8 k# c, d& b
v11=y*(-b1-m+r)+(1-y)*(-b1-b2+n);/ [. ^* p3 U' I: m3 G, o* H
v12=y*r-(1-y)*b2;, y( C: a8 Z0 O
v1=simplify(x*v11-(1-x)*v12)$ K; V% ]* M2 n2 Z# U' C9 D
fx=simplify(x*(v11-v1))
' f. J9 C4 |" r, @1 [) ]- w- {%2
$ a' R( l5 v6 sv21=x*(-c+p+m+g)+(1-x)*(-c+p+g);. C( n- b1 z& i0 P8 H
v22=x*(p-n)+(1-x)*p;
1 }+ {, y4 S! ], y0 j( M/ iv2=simplify(y*v21-(1-y)*v22)! w, y& H. ]  [& r
fy=simplify(y*(v21-v2))
: F+ F6 s% O! \( n. l  D# Y5 G%构建雅可比矩阵! ~3 B% B& `) O- U! ]6 ~
disp=(['雅可比矩阵']);' Z9 R+ C% [9 S6 h) X4 U/ h
A = [diff(fx,x) diff(fx,y);. ~. ~/ t% j6 a5 F7 ~3 {
    diff(fy,x) diff(fy,y)]6 r1 v% T$ S! d9 f$ j/ ?
%2求均衡点                                                                                ; ^7 {  o* i$ L% D9 I
equ=[fx==0,fy==0];                                                               
0 \& }* y4 n/ ^' [# W* H% k, ]answ=solve(equ,[x,y]);                                                               2 Z- ~! O, ~( F3 W  C
disp('均衡点:');A1=[answ.x,answ.y]
$ y9 C7 _% t, ?6 x, R+ P1 }%均衡点个数                                                                                 - F3 x& e7 d% j& i/ b
disp(['均衡点个数:']);length(answ.x)

9 c! }/ i1 l/ k- e# W/ ?' Y! L5 n
  • 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-8-24 06:22 , Processed in 0.109375 second(s), 23 queries , Gzip On.

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

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

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