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

请教大神,这个是为什么呢?

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x

. O- ~2 C. _5 S% ^; z( }
; }) r7 N5 f5 }' e  Z9 K! x请教大神,这个是为什么呢??ode45 无法执行赋值,因为左侧的大小为 2×4,右侧的大小为 2×5?非常感谢!
$ _8 p! @: f( {. u6 g! @1 F1 H  G. d
% u3 }+ B1 i. z. \; I4 Lfunction  f=boat1(t,x)% J2 a& ]  K$ f9 O
if(norm(x)>1e-5)   %限制条件,保证在x足够小的情况下可以结束计算6 n* S5 D  w5 D. H: k5 S
    f=[-2*x(1)/sqrt(x(1).^2+x(2).^2);1-2*x(2)/sqrt(x(1).^2+x(2).^2)]! C0 o" ^% h9 |. }% Q
else
* s7 h" a! L/ U( \+ B: \    f=[0,0]
0 r" r5 n) |9 V2 H1 I$ d5 p0 @0 M. n, @* m8 y
end! [" e: o$ y6 Q2 v) u0 z
% }' {, y; b" U

& l3 L. B9 }" {% ]- r2 ?, }y0=[1000 0];      %微分方程的初值7 R3 w0 ]9 ~: P9 D( ~7 |9 @
[t,x]=ode45(@boat1,[0,1000],y0)      %求解微分方程的数值解( t( S( J! i* F5 ]  y
subplot(1,2,1),plot(t,x);         %画出关于x(1)和x(2)的系统时间响应曲线
0 P9 k0 S  u$ t$ X; Iplot(x(:,1),x(:,2))          %画出关于x(1)和x(2)的系统相平面曲线5 Y$ |/ T2 E7 S& [
axis([-1 1000,0 1000])          %限定区间  d) c: [. ^1 G- O. d5 ^1 {" y7 U
text(10,15,'B')         %标注B点1 |/ F# A, U7 z' ^2 Y! Z- w
text(1000,0,'A')         %标注A点
! B0 Z& E* u, G" \3 o5 gxlabel('X')
+ r' C" ~5 S5 i( F; o0 Tylable(‘Y’); r; [7 @& _" q/ P; U

4 {" F3 D" e6 ~/ m! h, P) l- n3 I2 }; p( l, O5 x7 e; M: g$ j# f
无法执行赋值,因为左侧的大小为 2×4,右侧的大小为 2×5。
2 G, F* t$ i! H0 Z0 i* D2 g, N/ T/ c' k& ]7 e9 u8 N& Y
出错 ode45 (line 488)
! n& F: h, f7 D4 F% z6 `  w5 ~7 ?7 }/ e
        yout(:,idx) = yout_new;

该用户从未签到

2#
发表于 2021-2-5 20:39 | 只看该作者
你的代码有点问题。
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

推荐内容上一条 /1 下一条

EDA365公众号

关于我们|手机版|EDA365电子论坛网 ( 粤ICP备18020198号-1 )

GMT+8, 2025-10-31 14:57 , Processed in 0.125000 second(s), 23 queries , Gzip On.

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

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

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