unit3 matlab.ppt

3.4.3 符号积分（Integration） 积分有定积分和不定积分，运用函数int可以求得符号表达式的积分，即找出一个符号表达式F使得diff(F)=f，也可以说是求微分的逆运算。 int(f,’t’) %求符号变量t的不定积分 int(f,’t’,a,b) %求符号变量t的定积分 int(f,’t’,’m’,’n’) %求符号变量t的定积分 例：f=sym(cos(x)); int(f) int(f,0,pi/3) int(f,a,b) diff和int都可以对符号矩阵进行操作，例如diff(A),若A是一个 符号矩阵，则对A的每个元素（是一个符号表达式）进行微分运算 3.4.4 符号级数 1. symsum函数（Symbolic Summation） symsum(s,x,a,b) %计算表达式s的级数和， 前提是级数存在。 说明：x为自变量，x省略则默认为对自由变量 求和；s为符号表达式；[a,b]为参数x的取值范 围。 例：syms x k s1=symsum(1/k^2,1,10) %求级数的前10项和 s2=symsum (1/k^2,1,inf) %求级数和 s3=symsum(x^k,k,0,inf) %求k为自变量的级数和 2. taylor函数（Taylor Series） taylor(F,x,n) %求泰勒级数展开 说明：x为自变量，F为符号表达式；对F进行泰勒级数展开至n项，参数n省略则默认展开前5项。 例：syms x s1=taylor(exp(x),8) 3.5 符号积分变换 3.5.1傅里叶变换及其反变换 1. Fourier变换 F＝fourier(f,t,w)%求时域函数f(t)的fourier变换F 说明：返回结果F是符号变量w的函数，f为t的函数。若省略t则为自由变量，省略w则默认F是w的函数。 2. Fourier反变换 f=ifourier (F,w,t) 说明：ifourier函数的用法与fourier函数相同。 syms t w f=sym(’heaviside(t)) %Heaviside(t)为单位阶跃函数u(t) f = heaviside(t) F=fourier(f) F= pi*dirac(w)-i/w %dirac(w)为单位脉冲函数 ifourier(F) ans = heaviside(x) 注意比较fourier(f), fourier(f,’u’), fourier(f,’x’), fourier(f,’x’,’u’) 3.5.2拉普拉斯变换及其反变换 1. Laplace变换 F=laplace(f,t,s) %求时域函数f的Laplace变换F 说明：返回结果F为s的函数，当参数s省略，返回结果F默认为s的函数；f为t的函数，当参数t省略，默认自由变量为t。 2. Laplace反变换 f＝ilaplace(F,s,t) %求F的Laplace反变换f syms a t s F1=laplace(sin(a*t),t,s) %sinat的Laplace变换 F1 =a/(s^2+a^2) F=ilaplace(sym(1),s,t) %1的Laplace反变换 3.5.3 Z变换及其反变换 1. ztrans函数 F＝ztrans(f,n,z) %求时域序列f的Z变换F 说明：返回结果F是以符号变量z为自变量；当参数n省略，默认自变量为n；当参数z省略，返回结果默认为z的函数。 2. iztrans函数 f＝iztrans(F,z,n) %求F的z反变换f syms a n z t Fz3=ztrans(exp(-a*t),n,z) %e-at的Z变换 Fz3 =exp(-a*t)*z/(z-1) 比较ztrans(a^t,n,z)和ztrans(a^n,n,z) 3.6 符号方程的求解3.6.1 代数方程 MATLAB可以用solve命令给出方程的数值解。 solve(‘eq’,’v’)%求方程关于指定变量的解 solve(‘eq1’, ’eq2’,’v1’,’v2’,…) %求方程组关于指定变量的解 例如，解方程： solve(a*x^2+b*x+c) solve(a*x^2+b*x+c=0) solve(a*x^2+b*x+c=0,x) 上面三个命令的结果一样 【例3.22】求三元非线性方程组的解。 eq1=sym(x^2+2*x+1); eq2=sym(x+3*z=