mdsk.net
当前位置:首页 >> mAtlAB二阶微分方程数值解 >>

mAtlAB二阶微分方程数值解

用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]; %初值>> [t,y]= ode45(odefun,tspan,y0); >> A=[t,y];>> A运行结果

odefun=@(t,x)[x(2);3*x(2)-2*x(1)+1]; [t,y]=ode45(odefun,[0:0.01:2],[1 0]); plot(t,y) [t y] 结果 y(0.5000)=0.7896 y= dsolve('D2y-3*Dy+2*y=1','Dy(0)=0','y(0)=1');>> y y = exp(t) - exp(2*t)/2 + 1/2>> feval(@(t)exp(t) - exp(2*t)/2 + 1/2,0.5) ans = 0.7896

1、首先看一下matlab求解方程的方法,首先指明所解方程的变量,然后指明方程,未知数和限制条件,最后求解方程.2、在matlab命令行窗口中输入syms x [x,params,conds]=solve(sin(x)==1,'ReturnConditions', true) ,按回车键可以得到方程解.3、转换一下,可以看到sin(x)=1方程的解是如下图所示 .4、也可以求解下面的一个方程.5、输入syms a b c y x;[x,y]=solve([a*x^2+b*y+c==0,a*x+2*y==4],[x,y]).6、按回车键可以得到方程解.

1、把G=1/3*((5*Pp-2*P1)/(P1-2*Pp)-P1*c1^2/Pp*cp^2);改成G=1/3*((5*Pp-2*P1)/(P1+2*Pp)-P1*c1^2/(Pp*cp^2));表达式的错误导致出现高频振荡,积分步长非常小(大约在10^-8量级),所以很长时间算不出. 2、从结果看,仿真时间取50秒就足够了(已进入稳态). 3、在出现busy的情况下,可以考虑让 ode* 函数不返回参数,此种情况下,会自动绘制响应曲线,而且可以随时通过GUI停止运行,便于观察求解过程的具体情况.

x''+0.2x'-0.4x=0.2u(t), x(0)=x'(0)=0.1.你用Laplace Transform 做吧.L(x'')+0.2L(x')-0.4L(x)=0.2L(u(t))(p^2)*X(p)-p*x(0)-x'(0)+0.2(p*X(p)-x(0))-0.4*X(p)=0.2/p,(p^2+0.2*p-0.4)*X(p)=0.2/p,X(p)=0.2/(p^3+0.2

一个二阶微分方程:y''+y'+y=sin(t) 初始条件为y(0)=5,y'(0)=6.过程:先降阶为一阶微分方程组 y'=z z'=-z-y+sin(t) 编制如下函数m文件 function dy=weifen(t,x) dy=zeros(2,1);%y=x(1)%z=x(2) dy(1)=x(2); dy(2)=sin(t)-x(2)-x(1);******************* 然后用ode45解方程 [t,y]=ode45(@weifen,[0 20],[5 6]) plot(t,y)就可以画出y和y'的图像

用dsolve()函数,就可以解决.dsolve('3*D2x+500*Dx+2000*x','Dx(0)=2.5','x(0)=0.1')ans =(565^(1/2)*exp(t*((10*565^(1/2))/3 - 250/3))*(2*565^(1/2) + 65))/22600 + (565^(1/2)*(2*565^(1/2) - 65))/(22600*exp(t*((10*565^(1/2))/3 + 250/3))) %x(t)

y''-10y'+9y=e^{2x}令y=u*e^{2x},代入化简可得u''-6u'-7u-1=0令v=u+1/7得v''-6v'-7v=0解之得:v=A*e^{-x}+B*e^{7x}从而得:u=A*e^{-x}+B*e^{7x}-1/7y=(A*e^{-x}+B*e^{7x}-1/7)*e^{2x} =A*e^{x}+B*e^{9x}-(1/7)*e^{2x}y'=A*e^{x}+9B*e^

一、理论解,用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 中输入

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