原问题为利用ode45求解下面微分方程：
初值为
可以转为二阶微分方程
利用Matlab求解，下面为源代码
functionMyode45
tspan = [2,4]; %求解区间
y0 = [1 0]; %y0的初值形式为行向量
[t,x] = ode45(@odefun,tspan,y0); %使用ode45求解微分方程
plot(t,x(:,1),'-o',t,x(:,2),'-*') %作图
legend('y1','y‘')
title('y'' ''=-y-7y''y1^2+7')
xlabel('t')
ylabel('y')
function y =odefun(~,x) % 定义ode45函数的参数odefun
y = zeros(2,1); % 创建列向量
y(1) = x(2); % 表达式1：dx1/dt=x2 => y(1) =x1'=x(2)
y(2) = 7*(1-x(1)^2)*x(2)-x(1); % 表达式2：dx2/dt=7(1-x1^2)x2-x1 => y(2)=x2'=7*(1-x(1)^2)*x(2)-x(1)
end
end
原问题为利用ode45求解下面微分方程： 初值为 可以转为二阶微分方程 利用Matlab求解，下面为源代码 fun...
ode45 求解非刚性微分方程，中阶方法。 [TOUT,YOUT] = ode45(ODEFUN,TSPAN,Y0...
高阶微分方程 对于高阶微分方程 一般没有普遍的解法，求解告诫微分方程的基本思想是降阶，通过变量变换的方法将方程（2...
以下只是例题介绍，编程方式情参照例题，因为markdown用的不好，现在我都不知道怎样输入数学公式。所以私底下还得...
常微分方程 1、四阶经典步长Runge-Kutta法求解高阶常微分方程 2、用ode15s函数、ode23s函数和...
求解微分方程本来就是一个复杂艰难的问题，原因是其种类繁多，形式多样，所以增加了求解难度，完全掌握求解微分方程的办法...
在使用integrate.odeint求解微分方程时，有时会遇到如下错误： 我的代码出现的错误是odeint()的...
matlab 出现NaN是因为求解器不对
Problem 用Matlab软件包求解下列方程的全部实根： Methods roots函数 使用方法: r = ...
弦振动方程的解法 微分方程基础 本文所言的方程基本上是微分方程，而非中学阶段所言的代数方程，微分方程的求解目标一般...
本文标题：Matlab使用ode45求解高阶微分方程
本文链接：https://www.haomeiwen.com/subject/cueyhqtx.html
网友评论