应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 数学综合 » 求助，用MATLAB求解这个一阶微分方程，为啥一运行就出错，求指导
41/1返回列表
查看: 588  |  回复: 3
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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客户端
回复此楼

» 猜你喜欢
1楼 2017-08-17 21:37:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mdechuan

新虫 (小有名气)
y4（n+1）=1-y1（n）；

发自小木虫Android客户端
赞一下 回复此楼
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 的主题更新
41/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 295分求调剂 +6 ?要上岸? 2026-04-10 6/300 2026-04-10 21:13 by chemisry
[考研] 346，工科0854求调剂，专硕 +5 moser233 2026-04-10 6/300 2026-04-10 18:04 by 周邵阳
[考研] 085410-273求调剂 +4 X1999 2026-04-10 4/200 2026-04-10 15:29 by hemengdong
[考研] 085800 能源动力求调剂 +6 阿biu啊啊啊啊啊 2026-04-10 6/300 2026-04-10 15:03 by hemengdong
[考研] 调剂 +4 青灯不负 2026-04-09 4/200 2026-04-10 09:55 by LHGeng
[考研] 308求调剂 +21 倘若起风了呢 2026-04-05 21/1050 2026-04-10 08:13 by Sammy2
[考研] 085500求调剂材料 +8 易11122 2026-04-09 8/400 2026-04-09 23:15 by parmtree
[考研] 初试分332，一志愿报考西北工业大学， +11 故人?? 2026-04-09 11/550 2026-04-09 21:54 by JineShine
[考研] 085600材料与化工，求调剂 +11 won_qii 2026-04-07 11/550 2026-04-09 17:03 by luoyongfeng
[考研] 327求调剂 +10 Xxjc1107. 2026-04-06 11/550 2026-04-09 01:21 by lature00
[考研] 318求调剂 +13 ykyhsa 2026-04-05 15/750 2026-04-08 21:37 by wj165256
[考研] 298求调剂 +4 残荷新柳 2026-04-07 4/200 2026-04-07 23:02 by lbsjt
[考研] 363求调剂 +9 zh096 2026-04-04 9/450 2026-04-07 21:51 by 418490947
[考研] 22408 318分求调剂 +4 勤奋的小笼包 2026-04-06 6/300 2026-04-07 15:05 by 纸鹤555
[考研] 328求调剂 +4 ghhh88888 2026-04-06 5/250 2026-04-07 14:45 by ghhh88888
[考研] 372分材料与化工（085600）英二数二求调剂 +4 蓝笺片 2026-04-06 4/200 2026-04-07 12:30 by dongzh2009
[考研] 277工科求调剂 +10 1915668 2026-04-04 11/550 2026-04-06 23:53 by 心媛意北
[考研] 338求调剂 +4 我想上岸ii 2026-04-05 4/200 2026-04-06 21:04 by 木子君1218
[考研] 320求调剂 +3 一样圆 2026-04-04 3/150 2026-04-04 22:29 by 啵啵啵0119
[考研] 材料383求调剂 +5 郭阳阳阳成 2026-04-04 5/250 2026-04-04 19:06 by dongzh2009
信息提示
关闭
请填处理意见
关闭