knrt.net
当前位置:首页 >> 如何用mAtlAB求解二阶微分方程,以及程序实例 >>

如何用mAtlAB求解二阶微分方程,以及程序实例

1、对于解析值,你可以用dsolve()函数求解。如求微分方程x*y''+x﹡(y')^2-y'=0的解析解 >> syms y(x) >>Dy = diff(y);D2y = diff(y, 2); >>dsolve(x*D2y+x*(Dy)^2-Dy==0,'x') 运行结果 2、对于数值解,你可以查看二阶常微分方程这个例子,http...

如何用matlab求解二阶微分方程,对于一般的微分方程(组)可以用dsolve()函数求得其解析解或数值解,对于较复杂的微分方程(组)可以用ode45()函数求得其数值解。 例如:微分方程,用dsolve和ode45计算t=0.1,0.2时y的值。 y = dsolve('D2y+0...

用matlab的dsolve函数计算该二阶微分方程的解析值是有一定的难度。只能用数值分析的方法来求解其数值解。当然求其数值解,必须还得已知a、θ的具体数值。该问题可以用matlab的ode45函数求得。求解方法如下: 1、建立微分方程的自定义函数,odefun...

对于解析值,你可以用dsolve()函数求解。如求微分方程x*y''+x﹡(y')^2-y'=0的解析解 >> syms y(x) >>Dy = diff(y);D2y = diff(y, 2); >>dsolve(x*D2y+x*(Dy)^2-Dy==0,'x') 运行结果

一、理论解,用dsolve函数(在command window 中输入doc dsolve可以查看帮助) 示例: 代码: [x,y]=dsolve('D2x+2*Dx=x+2*y-exp(-t)','Dy=4*x+3*y+4*exp(-t)') 二、数值解,用ode45,或ode23, ode15s其他函数。(在command window 中输入doc od...

这个不是解方程问题,是求状态方程求解的问题。 你可以找本现代控制理论看下,里面有这种方程的求解公式。

参考代码: m=10600; k=31986497; c=8.9E4;dy = @(t,y)[y(2); (k*y(1)-c*y(2)/abs(y(2))*y(2)^2)/m];ode45(dy,[0 5],[0 3.66]) 其中,蓝线为x,绿线为x'。

用Matlab编程计算二阶常微分方程组的数值解的方法有很多种,最常用有ode45函数。 ode45函数的使用格式: 【t,y】=ode45(odefun,tspan,y0) 例如: >> odefun=@(t,y)[y(1)+2*y(2);3*y(1)+2*y(2)]; %自定义函数 >> tspan=[0 20]; %时间 >> y0=[1 1]...

首先z1z2z3的初值应该给出的, 否则无法求解,你少定解条件!!!我自己给了一个初始条件令他们都为0,取100个时间点为区间,然后进行了求解,结果如下: 你这个方程组好像结构动力学中的mck方程,结果类似正弦波,程序你可以私信我,建议你去再...

function dy=dong22(t,y) m=1; k=2; g= %注意输入g l0 l0= dy=[y(2); y(1)*y(4)^2+g*cos(y(3))-(k/m)*(y(1)-l0); y(4); -(2/y(1))*y(2)*y(4)-(g/y(1))*sin(y(3))]; 然后和一次微分方程的一样。

网站首页 | 网站地图
All rights reserved Powered by www.knrt.net
copyright ©right 2010-2021。
内容来自网络,如有侵犯请联系客服。zhit325@qq.com