| 查看: 3780 | 回复: 2 | ||||
[求助]
小白求助一个matlab初级问题“??? Undefined function or variable”
|
|
运行一个模型，偶尔能出结果图形，但多数情况下总是出现??? Undefined function or variable “t”或“y”。这应该是个初级问题，但本人纯粹新手，不知如何解决，抑郁了两天，向各位大侠求助！程序如下： function dy=ODE45_fun(t,y) a=5;b=4.5;c=3;d=3;e=4.9;f=4;g=3.5;h=5; dy(1)=y(1).*(1-y(1)).*((f-e+g-h).*y(2)-(g-h)); dy(2)=y(2).*(1-y(2)).*((a-c)-(a-c+d-b).*y(1)); dy=[dy(1);dy(2)]; ODE45_fun(t,y) %ODE45主函数，用以计算yh方程 %tspan为时间区间 tspan=[0,50]; %y0为初始值 y0=[0.45;0.76]; %调用函数进行计算 [t,y]=ode45('ODE45_fun',tspan,y0); %其他设置为缺省默认 data=[t,y]; save ODE45_data.txt data -ascii %数据太大，以ASCII文件保存 plot(t,y(:,1),'-',t,y(:,2),'m') grid on title('yh方程') |
回复此楼» 猜你喜欢
289求调剂
已经有5人回复
-
一志愿 南京航空航天大学大学 ，080500材料科学与工程学硕
已经有5人回复
-
一志愿西安交通大学 学硕 354求调剂
已经有3人回复
-
317求调剂
已经有11人回复
-
287求调剂
已经有6人回复
-
一志愿南昌大学，327分，材料与化工085600
已经有7人回复
-
材料080500调剂求收留
已经有8人回复
-
353求调剂
已经有3人回复
-
北科281学硕材料求调剂
已经有3人回复
-
A区线材料学调剂
已经有3人回复
1592203609
荣誉版主 (知名作家)
- 仿真EPI: 2
- 应助: 254 (大学生)
- 贵宾: 0.37
- 金币: 35191.9
- 散金: 468
- 红花: 164
- 沙发: 5
- 帖子: 5653
- 在线: 2195.1小时
- 虫号: 1766262
- 注册: 2012-04-19
- 性别: GG
- 专业: 有机合成
- 管辖: 仿真模拟
2楼2013-06-07 09:16:08
dbb627
荣誉版主 (著名写手)
-
专家经验: +4 - 仿真EPI: 6
- 应助: 289 (大学生)
- 贵宾: 0.589
- 金币: 24640.4
- 散金: 551
- 红花: 61
- 沙发: 1
- 帖子: 1246
- 在线: 1794.8小时
- 虫号: 149791
- 注册: 2005-12-29
- 性别: GG
- 专业: 污染控制化学
- 管辖: 计算模拟
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
1592203609: 金币+1, 谢谢回帖 2013-06-07 11:59:34
shirleyran: 金币+10, ★★★★★最佳答案, 你后面帮我修正的程序一下就运算出来了，激动得真的快掉眼泪了！金币有限，难以表达感激之情，谢谢！！！ 2013-06-07 15:48:38
感谢参与，应助指数 +1
1592203609: 金币+1, 谢谢回帖 2013-06-07 11:59:34
shirleyran: 金币+10, ★★★★★最佳答案, 你后面帮我修正的程序一下就运算出来了，激动得真的快掉眼泪了！金币有限，难以表达感激之情，谢谢！！！ 2013-06-07 15:48:38
|
估计你的ODE45_fun没有存m文件吧，你也可以这样写 a=5;b=4.5;c=3;d=3;e=4.9;f=4;g=3.5;h=5; ODE45_fun=@(t,y)[y(1).*(1-y(1)).*((f-e+g-h).*y(2)-(g-h));y(2).*(1-y(2)).*((a-c)-(a-c+d-b).*y(1))]; %ODE45主函数，用以计算yh方程 %tspan为时间区间 tspan=[0,50]; %y0为初始值 y0=[0.45;0.76]; %调用函数进行计算 [t,y]=ode45(ODE45_fun,tspan,y0); %其他设置为缺省默认 data=[t,y]; save ODE45_data.txt data -ascii %数据太大，以ASCII文件保存 plot(t,y(:,1),'-',t,y(:,2),'m') grid on title('yh方程') |
3楼2013-06-07 11:23:52
