第3讲矩阵、数组和符号运算2new.pptVIP

;第3讲矩阵、数组和符号运算;第3讲矩阵、数组和符号运算;第3讲矩阵、数组和符号运算;第3讲矩阵、数组和符号运算;2、符号表达式和矩阵的创建 a.字符串直接输入创建 符号表达式和符号方程对空格很敏感。因此，在创建符号表达式或符号方程时，不要在字符间任意加空格符； 符号计算中出现的数字也是当作符号处理的； f=a*x^2+b*x+c ?f = ?a*x^2+b*x+c f=a*x^2+b*x+c=0 f = ?a*x^2+b*x+c=0 ;第3讲矩阵、数组和符号运算;第3讲矩阵、数组和符号运算;第3讲矩阵、数组和符号运算;基本运算 四则运算 两个符号矩阵的大小相等方可进行加减运算，符号矩阵和符号标量的加减运算按照数组运算规则进行； 两个符号矩阵只有内积相等时才可以进行乘法运算； 符号的乘方运算 S^p，若 S 为符号表达式，p 可以为符号表达式或数值表达式；若 S 为符号矩阵，则 p 必须是整数。;第3讲矩阵、数组和符号运算;b符号简化 Matlab提供了多种化简符号表达式的方法： 降幂排列法（collect） 展开法（expand） 重叠法（horner） 因式分解法（factor） 一般化简法（simplify） 不定化简法（simple） 分子分母形式法（numden）;（1）降幂排列法?collect（又称合并同类项） 格式 R = collect(S) %对于多项式S中的每一函数，collect(S)按缺省变量x的次数合并系数。 R = collect(S,v)?? %对指定的变量v计算，操作同上。;;;（3）重叠法 格式 ?R = horner(P)?? %若P为一符号多项式的矩阵，该命令将矩阵的每一元素转换成嵌套形式的表达式R。;;;;（6）不定化简法（simple） simple (s)? ?????% s是矩阵或表达式，给出最简形式 ?[R,how]=simple (s)? %R为返回的最简形，how为简化过程中使用的主要方法。最笨拙，但综合前面5种化简方法的优点。用户可以根据所列出的结果进行比较和筛选。将表达式s的长度化到最短。若还想让表达式符合人的书写习惯，可使用函数Pretty。 ;第3讲矩阵、数组和符号运算;;;;; ;a. 符号极限 limit(F,x,a)：求当x--a时表达式F的极限值 limit(F,a)，变量为由findsym定义的默认变量 limit(F,x,a,’right’) /limit(F,x,a,’left’)（右/左极限值） Limit(F), a=0：求系统默认变量?0时表达式F的极限值 syms x a t h limit(sin(x)/x) ans = 1 limit((1+2*t/x)^(3*x),x,inf) ans = exp(6*t);;b.符号积分 int(S) 不定积分 int(S,v) int(S,a,b) 定积分 int(S,v,a,b);; ?r = symsum(s)?? ??%对符号表达式s中的符号变量k（由命令findsym(s)确定的）从0到k-1求和 r = symsum(s,v)? ??%对符号表达式s中指定的符号变量v从0到v-1求和 r = symsum(s,a,b)?? %对符号表达式s中的符号变量k（由命令findsym(s)确定的）从a到b求和 r = symsum(s,v,a,b)?? %对符号表达式s中指定的符号变量v从a到b求和;;1、积分变换 （1）傅立叶变换及其逆变换 B、对函数进行傅立叶逆变换的计算公式为： ;;;b、对函数L(s)进行拉普拉斯反变换，计算公式为：;A、 Laplace变换 L=laplace(F) x s，F为待进行拉普拉斯变换的代数表达式，默认变量为t。函数返回默认自变量为s。 L=laplace(F,z) x z ，指定返回结果为自变量z的函数。 L=laplace(F,w,z) w z，要求F对自变量w进行拉普拉斯变换，返回值的自变量为z。 ;试对函数 进行拉普拉斯变换，并将所得结果表示为变量v的函数;已知某控制系统的微分方程 其中T＝0.5s，K＝10.设初始条件为0，试求系统单位阶跃响应曲线;B、 Laplace反变换 F=ilaplace(L) L为待进行拉普拉斯反变换的代数表达式，默认变量为s。函数返回默认自变量为t。如果L=L(t),则返回结果为F(x). F=ilaplace(L,y) 指定返回F为自变