| 查看: 540 | 回复: 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客户端 |
回复此楼» 猜你喜欢
回收溶剂求助
已经有7人回复
-
职称评审没过，求安慰
已经有40人回复
-
硝基苯如何除去
已经有3人回复
-
A期刊撤稿
已经有4人回复
-
垃圾破二本职称评审标准
已经有17人回复
-
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿
已经有22人回复
-
申请26博士
已经有5人回复
-
EST投稿状态问题
已经有7人回复
-
毕业后当辅导员了，天天各种学生超烦
已经有4人回复
-
求助文献
已经有3人回复
happyday2016
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 366
- 散金: 70
- 帖子: 146
- 在线: 20.7小时
- 虫号: 4997281
- 注册: 2016-09-09
- 专业: 凝聚态物性 II ：电子结构
4楼2017-08-19 20:51:32
2楼2017-08-17 22:05:58
wurongjun
专家顾问 (职业作家)
-
专家经验: +831 - 数学EPI: 9
- 应助: 791 (博后)
- 贵宾: 0.308
- 金币: 24607.5
- 散金: 310
- 红花: 75
- 帖子: 3004
- 在线: 881.2小时
- 虫号: 1368482
- 注册: 2011-08-14
- 性别: GG
- 专业: 计算数学与科学工程计算
- 管辖: 数学
3楼2017-08-18 23:32:48