24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

容墨竹止

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 137
帖子: 2
在线: 4.9小时
虫号: 4279878
注册: 2015-12-10
专业: 动力学与控制

[求助] 【ode45】matlab 用积分出现关于积分公差的警告已有1人参与

这个问题一直解决不了，又很急，希望有懂的朋友帮帮我。

代码如下：

M文件:
function dx=F_LTV_0(t,x)
dx=zeros(6,1);
global m n b_0 w_c

r0=10;
re=13;

r=r0+(1-exp(-t))*(re-r0);

dr=exp(-t)*(re-r0);
ddr=-exp(-t)*(re-r0);

rr=sqrt(x(1)^2+x(2)^2+x(3)^2);

x10 = sqrt(1/6)*r;
x20 = 0;
x30 = sqrt(5/6)*r;
x40 = sqrt(1/6)*dr;
x50 = 0;
x60 = sqrt(5/6)*dr;
x70 = sqrt(1/6)*ddr;
x80 = 0;
x90 = sqrt(5/6)*ddr;

DX=[x(1)-x10
x(2)-x20
x(3)-x30
x(4)-x40
x(5)-x50
x(6)-x60];

% % %

p1=(-x(1)^2-x(2)^2+2*x(3)^2)*x(5)-3*x(2)*x(3)*x(6)+w_c*x(1)*(x(1)^2+x(2)^2-2*x(3)^2);
p2=(x(1)^2+x(2)^2-2*x(3)^2)*x(4)+3*x(1)*x(3)*x(6)+w_c*x(2)*(x(1)^2+x(2)^2-2*x(3)^2);
p3=3*x(3)*(x(2)*x(4)-x(1)*x(5)+w_c*(x(1)^2+x(2)^2));

P=[ 0;
0;
0;
p1;
p2;
p3]';

%
%
Q=[2*sign(x(1)-x10),2*sign(x(2)-x20),sign(x(3)-x30),2*n*x(5)+3*n^2*x(1)-x70,-2*n*x(4),-n^2*x30-x90];

AA=Q*DX;
BB=P*DX;

U=-m*rr^5/b_0*(AA/BB);

% % 状态方程------------------------------

dx(1)=x(4);
dx(2)=x(5);
dx(3)=x(6);
dx(4)=2*n*x(5)+3*n^2*x(1)+(U)*b_0/(m*rr^5)*p1;
dx(5)=-2*n*x(4)+(U)*b_0/(m*rr^5)*p2;
dx(6)=-n^2*x(3)+(U)*b_0/(m*rr^5)*p3;

%---------------------------------------------

dx=[dx(1);dx(2);dx(3);dx(4);dx(5);dx(6)];
end

主程序：
close all;
clear;
clc
global m n b_0 w_c

%% Constants
m=50;
n=2*pi/5400/2;
b_0=20;
w_c=500*pi;

r0 = 10;
x_1=sqrt(1/6)*r0;
x_2=0;
x_3=sqrt(5/6)*r0;
x_4=0;
x_5=0;
x_6=0;

%% Solve ODEs
xe= [x_1;x_2;x_3;x_4;x_5;x_6];
tspan = [0 10000];
options =odeset('RelTol',1e-9,'AbsTol',1e-9);
[t,x] = ode45(@F_LTV_0,tspan,xe,options);
y = x';

出现了如下警告并且ode不再继续运行。
警告: 在 t=1.596171e+00 处失败。在时间 t 处，若不将步长降至允许的最小值(3.552714e-15)以下，积分公差要求无法满足。

我把ode45换成其他的ode积分也不行。除非把绝对误差和相对误差调大到10^(-5)以上，但是这样出来的结果也不对。我怀疑是U的设计不对，但是反复检查后确实表达式是这样。

回复此楼