第4章高等数学计算VIP

1、符号表达式及其应用符号表达式及其应用微积分基本运算微积分基本运算taylor级数展开级数展开常微分方程求解常微分方程求解第第4章章 高等数学计算高等数学计算matlab的的符号计算符号计算符号计算相对于数值计算而存在符号计算相对于数值计算而存在,数值计算的特点是数值计算的特点是“近似近似”, 符号计算的特点是符号计算的特点是“准确准确”. 早期的早期的matlab不能实现符号计算不能实现符号计算,1993年年mathworks公司购买软件公司购买软件maple的使用权的使用权,并开发实现符号计算的工具箱并开发实现符号计算的工具箱symbolic math toolbox. matlab符号计算

2、能实现微积分计算、线性代数计算、符号计算能实现微积分计算、线性代数计算、方程求解、积分变换等。在符号计算中方程求解、积分变换等。在符号计算中, 符号表达式是符号表达式是主要操作对象主要操作对象.符号表达式符号表达式符号变量、运算符、函数、数字组成符号变量、运算符、函数、数字组成表达式中的符号变量必须事先创建表达式中的符号变量必须事先创建.创建的常方法是用创建的常方法是用关键词关键词syms 空格后跟空格后跟符号变量符号变量。2/16例例1. 用符号表达式定义用符号表达式定义 f = e 0.2x sin (0.5x)并绘图并绘图.syms x ; f = exp(-0.2*x)*sin(0.5

3、*x); ezplot(f,0,2*pi)3/16syms 符号变量符号变量1 符号变量符号变量2 f = exp(-1/5*x) * sin(1/2*x)ezplot( f )绘表达式绘表达式f=f(x)的图的图绘图区域绘图区域 -2*pix2*piezplot(f, a,b)绘图区域绘图区域：axa时的极限时的极限symsum(f,k,m,n) 求级数求级数7/16 nmkkf)(例例5.绘函数绘函数 的图形的图形,并分别计算并分别计算 一阶和二阶导数一阶和二阶导数syms xf=3*x2+6*x-1;g=x2+x-3;r=f/g;ezplot(r,-10,10)r1=diff(r,x);

4、simplify(r1)f1,g1=numden(r1)r2=diff(r,x,2);simplify(r2)f2,g2=numden(r2)8/163163)(22 xxxxxr222)3()17163()( xxxxxr3223)3()4151243(2)( xxxxxxrx=double(solve(f1,x);x = -1.4648 -3.8685例例6. 求例求例5中函数的极值点与拐点中函数的极值点与拐点x1=x(1);x2=x(2);y1=subs(r,x,x1);y2=subs(r,x,x2); ezplot(r,-10,10),hold onplot(x1,y1,ro,x2,y

5、2,ro)9/16极值点分布极值点分布xx=double(solve(f2)xx = -5.2635 -1.3682 - 0.8511i -1.3682 + 0.8511ix0=xx(1)y0=subs(r,x,x0)ezplot(r,-10,10),hold onplot(x0,y0,ro)求例求例5 5中函数的拐点中函数的拐点10/16拐点分布拐点分布例例7. 计算不定积分计算不定积分 dxbxeax)sin(syms a b xf=exp(a*x)*sin(b*x); int(f,x);pretty(ans)11/16 b exp(a x) cos(b x) a exp(a x) sin

6、(b x) - - + - 2 2 2 2 a + b a + bg=simplify(ans) g = exp(a*x)*(-b*cos(b*x)+a*sin(b*x)/(a2+b2)s2=symsum(k2,k,1,n); factor(s2) ans = 1/6*n*(n+1)*(2*n+1) nkks11例例8.计算级数计算级数 nkks122syms k ns1=symsum(k,k,1,n);factor(s1) ans = 1/2*n*(n+1)s3=symsum(1/k2,k,1,inf) s3= 1/6*pi2 1231kks12/16taylor(f,n,x) n-1次麦克

7、劳林多项式展开次麦克劳林多项式展开taylor(f,n,x,a) a点的点的n-1次泰勒多项式展开次泰勒多项式展开.)cos(451)(xxf 例例9. 麦克劳林展开函数麦克劳林展开函数syms xf=1/(5+4*cos(x)t=taylor(f,8)pretty(t) 2 4 49 6 1/9 + 2/81 x + 5/1458 x + - x 13122013/16命令格式命令格式：dsolve(eq1,con1,x)y的一阶导数的一阶导数 dy, y的二阶导数的二阶导数 d2y例例10.解微分方程解微分方程14/16y = dsolve(dy=1/(1+x2)-2*y2,y(0) =

8、0,x)y = 2*x/(2*x2+2)22211yxy 符号解符号解： y(x)= x / (1 + x 2) 练习题与思考题练习题与思考题15/161.用用matlab求极限求极限xxxx30sinsintanlim 2.用用matlab求导数求导数xxxfcos1sin1)( 3.用用matlab求定积分求定积分 022dxexx4.用用matlab求级数和求级数和 12)12(1nnn5.计算曲线段计算曲线段 f(x)=exp(a x)sin(b x), 绕绕x轴旋转的旋转曲面面积轴旋转的旋转曲面面积 20 x6. 编写一个函数文件，用于计算椭圆周长，要求输编写一个函数文件，用于计算椭圆周长，要求输入变