计算方法实验报告.docx

1、计算方法实验报告江 西 科 技 师 范 学 院实 验 报 告 课 程 计算方法 系 别 班 级 _ _ 学 号 _ _ 姓 名 _ 目录实验一 误差的传播与估计实验二 拉格朗日插值多项式实验三 变步长复合梯形求积公式实验四 解非线性方程二分法实验五 一元非线性方程的迭代解法实验六 列主元高斯消去法每次实验课必须带上此本子，以便教师检查预习情况和记录实验原始数据。实验时必须遵守实验规则。用正确的理论指导实践袁必须人人亲自动手实验，但反对盲目乱动，更不能无故损坏仪器设备。这是一份重要的不可多得的自我学习资料袁它将记录着你在大学生涯中的学习和学习成果。请你保留下来，若干年后再翻阅仍将感到十分新鲜，记

2、忆犹新。它将推动你在人生奋斗的道路上永往直前！实验一 误差的传播与估计一、 实验课程名称 数学实验二、 试验项目名称 误差的传播与估计三、 实验目的和要求理解误差在算术运算中的传播方式及如何在算术运算中控制误差的传播。四、 实验内容和原理内容：计算积分原理：算术运算中的舍入误差在运算过程中的传播常因算法而异，不同算法的数值稳定性往往差别很大。五、 主要仪器及耗材计算机与Windows 2000/XP系统；Turbo C/Win-TC/VB等软件。六、 操作方法与实验步骤步骤：（1）由分部积分法建立两种递推公式： （2）对这两种递推式进行程序实现，分别求出的值。 1）fprintf(当n=%1.

3、0f,E1=%1.15fn,1/exp(1) En=1/exp(1); for n=2:20 En=1-n*En; fprintf(当n=%1.0f,E%1.0f=%1.15fn,n,n,En) end结果：当n=1,E1=0.367879441171442当n=2,E2=0.264241117657115当n=3,E3=0.207276647028654当n=4,E4=0.170893411885385当n=5,E5=0.145532940573073当n=6,E6=0.126802356561559当n=7,E7=0.112383504069084当n=8,E8=0.10093196744

4、7330当.0916*当.0838*当.0773*当.0717*当n=13,E13=0.066947434300630当n=14,E14=0.062735919791180当n=15,E15=0.058961203132299当.0566*当.0374*当n=18,E18=0.325949464266159当n=19,E19=-5.193039821057027当n=20,E20=104.860796421140550 2）fprintf(当n=20,E1=%1.15fn,1/21) En=1/21; for n=20:-1:2 En=(1-En)/n; fprintf(当n=%1.0f,E%

5、1.0f=%1.15fn,n-1,n-1,En) end 结果：当.0476*当.0476*当n=18,E18=0.050125313283208当n=17,E17=0.052770815928711当.0557*当n=15,E15=0.059017539764444当n=14,E14=0.062732164015704当n=13,E13=0.066947702570307当.0717*当.0773*当.0838*当.0916*当n=8,E8=0.100931967445593当n=7,E7=0.112383504069301当n=6,E6=0.126802356561528当n=5,E5=0

6、.145532940573079当n=4,E4=0.170893411885384当n=3,E3=0.207276647028654当n=2,E2=0.264241117657115当n=1,E1=0.367879441171442七：分析：算术运算中的舍入误差在运算过程中的传播常因算法而异，不同的算法的数值稳定性往往差别很大实验二 拉格朗日插值多项式一、 实验课程名称 数学实验二、 实验项目名称 拉格朗日插值多项式三、 实验目的和要求理解插值函数的概念及插值多项式的唯一性，掌握代数插值中拉格朗日插值多项式的基本原理及程序实现。四、 实验内容和原理内容：从函数表：0.00.10.1950.30

7、.4010.50.398940.396950.391420.381380.368120.35206出发，用拉格朗日插值多项式计算的近似值。原理：利用个插值节点可构造一个最高次数不超过次的代数插值多项式，利用该代数插值多项式可做相应的函数逼近。五、 主要仪器及耗材计算机与Windows 2000/XP系统；Turbo C/Win-TC/VB等软件。六、 操作方法与实验步骤步骤：（1）根据实验内容，给出问题的程序代码； （2）上机输入和调试所编程序； （3）调试成功后，整理出实验报告。 function f=fun4(t) format long x=0.0 0.1 0.195 0.3 0.401

8、 0.5; y=0.39894 0.39695 0.39142 0.38138 0.36812 0.35206; s=0; for i=1:6 s1=1;s2=1; for j=1:i-1 s1=s1*(t-x(j); end for j=i+1:6 s2=s2*(t-x(j); end p=s1*s2; s3=1;s4=1; for j=1:i-1 s3=s3*(x(i)-x(j); end for j=i+1:6 s4=s4*(x(i)-x(j); end q=s3*s4; r=p/q; s=s+r*y(i);endf=s; 结果： fun4(0.15)ans = 0.3944728038

9、78061 fun4(0.31)ans = 0.380219062454732 fun4(0.47)ans = 0.357222112339485七、 实验结果与分析：结果： fun(o、15)=0.394472803878061 fun(0、31)=0.380219062454732 Fun(0、47)=0.357222112339485分析：利用n+1个节点可构造不超过n次袋鼠插值多项式实验三 变步长复合梯形求积公式一、 实验课程名称 数学实验二、 实验项目名称 变步长复合梯形求积公式三、 实验目的和要求理解积分数值计算方法的意义，明确数值积分中步长与精度的关系，掌握变步长复合梯形求积公式

10、的程序实现。四、 实验内容和原理内容：用变步长复合梯形求积公式计算积分的近似值，要求误差不超过。原理：将积分区间逐次分半，每分一次就用同一个复合积分公式算出相应的积分近似值，并利用前后两次计算结果来判断误差是否达到精度要求。五、 主要仪器及耗材计算机与Windows XP系统； C/VB编译环境。六、 操作方法与实验步骤步骤：（1）根据实验内容，给出问题的程序代码； （2）上机输入和调试所编程序； （3）调试成功后，整理出实验报告。 1）function f=fun2(x) f=4/(1+x2); 2）function f=t(h) a=0;b=1;n=(b-a)/h;s=0; for k=0

11、:n-1 s=s+fun2(a+k*h); end f=h/2*(fun2(a)+fun2(b)+2*s); 3）clc m=1;h=1; t1=t(h);t2=t(h/2); while abs(t1-t2)1/2*10(-5) m=m+1; t1=t2; t2=(4(m-1)*t(h/2)-t(h)/(4(m-1)-1); h=h/2;endmt2double(int(4/(1+x2),0,1)结果：m = 23t2 = 3.141593607264110ans = 3.141592653589793七、 实验结果与分析。结果与分析： m=23 t2=3.141593607264110 a

12、ns=3.141592653589793达到了精度要求。实验四 解非线性方程的二分法一、 实验课程名称 数学实验二、 实验项目名称 解非线性方程的二分法三、 实验目的和要求理解二分法的基本思想，掌握二分法的程序实现。四、 实验内容和原理内容：用二分法计算方程在区间内根的近似值，要求误差不超过。原理：将含根区间逐次分半，每分一次判断根所在的子区间，舍去无根区间，保留含根区间，依次类推，直至结果满足精度要求。五、 主要仪器及耗材计算机与Windows 2000/XP系统；Turbo C/Win-TC/VB等软件。六、 操作方法与实验步骤步骤：（1）根据实验内容，给出问题的程序代码； （2）上机输入

13、和调试所编程序； （3）调试成功后，整理出实验报告。 1、function f=f(x) f=x3-x2-1; 2、a=1;b=2; while abs(f(a)-f(b)10(-6) c=(a+b)/2; if f(a)*f(c)0.000001 x=x-y/z; y=x3-x-1;z=3*x2-1; endx结果：x = 1.324717957245858八、 实验结果与分析结果： x=1.324717957245858 ans=1.32471795724746分析：这种迭代法仅用于某个迭代公式反复校正根的近似值，使之逐步精确化，最后得到满足精确度要求的结果实验六 列主元高斯消去法一、 实

14、验课程名称 数学实验二、 实验项目名称 列主元高斯消去法三、 实验目的和要求理解线性方程组直接解法的基本思想，理解列主元高斯消去法的稳定性条件，掌握列主元高斯消去法的程序实现，体会模块化结构程序设计方法的优点。四、 实验内容和原理内容：用列主元消去法求解方程组 要求输出：计算解;原理：在高斯消去法的每一次消元计算之前，在系数矩阵或消元后的低阶矩阵中按照一定的范围选取绝对值最大的数做为主元素，以减少计算过程中舍入误差对解的影响。五、 主要仪器及耗材计算机与Windows 2000/XP系统；Turbo C/Win-TC/VB等软件。六、 操作方法与实验步骤步骤：（1）根据实验内容，给出问题的程序

15、代码； （2）上机输入和调试所编程序； （3）调试成功后，整理出实验报告。 a=1 1 1 1 1 1 1 7; 2 1 1 1 1 1 1 8; 3 2 1 1 1 1 1 10; 4 3 2 1 1 1 1 13; 5 4 3 2 1 1 1 17; 6 5 4 3 2 1 1 22; 7 6 5 4 3 2 1 28; m,n=size(a); x=zeros(m,1); for j=1:n-2 c=max(a(j:m,j); for i=j:m if a(i,j)=c t=a(j,:); a(j,:)=a(i,:); a(i,:)=t; endendfor i=j+1:m for p=

16、j+1:n a(i,p)=a(i,p)-a(i,j)*a(j,p)/a(j,j); endendfor i=j+1:m a(i,j)=0;endenda;x(m)=a(m,n)/a(m,m);for i=m-1:-1:1 s=0; for j=i+1:m s=s+a(i,j)*x(j); end x(i)=(a(i,n)-s)/a(i,i); end X A=1 1 1 1 1 1 1; 2 1 1 1 1 1 1; 3 2 1 1 1 1 1; 4 3 2 1 1 1 1; 5 4 3 2 1 1 1; 6 5 4 3 2 1 1; 7 6 5 4 3 2 1; B=7 8 10 13 17 22 28; X=AB七、 实验结果与分析结果与分析：结果：X = 1 1 1 2 1 3 1 4 1 5 1 6 1 7 1 8 1 9 1 10X = 1.000000000000000 1.000000000000000 1.000000000000000 1.000000000000002 0.999999999999996 1.000000000000003 0.999999999999999