| 查看: 760 | 回复: 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客户端 |
回复此楼» 猜你喜欢
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有4人回复
-
跳槽后在研项目怎么办？
已经有10人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有5人回复
-
什么时候开奖？
已经有10人回复
-
2026国自然函评费到账
已经有16人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有9人回复
-
让我中一个面上吧！
已经有14人回复
-
今天放榜吗？
已经有16人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有10人回复
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
40