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

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

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x
%19 O6 I3 h. t" d% `' e$ w3 \+ |
syms x y b b1 n m b2 p c g r;" n( m. {: Y# T0 A/ K9 f9 M/ P9 @/ t
v11=y*(-b1-m+r)+(1-y)*(-b1-b2+n);
( p0 B4 w$ _# J7 Ev12=y*r-(1-y)*b2;+ f) P$ j( U0 y, k! h9 i* \' F+ k0 P
v1=simplify(x*v11-(1-x)*v12)' S  _, Y5 ^# b- x
fx=simplify(x*(v11-v1)). t6 }/ b0 P6 K* w
%2/ x4 a& J# |( F& E* W% l
v21=x*(-c+p+m+g)+(1-x)*(-c+p+g);
4 ~" [" j% }) H; X9 N; T/ j1 hv22=x*(p-n)+(1-x)*p;- Z0 M7 S) c+ v
v2=simplify(y*v21-(1-y)*v22)
+ [- e# Y. W2 Tfy=simplify(y*(v21-v2)): }: S6 C; B! r
%构建雅可比矩阵
/ x* c. c8 h, X: @0 zdisp=(['雅可比矩阵']);
* m8 [; _% h0 X( w- U2 `8 SA = [diff(fx,x) diff(fx,y);
( g% T! ~' y, J0 _, e    diff(fy,x) diff(fy,y)]3 q1 |( Z! h) s- I" d* d9 k, H
%2求均衡点                                                                                
7 z* x/ N/ I  ^5 z* Y& ?9 {equ=[fx==0,fy==0];                                                               , ~+ y, j2 \4 p% D. r5 e: p1 |4 }
answ=solve(equ,[x,y]);                                                               / z- ^3 Y, T# _$ _6 R5 a$ P$ Y
disp('均衡点:');A1=[answ.x,answ.y]: r" o. [9 p' A' N, k
%均衡点个数                                                                                 
  t+ b6 n' l  G2 f7 L8 ]disp(['均衡点个数:']);length(answ.x)
+ n0 i8 G' u* ]! O
  • 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-11-23 15:16 , Processed in 0.156250 second(s), 24 queries , Gzip On.

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

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

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