| 查看: 6702 | 回复: 2 | |||
[求助]
【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的设计不对，但是反复检查后确实表达式是这样。 |
回复此楼» 猜你喜欢
哈尔滨理工大学物理系招收物理学考研调剂
已经有7人回复
-
0702一志愿吉大B区求调剂
已经有5人回复
-
物理学I论文润色/翻译怎么收费?
已经有53人回复
-
求调剂
已经有0人回复
-
0702一志愿吉大B区求调剂有论文
已经有0人回复
-
请问还有没有用Latex写文章的小伙伴们？
已经有0人回复
-
光学工程学硕调剂信息
已经有31人回复
-
欢迎加入课题组
已经有0人回复
-
散金币，求好运，祝面上顺利！
已经有68人回复
-
华东师范大学芯片设计徐珑真实水平咋样
已经有1人回复
-
调剂
已经有0人回复
chendequan
铁虫 (小有名气)
- 应助: 152 (高中生)
- 金币: 918.6
- 红花: 17
- 帖子: 262
- 在线: 208.6小时
- 虫号: 994474
- 注册: 2010-04-11
- 专业: 结构工程
2楼2017-04-14 10:46:35
3楼2017-04-14 20:53:34
