成都大学数值计算课程设计.docVIP

四川理工学院 《数值计算方法》课程设计 题目： 课程设计 专业： 信 计 班级： 2012级2班 姓名： 何 鑫 学号： 12071030205 目录 设计题一: 3 1.1问题分析及设计思路： 3 1.2程序清单： 3 1.3运行结果： 4 1.4结果分析与设计总结： 4 设计题二： 4 2.1问题分析及设计思路： 5 2.2程序清单： 5 2.3程序结果： 5 2.4结果分析与设计总结： 6 设计题三： 7 3.1问题分析及设计思路： 7 3.2程序清单： 7 3.3程序结果： 8 3.4结果分析与设计总结: 9 设计题四： 9 4.1问题分析及设计思路： 9 4.2程序清单： 10 4.3程序结果： 10 4.4结果分析与设计总结： 11 数值分析课程设计总结： 12 设计题一: 用复合梯形公式计算下面积分,取n=2,4,8,25,100: I= 并给出数值积分结果与精确值1之间的的误差. 1.1问题分析及设计思路： 复合梯形公式定理：将积分区间[a,b]分n等分，其中记分点xk=a+(k-1}h, (k=1,2,…,n)，其中 h=(b-a)/n，然后用f(x)在每个小区间上应用梯形公式即可。因此，在此题中将[0,]进行n等分，然后应用梯形公式依次计算每个分区间的积分，最后相加即为所计算的数值积分的结果，然后比较与精确值1之间的误差即可。程序设计思路：首先要在程序中编辑出精确值及积分的区间以便后边比较误差和等分区间，接着编写等分区间的程序与被积积分以及计算数值积分结果和精确值与计算结果差值（即为误差）的表示程序，先通过循环可依次计算n=2,4,8的结果，在分别单独计算n=25，100的结果，分别比较误差。 1.2程序清单： Inxact=1; a=0; b=pi/2; fprintf(\n 精确值Inxact=1\n); fprintf(\n n I Error\n); n=1; for k=1:3 n=2*n; h=(b-a)/n;i=1:n+1; x=a+(i-1)*h;f=sin(x); I=h*(sum(f)-(f(1)+f(length(f)))/2); fprintf( %3.0f %12.9f %12.9f\n, n, I, Inxact -I); end n=n+17; h=(b-a)/n;i=1:n+1; x=a+(i-1)*h;f=sin(x); I=h*(sum(f)-(f(1)+f(length(f)))/2); fprintf( %3.0f %12.9f %12.9f\n, n, I, Inxact -I); n=n*4; h=(b-a)/n;i=1:n+1; x=a+(i-1)*h;f=sin(x); I=h*(sum(f)-(f(1)+f(length(f)))/2); fprintf( %3.0f %12.9f %12.9f\n, n, I, Inxact -I); 1.3运行结果： 精确值Inxact=1 n I Error 2 0.948059449 0.051940551 4 0.987115801 0.012884199 8 0.996785172 0.003214828 25 0.999670992 0.000329008 100 0.999979438 0.000020562 1.4结果分析与设计总结： 分析：由程序的运行结果易得，当n的取值越小时，所得到的结果与精确值之间的误差越大，而当当n的取值越大时，所得到的结果与精确值之间的误差越小。由此可以得出，当区间等分越多，步长h的值越小，每个小区间计算的结果越精确，由此所计算出的结果误差越小。总结：在程序设计中，需要对定理有充分理解，灵活运用，将定理转换成matlab所能读取的程序；在编辑程序的过程中，要仔细，注意所保留的计算结果的有效数字