Undefined function or variable 'e'.出错求助

bryant

程序目的：求解微分方程
主程序：>> e=0.1,xi=0.01,Omega=0.5,omega=1;
               >> [t,y]=ode45('odey',[0:0.1:10],[0 0]);
子程序：function dy=odey(t,y)
               dy=zeros(2,1);
               dy(1)=y(2);
               dy(2)=e*sin(omega*t)-2*xi/Omega*y(2)-1/Omega^2*y(1)；
运行时显示：??? Undefined function or variable 'e'.
请问怎样正确定义参数及赋值

eight

原帖由 bryant 于 2008-5-5 14:39 发表

                               
登录/注册后可看大图

程序目的：求解微分方程
主程序：>> e=0.1,xi=0.01,Omega=0.5,omega=1;
               >> [t,y]=ode45('odey',[0:0.1:10],[0 0]);
子程序：function dy=odey(t,y)
               dy=zeros(2,1);
             ...

请阅读基础书吧，或者看看置顶区 for 新手系列 的帖子，再不行就搜索版面，这个问题讨论过无数次

ch_j1985

原帖由 bryant 于 2008-5-5 14:39 发表

                               
登录/注册后可看大图

程序目的：求解微分方程
主程序：>> e=0.1,xi=0.01,Omega=0.5,omega=1;
               >> [t,y]=ode45('odey',[0:0.1:10],[0 0]);
子程序：function dy=odey(t,y)
               dy=zeros(2,1);
             ...

帮你把代码修改了一下，见附件：

woshiaq

怎么下载不了附件呢

zhouyang664

回复 4 # woshiaq 的帖子

不要使用迅雷，右键-->目标另存为...
图片传不上来，自己摸索摸索！

zhouyang664

不要用迅雷下载，右键-->目标另存为...
图片传不上来，自己摸索摸索！

ChaChing

回复 3 # ch_j1985 的帖子

个人以为若程序不大, 直接贴上比较方便阅读, 也不必扣体能, 所以善作主张, 帮忙贴下
若主任们或老k认为不妥, 请告知并删帖!

%%% odey1_case.m %%%%
global e xi Omega omega
e=0.1;xi=0.01;Omega=0.5;omega=1; t=0:0.1:10;
[t,y]=ode45('odey1',t,[0 0]); plot(t,y(:,1),'r',t,y(:,2),'b')
function dy=odey1(t,y)

%%% odey1.m %%%
global e xi Omega omega
dy=[y(2);e*sin(omega*t)-2*xi/Omega*y(2)-1/Omega^2*y(1)];
end

woshiaq

耗掉了我好几个体能就是没下载下来
支持ChaChing的做法

ChaChing

回复 8 # woshiaq 的帖子

不客气! 论坛本是相互学习的地方, 大家都是有施有得, 不是吗!?
也希望ls常来相互学习
小建议, 有空也参与回覆讨论:@)