求用matlab 做数学实验报告! 1求微分方程特解. d^2y/dx^2+4dy/dx+29y=0 y(0)=0,y'(0)=152求微分方程的通解y''-2y'+5y=sin2x3解微分方程组y1'=y2y3y2'=-y1y3y3'=-0.51y1y2y1(0)=0,y2(0)=1,y3(0)=1
来源:学生作业帮助网 编辑:作业帮 时间:2024/07/01 02:23:02
![求用matlab 做数学实验报告! 1求微分方程特解. d^2y/dx^2+4dy/dx+29y=0 y(0)=0,y'(0)=152求微分方程的通解y''-2y'+5y=sin2x3解微分方程组y1'=y2y3y2'=-y1y3y3'=-0.51y1y2y1(0)=0,y2(0)=1,y3(0)=1](/uploads/image/z/6864952-40-2.jpg?t=%E6%B1%82%E7%94%A8matlab+%E5%81%9A%E6%95%B0%E5%AD%A6%E5%AE%9E%E9%AA%8C%E6%8A%A5%E5%91%8A%21+1%E6%B1%82%E5%BE%AE%E5%88%86%E6%96%B9%E7%A8%8B%E7%89%B9%E8%A7%A3.+d%5E2y%2Fdx%5E2%2B4dy%2Fdx%2B29y%3D0+y%280%29%3D0%2Cy%27%280%29%3D152%E6%B1%82%E5%BE%AE%E5%88%86%E6%96%B9%E7%A8%8B%E7%9A%84%E9%80%9A%E8%A7%A3y%27%27-2y%27%2B5y%3Dsin2x3%E8%A7%A3%E5%BE%AE%E5%88%86%E6%96%B9%E7%A8%8B%E7%BB%84y1%27%3Dy2y3y2%27%3D-y1y3y3%27%3D-0.51y1y2y1%280%29%3D0%2Cy2%280%29%3D1%2Cy3%280%29%3D1)
求用matlab 做数学实验报告! 1求微分方程特解. d^2y/dx^2+4dy/dx+29y=0 y(0)=0,y'(0)=152求微分方程的通解y''-2y'+5y=sin2x3解微分方程组y1'=y2y3y2'=-y1y3y3'=-0.51y1y2y1(0)=0,y2(0)=1,y3(0)=1
求用matlab 做数学实验报告! 1求微分方程特解. d^2y/dx^2+4dy/dx+29y=0 y(0)=0,y'(0)=15
2求微分方程的通解
y''-2y'+5y=sin2x
3解微分方程组
y1'=y2y3
y2'=-y1y3
y3'=-0.51y1y2
y1(0)=0,y2(0)=1,y3(0)=1
求用matlab 做数学实验报告! 1求微分方程特解. d^2y/dx^2+4dy/dx+29y=0 y(0)=0,y'(0)=152求微分方程的通解y''-2y'+5y=sin2x3解微分方程组y1'=y2y3y2'=-y1y3y3'=-0.51y1y2y1(0)=0,y2(0)=1,y3(0)=1
对于第一个问题可以直接采用函数求解的方法
y1=dsolve('D2y+4*Dy+29*y=0','Dy(0)=15',' y(0)=0','x')
y1=(3*sin(5*x))/exp(2*x)
y2=dsolve('D2y-2*Dy+5*y=sin(2*x)')
y2=sin(2*x)/5 + C5*cos(2*t)*exp(t) + C6*sin(2*t)*exp(t)
对于第三个问题,那么不能求出通解,所以只能借助于数值求解的方法来求解,数值求解采用ode45函数来求解的方法,具体过程如下:
先编写待求解的微分方程函数,打开编辑器:
%编写要求解的微分方程组函数表达式!
function dy = rigid(t,y)
dy = zeros(3,1); % 定义数组函数!
dy(1) = y(2) * y(3);%第一个微分方程;
dy(2) = -y(1) * y(3);%第二个微分方程;
dy(3) = -0.51 * y(1) * y(2);%第三个微分方程;
并以默认的文件名保存函数文件!
编写命令行,对微分方程求
options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4 1e-5]);%定义求解选项包含精度项!
[T,Y] = ode45(@rigid,[0 12],[0 1 1],options);%采用ode45求解方程组,并把求解结果保存到数组T,Y中!
plot(T,Y(:,1),'r-',T,Y(:,2),'b-.',T,Y(:,3),'k.')%作图!
grid minor%网格化!
运行上面的代码就得到函数的
附上图片!
如果有什么问题可以问我!