clzm.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...

1、对于求数值解的微分方程,你可以用ode45()函数求解。如求下列微分方程 func。m %自定义微分方程的函数 function z = func(t,y) z =[y(2);(1-y(1)^2)*y(2)-y(1)]; main。m %主程序 clear all;close all;clc y0 = [0.25;0]; h = 0.1; a = 0; b ...

>> clear >> syms a b c d e; >> y=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e') y = (2*a*e + b*d + d*(b^2 - 4*a*c)^(1/2))/(2*exp((t*(b - (b^2 - 4*a*c)^(1/2)))/(2*a))*(b^2 - 4*a*c)^(1/2)) - (2*a*e + b*d - d*(b^2 - 4*a*c)^(1/2))/(...

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

参考代码: 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'。

用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 + ...

dsolve('D2y=x','x') ans = x^3/6 + C2*x + C3 D2y的意思就是y的二阶微分项 不明白你的问题什么意思,要输入的话直接定义符号变量输入... syms D2x x D2x=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...

你没给初值,我假设一个[1 -3] clc f=inline('[y(2);(1-x)/x*y(2)+3/x*y(1)]','x','y');[x y]=ode45(f,[eps 10],[1 -3]);z=x.^2.*y(:,1).^2;plot(x,y(:,1),x,z)xlabel('x')legend('y','{x^2}{y^2}') i无穷大; iii无穷大

令x1 = x, x2 = x', x3=y, x4=y' 得 x1' = x2; x2' = -x2*sqrt(x2^2+x4^2); x3' = x4; x4' = -x4*sqrt(x2^2+x4^2) - 9.8; 程序为 f = @(t, x)[x(2); -x(2)*sqrt(x(2).^2+x(4).^2); x(4); -x(4)*sqrt(x(2).^2+x(4).^2)-9.8];[t, y ]= ode45(f, [0...

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