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

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

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x
%1
: t8 c. a8 k$ t' @3 A" bsyms x y b b1 n m b2 p c g r;
! i: T4 ?( k* D3 Iv11=y*(-b1-m+r)+(1-y)*(-b1-b2+n);
9 A2 n0 X% V" U6 S0 Tv12=y*r-(1-y)*b2;
! e! m6 ^" p; a  fv1=simplify(x*v11-(1-x)*v12)8 ^8 x6 T. w0 w' F6 G+ Y
fx=simplify(x*(v11-v1))% o- P8 e2 `( O- g, X" Y
%2% P* X0 |8 S. k0 q
v21=x*(-c+p+m+g)+(1-x)*(-c+p+g);
( I; k7 m2 ]* A9 o+ i; F$ _v22=x*(p-n)+(1-x)*p;
. c) f) J1 ?( cv2=simplify(y*v21-(1-y)*v22)3 R$ g: e: d; f/ b8 C4 W
fy=simplify(y*(v21-v2))2 \0 C# `& F/ n" x2 i
%构建雅可比矩阵
6 \! F# e4 g$ k; ~& Udisp=(['雅可比矩阵']);
) w% N$ f. F9 U8 H  B5 Q; {. qA = [diff(fx,x) diff(fx,y);
9 |2 r& ^5 ?" a' }: L" a( H    diff(fy,x) diff(fy,y)]
" z9 i- G6 ?. @, S. c: L3 g%2求均衡点                                                                                - n6 V. H2 d0 g0 P
equ=[fx==0,fy==0];                                                               
* \, A$ ~: x/ {9 }8 u$ Ransw=solve(equ,[x,y]);                                                               
6 X( w0 H7 W: M! A; jdisp('均衡点:');A1=[answ.x,answ.y]
8 t9 H$ E' S! [1 O%均衡点个数                                                                                 
6 K+ _& C8 R! i/ M, V+ t" ~2 _5 gdisp(['均衡点个数:']);length(answ.x)

/ w: S2 A0 E6 ^" j
  • 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 13:25 , Processed in 0.140625 second(s), 23 queries , Gzip On.

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

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

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