资源描述:
《工具箱三十多个工具箱大致可分为两类功能型》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、工具箱三十多个工具箱大致可分为两类:功能型工具箱和领域型工具箱。功能型工具箱主要用来扩充MATLAB的符号计算功能、图形建模仿真功能、文字处理功能以及与硬件实时交互功能,能用于多种学科。而领域型工具箱是专业性很强的。如控制工具箱(ControlToolbox)、信号处理工具箱(SignalProcessingToolbox)等。下面,将MATLAB工具箱内所包含的主要内容做简要介绍:Symbolic符号运算(1)Symbolic是应用SymbolicMathToolbox,以符号(如a,b,c,x,y,z)为对象的一种数学运算。适用于数学推理、可变精度运算等。在Matlab界面下调
2、用‘sym’、‘syms’等函数。例symstM=(1-t^2)*exp(-1/2*t^2);P=(1-t^2)*sech(t);ezplot(M);holdon;ezplot(P)运行结果Symbolic符号运算(2)代数运算、三角运算微、积分(diff、int)函数整理与化简(combine、simplify、simple)可变精度计算(vpa、digits)解方程(代数solve、微分方程dsolve)线性代数与矩阵运算(determ、linsolve)变换(laplace、fourier、ztrans)Symbolic符号运算(3)——函数整理与简化将复杂的表达式进一步简化例:验证
3、积分symsAttaowyf=int(A*exp(-i*w*t),t,-tao/2,tao/2)Yf=simple(yf)symsAttaowYf=2*A*sin(1/2*tao*w)/wSymbolic符号运算(4)——可变精度计算可以设置任意多的有效计算位数进行计算digits(16)p0=sym('1+sqrt(5)/2')p1=sym((1+sqrt(5))/2)e01=vpa(abs(p0-p1))%32位缺省精度下观察p0,p1间误差p0=1+sqrt(5)/2p1=7286977268806824*2^(-52)e01=.5000000000000001例:以n位相对精度计算
4、Symbolic符号运算(5)——解方程solve&dsolve例:图示微分方程的通解和奇解的关系例解:y=dsolve('y=x*Dy-(Dy)^2','x')clf,holdon,ezplot(y(2),[-6,6,-4,8],1)%画奇解cc=get(gca,‘Children’);%取奇解曲线的图柄set(cc,‘Color’,‘r’,‘LineWidth’,5)%把奇解画成粗红线fork=-2:0.5:2;ezplot(subs(y(1),'C1',k),[-6,6,-4,8],1);end%画解holdoff,title(‘fontname{隶书}fontsize{16}通
5、解和奇解’)运行结果注意Symbolic(符号运算)中的dsolve只能给出解析解,而不能给出数值解。例:当运行y=dsolve(‘D6y+a5*D5y+a4*D4y+a3*D3y+a2*D2y+a1*Dy+a0*y=b4*D4u+b3*D3u+b2*D2u+b1*Du+b0*u’,‘t’)时,由于没有解析解,运行结果会出现:Warning:Explicitsolutioncouldnotbefound.>InD:MATLAB6p5toolboxsymbolicdsolve.matline326Ind:MATLAB6p5workexample2.matline6Symboli
6、c符号运算(6)——线性代数与矩阵运算Determ(即det)、linsolvesymsa11a12a21a22;A=[a11,a12;a21,a22]DA=determ(A)IA=inv(A)DA=a11*a22-a12*a21IA=[a22/(a11*a22-a12*a21),-a12/(a11*a22-a12*a21)][-a21/(a11*a22-a12*a21),a11/(a11*a22-a12*a21)]例:求矩阵的行列式值,逆。Symbolic符号运算(7)——变换laplace、fourier、ztranssymsts;symsabpositiveDt=sym('Dirac
7、(t-a)');Ut=sym('Heaviside(t-b)');Mt=[Dt,Ut;exp(-a*t)*sin(b*t),t^2*exp(-t)];MS=laplace(Mt,t,s)MS=[exp(-a*s),exp(-b*s)/s][b/((s+a)^2+b^2),2/(s+1)^3]例:求的Laplace变换参考文献张志涌等编著,«精通Matlab6.5版»(第五章符号计算),北京航空航天大学出版社,2003.陈