研究生院matlab作业5.docVIP

MATLAB作业5 试求出下面线性微分方程的通解。 假设上述微分方程满足已知条件，试求出满足该条件的微分方程的解析解。 syms t y ; y=dsolve([D5y+13*D4y+64*D3y+152*D2y+176*Dy+80*y=,exp(-2*t)*(sin(2*t+pi/3)+cos(3*t))],y(0)=1,y(1)=3,y(pi)=2,Dy(0)=1,Dy(1)=2); vpa(y,20) ans = .20576131687242798354e-2*exp(-2.*t)*cos(3.*t)+.15538705805619602373e-1*exp(-2.*t)*sin(2.*t)+.76830587084294035587e-2*exp(-2.*t)*cos(2.*t)+98.159206062620455336*exp(-2.*t)*t+59.405044899367325899*exp(-2.*t)*t^3-106.24422608844727795*exp(-2.*t)*t^2-30.741892776456442810*exp(-2.*t)+.20576131687242798354e-2*exp(-2.*t)*sin(3.*t)+31.732152104579289128*exp(-5.*t) 试求解下面微分方程的通解以及满足条件下的解析解。 syms t x y; [x,y]=dsolve(D2x+5*Dx+4*x+3*y=exp(-6*t)*sin(4*t), 2*Dy+y+4*Dx+6*x=exp(-6*t)*cos(4*t),x(0)=1,x(pi)=2,y(0)=0);vpa(x,20) ans = .24816266536011758942e-1*exp(-6.*t)*cos(4.*t)-.16682998530132288094e-1*exp(-6.*t)*sin(4.*t)+.85861668591377882749e-1*exp(t)-.57658489325155296910e-1*exp(-5.1374586088176874243*t)+.94698055419776565523*exp(-1.3625413911823125757*t) vpa(y,20) ans = -.10607545320921117099*exp(-6.*t)*cos(4.*t)+.68348848603625673689e-1*exp(-6.*t)*sin(4.*t)-.28620556197125960916*exp(t)+.90450561852315609427e-1*exp(-5.1374586088176874243*t)+.30183045332815517077*exp(-1.3625413911823125757*t) 试求出微分方程的解析解通解，并求出满足边界条件的解析解。 syms x y; y=dsolve([D2y-(2-1/x)*Dy+(1-1/x)*y=, x^2*exp(-5*x)]) y = exp(t)*C2+exp((x-1)/x*t)*C1+exp(-5*x)/(x-1)*x^3 y=dsolve([D2y-(2-1/x)*Dy+(1-1/x)*y=, x^2*exp(-5*x)],y(1)=pi,y(pi)=1) y = -exp(t)/(-exp(1)*exp((x-1)/x*pi)+exp(1)*x*exp((x-1)/x*pi)-exp((x-1)/x)*exp(pi)*x+exp((x-1)/x)*exp(pi))*(-pi*x*exp((x-1)/x*pi)+pi*exp((x-1)/x*pi)+x^3*exp(-5*x)*exp((x-1)/x*pi)+exp((x-1)/x)*x-exp((x-1)/x)-exp((x-1)/x)*x^3*exp(-5*x))-exp((x-1)/x*t)*(-exp(pi)*x^3*exp(-5*x)+x^3*exp(-5*x)*exp(1)+exp(pi)*pi*x-exp(1)*x-exp(pi)*pi+exp(1))/(x-1)/(exp(1)*exp((x-1)/x*pi)-exp((x-1)/x)*exp(pi))+exp(-5*x)/(x-1)*x^3 Lotka-Volterra扑食模型方程为，且初值为，试求解该微分方程，并绘制相应的曲线。 syms x y t; f=inline([4*x(1)-2*x(1)*x(2); x(1)*x(2)-3*x(2)