MATLAB中用ODE45解两个2阶非线性常微分方程组,Command Window出现警告：写在补充里面>> [t,y]=ode45(@RK4_shiyan01,[0:0.001:1000],[0,0,0,0])Warning: Failure at t=1.768731e-001. Unable to meet integration tolerances without reducin

来源：学生作业帮助网 编辑：六六作业网 时间：2024/11/27 09:00:04
MATLAB中用ODE45解两个2阶非线性常微分方程组,CommandWindow出现警告：写在补充里面>>[t,y]=ode45(@RK4_shiyan01,[0:0.001:1000],[0,0,

MATLAB中用ODE45解两个2阶非线性常微分方程组,Command Window出现警告：写在补充里面>> [t,y]=ode45(@RK4_shiyan01,[0:0.001:1000],[0,0,0,0])Warning: Failure at t=1.768731e-001. Unable to meet integration tolerances without reducin
MATLAB中用ODE45解两个2阶非线性常微分方程组,Command Window出现警告：写在补充里面
>> [t,y]=ode45(@RK4_shiyan01,[0:0.001:1000],[0,0,0,0])
Warning: Failure at t=1.768731e-001. Unable to meet integration tolerances without reducing the step
size below the smallest value allowed (4.440892e-016) at time t.
> In ode45 at 371
方程如下：
function dx=RK4_shiyan01(t,x)
a1=4.73001;a2=7.85320;b1=0.005;b2=0.2;b3=10;b4=50;b5=7.6;b6=20;b7=0.05;b8=20;b12=0.8355;b21=-0.8355;e=-7.5451;g11=-12.3025;g22=-46.0525;
h11=-6.1513;h12=-8.3806;h21=-6.7096;h22=-23.0263;d=0.083;
dx=[x(3);x(4);
-(a1^4+(b5^2-0.5*b4)*g11)*x(1)-b4*e*x(2)-b1*a1^4*x(3)-2*b2^0.5*b5*b12*x(4)+d*b6^2*sin(b6*t)*(sinh(a1)-sin(a1)-((cosh(a1)-cos(a1))/(sinh(a1)-sin(a1)))*(cos(a1)+cosh(a1)-2)/a1-b7*g11*x(1)*(2*b5*sin(b8*t)+b2^0.5*b8*cos(b8*t))+b2^0.5*b8*b7*(h11*x(1)+h12*x(2))*cos(b8*t)-2*b2^0.5*b7*b12*x(4)*sin(b8*t)-b3*((x(1)^2+2*b1*x(1)*x(3))*g11+((x(2)^2+2*b1*x(2)*x(4))*g22)*x(1)*g11));
-(a2^4+(b5^2-0.5*b4)*g22)*x(2)-b4*e*x(1)-b1*a2^4*x(4)-2*b2^0.5*b5*b21*x(3)+d*b6^2*sin(b6*t)*(sinh(a2)-sin(a2)-((cosh(a2)-cos(a2))/(sinh(a2)-sin(a2)))*(cos(a2)+cosh(a2)-2)/a2-b7*g22*x(2)*(2*b5*sin(b8*t)+b2^0.5*b8*cos(b8*t))+b2^0.5*b8*b7*(h21*x(1)+h22*x(2))*cos(b8*t)-2*b2^0.5*b7*b21*x(3)*sin(b8*t)-b3*((x(1)^2+2*b1*x(1)*x(3))*g11+((x(2)^2+2*b1*x(2)*x(4))*g22)*x(2)*g22))];

MATLAB中用ODE45解两个2阶非线性常微分方程组,Command Window出现警告：写在补充里面>> [t,y]=ode45(@RK4_shiyan01,[0:0.001:1000],[0,0,0,0])Warning: Failure at t=1.768731e-001. Unable to meet integration tolerances without reducin
警告的翻译：
无法满足积分条件,因此在时间t=1.768731e-001s处无法积分
帮你运行了一下,原因是在1.768731e-001s处,有个值趋于负的无穷大,因此无法积分
可能的原因：
1.方程出错；
2.时间t的范围；
3.初值

在matlab中用 ode45解方程 y'=x-y,y(0)=1,0 MATLAB中用ODE45解两个2阶非线性常微分方程组,Command Window出现警告：写在补充里面>> [t,y]=ode45(@RK4_shiyan01,[0:0.001:1000],[0,0,0,0])Warning: Failure at t=1.768731e-001. Unable to meet integration tolerances without reducin matlab中用ode45求微分方程的数值解问题[T,Y]=ode45('function1',[0 3000],[2 0]);其中[0 3000]是起始时间和结束时间.那步长是在哪里设置的,怎么没有步长呢? 用matlab中ode45解方程组 x1'=-2*x1+x2 x2'=-3*x1+6*x2 初始值自定 求教matlab高手ode45的用法!dsolve('Dt=m/(p/v-k*v^2)','t(0)=0','v')我懂得用这个函数解出微分方程的解析解,但不知道怎么用数值解（ode45）解这个微分方程,要怎么用ode45? 应用Matlab的ode45解含有约束条件的多参数微分方程 matlab ode45求解微分方程 D2y-0.01*Dy.^2+2*y=sin(t),初始值0 matlab中,ode45能解决含有（dy/dt）^2和d2y/d2t的方程么 matlab ode45 可以解二阶微分方程吗 怎样用matlab中用ode45求解带有分段参数的微分方程组有如下微分方程组：dx(1)=x(2)*dc(t)/dt; 其中dc(t)/dt是对c(t)求t的导数dx(2)=x(1)+a;c(t)=1/e(t);e(t)的方程式是知道的但是很复杂,通过手工算出dc(t)/dt matlab中的ode45怎样设置定步长解微分方程?matlab中的ode45是变步长解法,定步长在simulink里有ode4和ode5,但是在matlab里面好像没有现成的函数,听说ode45可以实现定步长的输出,如何在matlab中（不是simu matlab中ode45 及odeset的使用.在蝴蝶效应的解方程组中,方程是%x'=a*(X(2)-X(1))%y'=X(1)*(c-X(3))-X(2)%z'=X(1)*X(2)-b*X(3)使用ode45解方程.语句为lorenz=@(t,X) [a*(X(2)-X(1));X(1)*(c-X(3))-X(2);X(1)*X(2)-b*X(3)];[t X]=ode45(lorenz matlab解微分方程组求微分方程组的数值解并作图.要求是用ode45的方法来解 用MATLAB求微分方程分别用Euler法和ode45解下列常微分方程并与解析解比较 用matlab解二阶微分方程y''(t)= -3 cos(2t) +2sin(t)+t-3.8使用ode23或ode45 matlab用ode45 解微分方程如何用ODE45解微分方程 其中初值给出方程如下,请帮忙写出方程在程序中如何表达,我用function dy=func(x,y,r)dy = zeros(size(y));dy(1) = y(1)dy(2) = sqrt((1+r^2)/r^2*log((1+r^2)/(1+y(1)^2))+y(1 Matlab使用ode45解微分方程组的基本问题,下面有一个简单的微分方程组求代码d(X1)/dt=4X1+X2 ^2d(X2)/dt=2*d(X1)/dt+5*X2初始条件是X1=X2=0要用ode45解这个微分方程组该怎么写代码呢? matlab ode45 求方程（1+x^2)y''=2xy',y(0)=1,y'(0)=3 的数值解,我的运行不行,我编的函数文件：function dy=sdre(x,y) dy=[y(2);2*x/(1+x^2)*y(2)]; 程序 [x,y]=ode45('sdre',[0 10],[1 3]);plot(x,y)