计算方法实验报告

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

实验一误差的传播与估计实验课程名称数学实验试验项目名称误差的传播与估计实验目的和要求理解误差在算术运算中的传播方式及如何在算术运算中控制误差的传播。实验内容和原理内容：计算积分原理：算术运算中的舍入误差在运算过程中的传播常因算法而异，不同算法的数值稳定性往往差别很大。主要仪器及耗材计算机与Windows2000/XP系统；TurboC/Win-TC/VB等软件。操作方法与实验步骤步骤：（1）由分部积分法建立两种递推公式：（2）对这两种递推式进行程序实现，分别求出的值。1）fprintf('当n=%1.0f,E1=%1.15f\n',1/exp(1))En=1/exp(1);forn=2:20En=1-n*En;fprintf('当n=%1.0f,E%1.0f=%1.15f\n',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.100931967447330当n=9,E9=0.091612292974027当n=10,E10=0.083877070259732当n=11,E11=0.077352227142953当n=12,E12=0.071773274284567当n=13,E13=0.066947434300630当n=14,E14=0.062735919791180当n=15,E15=0.058961203132299当n=16,E16=0.056620749883223当n=17,E17=0.037447251985213当n=18,E18=0.325949464266159当n=19,E19=-5.193039821057027当n=20,E20=104.8607964211405502）fprintf('当n=20,E1=%1.15f\n',1/21)En=1/21;forn=20:-1:2En=(1-En)/n;fprintf('当n=%1.0f,E%1.0f=%1.15f\n',n-1,n-1,En)end结果：当n=20,E1=0.047619047619048当n=19,E19=0.047619047619048当n=18,E18=0.050125313283208当n=17,E17=0.052770815928711当n=16,E16=0.055719363768899当n=15,E15=0.059017539764444当n=14,E14=0.062732164015704当n=13,E13=0.066947702570307当n=12,E12=0.071773253648438当n=11,E11=0.077352228862630当n=10,E10=0.083877070103397当n=9,E9=0.091612292989660当n=8,E8=0.100931967445593当n=7,E7=0.112383504069301当n=6,E6=0.126802356561528当n=5,E5=0.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.4010.50.398940.396950.391420.381380.368120.35206出发，用拉格朗日插值多项式计算的近似值。原理：利用个插值节点可构造一个最高次数不超过次的代数插值多项式，利用该代数插值多项式可做相应的函数逼近。主要仪器及耗材计算机与Windows2000/XP系统；TurboC/Win-TC/VB等软件。操作方法与实验步骤步骤：（1）根据实验内容，给出问题的程序代码；（2）上机输入和调试所编程序；（3）调试成功后，整理出实验报告。functionf=fun4(t)formatlongx=[0.00.10.1950.30.4010.5];y=[0.398940.396950.391420.381380.368120.35206];s=0;fori=1:6s1=1;s2=1;forj=1:i-1s1=s1*(t-x(j));endforj=i+1:6s2=s2*(t-x(j));endp=s1*s2;s3=1;s4=1;forj=1:i-1s3=s3*(x(i)-x(j));endforj=i+1:6s4=s4*(x(i)-x(j));endq=s3*s4;r=p/q;s=s+r*y(i);endf=s;结果：fun4(0.15)ans=0.394472803878061fun4(0.31)ans=0.380219062454732fun4(0.47)ans=0.357222112339485实验结果与分析：结果：fun(o、15)=0.394472803878061fun(0、31)=0.380219062454732Fun(0、47)=0.357222112339485分析：利用n+1个节点可构造不超过n次袋鼠插值多项式实验三变步长复合梯形求积公式实验课程名称数学实验实验项目名称变步长复合梯形求积公式实验目的和要求理解积分数值计算方法的意义，明确数值积分中步长与精度的关系，掌握变步长复合梯形求积公式的程序实现。实验内容和原理内容：用变步长复合梯形求积公式计算积分的近似值，要求误差不超过。原理：将积分区间逐次分半，每分一次就用同一个复合积分公式算出相应的积分近似值，并利用前后两次计算结果来判断误差是否达到精度要求。主要仪器及耗材计算机与WindowsXP系统； C/VB编译环境。操作方法与实验步骤步骤：（1）根据实验内容，给出问题的程序代码；（2）上机输入和调试所编程序；（3）调试成功后，整理出实验报告。1）functionf=fun2(x)f=4/(1+x^2);2）functionf=t(h)a=0;b=1;n=(b-a)/h;s=0;fork=0:n-1s=s+fun2(a+k*h);endf=h/2*(fun2(a)+fun2(b)+2*s);3）clcm=1;h=1;t1=t(h);t2=t(h/2);whileabs(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+x^2)',0,1))结果：m=23t2=3.141593607264110ans=3.141592653589793实验结果与分析。结果与分析：m=23t2=3.141593607264110ans=3.141592653589793达到了精度要求。实验四解非线性方程的二分法实验课程名称数学实验实验项目名称解非线性方程的二分法实验目的和要求理解二分法的基本思想，掌握二分法的程序实现。实验内容和原理内容：用二分法计算方程在区间内根的近似值，要求误差不超过。原理：将含根区间逐次分半，每分一次判断根所在的子区间，舍去无根区间，保留含根区间，依次类推，直至结果满足精度要求。主要仪器及耗材计算机与Windows2000/XP系统；TurboC/Win-TC/VB等软件。操作方法与实验步骤步骤：（1）根据实验内容，给出问题的程序代码；（2）上机输入和调试所编程序；（3）调试成功后，整理出实验报告。1、functionf=f(x)f=x^3-x^2-1;2、a=1;b=2;whileabs(f(a)-f(b))>10^(-6)c=(a+b)/2;iff(a)*f(c)<0a=a;b=c;elsea=c;b=b;endendalpha=(a+b)/2double(solve('x^3-x^2-1'))结果：alpha=1.465571284294128ans=1.465571231876768-0.232785615938384+0.792551992515448i-0.232785615938384-0.792551992515448i实验结果与分析结果与分析：alpha=1.465571284294128-0.2327856159383840.792551992515448i-0.2327856159383840.792551992515448实验五一元非线性方程的迭代解法实验课程名称数学实验实验项目名称一元非线性方程的迭代解法实验目的和要求通过该综合性实验的基本训练，使学生进一步掌握科学计算方法的基本知识；掌握使用常用计算机语言编程进行数值计算的方法；巩固学生对一元非线性方程的迭代解法的理解及对迭代解法收敛性及收敛速度的理解；提高学生对数值计算中误差分析的能力。实验内容和原理内容：用单点迭代法及牛顿迭代法解一元非线性方程在[1，3]内的根，要求误差不超过。原理：迭代法是一种逐次逼近法，这种方法使用某个迭代公式反复校正根的近似值，使之逐步精确化，最后得到满足精度要求的结果。主要仪器及耗材计算机与Windows2000/XP系统；TurboC/Win-TC/VB等软件。操作方法与实验步骤步骤：（1）建立迭代公式；（2）选择合适的初值；（3）编程计算；（4）上机输入和调试所编程序；（5）进行误差分析及收敛性分析；（6）整理出实验报告。a=1;b=3;x=(a+b)/2;y=x^3-x-1;z=3*x^2-1;whileabs(y/z)>0.000001x=x-y/z;y=x^3-x-1;z=3*x^2-1;endx结果：x=1.324717957245858实验结果与分析结果：x=1.324717957245858ans=1.32471795724746分析：这种迭代法仅用于某个迭代公式反复校正根的近似值，使之逐步精确化，最后得到满足精确度要求的结果实验六列主元高斯消去法实验课程名称数学实验实验项目名称列主元高斯消去法实验目的和要求理解线性方程组直接解法的基本思想，理解列主元高斯消去法的稳定性条件，掌握列主元高斯消去法的程序实现，体会模块化结构程序设计方法的优点。实验内容和原理内容：用列主元消去法求解方程组要求输出：计算解;原理：在高斯消去法的每一次消元计算之前，在系数矩阵或消元后的低阶矩阵中按照一定的范围选取绝对值最大的数做为主元素，以减少计算过程中舍入误差对解的影响。主要仪器及耗材计算机与Windows2000/XP系统；TurboC/Win-TC/VB等软件。操作方法与实验步骤步骤：（1）根据实验内容，给出问题的程序代码；（2）上机输入和调试所编程序；（3）调试成功后，整理出实验报告。a=[11111117;21111118;321111110;432111113;543211117;654321122;765432128];[m,n]=size(a);x=zeros(m,1);forj=1:n-2c=max(a(j:m,j));fori=j:mifa(i,j)==ct=a(j,:);a(j,:)=a(i,:);a(i,:)=t;endendfori=j+1:mforp=j+1:na(i,p)=a(i,p)-a(i,j)*a(j,p