用 四阶龙格库塔 解三阶微分方程 公式和程序0.001y'''+0.11y''+y'+10y=10初值X=0,y=0要三阶微分方程的推导公式,如果有程序可以加分.
来源:学生作业帮助网 编辑:作业帮 时间:2024/06/01 17:03:29
用 四阶龙格库塔 解三阶微分方程 公式和程序
0.001y'''+0.11y''+y'+10y=10
初值X=0,y=0
要三阶微分方程的推导公式,如果有程序可以加分.
初值给的不够啊.
不仅要给y,还要给y的一阶导和二阶导.否则数值解没法弄.
理论解的话,可能还可以含有一些系数.
在Matlab下输入:edit,然后将下面两行百分号之间的内容,复制进去,保存
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function DYDt=zhidao_amanm(t,Y)
y=Y(1);
dydt=Y(2);
d2ydt2=Y(3);
DYDt=[Y(2);Y(3);(10-10*y-dydt-0.11d2ydt2)/0.001];
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
在Matlab命令行下面输入:
t_start=0;
t_end=1;
yy0=[0.1;0;0]; %初值
[t,y]=ode45('zhidao_amanm',[0,t_end],yy0);
plot(t,y(:,1));
xlabel('t');
ylabel('y')
y(end,1)
得到的结果:
ans =
1.0056
用昨天给的C代码计算出来的结果是:
1.005630
用 四阶龙格库塔 解三阶微分方程 公式和程序0.001y'''+0.11y''+y'+10y=10初值X=0,y=0要三阶微分方程的推导公式,如果有程序可以加分.
用公式法求微分方程
微分方程和偏微分方程入门用什么书好
这个方程可以用一阶微分方程公式吗?详解
直接用一阶微分方程的公式怎么算不出来
matlab 用eig和diag解微分方程
偏微分方程和常微分方程的区别?
怎么区别线性微分方程和非线性微分方程?
怎么区别常微分方程和偏微分方程?
求一阶线性微分方程为什么用常数变易法,不直接用通解公式
用公式法求微分方程xy'+y=lnx/x通解
一阶线性微分方程怎么求通解 找不到思路 想直接用公式又配不好
微分方程 有没有一个万用的公式?
微分方程图中的公式怎么得到的
怎么用simulink解微分方程要有例子和解答,
怎样用公式法求微分方程:xy'+y=x^2+3x+2的通解和特解
用MATLAB求微分方程分别用Euler法和ode45解下列常微分方程并与解析解比较
用mathematica解决微分方程