工数部分题编程.doc

工科数学分析编程作业 制作人： 李光 学 号： 201181624 班 级：电医1101班 上一次打错啦！我学工数，不学微积分，更新一下！ 作业一（matlab版） 绘出sinx/x及(1+x)^(1/x)图像,并观察两个极限 思路及算法 利用matlab中的函数计算及二维图像绘制功能可绘出图像 算法简述：限定定义域——函数——取屏幕——图像及类型 源码 x=0.1: 0.1:10; y1=sin(x)./x; figure,plot(x,y1); y2=(1+x).^(1./x); figure,plot(x,y2) 图像 结论 由图可知，当x无限趋近于0，sinx/x趋近于1； 当x无限趋近于0，（1+x）^(1/x)趋近于2.7; 由此可推断，当x无限趋近于0，sinx/x极限为1； 当x无限趋近于0，（1+x）^(1/x) 极限为e； 作业三（c 语言版） 计算 , 的函数值。计算结果取7位有效数字。 思路及算法 采用数值积分的方法，使用for语句对所有的矩形面积求和，对积分部分进行近似计算：将0到X分割，把每小块面积求出后作和，规定步长为0可认为已取极限。当每次取不同的X值，计算出一个近似值带入函数得到一个函数值，输出即可。 算法简述（for循环细节省略）： //伪代码 for(k=1;k=30;k=k+1) { f___0.0; for(i=0;i=k*1.0/10;i=i+0 f__f+pow(e,-i*i/2)*0; f__0.5+(1.0/sqrt(2*pi))*f ; printf //伪代码 源码 #includestdio.h #includemath.h #define e 2#define pi 3.1415926 void main() { double i,f; int k; for(k=1;k=30;k=k+1) { f=0.0; for(i=0;i=k*1.0/10;i=i+0 f=f+pow(e,-i*i/2)*0; f=0.5+(1.0/sqrt(2*pi))*f ; printf(k=%2d时，函数值为%.7f\n,k,f); } 数据 结论 数值积分适合对一个很难利用牛顿莱布尼兹公式求解的问题求解，在这个过程中，随着步长的减小，计算精度越高，但是运算量也会增大，因此采用计算机编程计算是解决这类问题的方法。 由于步长较小，本人程序计算较慢，每x耗时5s左右。 计算结果如上图所示，其可精确到0.0000001. 附录： 测试机房：412#25 本报告为李光原创，并以实名上传至百度文库（文档名——工数部分题编程）。 如有雷同，纯属巧合！ 上一次打错啦！我学工数，不学微积分 BEGIN 宏定义常量e及pi for循环:控制x范围及划分 for循环: :控制x划分（取极限） 累加面积（e^(t^2/2)的） 计算F（X） 变量I,f,k定义 输 出 STOP