




版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、课程名称: 计算机数值方法 实验项目:方程求根,线性方程组的直接解法与迭代 解法,代数插值,最小二乘法拟合多项式实验地点: 逸夫楼402 专业班级: 学号: 学生姓名: 指导教师: 于亚男 2012年4月26日太原理工大学学生实验报告学院名称软件学院专业班级学号学生姓名 实验日期4月26日成绩课程名称计算机数值方法实验题目 方程求根1、 实验目的和要求(1)了解非线性方程求根的常见方法,如二分法、牛顿法、割线法。(2)加深对方程求根方法的认识,掌握算法。二、实验内容和原理熟悉使用二分法、迭代法、牛顿法、割线法等方法对给定的方程进行根的求解。选择上述方法中的两种方法求方程:f(x)=x3+4x2
2、-10=0在1,2内的一个实根,且要求满足精度|x*-xn|<0.5×10-53、 主要仪器设备 HP笔记本,VC+6.04、 操作方法与实验步骤 1.二分法:#include <stdio.h>#include <math.h>double f(double x) double y; y=pow(x,3)+4*pow(x,2)-10; return (y);main() double a,b,yl,yr,ym,x1,xr,xm;a=1;b=2; yl=f(a); yr=f(b); while(yl*yr>0) yl=f(a); yr=f(b);
3、xr=a; x1=b; while(fabs(xr-x1)>0.000001) xm=(x1+xr)/2; ym=f(xm); yl=f(x1); yr=f(xr); if(yl*ym<0&&yr*ym>0) xr=xm; else x1=xm;printf("%ft%ft%fn",x1,xr,xm); printf("n所求根为x=%fn",xm); return 0;2.牛顿法:#include <stdio.h>#include <math.h>double f(double x) retu
4、rn (pow(x,3)+4*pow(x,2)-10);double f1(double x) return (3*pow(x,2)+10*x);int main() double x,x1,y1,y2; printf("请输入一个任意实数:X="); scanf("%lf",&x); printf("方程的解为:n"); do x1=x; y1=f(x); y2=f1(x1); x=x1-y1/y2; while (fabs(x-x1)>=5e-6); printf(" %lfn",x1);3.追赶
5、法5、 实验数据记录和处理1. 二分法2. .牛顿法6、 实验结果与分析 通过这个两个程序可看出,二分法的计算量更大一些。7、 讨论、心得 通过这个实验,我了解了线性方程的一些求根方法,对于方程近似值的求解有了更多的理解。实验地点 指导教师于亚男太原理工大学学生实验报告学院名称专业班级学号学生姓名 实验日期4月26日成绩课程名称计算机数值方法实验题目 求解线性方程组一、实验目的和要求(1) 了解直接法解线性方程组(2) 对Gauss消元法、LU分解法、追赶法有更深刻的理解二、实验内容和原理合理利用Gauss消元法、LU分解法、追赶法求解下列方程组: (n=5,10,100,)三主要仪器设备 H
6、P笔记本,VC+6.0四操作方法与实验步骤 1.高斯:#include<stdio.h> void main() int i,j,f,n; float a55,b5,x5,l,k; printf("Input n:"); scanf("%d",&n); printf("Input a:"); for(i=0;i<n;i+) for(f=0;f<n;f+) scanf("%f",&aif); printf("Input b:"); for(i=0;i<
7、n;i+) scanf("%f",&bi); for(i=0;i<n-1;i+) for(j=i+1;j<n;j+) l=-aji/aii; bj=bj+bi*l; for(f=i;f<n;f+) ajf=ajf+aif*l; xn-1=bn-1/an-1n-1; for(i=n-2;i>=0;i-) xi=bi; k=0; for(f=i+1;f<n;f+) k=k-aif*xf; xi=(k+bi)/aii; printf("The answer is "); for(i=0;i<n;i+) printf(
8、"%.4f ",xi);2.LU分解:#include <stdio.h> #include <math.h> #define L 30 double aLL,bL,lLL,uLL,xL,yL; int main() int n,i,j,k,r; printf("n=n"); scanf("%d",&n); printf("ann=n"); for(i=1;i<=n;+i) for(j=1;j<=n;+j) scanf("%lf",&aij);
9、 printf("bn=n"); for(i=1;i<=n;+i) scanf("%lf",&bi); for(i=1;i<=n;+i) for(j=1;j<=n;+j) lij=0; uij=0.0; for(k=1;k<=n;+k) for(j=k;j<=n;+j) ukj=akj; for(r=1;r<k;+r) ukj-=lkr*urj; for(i=k+1;i<=n;+i) lik=aik; for(r=1;r<k;+r) lik-=lir*urk; lik/= ukk; lkk=1.0;
10、 for(i=1;i<=n;+i) yi = bi; for(j=1;j<i;+j) yi-=lij*yj; for(i=n;i>0;-i) xi = yi; for(j=i+1;j<=n;+j) xi-=uij*xj; xi/= uii; for(i=1;i<=n;+i) printf("%0.2lfn",xi); return 0; 五实验数据记录和处理1.高斯2.LU分解六、实验结果与分析 本次实验数据较多,在输入上要多费点功夫,一不小心就全部都错了。在今后编程过程中,一定要小心谨慎。七、讨论、心得 通过本次实验,我深刻理解了直接法在计算
11、机上解线性方程组的有效性,对于Gauss消元法、LU分解法也有了深刻的理解。实验地点 逸夫楼402指导教师于亚男太原理工大学学生实验报告学院名称软件学院专业班级学号学生姓名 实验日期4月26日成绩课程名称计算机数值方法实验题目 线性方程组的迭代解法1、 实验目的和要求掌握雅可比迭代法和高斯-赛德尔迭代法对方程组进行求解二、实验内容和原理使用雅可比迭代法或高斯-赛德尔迭代法对下列方程组进行求解。三、主要仪器设备 HP笔记本,VC+6.0四、操作方法与实验步骤 雅可比迭代法:#include <stdio.h> #include <math.h> int main() do
12、uble a33=10,-1,-2,-1,10,-2,-1,-1,5,b3=7.2,8.3,4.2;/定义方程组 float x3=0,0,0,sum; int i,j,k,n=3; printf("tt X1tt X2tt X3n"); for(k=0;k<8;k+) for(i=0;i<3;i+) sum=0; for(j=0;j<n;j+) if(i=j)continue; sum=sum+aij*xj; xi=(bi-sum)/aii; printf("第%d次迭代:t",k+1); for(i=0;i<n;i+) pri
13、ntf("%ft",xi); printf("n"); 五、实验数据记录和处理六、实验结果与分析 在本次实验中,编程不太容易,对c/c+的学习应该更进一步。七、讨论、心得 通过这次实验,我明白了雅克比迭代的一般性解法,对于编程的应用也有了更深刻的理解。实验地点 逸夫楼402指导教师于亚男太原理工大学学生实验报告学院名称软件学院专业班级学号学生姓名 实验日期4月26日成绩课程名称计算机数值方法实验题目 代数插值1、 实验目的和要求 掌握拉格朗日插值法和牛顿插值法求近似解二、实验内容和原理使用拉格朗日插值法或牛顿插值法求解:已知f(x)在6个点的函数值如下表
14、所示,运用插值方法,求f(0.596)的近似值。X0.400.550.650.800.901.05f(x)0.410750.578150.696750.888111.026521.25386三、主要仪器设备 HP笔记本,VC+6.0四、操作方法与实验步骤#include "stdafx.h"#include <stdio.h> #include <stdlib.h> #include <conio.h> #include <malloc.h>void difference(float *x,float *y,int n) fl
15、oat *f; int k,i; f=(float *) malloc (n*sizeof(float); for(k=1;k<=n;k+) f0=yk; for(i=0;i<k;i+)fi+1=(fi-yi)/(xk-xi); yk=fk; return; int main() int i,n; float x10,y10,xx,yy; printf("输?入?结á点?个?数簓n:阰"); scanf("%d",&n);printf("n"); for(i=0;i<=n-1;i+) printf(
16、"x%d=",i); scanf("%f",&xi);printf("y%d=",i);scanf("%f",&yi);printf("n"); difference(x,(float *)y,n); printf("所ù求ó插?值X:阰"); scanf("%f",&xx); yy=y20; for(i=n-1;i>=0;i-)yy=yy*(xx-xi)+yi; printf("n近ü
17、似?值为a:阰F(%f)=%fn",xx,yy); 五、实验数据记录和处理六、实验结果与分析 拉格朗日插值的优点是插值多项式特别容易建立,缺点是增加节点是原有多项式不能利用,必须重新建立,即所有基函数都要重新计算,这就造成计算量的浪费。所以要用到牛顿插值多项式,两种方法相结合才能更好的解决多项式的问题。实验地点 逸夫楼402指导教师于亚男太原理工大学学生实验报告学院名称软件学院专业班级学号学生姓名实验日期4月26日成绩课程名称计算机数值方法实验题目 最小二乘法拟合多项式一、实验目的和要求掌握用最小二乘法拟合数据的多项式,并求平方误差二、实验内容和原理给定数据点(xi ,yi),用最小
18、二乘法拟合数据的多项式,并求平方误差。xi00.50.60.70.80.91.0yi11.751.962.192.442.713.00三、主要仪器设备 HP笔记本,VC+6.0四、操作方法与实验步骤#include<iostream.h>#include<fstream.h>#define N 15double power(double &a,int n)double b=1;for(int i=0;i<n;i+)b*=a;return b;void Gauss();double XN,YN,sumXN,sumYN,aNN,bN,lNN,xN;void m
19、ain()ofstream outdata;ifstream indata;double s;int i,j,k,n,index;cout<<"请输入已知点的个数n="cin>>n;cout<<endl;cout<<"请输入X和Y:"<<endl; /输入给定数据for(i=0;i<n;i+)cout<<"X"<<i<<"="cin>>Xi;sumX1+=Xi;cout<<"Y&q
20、uot;<<i<<"="cin>>Yi;sumY1+=Yi;cout<<endl;cout<<"sumX1="<<sumX1<<"t"<<"sumY1="<<sumY1<<endl;cout<<"请输入拟合次数index="cin>>index;cout<<endl;i=n;sumX0=i;for(i=2;i<=2*index;i+)sumXi=0;for(j=0;j<n;j+)sumXi+=power(Xj,i);cout<<"sumX"<<i<<"="<<sumXi<<endl;for(i=2;i<=index+1;i+)sumYi=0;for(j=0;j<n;j+)sumYi+=power(Xj,i-1)*Yj;cout<<"sumY"<<i<<"="<
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 适应性教学策略与实施计划
- 四年级下数学教案-空间与图形-苏教版
- 全国泰山版初中信息技术八年级上册第二章第一节《数字化视频播放》教学设计
- 视网膜坏死的护理
- 人才继任计划的制定与落实
- 工业园区安全隐患的评估与治理计划
- 2025年无机矿物填充塑料合作协议书
- 言语治疗辨音训练
- 九年级道德与法治下册 第一单元 我们共同的世界 第二课 构建人类命运共同体教学实录设计(pdf) 新人教版
- 南京会展工作简历模板
- 气瓶充装安全及培训课件PPT幻灯片
- (高清版)辐射供暖供冷技术规程JGJ142-2012
- JTT 1295—2019道路大型物件运输规范_(高清-最新)
- 电缆原材料检验规范
- 土壤固化土施工技术导则
- VAR模型Johansen协整检验在eviews中的具体操作步骤及结果解释
- 冷冻面团项目市场分析
- 北工大电子课程设计报告数电自行车里程表.
- 加油站法律法规符合性评价
- 5外科--丹毒下肢丹毒中医诊疗方案2017年版
- 数学物理方法第十一章PPT课件
评论
0/150
提交评论