实验二MATLAB中极限和微分积分运算.pptVIP

实验二 MATLAB中的极限、微分和积分运算 ;一、实验目的;二、相关知识;MATLAB提供的命令函数limit()可以完成极限运算，其调用格式如下： limit(F,x,a,’left’) 该命令对表达式F求极限，独立变量x从左边趋于a，函数中除F外的参数均可省略，’left’可换成’right’。举例如下： ;例1：求极限 。 解：可用以下程序完成： clear F=sym(‘(1+a/x)^x’) limit(F,’x’,inf,’left’);结果为exp(a)。其中，语句F=sym(‘(1+a/x)^x’)表示定义符号表达式 也可用以下的语句来完成：;clear; syms x %这里是把x先说明成符号。 F=(1+a/x)^x %这里的定义形式和前面不同。 limit(F, x, inf, ’left’) %这里的x本身就是符号，因 此不需要单引号。;MATLAB提供的函数diff()可以完成对给定函数求导函数的运算，其调用格式如下： diff(fun,x,n) 其意义是求函数fun关于变量x的n阶导数，n为1时可省略。这里的fun用上例的后一种方式来定义较为妥当。我们看下面的例： ;;MATLAB提供的命令函数int()可以完成积分运算，其调用格式有如下几种： int(fun) 计算函数fun关于默认变量的不定积分 int(fun,x) 计算函数fun关于变量x的不定积分 int(fun,x, a, b)计算函数fun关于变量x从a到b的定积分 我们通过例子来学习具体的用法： ; ;;例3：计算二重积分 ，其中D为曲 线 和 所围成的区域。 ;例4：被积曲面S为球面 在第一卦限部分的外则，计算曲面积分 。;然后计算该二次积分，程序如下： ; 这里我们看到，所有的积分计算都是利用函数int完成的，当我们遇到二重积分、三重积分和曲线、曲面积分时需要先化为相应的累次积分，再用int来完成积分的计算。 ;综合实验;a=0;b=1;m=1000; s=0;H=exp(1);%s设置为落在曲边梯形内的点数 for i=1:m xi=rand();yi=H*rand(); if yiexp(xi) s=s+1; end; %如果随机点落在曲边梯形内,s增加1 end; fprintf(%s%g\n,exp(x)在[0,1]上的积分约等于,H*(b-a)*s/m) ;2.Matlab近似计算定积分rsums(f,a,b); rsums exp(-5*x^2);3.泰勒级数逼近计算器taylortool ;4.最大公约数问题;现在考虑这样的一个问题：任取两个正整数，它们互质的可能性有多大？;实际上上述所求概率的精确值应该是;5.勾股数;l=0; m=[]; for a=1:100 for c=a+1:100 b=sqrt(c^2-a^2); if (b==floor(b))(ba) l=l+1; m(:,l)=[a,b,c]; end end end m;考察c-b=1时的勾股数，将前面程序略加修改可得到下表：;c-b=2的情况留着自己讨论，现在给出c-b=3时的勾股数，可以得到下表