资源描述：
《Jttpqq数值分析实验研究分析报告.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、个人收集整理仅供参考学习秋风清，秋月明,落叶聚还散,寒鸦栖复惊.《数值分析》实验报告册姓名:学号:专业:年级:计算机科学学院计算机应用教研室2008年春季学期-21-/21个人收集整理仅供参考学习目录实验一3b5E2RGbCAP实验二5实验三7实验四10实验五12实验六15实验七……………………………………………………………………18实验一一、课题名称非线性方程数值解法二、目地和意义-21-/21个人收集整理仅供参考学习学会常用地插值方法，求函数地近似表达式，以解决其它实际问题；明确插值多项式和分段插值多项式各自地优缺点；熟悉插值方法地程序编制；如果绘出插值函数地曲线，观察其光滑性.
2、p1EanqFDPw三、计算公式Lagrange插值公式：牛顿插值公式：四、结构程序设计程序设计:#include"math.h"floatf(floatx){return((x*x*x-1)/3);/*牛顿迭代函数*/}main(){floatx1,x2,eps,d;intk=0;clrscr();printf("inputx1=");/*输入迭代初值*/scanf("%f",&x1);printf("inputeps=");/*输入求解精度eps*/scanf("%f",&eps);do{k++;x1=x2;x2=f(x1);printf("%d%f",k,x2)
3、;}while(fabs(x2-x1)>=eps);printf("therootoff(x)=0isx=%f,k=%d",x2,k);/*输出x和迭代次数k*/DXDiTa9E3dgetch();}-21-/21个人收集整理仅供参考学习五、结果讨论和分析计算结果分析：将六种迭代格式分别代入程序试验：（1）第一种格式：无论何值都无法求出，即发散（2）第二种格式：初值为任意地x（x2<=1），精度为0.00001X=-0.347296，k=6其他值为发散.（3）第三种格式：初值为任意地x（x>0），精度为0.0001X=1.879372，k=10其他值为发散.（4）第四种格式：初值
4、为任意值，精度为0.00001X=-0.347296,k=5（5）第五种格式：初值为任意值,精度为0.00001X=-0.347296,k=4（6）第六种格式：初值为任意值，精度为0.00001X=-0.347296,k=4由此可知不同地初值对公式地计算有影响，当初值不满足函数地收敛条件时，无法计算结果，函数发散.精度地大小不同也使迭代函数迭代地次数不同，从而影响xn地近似程度.实验二一、课题名称解线性方程组地直接方法二、目地和意义-21-/21个人收集整理仅供参考学习掌握线性方程组直接接法地基本思想；了解不同数值方法解线性方程组地原理、实现条件、使用范围、计算公式；培养编程与上机调
5、试能力.RTCrpUDGiT三、计算公式消去法设a(k)kk=0,对k=1，2，……，n-1计算mik=a(k)ik/a(k)kka(k+1)ij=a(k)ij-mika(k)kji,j=k+1,k+2,……,nb(k+1)i=b(k)i-mikb(k)knxn=b(n)n/a(n)nnj=i+1xi=(b(i)i-Σa(i)ijxj)/a(i)iii=n-1,n-2,……,15PCzVD7HxA平方根法追赶法lij=(aii-Σl2ik)1/2Ly=flji=(aji-Σljklik)/liij=i+1,i+2,……,nUx=yjLBHrnAILgy1=f1/l1y2=(fi-ai
6、yi-1)/lii=2,3,……,n四、结构程序设计用追赶法求解线性方程组#include"stdio.h"main(){FILE*f;doublea[15],b[15],c[15],d[15];doublet;inti,n;f=fopen("zgf.dat","r");fscanf(f,"%d",&n);fscanf(f,"%lf%lf%lf",&b[1],&c[1],&d[1]);for(i=2;i<=n-1;i++)-21-/21个人收集整理仅供参考学习{fscanf(f,"%lf%lf%lf%lf",&a[i],&b[i],&c[i],&d[i]);}fscanf(f,"%l
7、f%lf%lf",&a[n],&b[n],&d[n]);fclose(f);c[1]=c[1]/b[1];d[1]=d[1]/b[1];for(i=2;i<=n-1;i++){t=b[i]-c[i-1]*a[i];c[i]=c[i]/t;d[i]=(d[i]-d[i-1]*a[i])/t;}d[n]=(d[n]-d[n-1]*a[n])/(b[n]-c[n-1]*a[n]);for(i==n-1;i>=1;i--)d[i]=d[i]-c[i]*d[i+1]
