(MATLAB学习)第四节 MATLAB符号计算.ppt

第四节 MATLAB符号计算 y=sym(4); z=y/x1 z = 4/3 double(z) ans = 1.3333 4.分式通分： numden 例：求极限 syms a b c x f=a*x^2+b*x+c; diff(f,x) ans = 2*a*x+b diff(f,x,2) ans = 2*a diff(f,a) ans = x^2 syms x y f=x*y+y; int(f,x) ans = 1/2*x^2*y+x*y int(f,y) ans = 1/2*x*y^2+1/2*y^2 syms x a b f=x^2+1; int(f,x,a,b) ans = 1/3*b^3-1/3*a^3+b-a 多重积分: syms x y f=x^y; int(int(f,x,0,1),y,1,2) ans = log(3)-log(2) double(ans) ans = 0.40546510810816 symsum 函数 符号求和由函数symsum来实现。 该函数的一般调用格式为： symsum(s)： 按照默认自变量，对s离散求和 symsum(s,v)： 以v为自变量，对s离散求和 symsum(s,a,b)： 按照默认自变量， a,b分别表示下限和上限 symsum(s,v,a,b)：以v为默认自变量， a,b分别表示下限和上限 例： syms k s=symsum(1/k^2,1,inf) s = 1/6*pi^2 syms k x s=symsum(x^k,k,0,inf) s = -1/(x-1) 作业： syms a b x X Y k=sym(’3’); z=sym(c*sqrt(delta)+y*sin(theta)); f=a*z*X+(b*x^2+k)*Y; syms a b t u v x y; A=[a+b*x,sin(t)+u; x*exp(-t),log(y)+v]; syms a x; f=a*sin(x)+5; f1=subs(f,sin(x),sym(y)) f2=subs(f,{a,x},{2,sym(pi/3)}) f3=subs(f,{a,x},{2,pi/3}) f4=subs(subs(f,a,2),x,0:pi/6:pi) f5=subs(f,{a,x},{0:6,0:pi/6:pi}) 1.化简 2. 试写出下面命令的结果。 findsym(f) findsym(f,1) findsym(f,2) findsym(f,3) 3.设A是一个符号矩阵（定义如下），试指出findsym(A,1) 的输出结果，并由此得出结论： findsym 确定自由变量时，是对整个矩阵进行的，还是对各个矩阵元素分别进行？ 4.(1)设 求复合函数 f(g(z)); (2)设 求复合函数 f(g(z))。 5.替换函数subs的使用。 (1)试指出下面程序中的f1、f2、f3、f4、f5的值。 syms x y; f=2*x+y; subs(f,{x,y},{x+y,x-y}) (2)试指出下面替换的结果，并说明理由。 4.4 常见积分变换函数 常见的积分变换有傅立叶变换、拉普拉斯变换和Z变换。 1．傅立叶(Fourier)变换 在MATLAB中，进行傅立叶变换的函数是： fourier(f):默认得到关于w的函数 fourier(f,v):指定函数为v fourier(f,u,v)：f有多个变量 2．傅立叶(Fourier)反变换 f = ifourier(F) f = ifourier(F,u) f = ifourier(F,v,u) 例题： 2．拉普拉斯(Laplace)变换 在MATLAB中，进行拉普拉斯变换的函数是： laplace(fx,x,y)： 求函数f(x)的拉普拉斯像函数F(y)， 默认为F(s)。 ilaplace(Fw,s,y)：求拉普拉斯像函数F(s)的原函数f(y)， 默认为f(x) 例9-5 计算y=x3a的拉普拉斯变换及 逆变换。 * * 4.1 数值运算与符号运算 4.2 符号变量和符号表达式 4.3 符号表示式的运算 4.4 微积分 4.5 方程求解 学习重点： 1、理解符号运算的有关概念 2、掌握使用符号运算解决符号推导、微积分