第4讲　matlab符号运算.ppt

例4-7 1. 求出 的分子、分母 2. 求出符号矩阵各元素的分子、分母 eg4-7 (4)．符号表达式的替换 [Y,SIGMA]=subexpr(S,SIGMA)：用SIGMA的值代替符号表达式S中重复出现的字符串 subs (S)：用MATLAB工作空间中的变量替换S符号表达式中的所有变量 subs (S,NEW)：用变量NEW替换符号表达式S中的自由变量 subs (S,OLD,NEW)：用变量NEW替换符号表达式S中的变量OLD 例4-8： syms a b x; f=a*sin(x)+b; f1=subs(f,a,sym(2)) %以符号常量2替换符号变量a f2=subs(f,{a,b},{sym(2),sym(3)}) %分别以符号常量2、3替换符号变量a、b f3=subs(f,{a,b},{3,2}) %分别以标量3、2替换符号变量a、b f4=subs(f,sin(x),sym(y)) %以符号变量y替换符号表达式sin(x) eg4-8 4. 符号微积分与积分变换 符号极限函数limit的调用格式为： limit(f,x,a)：计算当变量x趋近于a时，f(x)函数的极限值 (2) limit(f,a)：求f(x)的极限值，此时自变量为默认自变量 (3) limit(f)：求f(x)的极限值，此时默认变量趋近于0 (4) limit(f,x,a,‘right’)：求f的极限值，right表示变量x从右边趋近于a (5) limit(f,x,a,‘left’)：求f的极限值，‘left’表示变量x从左边趋近于a 例4-9：求下列极限 eg4-9 diff(f) — 对缺省变量求微分 diff(f,v) — 对指定变量v求微分 diff(f,v,n) —对指定变量v求n阶微分 int(f) — 对f表达式的缺省变量求积分 int(f,v) — 对f表达式的v变量求积分 int(f,v,a,b) — 对f表达式的v变量在（a,b）区 间求定积分 5.符号微分和积分 不定 积分 例4-10 计算二重不定积分 syms x y F=int(int(x*exp(-x*y),x),y) F = 1/y*exp(-x*y) 级数求和运算函数： symsum(s,x,a,b) 功能：计算符号表达式s的级数和 其中s为通项式，x自变量，x省略则默认为对自由变量求和；[a,b]为参数x的取值范围 级数求和 eg4-11 taylor(f,n) —— 泰勒级数展开 ztrans(f) —— Z变换 iztrans(f) —— 反Z变换 laplace(f) —— 拉氏变换 ilaplace(f) —— 反拉氏变换 fourier(f) —— 傅里叶变换 ifourier(f) —— 反傅里叶变换 积分变换 例4-12 计算指数函数eAt 利用拉氏反变换法计算eAt的公式为： eAt = L-1[(SI-A)-1] 系统矩阵A= eAt = 结果： a=[0 1;-2 -3]; syms s B=(s*eye(2)-a) B = [s, -1] [2, s+3] b=inv(B) b=[ (s+3)/(s^2+3*s+2), 1/(s^2+3*s+2)] [ -2/(s^2+3*s+2), s/(s^2+3*s+2)] b11=ilaplace(sym(b,1,1));b(1,1)=b11; b12=ilaplace(sym(b,1,2));b(1,2)=b12; b21=ilaplace(sym(b,2,1));b(2,1)=b21; b22=ilaplace(sym(b,2,2));b(2,2)=b22; b b = [ -exp(-2*t)+2*exp(-t), exp(-t)-exp(-2*t)] [ -2*exp(-t)+2*exp(-2*t), 2*exp(-2*t)-exp(-t)] 6.符号代数方程求解 matlab符号运算能够解一般的线性方程、非线性方程及一般的代数方程、代数方程组。当方程组不存在符号解时，又无其他自由参数，则给出数值解。 命令格式： solve(f) —— 求一个方程的解 solve(f1,f2, …fn) —— 求n个方程的解 例4-13. 求解方程ax2+bx+c＝0 f=a*x^2+b*x+c;