第三讲 矩阵、数组与符号运算2new.ppt

第3讲 矩阵、数组和符号运算2 第3讲 Matlab流程控制结构3 （1）降幂排列法?collect（又称合并同类项） 格式 R = collect(S) %对于多项式S中的每一函数，collect(S)按缺省变量x的次数合并系数。 R = collect(S,v)?? %对指定的变量v计算，操作同上。 已知表达式 ，将f按x 进行降幂排列，将g按 降幂排列 （2）展开法 格式? R = expand(S)?? %对符号表达式S中每个因式的乘积进行展开计算。该命令通常用于计算多项式函数、三角函数、指数函数与对数函数等表达式的展开式。 （3）重叠法 格式 ?R = horner(P)?? %若P为一符号多项式的矩阵，该命令将矩阵的每一元素转换成嵌套形式的表达式R。 （4）因式分解法（factor） 格式? factor(s)??? %符号表达式s的因式分解函数，说明： S为符号矩阵或符号表达式，常用于多项式的因式分解。 （5）一般化简法（simplify） 格式? R = simplify(S) 说明? ：充分考虑符号表达式的各种运算、特殊函数的运算性质，经计算机比较后，给出认为表达式相对简单的一种化简方法。 试编制程序，判断以下等式是否相等 （6）不定化简法（simple） simple (s)? ?????% s是矩阵或表达式，给出最简形式 ?[R,how]=simple (s)? %R为返回的最简形，how为简化过程中使用的主要方法。最笨拙，但综合前面5种化简方法的优点。用户可以根据所列出的结果进行比较和筛选。将表达式s的长度化到最短。若还想让表达式符合人的书写习惯，可使用函数Pretty。 （7）分子分母法（numden） ?[N,D] = numden(A) ? 将符号或数值矩阵A中的每一元素转换成整系数多项式的有理式形式，其中分子与分母是相对互素的。输出的参量N为分子的符号矩阵，输出的参量D为分母的符号矩阵。 C.符号变量（表达式）的提取与代入 （1）subexpr 格式 [Y,SIGMA] = subexpr(X,SIGMA) 说明? 找出符号表达式 X中相同的子表达式，将X中相同的、比较复杂的子字符串用符号SIGMA代替。 t = solve(a*x^3+b*x^2+c*x+d = 0); [r,s] = subexpr(t,s); pretty(t) pretty(r) 输入命令，并观察相同因子的提取 （2）符号变量的代入 格式subs(S,old,new) 用符号new 代替表达式S 中的符号old Syms x y subs(x+y,x,4) ans =4+y syms a b f=subs(cos(a)+sin(b),[a,b],[sym(alpha),2]) f =cos(alpha)+sin(2) 试将以下表达式中的变量x分别用y和m＋nt代替。 f=a*x^2+b*x+c 练习：求下列函数极限 练习：计算以下积分结果 练习：试求以下级数和 7、高级符号运算功能 1、积分变换 （1）傅立叶变换及其逆变换 A、对函数进行傅立叶变换的计算公式为： A、傅立叶变换fourier 格式：F=fourier(f):对自变量为x的表达式f(x)进行变换，返回值为F(ω)。 F=fourier(f,v):指定变换结果为变量v的函数。 F=fourier(f,u,v)：指定对函数表达式作关于变量u的傅立叶变换。 例：对以下函数求傅立叶变换 (1)拉普拉斯变换及其反变换 A、对函数f(x)进行拉普拉斯变换，计算公式为： 试求以下函数的拉普拉斯反变换，并把L3、L4表示为变量x的函数。 练习 线性系统冲击响应 编写求任意高阶连续常系数线性系统脉冲响应的程序。 解:任意阶次的连续线性系统可用下列微分方程表述： 写成传递函数形式为 其特性可用系统传递函数的分子分母系数向量b和a来表示。 传递函数反变换的求法 如果分母系数多项式没有重根，则可以把两个多项式之比分解成n个一阶部分分式之和。即： 其中p1,p2,…,pn是分母多项式的n个根，而r1,r2,…,rn是则对应于这n个根的留数。一阶分式的反变换可以查表得到，容易写出冲击响应的公式如下： 可见只要求出根向量p和留数向量r，线性方程的解就得到了。求根是代数问题，当阶次很高时，代数方程没有解析解。可喜的是MATLAB提供了用数值方法求根和留数的函数residue.m，它的调用格式为： [r,p]=residue(b,a) 课堂练习 给出系统的传递函数为 求系统的脉冲响应h(t)，如右图所示。