第七讲matlab符号运算2.pptVIP

第七讲 符号运算 1、符号计算基础 2、微积分 3、简化方程表达式 4、符号积分变换 5、解方程 6、图示化符号函数计算器 其他函数： 符号表达式的因式分解与展开的函数，调用格式为： factor(s)：对符号表达式s分解因式。 expand(s)：对符号表达式s进行展开。 collect(s)：对符号表达式s合并同类项。 collect(s,v)：对符号表达式s按变量v合并同类项。 符号表达式的化简 simplify(s)：应用函数规则对s进行化简。 simple(s)：调用MATLAB的其他函数对表达式进行综合化简，并显示化简过程。 * gjy112@ * 江苏科技大学数理学院 将符号多项式s用嵌套形式表示，即用多层括号的形式表示。Horner函数可以实现此功能。该函数的格式为： horner(s) 【例10】将表达式x^3-6*x^2+11*x-6用嵌套形式表示。 syms x horner(x^3-6*x^2+11*x-6) ans = -6+(11+(-6+x)*x)*x 1、Fourier变换及其反变换 例如：求 的Fourier变换。本例演示三个重要内容：单位阶跃函数和 单位脉冲函数的符号表示；fourier指令的使用；simple指 令在MATLAB不同版本中的表现差异。 （1）求Fourier变换 syms t w;ut=sym(Heaviside(t)); %定义0时刻起跳的单位阶跃函数 UT=fourier(ut) %实施Fourier变换，给出与理论一致的结果 UTC=maple(convert,UT,piecewise,w) %计算结果起指示作用 UTS=simple(UT) %在此是5.3版的运算结果，简化导致漏项！ UT =pi*Dirac(w)-i/w UTC = PIECEWISE([undefined, w = 0],[0, otherwise]) UTS =-i/w HEAVISIDE(X) is 0 for X 0, 1 for X 0, and NaN for X == 0. diff(heaviside(x),x) = dirac(x) （2）求Fourier反变换进行验算 Ut=ifourier(UT,w,t) %结果与原函数相等 Uts=ifourier(UTS,w,t) %结果与原函数不等 Ut = 1/2+1/2*Heaviside(t)-1/2*Heaviside(-t) Uts = 1/2*Heaviside(t)-1/2*Heaviside(-t) 求 的Laplace变换。 syms t s;syms a b positive %对常数进行“限定性”设置 Dt=sym(Dirac(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/(1+s)^3] 2、Laplace变换及其反变换 3、Z变换及其反变换 求序列 的Z变换，并用反变换验算 （1）单位函数及性质验证 syms n Delta=sym(charfcn[0](n));%定义单位函数 D0=subs(Delta,n,0); %计算 D15=subs(Delta,n,15);%计算 disp([D0,D15]);disp([D0,D15]) [D0,D15] [ 1, 0] （2）求序列 syms z;fn=2*Delta+6*(1-(1/2)^n); FZ=simple(ztrans(fn,n,z)); disp(FZ = ); pretty(FZ),FZ_n=iztrans(FZ,z,n) FZ = 2 4 z + 2 -------------- 2 2 z - 3 z + 1 FZ_n = 2*charfcn[0](n)+6-6*(1/2)^n 的Z变换 解方程函数的格式为： solve(expr1,expr2,...,exprN,var1,var2,...varN