MATLAB(符号运算)解析.pptVIP

或 或 或 y的一阶导数—— Dy y的二阶导数—— D2y y的 n 阶导数—— Dny [y1,y2…]=dsolve(x1,x2,…xn) —— 返回 微分方程的解 一阶微分方程 dsolve(Dx=y,Dy=x,x(0)=0,y(0)=1) ans = x(t) = sin(t), y(t) = cos(t) 二阶微分方程 dsolve(D2y=-a^2*y,y(0)=1,Dy(pi/a)=0) ans = cos(a*x) 例. y=dsolve(D2y+2*Dy+2*y=0,y(0)=1,Dy(0)=0) ans = exp(-x)*cos(x)+exp(-x)*sin(x) ezplot(y) —— 方程解y(t)的时间曲线图 求该方程的解 例：直线上的运动体受两个外力：一个正比于位移并与位移方向相反；另一个正比于速度并与速度相反，研究它的运动。 数学模型 令： 得： clear syms x w b t x=dsolve(D2x+2*b*Dx+w^2*x=0,Dx(0)=4,x(0)=0) dx=diff(x,t) w=2 b=0.1 x1=subs(x) dx1=subs(dx) t=0:0.1:100; xt=subs(x1,t); dxt=subs(dx1,t); subplot(22