第三章-matlab语言的符号运算.ppt

例题：计算拉普拉斯反变换例题：计算Z变换计算反Z变换。4.符号代数方程求解matlab符号运算能够解一般的线性方程、非线性方程及一般的代数方程、代数方程组。当方程组不存在符号解时，又无其他自由参数，则给出数值解。命令格式：solve(f)——求一个方程的解solve(f1,f2,…fn)——求n个方程的解例1.f=ax2+bx+c求解f=a*x^2+b*x+c;solve(f)——对缺省变量x求解ans=[1/2/a*(-b+(b^2-4*a*c)^(1/2))][1/2/a*(-b-(b^2-4*a*c)^(1/2))]计算机格式一般格式例2.符号方程cos(x)=sin(x)tan(2*x)=sin(x)求解 f1=solve(cos(x)=sin(x)), f1= 1/4*pisolve(f,b)——对指定变量b求解ans=-(a*x^2+c)/xf3= matlab6.2的解 [0][pi][atan(1/2*(-2*3^(1/2))^(1/2),1/2+1/2*3^(1/2))][atan(-1/2*(-2*3^(1/2))^(1/2),1/2+1/2*3^(1/2))][atan(1/2*2^(1/2)*3^(1/4)/(1/2-1/2*3^(1/2)))+pi][-atan(1/2*2^(1/2)*3^(1/4)/(1/2-1/2*3^(1/2)))-pi]f2=solve(tan(2*x)=sin(x))f2= matlab4.2的解[0][acos(1/2+1/2*3^(1/2))][acos(1/2-1/2*3^(1/2))]例.解方程组x+y+z=1x-y+z=22x-y-z=1g1=x+y+z=1,g2=x-y+z=2,g3=2*x-y-z=1f=solve(g1,g2,g3)f=solve(x+y+z=1,x-y+z=2,2*x-y-z=1)f=z=5/6,y=-1/2,x=2/3f=solve(x+y+z=1,x-y+z=2,2*x-y-z=1)f=x:[1x1sym] f.x ans=2/3 y:[1x1sym] f.y ans=-1/2z:[1x1sym] f.z ans=5/6[x,y,z]=solve(x+y+z=1,x-y+z=2,2*x-y-z=1)x=2/3y=-1/2z=5/65.符号微分方程求解——用一个函数可以方便地得到微分方程的符号解符号微分方程求解指令：dsolve命令格式：dsolve(f,g)f——微分方程，可多至12个微分方程的求解；g为初始条件默认自变量为x,可任意指定自变量t,u等微分方程的各阶导数项以大写字母D表示或或或y的一阶导数——Dyy的二阶导数——D2yy的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)的时间曲线图求该方程的解 对于教材中提到的状态方程求解、差分方程求解和约当标准型我们就不作要求了