| 查看: 6723 | 回复: 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的设计不对，但是反复检查后确实表达式是这样。 |
回复此楼» 猜你喜欢
**
已经有1人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有3人回复
-
物理学I论文润色/翻译怎么收费?
已经有236人回复
-
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有2人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有11人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有13人回复
-
桂林理工大学物理学专业招收调剂
已经有13人回复
-
VASP 的一组 GPU / CPU 基准测试记录
已经有0人回复
-
美国顶级物理期刊《应用物理快报》（APL）的编辑欺骗和歧视作者及AIP的官僚主义傲慢
已经有6人回复
-
津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生一名
已经有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
5