数值分析c语言代码

《数值分析c语言代码》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、UniversityofElectronicScienceandTechnologyofChina信息与软件工程学院标准实验报告（实验）课程名称计算数学基础电子科技大学教务处制表电孑科技丈学实验报告学生姓名：何安迪学号：2013220701020指导教师：顾小丰实验地点：实验中心306实验时间：2015.11.28一、实验室名称：实验中心306二.实验项目名称：计算数学(数值分析)的核心算法实现三、实验原理通过《计算数学基础》的课堂学习，参照《数值分析简明教程》的内容，熟练掌握插值方法、数值积分、方程求根的迭代法、线性方程组

2、的迭代法、线性方程组的直接法等解题技巧。使用算法对例题进行求解，并且通过在VisualStudio中编写C程序实现解题的算法。四、实验目的基于课程的学习，使用VisualStudio编写C程序实现计算数学(数值分析)的核心算法，熟练掌握课堂上学习的各种算法，充分理解课本上的知识，并将之用于实例求解。五、实验内容1、最小二乘法：在卜1,1］区间上取N=21个等距结点，计算出以相应节点上於的值做为数据样木，作出m二3,5,7,9次的最小二乘多项式。并计算出不同阶最小二乘多项式给出的最小偏差20a(in)=^(yill(xi)-y

3、i)2<=0六、实验环境OS:MicrosoftWindows764bit环境：MicrosoftVisualStudio2012七、实验步骤1.实验一：最小二乘法，同样用到上面构造的高斯函数，根据最小二乘法的分划和多项式拟合公式构建实验二的函数voidproblem2(),代码如下。voidproblem2()//最小二乘法{doublex[21]={0}doublea[4][4]={0doubleb[6][6]={0doublec[8][8]={0y[21]={0}；},result_a[4]={0},result_b[

4、6]={0}，result_c[8]={0｝；｝；｝；doubled[10][10]=0},result_d[10]={0};doubleright[10]={0};double*result[4]={resuresuIt_b»result_c,result_d};for(inti=0;i<21;i++){x[i]=-1+i*2/20.0;y[i]=exp(x[i]);}for(inti=0;i<10;i++){for(intj=0;j<10;j++){for(intk=0;k<21;k++){d[i][j]+=pow(x[

5、k],i+j);}if(i<4&&j<4){a[i][j]=d[i][j];b[i][j]=d[i][j];c[i][j]=d[i][j];}elseif(i<6&&j<6){b[i][j]=d[i][j];c[i][j]=d[i][j];}elseif(i<8&&j<8){c[i][j]=d[i][j];for(inti=0;i<10;i++){for(intj=0；j<21;j++){right[i]+=pow(x[j],i)*y[j];}}right,resu4)；right,result_b.6)；right,res

6、ult_c.8)；right,result_d.10)；Gauss(*aJGauss(*bJGauss(*cJGauss(*dJfor(inti=0；i<4;i++){intk=i*2+3;printfC内"x的％d阶最小二乘多项式为:k);for(intj=0;j

7、(nn);}}八.实验数据及结果分析1.主界面，供选择模拟哪个实验的算法。棗Sfi1i组组三最©^线111213I4S16t±—匸匸UL〔r」*I•图1主界面2.录入2,选择实验一，计算得3、5、6、9阶的最小二乘多项式结果如下,根据验证结果正确。三次样乗航值得得到根号S的近似值为：2.22865014甲再次选择输入，2父近2,最小二乘法，在【7.11内的3阶最〃匸乘多顶式^.0.9955«346«x*0<0.99759453>xAl<0.54018639«xA2柱―】内「X的5阶最”匸乘多顶式为，1.9000383

8、0^^0<1.00001996^1<0.49925038^^2«0.16649757^3•0.M37484l-x*4KI-1U】内『浜的7阶最"匸乘多顶式为.0.99999983«x^0*0.99999992^*1<0e5««0577«x*2*0.16666786*x*3•8・04*636