数学实验高等数学计算.pptVIP

级数求和运算 S=symsum(f，n，a，b) 例11.计算级数 S2=symsum(k^2,k,1,n); factor(S2) S2= 1/6*n*(n+1)*(2*n+1) syms k n S1=symsum(k,k,1,n);factor(S1) S1 = 1/2*n*(n+1) S3=symsum(1/k^2,k,1,inf) S3= 1/6*pi^2 12/16 命令格式：dsolve(‘eq1’,···,’con1’,···,’x’) y的一阶导数—— Dy, y的二阶导数—— D2y 例10.解微分方程 14/16 y = dsolve(Dy=1/(1+x^2)-2*y^2,y(0) = 0,x) y = 2*x/(2*x^2+2) 符号解： y(x)= x / (1 + x 2) 旋转曲面绘制方法 非负函数 y =f(x)在有限区间上的图形为上半平面的一条曲线，曲线绕x轴旋转时，产生以x为对称轴的旋转曲面，方程为 绘制网面需创建三维坐标矩阵,对某一确定的x=t,旋转曲面上对应于过点x=t垂直于X轴的圆,该圆周上所有点的X坐标不变，y和z的坐标则满足圆的方程 f=inline(exp(-0.2*x).*sin(0.5*x)); t=(0:20)*pi/10; theta=t;r=f(t); x=t*ones(size(t)); y=r*cos(theta); z=r*sin(theta); mesh(x,y,z) colormap([0 0 0]) axis off view(-17,54) 曲线 绕X轴旋转图形绘制 练习题与思考题 15/16 1.用MATLAB求极限 2.用MATLAB求导数 3.用MATLAB求定积分 4.用MATLAB求级数和 5.计算曲线段 f(x)=exp(a x)sin(b x), 绕X轴旋转的旋转曲面面积 6. 编写一个函数文件，用于计算椭圆周长，要求输入变量是椭圆的长半轴和短半轴。 7. 符号变量与数值型变量有何区别？ 8. 用两个窗口绘下列函数图形,分析函数性质 16/16 思考题与练习题 用syms x 定义了符号变量,表达式 Y=exp(-0.2*x)*sin(0.5*x)与一般表达式有何不同 2.定积分符号计算与数值计算有何不同？ 3.旋转曲面的面积计算公式如何构造？ 4.写出曲线 y=f(x)绕y轴旋转的旋转曲面方程 5.下面两个曲面是由同一个平面曲线旋转产生的,这个平面曲线的方程是什么? 符号表达式及其应用 微积分基本运算 Taylor级数展开 常微分方程求解 ? ? ? ? 第3章 高等数学计算 MATLAB的符号计算 在数学、物理及各类工程应用中, 除了数值计算外, 还常常需用符号计算. 一般的计算机语言平台只能 实现数值计算, 早期的MATLAB也不能实现符号计 算。1993年MathWorks公司购买了MAPLE的使用 权,并开发出实现符号计算的工具箱——Symbolic Math Toolbox. 符号计算需要使用专门函数，功能包括有微积分、 线性代数、方程求解、积分变换等。在MATLAB数 值计算和字符串操作中,变量通过赋值语句创建.但 是符号计算中的变量在使用前，必须要用关键词 syms (或其它函数)创建。 2/16 例1. 用符号表达式定义 f = e – 0.2x sin (0.5x)并绘图. syms x ; f = exp(-0.2*x)*sin(0.5*x); ezplot(f,[0,2*pi]) 3/16 syms 符号变量1 符号变量2 … f = exp(-1/5*x) * sin(1/2*x) ezplot( f ) 绘表达式f=f(x)的图 绘图区域 -2*pix2*pi ezplot(f, [a,b]) 绘图区域：axb S1=subs(S, ‘old’, ‘new’)修改表达式 例2. 输入不同的参数a,b绘制函数的图形 f(x)=exp(a x)sin(b x) function mlab32(a1,b1) syms a b x f=exp(a*x)*sin(b*x); f1=subs(f,a,a1); f2=subs(f1,b,b1); ezplot(f2,[0,2*pi]) 符号表达式中变量替换 4/16 mlab32(-0.6,0.5)