五讲符号计算.pptVIP

* f4=int(int(f,a),x) %双重积分 f2=int(f,x,0,2) f2 =8/3*a+2*b+2*c f3=int(f,a) %对a积分 f3=1/2*a^2*x^2+b*x*a+c*a f4 =1/6*a^2*x^3+1/2*b*x^2*a+c*a*x * 例 求椭球的体积。 命令如下： syms a b c z; f=pi*a*b*(c^2-z^2)/c^2; V=int(f,z,-c,c) V = 4/3*pi*a*b*c 例 轴的长度为10米，若该轴的线性密度计算公式是f(x)=6+0.3x千克/米(其中x为距轴的端点距离)，求轴的质量。 (1)符号函数积分。在MATLAB命令窗口，输入命令： sym x; f=6+0.3*x; m=int(f,0,10) (2)数值积分。 先建立一个函数文件fx.m： function fx=fx(x) fx=6+0.3*x; 再在MATLAB命令窗口，输入命令： m=quad(fx,0,10,1e-6) * 3. 符号积分变换 数学变换中通过积分运算将一类函数转变为另一种函数，使得计算分析更为简洁 三个主要积分变换： 傅里叶变换 拉普拉斯变换 Z 变换 * （1）. 傅立叶(Fourier)变换 在MATLAB中，进行傅立叶变换的函数是： fourier(fx,x,t) 求函数f(x)的傅立叶像函数F(t)。 ifourier(Fw,t,x) 求傅立叶像函数F(t)的原函数f(x)。 * 例 求函数的傅立叶变换及其逆变换。 命令如下： syms x t; y=abs(x); Ft=fourier(y,x,t) %求y的傅立叶变换 fx=ifourier(Ft,t,x) %求Ft的傅立叶逆变换 Ft =-2/t^2 fx = x*(2*heaviside(x)-1) %heaviside阶跃函数 * （2）. 拉普拉斯(Laplace)变换 在MATLAB中，进行拉普拉斯变换的函数是： laplace(fx,x,t) 求函数f(x)的拉普拉斯像函数F(t)。 ilaplace(Fw,t,x) 求拉普拉斯像函数F(t)的原函数f(x)。 * 例 计算y=x2的拉普拉斯变换及其逆变换. 命令如下： x=sym(x);y=x^2; Ft=laplace(y,x,t) %对函数y进行拉普拉斯变换 fx=ilaplace(Ft,t,x) %对函数Ft进行拉普拉斯逆变换 Ft = 2/t^3 fx =x^2 * （3）. Z变换 对数列f(n)进行z变换的MATLAB函数是： ztrans(fn,n,z) 求fn的Z变换像函数F(z) iztrans(Fz,z,n) 求Fz的z变换原函数f(n) * 例 求数列 fn=e-n的Z变换及其逆变换。 命令如下： syms n z fn=exp(-n); Fz=ztrans(fn,n,z) %求fn的Z变换 f=iztrans(Fz,z,n) %求Fz的逆Z变换 Fz =z/exp(-1)/(z/exp(-1)-1) f =exp(-1)^n * 三、符号方程和符号微分方程 * 1. 符号代数方程的符号求解 求解函数是solve，调用格式为： solve(eqn1,eqn2,…,eqnN,var1,var2,…,varN) 其中eqn1,eqn2,eqnN等为符号表达式或不带符号的字符串；var1,var2,…,varN 为自变量 * 例 解方程。 命令如下： x1=solve(1/(x+2)+4*x/(x^2-4)=1+2/(x-2),x) f=sym(x-(x^3-4*x-7)^(1/3)=1); x2=solve(f) x3=solve(2*sin(3*x-pi/4)=1) x4=solve(‘x+x*exp(x)-10’,‘x’) x1 =1 x2 =3 x3 =5/36*pi x4 =1.633506170155846384193165178978923428636 * 例 求方程组 的解。 命令如下： [x y]=solve(1/x^3+1/y^3=28,1/x+1/y=4,x,y) x = 1 1/3 y = 1/3 1 * 解微分方程的基本操作命令见下表 在求解微分方程中，用Dy表示 ，