[求助]大家帮我看看这个程序的一个warning，怎么解决

lujixx · 发表于 2005-12-26 13:27

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

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

x

m文件： function dx=rotor_fun(t,x); %canshu g=9.8; c1=1.25*10^(-5); delta=4.5852*10^(-2); beta=1.1404*10^(-2); omega=50; omegag=omega*g^(1/2)/c1^(1/2); dx=zeros(4,1); %fangcheng dx(1,1)=x(3); dx(2,1)=x(4); dx(3,1)=beta*cos(t)-delta/omega*(((2*x(1)^2+2*x(2)^2-4*x(1)*x(4)+4*x(2)*x (3))/(1-x(1)^2-x(2)^2)^(1/2)/(2+x(1)^2+x(2)^2)+(pi-16/pi/(2+x(1)^2+x(2)^2))* (x(1)*x(3)+x(2)*x(4))/(x(1)^2+x(2)^2)^(1/2)/(1-x(1)^2-x(2)^2)^(3/2))*cos(t)+ (pi*(x(1)^2+x(2)^2)^(1/2)*(1-2*(x(1)*x(4)-x(2)*x(3))/(x(1)^2+x(2)^2))/(1-x(1 )^2-x(2)^2)^(1/2)/(2+x(1)^2+x(2)^2)+4*(x(1)*x(3)+x(2)*x(4))/(2+x(1)^2+x(2)^2 )/(1-x(1)^2-x(2)^2))*sin(t)); dx(4,1)=beta*sin(t)-1/omega^2 -delta/omega*(((2*x(1)^2+2*x(2)^2-4*x(1)*x( 4)+4*x(2)*x(3))/(1-x(1)^2-x(2)^2)^(1/2)/(2+x(1)^2+x(2)^2)+(pi-16/pi/(2+x(1)^ 2+x(2)^2))*(x(1)*x(3)+x(2)*x(4))/(x(1)^2+x(2)^2)^(1/2)/(1-x(1)^2-x(2)^2)^(3/ 2))*sin(t)-(pi*(x(1)^2+x(2)^2)^(1/2)*(1-2*(x(1)*x(4)-x(2)*x(3))/(x(1)^2+x(2) ^2))/(1-x(1)^2-x(2)^2)^(1/2)/(2+x(1)^2+x(2)^2)+4*(x(1)*x(3)+x(2)*x(4))/(2+x( 1)^2+x(2)^2)/(1-x(1)^2-x(2)^2))*cos(t)); m文件： clear all; clc; x0=[0.001; 0; 0; 0]; t0=[0:10000]; [t,x]=ode23s('rotor_fun',t0,x0); x plot(t,x(:,1)); 。Warning: Failure at t=1.251361e+002. Unable to meet integration tolerance s without reducing the step size below the smallest value allowed (4.445726e -013) at time t. (Type "warning off MATLAB:ode23s:IntegrationTolNotMet" to suppress this warn ing.) > In D:\MATLAB\toolbox\matlab\funfun\ode23s.m at line 431 In D:\MATLAB\work\rotor.m at line 5

xueyi · 发表于 2005-12-26 16:36

warning的意思是说：程序运行到t=1.251361e+002时出错。 在时间t时，如果不减小步长，就不能满足积分误差最低容许值4.445726e-013，也就是说误差在这一点大于4.445726e-013。 它给的解决办法是编辑ode。m使其误差大于默认值时不报错， 我建议尽量不要改变其源程序，如果你的数学水平不是很高的话。 问题可以通过减小积分步长来解决

lujixx · 发表于 2005-12-26 18:01

在不改变源程序的情况下，怎么改变步长啊？ options？

xueyi · 发表于 2005-12-27 16:38

t0=[0:10000]; 设置是步长，默认为1 你可以这样设t0=[0:0.1：10000];此时步长就是0.1

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