[求助]关于微分方程数值解

无水1324 · 发表于 2006-3-24 11:40

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

请多指教！ 假如我用ODE45求解一个微分方程： x''+c(t)x'+k(t)f(x)=F(t) f(x)=x-1,x>1 f(x)=0,abs(x)<1 f(x)=x+1,x<-1 怎么样编程请多指教 谢谢！！

happy · 发表于 2006-3-25 19:03

在odefun中 f(x>1)=x-1; f(abs(x)<1)=0; f(x<-1)=x+1; 后面的正常写就行了

无水1324 · 发表于 2006-3-25 19:09

谢谢，教授！ function xx=x1(t,x) if x(1)>1; f=x(1)-1 elseif x(1)<-1; f=x(1)+1; else f=0 end w=0.5;k=0.2; xx=[x(2);-0.4*x(2)-f*(1+k*sin(w*t))+0.2+0.5*w^2*sin(w*t)]; %%变参数的微分方程求解 [t,x]=ode45(@x1,[0,100],[1.2,0.05]); subplot(121);plot(t,x);subplot(122);plot(x(:,1),x(:,2)); 是不是这样

happy · 发表于 2006-3-25 19:15

<DIV class=quote>以下是引用无水1324在2006-3-25 19:09:57的发言： 谢谢，教授！ function xx=x1(t,x) if x(1)>1; f=x(1)-1 elseif x(1)<-1; f=x(1)+1; else f=0 end w=0.5;k=0.2; xx=[x(2);-0.4*x(2)-f*(1+k*sin(w*t))+0.2+0.5*w^2*sin(w*t)]; %%变参数的微分方程求解 [t,x]=ode45(@x1,[0,100],[1.2,0.05]); subplot(121);plot(t,x);subplot(122);plot(x(:,1),x(:,2)); 是不是这样 </DIV>
 不用这么麻烦，红色标记部分改成 f(x>1)=x-1; f(abs(x)<1)=0; f(x<-1)=x+1;

happy · 发表于 2006-3-25 19:16

另外注意，你这里少了abs(x)=1的状态

无水1324 · 发表于 2006-3-25 19:28

哦1!!这个运行的时候有问题 错误信息如下： ??? In an assignment A(I) = B, the number of elements in B and I must be the same.
Error in ==> C:\MATLAB6p5\work\x2.m On line 2 ==> f(x>1)=x-1;
Error in ==> C:\MATLAB6p5\toolbox\matlab\funfun\private\odearguments.m On line 104 ==> f0 = feval(ode,t0,y0,args{:});
Error in ==> C:\MATLAB6p5\toolbox\matlab\funfun\ode45.m On line 155 ==> [neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, args, ...

happy · 发表于 2006-3-26 19:04

<DIV class=quote>以下是引用happy在2006-3-25 19:15:38的发言： 不用这么麻烦，红色标记部分改成 f(x>1)=x-1; f(abs(x)<1)=0; f(x<-1)=x+1;</DIV>
 上面的x改成x(1)就行了

账号		自动登录	找回密码
密码			我要加入

[求助]关于微分方程数值解

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

回复：（无水1324）[求助]关于微分方程数值解

回复：（无水1324）谢谢，教授！function xx=x1(t,x)i...

回复：（无水1324）[求助]关于微分方程数值解

回复：（happy）回复：（无水1324）谢谢，教授！funct...