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

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

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x
%1
/ t9 X' M/ f3 M, ?  O* gsyms x y b b1 n m b2 p c g r;8 k- \, c) P& r
v11=y*(-b1-m+r)+(1-y)*(-b1-b2+n);
$ \& n) m# G  p, [0 l. Rv12=y*r-(1-y)*b2;
5 \  R( d- L  _1 V* Z) [; J8 `9 fv1=simplify(x*v11-(1-x)*v12)+ b; Y, L; i$ `/ W
fx=simplify(x*(v11-v1))
) I- b* d6 M3 x  q2 ~%2
+ }/ C* R# C- Q" bv21=x*(-c+p+m+g)+(1-x)*(-c+p+g);$ W3 D9 n& V: y7 F) ~& T  f' b
v22=x*(p-n)+(1-x)*p;
7 E% ]9 M/ U. S5 O. O7 |v2=simplify(y*v21-(1-y)*v22)
  T5 L' L, u8 A% b8 ?. _, tfy=simplify(y*(v21-v2)), Q% s1 ]$ M1 H) o2 z. B# k
%构建雅可比矩阵8 S4 }& X" d" C  `  G# o; q
disp=(['雅可比矩阵']);
$ f% o6 {6 J) L6 I" y' g2 ]A = [diff(fx,x) diff(fx,y);
" [6 v; e1 G* ^+ r* f    diff(fy,x) diff(fy,y)]9 F0 q  t# c- l. A' Y: [
%2求均衡点                                                                                + O0 q8 o( \5 S3 w3 R) l7 g
equ=[fx==0,fy==0];                                                               " T3 k3 E: s/ {. e, h  R; r0 i* t
answ=solve(equ,[x,y]);                                                               
1 H0 F$ C! V8 w9 l- E: |! Ddisp('均衡点:');A1=[answ.x,answ.y]) y' v* {  A7 F9 s7 y9 \
%均衡点个数                                                                                 
+ Y) Z" j- `3 v4 ydisp(['均衡点个数:']);length(answ.x)
, _: @% W9 C1 S# S0 t( F, P! |5 ?
  • 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 16:02 , Processed in 0.156250 second(s), 23 queries , Gzip On.

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

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

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