| 查看: 6851 | 回复: 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的设计不对，但是反复检查后确实表达式是这样。 |
回复此楼» 猜你喜欢
各位大佬，求一份最新Wien2k版本作为学习用途
已经有0人回复
-
南科大活性流体和软物质课题组诚招2027级博士、硕士研究生和博后
已经有523人回复
-
物理学I论文润色/翻译怎么收费?
已经有266人回复
-
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向
已经有130人回复
-
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
PVA溶解
已经有2人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
3楼2017-04-14 20:53:34
chendequan
铁虫 (小有名气)
- 应助: 152 (高中生)
- 金币: 918.6
- 红花: 17
- 帖子: 262
- 在线: 208.7小时
- 虫号: 994474
- 注册: 2010-04-11
- 专业: 结构工程
2楼2017-04-14 10:46:35
