找回密码
 注册
查看: 519|回复: 1
打印 上一主题 下一主题

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

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x
%1
8 o9 q5 \5 L8 tsyms x y b b1 n m b2 p c g r;
  A  Z4 d& `9 D/ pv11=y*(-b1-m+r)+(1-y)*(-b1-b2+n);) `2 q1 e' D5 }; |
v12=y*r-(1-y)*b2;! ^- ?8 d$ w0 f4 N
v1=simplify(x*v11-(1-x)*v12)  M  D1 [4 J$ N, ?
fx=simplify(x*(v11-v1))
. I7 a+ w& f6 y( s2 ?%2
; [. F, D2 K! p- sv21=x*(-c+p+m+g)+(1-x)*(-c+p+g);$ r7 x' n4 V! t! N; R
v22=x*(p-n)+(1-x)*p;
# m4 w( F! q4 D1 A7 P9 n4 r4 yv2=simplify(y*v21-(1-y)*v22)
8 d! ?5 w& H$ z& R( b3 Lfy=simplify(y*(v21-v2))0 O0 d9 d! l8 f/ u3 s7 R
%构建雅可比矩阵
* X2 i$ T; p2 O3 o& Q# {6 ^9 D& |  }disp=(['雅可比矩阵']);
3 n5 S  {- i$ m; e8 F& |A = [diff(fx,x) diff(fx,y);
& w0 w! S" r, F$ l; ^8 }    diff(fy,x) diff(fy,y)]
" r+ l, b! P) Q2 ~1 W7 i4 j% Q# Y  F! X%2求均衡点                                                                                ! t, ^! f6 x2 D/ t9 t
equ=[fx==0,fy==0];                                                               
/ K% e- b- w* }- bansw=solve(equ,[x,y]);                                                               3 d% H, q- Z8 k- \
disp('均衡点:');A1=[answ.x,answ.y]( I7 ^! R8 k3 e# S
%均衡点个数                                                                                 2 E) P* k( d% @5 F* S
disp(['均衡点个数:']);length(answ.x)
0 N7 |1 q% k6 t
  • 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-5-23 22:44 , Processed in 0.093750 second(s), 24 queries , Gzip On.

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

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

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