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

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

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x
%1. S4 H) q4 s) C/ @! A9 K) Z
syms x y b b1 n m b2 p c g r;- x; O6 w" [0 K: r) Y7 E
v11=y*(-b1-m+r)+(1-y)*(-b1-b2+n);9 f/ Q! z) T) Q1 n
v12=y*r-(1-y)*b2;
1 J: O- Z! S8 d7 a! Mv1=simplify(x*v11-(1-x)*v12)
% ~1 {, j) h2 R* ^( qfx=simplify(x*(v11-v1))4 J$ g5 h. x+ C* ^) j2 n' f1 m) \
%2, @3 t) n. g( f
v21=x*(-c+p+m+g)+(1-x)*(-c+p+g);0 n$ b8 ^2 E5 |/ e& h- O
v22=x*(p-n)+(1-x)*p;
; D4 d! N' N3 Y/ X/ e2 v& L1 hv2=simplify(y*v21-(1-y)*v22)
% z% z) [+ |, i$ `fy=simplify(y*(v21-v2))
  H6 U5 R2 I0 S& X" U$ ~%构建雅可比矩阵
3 a6 p7 s) _9 w1 |, k5 Mdisp=(['雅可比矩阵']);/ S- d1 {, D+ S4 i" p  ]1 s5 K
A = [diff(fx,x) diff(fx,y);, W8 y0 L2 [4 u
    diff(fy,x) diff(fy,y)]
2 l$ K* j& J1 q%2求均衡点                                                                                ; y* }( r0 {2 j, w
equ=[fx==0,fy==0];                                                               4 T& F* ^2 U, U% L) W% h' E# |
answ=solve(equ,[x,y]);                                                               8 R" `/ s* i" k8 E8 F
disp('均衡点:');A1=[answ.x,answ.y]* {/ B  a' A- V& [+ U4 w
%均衡点个数                                                                                 
; S9 a$ ~2 T# n; y& U/ W+ ~disp(['均衡点个数:']);length(answ.x)
% C9 C. ?, U& e' B
  • 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-10-4 13:16 , Processed in 0.125000 second(s), 24 queries , Gzip On.

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

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

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