MATLAB实验一《MATLAB在符号计算方面的应用》.doc

实验一 MATLAB在符号计算方面的应用 一、实验目的 1、掌握符号对象的创建和使用。 2、掌握符号数字及符号表达式的操作。 3、熟悉MATLAB在符号微积分、微分方程的符号解法、符号变换和符号卷积、线性代数方程的符号解等方面的应用。 了解Maple的深层符号计算资源。二、实验主要仪器与设备 装配有MATLAB.6软件的计算机三、预习要求 做实验前必须认真复习第二章MATLAB的符号计算功能。四、实验内容及实验步骤 指令窗中运行以下指令，写出运行结果，熟悉符号类数字与数值类数字之间的差异。本例演示：它们在创建方式、显示方式、类别判断和具体数值上的不同。 a=pi+sqrt(5) %创建一个数值类常数 sa=sym(pi+sqrt(5)) %创建一个符号类常数 Ca=class(a) %判断a的数据类别 Csa= class(sa) %判断sa的数据类别 vpa(sa-a) %在32位精度意义上计算两类数值之间的差 运行结果： a = 5.3777 sa = pi+sqrt(5) Ca = double Csa = sym ans = .138223758410852e-16 2、分两种情况用符号计算研究方程uz2+vz+w=0的解，写出输入的指令及运行结果。①不指定变量的情况；②把z指定为变量的情况；本例演示：MATLAB符号计算自动识别变量的能力；指定变量不同，所得解也不同；熟悉符号法解题的基本指令；熟悉在符号计算表达式中识别符号变量的专用指令findsym。syms u v w z %定义符号参数和变量u,v,w,z Eq= u*z^2+v*z+w ; %构成符号表达式 result_1=solve(Eq) % solve(Eq)意味着采用自动识别变量解方程u*z^2+v*z+w=0。 %MATLAB在符号计算中，x是首选符号变量，其后的次序排列规则是：与x的ASCII码值之差的绝对值小的字母优先；差绝对值相同时，ASCII码值大的字母优先。 findsym(Eq,1) %自动识别符号变量 运行结果： result_1 = -u*z^2-v*z ans = w ②result_2=solve(Eq,z) %对于指定变量z解方程 运行结果： result_2 = -1/2*(v-(v^2-4*u*w)^(1/2))/u -1/2*(v+(v^2-4*u*w)^(1/2))/u 3、输入指令简化，并写出运行结果。本例演示：simple的用法。syms x %创建一个符号对象x f= (1/x^3+6/x^2+12/x+8)^(1/3) ; %创建一个符号类表达式f g1=simple(f) %对表达式f进行简化 g2= simple(g1) %对表达式g1进行简化 运行结果： g1 = (2*x+1)/x g2 = 2+1/x 4、求。：clearsyms x %定义一个符号变量 f= 1/x*sqrt((1+x)/x) %创建一个符号表达式 s= int(f) %对f=f(x)以x为自变量进行积分 s1=simple(s) %对积分结果进行简化 s2= simple(s1) %再次对积分结果进行简化 运行结果： f = 1/x*((1+x)/x)^(1/2) s = ((1+x)/x)^(1/2)/x*(-2*(x^2+x)^(3/2)+2*(x^2+x)^(1/2)*x^2+log(1/2+x+(x^2+x)^(1/2))*x^2)/(x*(1+x))^(1/2) s1 = 1/x*(-2*(x*(1+x))^(1/2)+log(1/2+x+(x*(1+x))^(1/2))*x) 5、求的Fourier变换，在此x是参数，t是时间变量。本例演示：fourier的默认调用格式的使用要十分谨慎；在被变换函数中包含多个符号变量的情况下，对被变换的自变量给予指明，可保证计算结果的正确。syms t x w; %定义符号变量t,x,w ft= exp(-(t-x)) ; %创建符号表达式f=f(t) F1= simple(fourier(ft,t,w)) %求f(t)以t作为自变量的傅里叶变换， %傅里叶变换的结果以w作为自变量 F2=simple(fourier(ft))