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

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

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x
%1
2 `  D3 |7 n! k* C: M" {syms x y b b1 n m b2 p c g r;
( J3 |# W5 ?" R6 pv11=y*(-b1-m+r)+(1-y)*(-b1-b2+n);  p8 K( |: @+ D0 z- R
v12=y*r-(1-y)*b2;
# G9 O/ i( w: R0 @' Nv1=simplify(x*v11-(1-x)*v12)6 V5 |7 V  h7 R2 S- _: R. K8 a8 l
fx=simplify(x*(v11-v1))
' K0 f+ Y) X, y' E# t6 p%20 f8 v) r1 ^) D
v21=x*(-c+p+m+g)+(1-x)*(-c+p+g);
4 C+ L' l+ @3 L& u1 z: vv22=x*(p-n)+(1-x)*p;7 |$ ]' Y, O! F: i) s
v2=simplify(y*v21-(1-y)*v22)
# w- M2 u8 q0 T' n. lfy=simplify(y*(v21-v2))) a# D6 S  q) F' q) g9 ?1 P, B
%构建雅可比矩阵; M" O2 d7 P) m; @# i0 I- X: o
disp=(['雅可比矩阵']);( M: `' ?, n* H- t
A = [diff(fx,x) diff(fx,y);4 ~2 L. A: o# w* c  u4 R
    diff(fy,x) diff(fy,y)]
1 U' F* x7 M( C% A%2求均衡点                                                                                ! T( @; v9 A( z, S# q
equ=[fx==0,fy==0];                                                               
7 ?$ j% G7 ^& t" C& Ransw=solve(equ,[x,y]);                                                               : k  c7 w# I3 U7 s& k
disp('均衡点:');A1=[answ.x,answ.y]- @* k1 W0 ~) V/ R3 F+ \. [& [
%均衡点个数                                                                                 
& s# G8 d; m& L) O: Y; z. Zdisp(['均衡点个数:']);length(answ.x)
+ |& T) l7 f0 ]) l
  • 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-8-1 12:31 , Processed in 0.125000 second(s), 23 queries , Gzip On.

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

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

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