[工学]第四章MATLAB符号计算.pptVIP

第四章 MATLAB符号计算 第四章   MATLAB符号运算 pretty 对特别复杂的表达式化简为符合习惯的形式 例： pretty(cos(x)^2-sin(x)^2) 例： syms x y z a b c f=(x+y)*(a+b+c)^z/(x+a)^2-1/(x-1); pretty(f) （1） 符号极限 limit(F,x,a) 计算符号表达式F在x→a条件下的极限； limit(F,a) 计算符号表达式F中默认自变量→a条件下的极限； limit(F) 计算符号表达式F在默认自变量→ 0时的极限； limit(F,x,a，‘right’)和limit(F,x,a,’left’) 计算符号表达式 F在x→a条件下的右极限和左极限。 例1：syms x a t h limit(sin(x)/x) %求sin(x)/x在x→0时的极限 例2：limit((1+2*t/x)^(3*x),x,inf) %求x→∞时的极限 例3：limit(1/x,x,0,right) %求1/x在x→0时的右极限 limit(1/x,x,0,left) 矩阵的符号积分 例： syms x y t; A=[cos(x*t),sin(x*t);-sin(x*t),cos(x*t)] int(A,t) %各元素对t的不定积分 int(y*log(1+y),0,1) %对默认自变量y从0~1的定积分 (3)符号求和 symsum(S) 求符号表达式S对于默认自变量的不定和。 symsum(S，v) 求符号表达式S对于自变量v的不定和。 symsum(S,a,b) 求符号表达式S对于默认自变量从a到b的有限和。 syms x t diff(sin(x^2)) ans = 2*cos(x^2)*x diff(t^6,6) ans = 720 clear [x,y]=meshgrid(-2:0.2:2,-2:0.2:2); z=x.*exp(-x.^2-y.^2); [px,py]=gradient(z,0.2,0.2); contour(z) %等高线 hold on quiver(px,py) %画矢量图 hold off 傅立叶逆变换调用格式为： ◆f=ifourier(F)：符号表达式对象F的傅立叶逆变换。缺省的 自变量为ω，缺省返回是关于 x 的函数。 ◆f = ifourier(F,u)：返回函数 f 是关于符号表达式对象 u 的 函数，而不是缺省的 x的函数。 ◆f=ifourier(F,v,u)：对关于 v 的函数 F 进行变换，返回关 于 u 的函数 f。 （7） Taylor 级数展开 在符号数学工具中，表达式的Taylor级数展开由函数taylor实现，其调用格式为： taylor(f) 计算符号表达式f在默认自变量等于0处的5阶 Taylor级数展开式。 taylor(f,n,v) 计算符号表达式f在自变量v=0处的n-1阶Taylor 级数展开式。 taylor(f,n,v,a) 计算符号表达式f在自变量v=a处的n-1阶 Taylor级数展开式。 例： 分别计算表达式 的5阶Taylor级数展开式和 的5阶及12阶Taylor级数展开式。 在命令窗口创建符号变量x，分别计算上面各表达式： syms x f1=1/(5+cos(x)); r1=taylor(f1) f2=exp(x*sin(x)); r2=taylor(f2) r3=taylor(f2,12) 求非线性方程组 的解。 syms x y z f=x^2-y^2+z-10; g=x+y-5*z; h=2*x-4*y+z; [x,y,z]=solve(f,g,h) S=solve(f,g,h); %S为结构变量 [S.x,S.y,S.z