第3篇__MATLAB符号计算.ppt

3.1 符号计算基础 3.1 符号计算基础 3.1.1 符号对象 例3.1考察符号变量和数值变量的差别。 例3.2比较符号常数与数值在代数运算时的差别。 (2)syms函数 2. 建立符号表达式 例3.4建立x,y的一般二元函数。 3、默认自变量情况 3.1.2 基本的符号运算 (2)因式分解与展开 例3.7 对符号矩阵A的每个元素分解因式。 例3.8 计算表达式S的值。 (3)表达式化简 2. 符号矩阵运算 3.1.3 符号表达式中变量的确定 3.2 符号导数及其应用 3.2.1函数的极限 例3.10 求极限。 3.2.2 符号函数求导及其应用 例3.11 求函数的导数。 例3.12在曲线y=x3+3x-2上哪一点的切线与直线y=4x-1平行。 3.3 符号积分 3.3.1不定积分 例3.13 求不定积分 3.3.2 符号函数的定积分 例3.15 求椭球的体积 例3.16 轴的长度为10米，若该轴的线性密度计算公式是f(x)=3+0.3x千克/米(其中x为距轴的端点距离)，求轴的质量。 3.3.3积分变换 2. 拉普拉斯(Laplace)变换 3. Z变换 3.4 级数 3.4.1 级数的符号求和 3.4.2 函数的泰勒级数 例3.22将多项式表示成x+1的幂的多项式。 3.5 代数方程符号求解 3.5.1线性方程组的符号求解 3.5.2 非线性方程组的符号求解 3.6常微分方程的符号求解 1) 函数dsolve; 2)调用格式：dsolve(eqn1,condition,var) 方法：在函数dsolve所包含的equation中，用字母D表示求微分，D的数字表示几重微分，D后的变量为因变量。如Dy表示一阶微分项y’;D2y表示二阶微分项y’’;并且默认所有这些变量都是对自变量t求导。 dsolve(’D2y=1+Dy’) %求y’’=1+y’的通解 ans= -t+C1+C2*exp(t) dsolve(’D2y=1+Dy’,’y(0)=1’,’Dy(0)=0’) %求y’’=1+y’的解,加初始条件 ans= -t+exp(t) 例3.27 求微分方程的通解。 * 常微分方程组求解 1、求极限 和 1)ztrans(fn,n,z) 求fn的Z变换像函数F(z) 2)iztrans(Fz,z,n) 求Fz的z变换原函数f(n) 例3.19求数列 fn=e-n的Z变换及其逆变换。 syms n z fn=exp(-n); Fz=ztrans(fn,n,z) %求fn的Z变换 f=iztrans(Fz,z,n) %求Fz的逆Z变换 symsum (s,v,a,b) : %自变量v在[a,b]之间取值 时，对通项s求和。 例3.20 分别求级数1+1/2+1/3+…+1/k+…和级数1/2+1/(2×3)+1/(3×4)+…+1/(k(k+1))+…之和。 syms k symsum(1/k,k,1,inf) % +1/2+1/3+…+1/k+… symsum(1/(k*(k+1)),k,1,inf) % 1/2+1/(2×3)+1/(3×4)+…+1/(k(k+1))+… taylor(F,v,n): %求F对自变量v的taylor级数展 开至n级。 例3.21求函数在指定点的泰勒展开式。 x=sym(x); f1=(1+x+x^2)/(1-x+x^2); f2=sqrt(1-2*x+x^3)-(1-3*x+x^2)^(1/3); taylor(f1,x,5) %求展开到x的4次幂时应择n=5 taylor(f2,3) x=sym(x); p=1+3*x+5*x^2-2*x^3; f=taylor(p,x,-1,4) 例3.23求sin(x)的幂次方小于10次方的展开式。 syms x taylor(sin(x),10) %求sin(x)taylor级数展开 Ans=? 1)函数linsolve; 2)调用格式为：linsolve(A,b) 例3.24 求线性方程组AX=b的解。 解方程组(1)的命令如下： A=[34,8,4;3,34,3;3,3,8]; b=[4;3;2]; X=linsolve(A,b) %调用linsolve函数求(1)的解 A\b