[求助]关于微分方程数值解
请多指教！<BR>假如我用ODE45求解一个微分方程：<BR>x''+c(t)x'+k(t)f(x)=F(t)<BR><BR>f(x)=x-1,x>1<BR><BR>f(x)=0,abs(x)<1<BR><BR>f(x)=x+1,x<-1<BR><BR>怎么样编程请多指教<BR>谢谢！！<BR><BR><BR>回复：（无水1324）[求助]关于微分方程数值解
在odefun中<BR><BR>f(x>1)=x-1;<BR>f(abs(x)<1)=0;<BR>f(x<-1)=x+1;<BR><BR>后面的正常写就行了 谢谢，教授！<BR>function xx=x1(t,x)<BR>if x(1)>1;<BR> f=x(1)-1<BR>elseif x(1)<-1;<BR> f=x(1)+1;<BR>else<BR> f=0<BR>end<BR>w=0.5;k=0.2;<BR>xx=;<BR><BR>%%变参数的微分方程求解<BR>=ode45(@x1,,);<BR>subplot(121);plot(t,x);subplot(122);plot(x(:,1),x(:,2));<BR>是不是这样<BR>回复：（无水1324）谢谢，教授！function xx=x1(t,x)i...
<DIV class=quote><B>以下是引用<I>无水1324</I>在2006-3-25 19:09:57的发言：</B><BR>谢谢，教授！<BR>function xx=x1(t,x)<BR><FONT color=#ff0000>if x(1)>1;<BR> f=x(1)-1<BR>elseif x(1)<-1;<BR> f=x(1)+1;<BR>else<BR> f=0<BR>end<BR></FONT>w=0.5;k=0.2;<BR>xx=;<BR><BR>%%变参数的微分方程求解<BR>=ode45(@x1,,);<BR>subplot(121);plot(t,x);subplot(122);plot(x(:,1),x(:,2));<BR>是不是这样<BR> </DIV><br>不用这么麻烦，红色标记部分改成<BR>f(x>1)=x-1;<BR>f(abs(x)<1)=0;<BR>f(x<-1)=x+1;
回复：（无水1324）[求助]关于微分方程数值解
另外注意，你这里少了abs(x)=1的状态 <P>哦1!!这个运行的时候有问题<BR>错误信息如下：<BR>???In an assignmentA(I) = B, the number of elements in B and<BR> I must be the same.</P><P>Error in ==> C:\MATLAB6p5\work\x2.m<BR>On line 2==> f(x>1)=x-1;</P>
<P>Error in ==> C:\MATLAB6p5\toolbox\matlab\funfun\private\odearguments.m<BR>On line 104==> f0 = feval(ode,t0,y0,args{:});</P>
<P>Error in ==> C:\MATLAB6p5\toolbox\matlab\funfun\ode45.m<BR>On line 155==> [neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, args, ...<BR></P>
回复：（happy）回复：（无水1324）谢谢，教授！funct...
<DIV class=quote><B>以下是引用<I>happy</I>在2006-3-25 19:15:38的发言：</B><BR><BR><BR>不用这么麻烦，红色标记部分改成<BR>f(x>1)=x-1;<BR>f(abs(x)<1)=0;<BR>f(x<-1)=x+1;</DIV><br>上面的x改成x(1)就行了
页:
[1]