24小时热门版块排行榜    

Znn3bq.jpeg
查看: 587  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿211,化学学硕,310分,本科重点双非,求调剂 +17 努力奋斗112 2026-04-06 20/1000 2026-04-11 00:31 by wangjihu
[考研] 368求调剂 +3 17385968172 2026-04-10 3/150 2026-04-10 20:12 by 电子奥德彪
[考研] 347材料专硕求调剂 +19 zj8215216 2026-04-06 19/950 2026-04-10 09:36 by 690616278
[考研] 一志愿 江南大学 085602 化工专硕 338分求调剂 +16 路痴小琪 2026-04-05 16/800 2026-04-10 08:08 by kangsm
[考研] 269求调剂 +7 跪求收留。 2026-04-04 7/350 2026-04-09 19:06 by 探123
[考研] 考研调剂 +13 冰冰,,, 2026-04-07 13/650 2026-04-09 17:01 by Lilly_Li
[考研] 083200 初试305分 求调剂 暂不考虑跨专业 +15 Claireyyyy 2026-04-09 15/750 2026-04-09 16:11 by zhuimr
[考研] 求调剂 +8 吃口冰激凌 2026-04-07 8/400 2026-04-09 08:03 by 5268321
[考研] 085404,334分,求调剂 +5 sunjie8888 2026-04-08 8/400 2026-04-09 07:26 by sunjie8888
[考研] 化学308分求调剂 +21 你好明天你好 2026-04-07 23/1150 2026-04-08 22:32 by 凯凯要变帅
[考博] 材料方向考博,求推荐 +3 言语aaa 2026-04-05 4/200 2026-04-08 22:22 by nxgogo
[考研] 287求调剂 +6 Fnhc 2026-04-07 6/300 2026-04-08 10:05 by xingguangj
[考研] 071000生物学,一志愿深圳大学296分,求调剂 +12 TIckLw 2026-04-06 13/650 2026-04-07 20:34 by lijunpoly
[考研] 325 调剂 +6 QQ小虾 2026-04-07 6/300 2026-04-07 15:17 by Ccclqqq
[考研] 081700学硕,323分,一志愿中国海洋大学求调剂学校 +19 披星河 2026-04-04 19/950 2026-04-07 15:00 by 上岸快快
[考研] 考研调剂 +5 美丽的youth_ 2026-04-04 6/300 2026-04-06 06:57 by houyaoxu
[考研] 296求调剂 +3 汪!?! 2026-04-05 4/200 2026-04-05 20:13 by 啵啵啵0119
[考研] 308求调剂 +3 终不似从前 2026-04-05 3/150 2026-04-05 20:07 by 啵啵啵0119
[考研] 326求调剂 +3 顾若浮生 2026-04-05 3/150 2026-04-05 18:32 by 蓝云思雨
[考研] 313求调剂 +5 海日海日 2026-04-04 7/350 2026-04-05 13:58 by imissbao
信息提示
请填处理意见