资源描述：
《一元函数积分学与空间图形的画法》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、项目二一元函数积分学与空间图形的画法实验1一元函数积分学(基础实验)实验目的掌握用Mathematica计算不定积分与定积分的方法.通过作图和观察,深入理解定积分的概念和思想方法.初步了解定积分的近似计算方法.理解变上限积分的概念.提高应用定积分解决各种问题的能力.基本命令1.计算不定积分与定积分的命令Integrate求不定积分时,其基本格式为Integrate[f[x],x]如输入Integrate[x^2+a,x]则输出其中a是常数.注意积分常数C被省略.求定积分时,其基本格式为Integ
2、rate[f[x],{x,a,b}]其中a是积分下限,b是积分上限.如输入Integrate[Sin[x],{x,0,Pi/2}]则输出1注:Mathematica有很多的命令可以用相应的运算符号来代替.例如,命令Integrate可用积分号代替,命令Sum可以用连加号代替,命令Product可用连乘号代替.因此只要调出这些运算符号,就可以代替通过键盘输入命令.调用这些命令,只要打开左上角的File菜单,点击Palettes中的BasicCalculations,再点击Calculus就可以得到
3、不定积分号、定积分号、求和号、求偏导数号等等.为了行文方便,下面仍然使用键盘输入命令,但读者也可以尝试用这些数学符号直接计算.2.数值积分命令NIntegrate用于求定积分的近似值.其基本格式为NIntegrate[f[x],{x,a,b}]如输入NIntegrate[Sin[x^2],{x,0,1}]则输出0.3102683.循环语句For循环语句的基本形式是For[循环变量的起始值,测试条件,增量,运算对象]运行此命令时,将多次对后面的对象进行运算,直到循环变量不满足测试条件时为止.这里必
4、须用三个逗号分开这四个部分.如果运算对象由多个命令组成,命令之间用分号隔开.例如,输入t=0;For[j=1,j<=10,j++,t=t+j];t则循环变量j从取值1开始,到10结束.每次增加1.执行结果,输出变量t的最终值1+2+…+10=55.注:For语句中的j++实际表示j=j+1.74实验举例用定义计算定积分当在上连续时,有因此可将与作为的近似值.为了下面计算的方便,在例1.1中定义这两个近似值为和n的函数.例1.1(教材例1.1)计算的近似值.输入s1[f_,{a_,b_},n_]:
5、=N[(b-a)/n*Sum[f[a+k*(b-a)/n],{k,0,n-1}]];s2[f_,{a_,b_},n_]:=N[(b-a)/n*Sum[f[a+k*(b-a)/n],{k,1,n}]];再输入Clear[f];f[x_]=x^2;js1=Table[{2^n,s1[f,{0,1},2^n],s2[f,{0,1},2^n]},{n,1,10}];TableForm[js1,TableHeadings->{None,{"n","s1","s2"}}]则输出ns1s220.1250.62
6、540.218750.4687580.2734380.398438160.3027340.365234320.3178710.349121640.3255620.3411871280.3294370.337252560.3313830.3352895120.3323570.33431110240.3328450.333822这是的一系列近似值.且有例1.2计算的近似值.输入Clear[g];g[x_]=Sin[x]/x;js2=Table[{n,s2[g,{0,1},n]},{n,3,50}]则
7、得到定积分的一系列近似值:{{3,0.91687},{4,0.924697},{5,0.929226},…,{48,0.944421},{49,0.944455},{50,0.944488}}注：用这种方法(矩形法)得到的定积分的近似值随n收敛很慢.74可以用梯形法或抛物线法改进收敛速度(见教材中的有关章节).如果用Nintegrate命令可以得到本题的比较精确的近似值为0.946083.例1.3用定义求定积分的动画演示.输入Clear[f,x,a,b];f[x_]=x^2;a=0;b=1.5;
8、m=0;g1=Plot[f[x],{x,a,b},PlotStyle->{RGBColor[1,0,0]},DisplayFunction->Identity];For[j=3,j<=50,j+=2,m=j;tt1={};tt2={};For[i=0,i
