24小时热门版块排行榜    

查看: 576  |  回复: 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客户端
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mdechuan

新虫 (小有名气)

2楼2017-08-17 22:05:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wurongjun

专家顾问 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
把y4(n)=1-y1(n);
放到
y2(n+1)=(1+0.01*(b-c)).*y2(n)-a.*(y1(n)-y4(n)).*0.01;
之前就行了!
善恶到头终有报,人间正道是沧桑.
3楼2017-08-18 23:32:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

happyday2016

新虫 (小有名气)

引用回帖:
3楼: Originally posted by wurongjun at 2017-08-18 23:32:48
把y4(n)=1-y1(n);
放到
y2(n+1)=(1+0.01*(b-c)).*y2(n)-a.*(y1(n)-y4(n)).*0.01;
之前就行了!

谢谢!

发自小木虫Android客户端
4楼2017-08-19 20:51:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 happyday2016 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 梁成伟老师课题组欢迎你的加入 +9 一鸭鸭哟 2026-03-14 11/550 2026-03-19 17:22 by !本暗一次!
[考研] 085600材料与化工调剂 324分 +10 llllkkkhh 2026-03-18 12/600 2026-03-19 14:33 by llllkkkhh
[考研] 求调剂 +3 Mqqqqqq 2026-03-19 3/150 2026-03-19 14:11 by peike
[考研] 化学求调剂 +3 临泽境llllll 2026-03-17 4/200 2026-03-19 13:59 by houyaoxu
[考研] 281求调剂(0805) +9 烟汐忆海 2026-03-16 19/950 2026-03-19 11:42 by laoshidan
[考研] 328求调剂,英语六级551,有科研经历 +3 生物工程调剂 2026-03-17 7/350 2026-03-18 20:41 by Wangjingyue
[考研] 化学工程321分求调剂 +15 大米饭! 2026-03-15 18/900 2026-03-18 14:52 by haxia
[考研] 070300化学319求调剂 +6 锦鲤0909 2026-03-17 6/300 2026-03-18 13:22 by Iveryant
[考研] 278求调剂 +5 烟火先于春 2026-03-17 5/250 2026-03-18 08:43 by 星空星月
[考研] 考研化学学硕调剂,一志愿985 +4 张vvvv 2026-03-15 6/300 2026-03-17 17:15 by ruiyingmiao
[考研] 一志愿南京大学,080500材料科学与工程,调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
[考研] 一志愿,福州大学材料专硕339分求调剂 +3 木子momo青争 2026-03-15 3/150 2026-03-17 07:52 by laoshidan
[考研] 274求调剂 +5 时间点 2026-03-13 5/250 2026-03-17 07:34 by 热情沙漠
[考研] 0703化学调剂 290分有科研经历,论文在投 +7 腻腻gk 2026-03-14 7/350 2026-03-16 10:12 by houyaoxu
[考研] 288求调剂 +4 奇点0314 2026-03-14 4/200 2026-03-14 23:04 by JourneyLucky
[考研] 中科大材料与化工319求调剂 +3 孟鑫材料 2026-03-14 3/150 2026-03-14 20:10 by ms629
[基金申请] 现在如何回避去年的某一个专家,不知道名字 +3 zk200107 2026-03-12 6/300 2026-03-14 17:13 by zk200107
[考研] 复试调剂 +3 呼呼?~+123456 2026-03-14 3/150 2026-03-14 16:53 by WTUChen
[考研] 297求调剂 +4 学海漂泊 2026-03-13 4/200 2026-03-14 11:51 by 热情沙漠
[考研] 304求调剂 +7 7712b 2026-03-13 7/350 2026-03-13 21:42 by peike
信息提示
请填处理意见