| 查看: 6896 | 回复: 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的设计不对，但是反复检查后确实表达式是这样。 |
回复此楼» 猜你喜欢
求助中医院临床中医中标的青年标书学习，谢谢
已经有6人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有259人回复
-
物理学I论文润色/翻译怎么收费?
已经有198人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
课题组招1-2名硕士
已经有1人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有1人回复
-
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
apl投稿状态
已经有0人回复
chendequan
铁虫 (小有名气)
- 应助: 152 (高中生)
- 金币: 918.6
- 红花: 17
- 帖子: 262
- 在线: 209小时
- 虫号: 994474
- 注册: 2010-04-11
- 专业: 结构工程
2楼2017-04-14 10:46:35
3楼2017-04-14 20:53:34
