| 查看: 724 | 回复: 3 | |||
happyday2016新虫 (小有名气)
|
[求助]
求助，用MATLAB求解这个一阶微分方程，为啥一运行就出错，求指导 已有1人参与
|
|
求大神指导，为啥一运行MATLAB这个程序就出现这个错误 Index exceeds matrix dimensions.Error in Untitled7 (line 18) t=0:0.01:10; omg0=2.87; %单位 MHz omgm=0; %单位 MHz omg=2.5e6 %单位 rad/s y1(1)=0; y2(1)=0; y3(1)=0; y4(1)=0; a=(omg/2)*i; b=omg0*i; c=omgm*i; for n=1:1000; y1(n+1)=y1(n)+a.*(y3(n)-y2(n)).*0.01; y2(n+1)=(1+0.01*(b-c)).*y2(n)-a.*(y1(n)-y4(n)).*0.01; y3(n+1)=(1+0.01*(c-b)).*y3(n)+a.*(y1(n)-y4(n)).*0.01; y4(n)=1-y1(n); end t=0:0.01:10; hold on subplot(2,2,1); plot(t,y1,'r') subplot(2,2,2); plot(t,y2,'b') subplot(2,2,3); plot(t,y3,'k') subplot(2,2,4); plot(t,y4,'g') 一运行就出现这个错误 Index exceeds matrix dimensions. Error in Untitled7 (line 18) y2(n+1)=(1+0.01*(b-c)).*y2(n)-a.*(y1(n)-y4(n)).*0.01; 发自小木虫Android客户端 |
回复此楼» 猜你喜欢
微信指数没变化，科研之友没阅读
已经有18人回复
-
基础研究怎么拉横向，学校到款任务越来越多，难以完成 拉横向，都有哪些途径啊
已经有8人回复
-
HXDI做水性聚氨酯乳液，是不是特别容易出渣
已经有3人回复
-
小木虫看见有人已经知道结果了
已经有19人回复
-
2026年国自然面上资助率
已经有8人回复
-
系统今天又提示维护了，估计离放榜不远了
已经有15人回复
-
申请2027年材料/化学类博士
已经有5人回复
-
你们的时间戳变了吗
已经有4人回复
-
今年的WR进展到哪一步了？
已经有9人回复
-
产物和副产物价值比较
已经有5人回复
2楼2017-08-17 22:05:58
wurongjun
专家顾问 (职业作家)
-
专家经验: +831 - 数学EPI: 9
- 应助: 791 (博后)
- 贵宾: 0.308
- 金币: 24609
- 散金: 310
- 红花: 75
- 帖子: 3004
- 在线: 881.4小时
- 虫号: 1368482
- 注册: 2011-08-14
- 性别: GG
- 专业: 计算数学与科学工程计算
- 管辖: 数学
3楼2017-08-18 23:32:48
happyday2016
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 366
- 散金: 70
- 帖子: 146
- 在线: 20.7小时
- 虫号: 4997281
- 注册: 2016-09-09
- 专业: 凝聚态物性 II ：电子结构
4楼2017-08-19 20:51:32