计算方法与实习上机实验报告
计算方法上机实验报告
《计算方法》上机实验报告班级：XＸＸXXX小组成员:XXＸXXＸＸXXXＸXXXXXＸＸXXXＸXXXXXＸ任课教师：XXX二〇一八年五月二十五日前言通过进行多次得上机实验,我们结合课本上得内容以及老师对我们得指导,能够较为熟练地掌握Ｎewton迭代法、Ｊacｏｂi迭代法、Gausｓ－Seｉｄｅl 迭代法、Ｎewｔon 插值法、Ｌａgｒanｇe 插值法与Ｇaｕｓs 求积公式等六种算法得原理与使用方法，并参考课本例题进行了MATＬAB程序得编写。
以下为本次上机实验报告,按照实验内容共分为六部分.实验一：一、实验名称及题目:Newｔon迭代法例2、7(P3８）:应用Newｔon迭代法求在附近得数值解，并使其满足、二、解题思路：设就是得根,选取作为初始近似值，过点做曲线得切线,得方程为,求出与轴交点得横坐标,称为得一次近似值,过点做曲线得切线,求该切线与轴得横坐标称为得二次近似值，重复以上过程，得得近似值序列，把称为得次近似值，这种求解方法就就是牛顿迭代法。
三、Matlab程序代码：fｕncｔion newtｏn_itｅrａｔioｎ(x０,tol）syms z %定义自变量fｏrmaｔlｏnｇ％定义精度ｆ=ｚ*z＊z-z－1;f1=diff(f）；％求导y=suｂs(ｆ，ｚ,x0);y1=sｕbs(f１,ｚ，x0)；%向函数中代值x1=x0－y/y1; k=1；wｈilｅａｂｓ(ｘ1—x0)〉=ｔolx０=x1;y=subｓ（f,z,ｘ0)；y1=subｓ(ｆ1,ｚ,ｘ0);x1＝x0-ｙ/y1；k=k+1;ｅndｘ=douｂle（x１)Ｋ四、运行结果：实验二:一、实验名称及题目：Jacｏbｉ迭代法例３、7(P74）:试利用Jacobi迭代公式求解方程组要求数值解为方程组得精确解、二、解题思路：首先将方程组中得系数矩阵分解成三部分,即:，为对角阵,为下三角矩阵,为上三角矩阵。
之后确定迭代格式,,（, 即迭代次数),称为迭代矩阵。
上机实验报告(精选11篇)
上机实验报告篇1用户名se××××学号姓名学院①实验名称：②实验目的：③算法描述（可用文字描述，也可用流程图）：④源代码：（.c的文件）⑤用户屏幕（即程序运行时出现在机器上的画面）：2．对c文件的要求：程序应具有以下特点：a可读性：有注释。
b交互性：有输入提示。
c结构化程序设计风格：分层缩进、隔行书写。
3.上交时间：12月26日下午1点-6点，工程设计中心三楼教学组。
请注意：过时不候哟！四、实验报告内容0．顺序表的插入。
1.顺序表的删除。
2．带头结点的单链表的\'插入。
3.带头结点的单链表的删除。
注意：1.每个人只需在实验报告中完成上述4个项目中的一个，具体安排为：将自己的序号对4求余，得到的数即为应完成的项目的序号。
例如：序号为85的同学，85%4=1，即在实验报告中应完成顺序表的删除。
2.实验报告中的源代码应是通过编译链接即可运行的。
3.提交到个人空间中的内容应是上机实验中的全部内容。
上机实验报告篇2一、《软件技术基础》上机实验内容1．顺序表的建立、插入、删除。
2．带头结点的单链表的建立（用尾插法）、插入、删除。
二、提交到个人10m硬盘空间的内容及截止时间1．分别建立二个文件夹，取名为顺序表和单链表。
2．在这二个文件夹中，分别存放上述二个实验的相关文件。
每个文件夹中应有三个文件(.c文件、.obj文件和.exe文件)。
3. 截止时间：12月28日（18周周日）晚上关机时为止，届时服务器将关闭。
三、实验报告要求及上交时间（用a4纸打印）1．格式：《计算机软件技术基础》上机实验报告用户名se××××学号姓名学院①实验名称：②实验目的：③算法描述（可用文字描述，也可用流程图）：④源代码：（.c的文件）⑤用户屏幕（即程序运行时出现在机器上的画面）：2．对c文件的要求：程序应具有以下特点：a 可读性：有注释。
b 交互性：有输入提示。
(完整word版)计算方法A上机实验报告
计算方法A上机实验报告姓名：苏福班级：硕4020 学号：3114161019一、上机练习目的1）复习和巩固数值计算方法的基本数学模型，全面掌握运用计算机进行数值计算的具体过程及相关问题。
2）利用计算机语言独立编写、调试数值计算方法程序，培养学生利用计算机和所学理论知识分析解决实际问题的能力。
二、上机练习任务1）利用计算机语言编写并调试一系列数值方法计算通用程序，并能正确计算给定题目，掌握调试技能。
2）掌握文件使用编程技能，如文件的各类操作，数据格式设计、通用程序运行过程中文件输入输出运行方式设计等。
3）写出上机练习报告。
三、上机题目1. 共轭梯度法求解线性方程组。
（第三章）2. 三次样条插值（第四章）3. 龙贝格积分（第六章）4. 四阶龙格－库塔法求解常微分方程的初值问题四、上机报告题目1：共轭梯度法求解线性方程组1．算法原理共轭梯度法是把求解线性方程组的问题转化为求解一个与之等价的二次函数极小值的问题。
从任意给定的初始点出发，沿一组关于矩阵A共轭的方向进行线性搜索，在无舍入误差的假定下，最多迭代n 次（其中n 为矩阵A 的阶数），就可求得二次函数的极小值，也就求得了线性方程组Ax b =的解。
定理：设A 是n 阶对称正定矩阵，则x *是方程组Ax b =的解得充分必要条件是x *是二次函数1()2TT f x x Ax b x =-的极小点，即 ()()min nx R Ax b f x f x **∈=⇔=共轭梯度法的计算公式：(0)(0)(0)()()()()(1)()()(1)(1)(1)()()()(1)(1)()k T k k k T k k k k k k k k T k k k T k k k k k d r b Ax r d d Ad xx d r b Ax r Ad d Ad d r d ααββ++++++⎧==-⎪⎪=⎪⎪=+⎪⎨=-⎪⎪⎪=-⎪⎪=+⎩2. 程序框图（1）编写共轭梯度法求解对称正定矩阵的线性方程组见附录（myge.m）：function x=myge(A,b)输入对称正定矩阵及对应的列向量，初始向量设为0，精度取为810 。
计算方法上机实验
1.拉格朗日插值多项式，用于离散数据的拟合#include <stdio.h>#include <conio.h>#include <alloc.h>float lagrange(float *x,float *y,float xx,int n) /*拉格朗日插值算法*/{ int i,j;float *a,yy=0.0; /*a作为临时变量，记录拉格朗日插值多项式*/a=(float *)malloc(n*sizeof(float));for(i=0;i<=n-1;i++){ a[i]=y[i];for(j=0;j<=n-1;j++)if(j!=i) a[i]*=(xx-x[j])/(x[i]-x[j]);yy+=a[i];}free(a);return yy;}main(){ int i,n;float x[20],y[20],xx,yy;printf("Input n:");scanf("%d",&n);if(n>=20) {printf("Error!The value of n must in (0,20)."); getch();return 1;} if(n<=0) {printf("Error! The value of n must in (0,20)."); getch(); return 1;} for(i=0;i<=n-1;i++){ printf("x[%d]:",i);scanf("%f",&x[i]);}printf("\n");for(i=0;i<=n-1;i++){ printf("y[%d]:",i);scanf("%f",&y[i]);}printf("\n");printf("Input xx:");scanf("%f",&xx);yy=lagrange(x,y,xx,n);printf("x=%f,y=%f\n",xx,yy);getch();}2.牛顿插值多项式，用于离散数据的拟合#include <stdio.h>#include <conio.h>#include <alloc.h>void difference(float *x,float *y,int n){ float *f;int k,i;f=(float *)malloc(n*sizeof(float));for(k=1;k<=n;k++){ f[0]=y[k];for(i=0;i<k;i++)f[i+1]=(f[i]-y[i])/(x[k]-x[i]);y[k]=f[k];}return;}main(){ int i,n;float x[20],y[20],xx,yy;printf("Input n:");scanf("%d",&n);if(n>=20) {printf("Error! The value of n must in (0,20)."); getch(); return 1;} if(n<=0) {printf("Error! The value of n must in (0,20).");getch(); return 1;} for(i=0;i<=n-1;i++){ printf("x[%d]:",i);scanf("%f",&x[i]);}printf("\n");for(i=0;i<=n-1;i++){ printf("y[%d]:",i);scanf("%f",&y[i]);}printf("\n");difference(x,(float *)y,n);printf("Input xx:");scanf("%f",&xx);yy=y[20];for(i=n-1;i>=0;i--) yy=yy*(xx-x[i])+y[i];printf("NewtonInter(%f)=%f",xx,yy);getch();}3.高斯列主元消去法,求解其次线性方程组第一种#include<stdio.h>#include <math.h>#define N 20int main(){ int n,i,j,k;int mi,tmp,mx;float a[N][N],b[N],x[N];printf("\nInput n:");scanf("%d",&n);if(n>N){ printf("The input n should in(0,N)!\n");getch();return 1;}if(n<=0){ printf("The input n should in(0,N)!\n");getch();return 1;}printf("Now input a(i,j),i,j=0...%d:\n",n-1); for(i=0;i<n;i++){ for(j=0;j<n;j++)scanf("%f",&a[i][j]);}printf("Now input b(i),i,j=0...%d:\n",n-1);for(i=0;i<n;i++)scanf("%f",&b[i]);for(i=0;i<n-2;i++){ for(j=i+1,mi=i,mx=fabs(a[i][j]);j<n-1;j++) if(fabs(a[j][i])>mx){ mi=j;mx=fabs(a[j][i]);}if(i<mi){ tmp=b[i];b[i]=b[mi];b[mi]=tmp;for(j=i;j<n;j++){ tmp=a[i][j];a[i][j]=a[mi][j];a[mi][j]=tmp;}}for(j=i+1;j<n;j++){ tmp=-a[j][i]/a[i][i];b[j]+=b[i]*tmp;for(k=i;k<n;k++)a[j][k]+=a[i][k]*tmp;}}x[n-1]=b[n-1]/a[n-1][n-1];for(i=n-2;i>=0;i--){ x[i]=b[i];for(j=i+1;j<n;j++)x[i]-=a[i][j]*x[j];x[i]/=a[i][i];}for(i=0;i<n;i++)printf("Answer:\n x[%d]=%f\n",i,x[i]); getch();return 0;}第二种#include<math.h>#include<stdio.h>#define NUMBER 20#define Esc 0x1b#define Enter 0x0dfloat A[NUMBER][NUMBER+1] ,ark;int flag,n;exchange(int r,int k);float max(int k);message();main(){float x[NUMBER];int r,k,i,j;char celect;clrscr();printf("\n\nUse Gauss.");printf("\n\n1.Jie please press Enter."); printf("\n\n2.Exit press Esc.");celect=getch();if(celect==Esc)exit(0);printf("\n\n input n=");scanf("%d",&n);printf(" \n\nInput matrix A and B:"); for(i=1;i<=n;i++){printf("\n\nInput a%d1--a%d%d and b%d:",i,i,n,i);for(j=1;j<=n+1;j++) scanf("%f",&A[i][j]);}for(k=1;k<=n-1;k++){ark=max(k);if(ark==0){printf("\n\nIt's wrong!");message();}else if(flag!=k)exchange(flag,k);for(i=k+1;i<=n;i++)for(j=k+1;j<=n+1;j++)A[i][j]=A[i][j]-A[k][j]*A[i][k]/A[k][k];}x[n]=A[n][n+1]/A[n][n];for( k=n-1;k>=1;k--){float me=0;for(j=k+1;j<=n;j++){me=me+A[k][j]*x[j];}x[k]=(A[k][n+1]-me)/A[k][k];}for(i=1;i<=n;i++){printf(" \n\nx%d=%f",i,x[i]);}message();}exchange(int r,int k){int i;for(i=1;i<=n+1;i++)A[0][i]=A[r][i];for(i=1;i<=n+1;i++)A[r][i]=A[k][i];for(i=1;i<=n+1;i++)A[k][i]=A[0][i];}float max(int k){int i;float temp=0;for(i=k;i<=n;i++)if(fabs(A[i][k])>temp){temp=fabs(A[i][k]);flag=i;}return temp;}message(){printf("\n\n Go on Enter ,Exit press Esc!");switch(getch()){case Enter: main();case Esc: exit(0);default:{printf("\n\nInput error!");message();} }}4.龙贝格求积公式，求解定积分#include<stdio.h>#include<math.h>#define f(x) (sin(x)/x)#define N 20#define MAX 20#define a 2#define b 4#define e 0.00001float LBG(float p,float q,int n){ int i;float sum=0,h=(q-p)/n;for (i=1;i<n;i++)sum+=f(p+i*h);sum+=(f(p)+f(q))/2;return(h*sum);}void main(){ int i;int n=N,m=0;float T[MAX+1][2];T[0][1]=LBG(a,b,n);n*=2;for(m=1;m<MAX;m++){ for(i=0;i<m;i++)T[i][0]=T[i][1];T[0][1]=LBG(a,b,n);n*=2;for(i=1;i<=m;i++)T[i][1]=T[i-1][1]+(T[i-1][1]-T[i-1][0])/(pow(2,2*m)-1);if((T[m-1][1]<T[m][1]+e)&&(T[m-1][1]>T[m][1]-e)){ printf("Answer=%f\n",T[m][1]); getch();return ;}}}5.牛顿迭代公式，求方程的近似解#include<stdio.h>#include<math.h>#include<conio.h>#define N 100#define PS 1e-5#define TA 1e-5float Newton(float (*f)(float),float(*f1)(float),float x0 ) { float x1,d=0;int k=0;do{ x1= x0-f(x0)/f1(x0);if((k++>N)||(fabs(f1(x1))<PS)){ printf("\nFailed!");getch();exit();}d=(fabs(x1)<1?x1-x0:(x1-x0)/x1);x0=x1;printf("x(%d)=%f\n",k,x0);}while((fabs(d))>PS&&fabs(f(x1))>TA) ;return x1;}float f(float x){ return x*x*x+x*x-3*x-3; }float f1(float x){ return 3.0*x*x+2*x-3; }void main(){ float f(float);float f1(float);float x0,y0;printf("Input x0: ");scanf("%f",&x0);printf("x(0)=%f\n",x0);y0=Newton(f,f1,x0);printf("\nThe root is x=%f\n",y0); getch();}。
《数值计算方法》上机实验报告
《数值计算方法》上机实验报告华北电力大学实验名称数值il•算方法》上机实验课程名称数值计算方法专业班级：电力实08学生姓名：李超然学号:200801001008 成绩: 指导教师:郝育黔老师实验日期：2010年04月华北电力大学实验报告数值计算方法上机实验报吿一.各算法的算法原理及计算机程序框图1、牛顿法求解非线性方程*对于非线性方程，若已知根的一个近似值，将在处展开成一阶xxfx ()0, fx ()xkk泰勒公式"f 0 / 2 八八,fxfxfxxxxx 0 0 0 0 0 kkkk2!忽略高次项，有，fxfxfxxx 0 ()()()，，, kkk右端是直线方程，用这个直线方程来近似非线性方程。
将非线性方程的**根代入，即fx ()0, X ，* fxfxxx 0 0 0 0, ,, kkkfx 0 fx 0 0,解出fX 0 *k XX,, k' fx 0 k水将右端取为，则是比更接近于的近似值，即xxxxk, Ik, Ikfx ()k 八XX, Ikk* fx()k这就是牛顿迭代公式。
,2,计算机程序框图:，见,,3,输入变量、输出变量说明:X输入变量:迭代初值，迭代精度，迭代最大次数，\0输出变量:当前迭代次数，当前迭代值xkl,4,具体算例及求解结果:2/16华北电力大学实验报吿开始读入l>k/fx()0?,0fx 0 Oxx，，01* fx ()0XX,,，?10kk, ,1，kN, ?xx, 10输出迭代输出X输出奇异标志1失败标志,3,输入变量、输出变量说明: 结束例：导出计算的牛顿迭代公式，并il •算。
（课本P39例2-16） 115cc （0）, 求解结果:10. 75000010.72383710. 72380510. 7238052、列主元素消去法求解线性方程组，1,算法原理:高斯消去法是利用现行方程组初等变换中的一种变换，即用一个不为零的数乘 -个 方程后加只另一个方程，使方程组变成同解的上三角方程组，然后再自下而上 对上三角3/16华北电力大学实验报告方程组求解。
东南大学计算方法上机报告实验报告完整版
实习题11. 用两种不同的顺序计算644834.11000012≈∑=-n n，试分析其误差的变化解：从n=1开始累加，n 逐步增大，直到n=10000；从n=10000开始累加，n 逐步减小，直至1。
算法1的C 语言程序如下： #include<stdio.h> #include<math.h> void main() { float n=0.0; int i; for(i=1;i<=10000;i++) { n=n+1.0/(i*i); } printf("%-100f",n); printf("\n"); float m=0.0; int j; for(j=10000;j>=1;j--) { m=m+1.0/(j*j); } printf("%-7f",m); printf("\n"); }运行后结果如下：结论： 4.设∑=-=Nj N j S 2211，已知其精确值为)11123(21+--N N 。
1）编制按从大到小的顺序计算N S 的程序； 2）编制按从小到大的顺序计算N S 的程序；3）按2种顺序分别计算30000100001000,,S S S ，并指出有效位数。
解：1）从大到小的C语言算法如下：#include<stdio.h>#include<math.h>#include<iostream>using namespace std;void main(){float n=0.0;int i;int N;cout<<"Please input N"<<endl;cin>>N;for(i=N;i>1;i--){n=n+1.0/(i*i-1);N=N-1;}printf("%-100f",n);printf("\n");}执行后结果为：N=2时，运行结果为：N=3时，运行结果为：N=100时，运行结果为：N=4000时，运行结果为：2）从小到大的C语言算法如下：#include<stdio.h>#include<math.h>#include<iostream>using namespace std;void main(){float n=0.0;int i;int N;cout<<"Please input N"<<endl;cin>>N;for(i=2;i<=N;i++){n=n+1.0/(i*i-1);}printf("%-100f",n);printf("\n");}执行后结果为：N=2时，运行结果为：N=3时，运行结果为：N=100时，运行结果为：N=4000时，运行结果为：结论：通过比较可知：N 的值较小时两种算法的运算结果相差不大，但随着N 的逐渐增大，两种算法的运行结果相差越来越大。
精选计算实习报告（精选15篇）
精选计算实习报告（精选15篇）精选计算篇1专业：07计算机科学与技术学号：姓名：我到桂林市防雷中心参加社会实践活动，到8月30日实践结束。
在单位领导及同事的指导和帮助下，我慢慢了解了单位的组织机构及管理体制。
并学到了很多学校里学不到的实用知识，丰富了社会阅历。
单位给我的工作定位是从事网络基础工作以及一些计算机组装和日常维护等等，在此思想的指导下，我承担了单位有关计算机方面的一些基础工作，包括线路检修，交换设备更换，新用户入户等等这方面工作并很快掌握了其中的原理和方法，从学校的理论知识到过去真正的实践工作，其中的角色转变离不开单位各级领导以及同事们的帮助和指导。
网站建设方面，对此进行了改造，使用全新的后台，使单位网站的功能更加的强大，管理更加的有条理和方便。
取得了不错的社会效应。
我积极参加了单位组织的各种活动，包括乒乓球、羽毛球、篮球比赛。
在实习期间，我深刻的感受到，单位是一个团结、上进、充满活力的集体，每天大家都是笑脸相迎，即使面临很大的压力，办公室里仍然会听到笑声，面对同事，大家总是热情真诚，面对工作上的困难，大家总是互相帮助，直至解决困难。
整个集体和睦相处，气氛很和谐、很亲切。
领导与下属之间、同事与同事之间非常的团结、凝聚，气氛非常的和谐，就像一个温馨的大家庭，而领导就像这个家庭的家长，给每个人很大的空间自由发挥。
从他们身上，我真正体会到什么叫爱岗敬业，体会到了事业单位的严格管理制度，知道了什么求真务实。
特别令我感动的是，每当我遇到困难向大家求助时，谁都会无私的告诉我。
我对能到这样的单位实习感到骄傲、感到自豪。
我很庆幸自己能在这样有限的时间里，在这么和谐的气氛中工作，学习，和同事们一起分享快乐，分担工作，所以我努力向同事学习，不懂就问，认真完成领导和同事交给我的每一项工作，领导和同事也都尽力帮助我，给我讲授和工作有关的知识，耐心解答我的疑难困惑，并给我制订了一系列的，帮助我达到实习的目的。
在这一段时间里我不仅很好的运用了所学的专业知识，而且还学到了很多在学校里学不到的实用的待人处世之道，扩大了知识面，也丰富了社会实践经历。
计算方法与实习上机实验报告
计算方法与实习上机实验报告一、引言本文旨在介绍和展示我们在“计算方法”课程中的实习上机实验环节所完成的一些关键任务和所取得的成果。
该实验课程的目标是让我们更深入地理解和应用各种计算方法，并在实际操作中提高我们的编程和问题解决能力。
二、实验内容与目标实验的主要内容是利用各种计算方法解决实际数学问题。
我们被要求使用编程语言（如Python或Java）来实现和解决这些问题。
这些问题包括使用牛顿法求解平方根，使用蒙特卡洛方法计算圆周率，以及使用最优化方法求解函数的最小值等。
实验的目标不仅是让我们掌握计算方法的基本理论，更是要让我们能够在实际操作中运用这些方法。
我们需要在实习过程中，通过与同伴们合作，共同解决问题，提高我们的团队合作能力和问题解决能力。
三、实验过程与问题解决策略在实验过程中，我们遇到了许多问题，如编程错误、理解困难和时间压力等。
我们通过相互讨论、查阅资料和寻求教师帮助等方式，成功地解决了这些问题。
例如，在实现牛顿法求解平方根时，我们一开始对导数的计算和理解出现了一些错误。
但我们通过查阅相关资料和讨论，最终理解了导数的正确计算方法，并成功地实现了牛顿法。
四、实验结果与结论通过这次实习上机实验，我们不仅深入理解了计算方法的基本理论，还在实际操作中提高了我们的编程和问题解决能力。
我们的成果包括编写出了能有效求解平方根、计算圆周率和求解函数最小值的程序。
这次实习上机实验非常成功。
我们的团队不仅在理论学习和实践操作上取得了显著的进步，还在团队合作和问题解决方面积累了宝贵的经验。
这次实验使我们对计算方法有了更深的理解和认识，也提高了我们的编程技能和解决问题的能力。
五、反思与展望回顾这次实验，我们意识到在实验过程中，我们需要更好地管理我们的时间和压力。
在解决问题时，我们需要更有效地利用我们的知识和资源。
在未来，我们希望能够更加熟练地运用计算方法，并能够更有效地解决问题。
我们也希望能够将所学的计算方法应用到更广泛的领域中，如数据分析、科学研究和工业生产等。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
论文报告 PDF 文件使用 "pdfFactory Pro" 试用版本创建
2
计算方法
实验一 舍入误差与数值稳定性
实习题 1：用 2 种不同的顺序计算∑ 差的变化。 程序一（顺序） ：
#include<iostream.h> #include<math.h> #include<iomanip.h> void main() { double y1=1,y2; int n=1; cout<<"y[1]="<<setw(8)<<y1; while(1) { y2=y1+1/pow(n+1,2); cout<<"y["<<n+1<<"]="<<setw(8)<<y2; if(n>=9999) break; y1=y2; n++; if(n%3==0) cout<<'\n'; } }
9
计算方法
float Bisection(float a,float b,float(*f)(float)) { float c,fc,fa=(*f)(a),fb=(*f)(b); int n=1; cout<<setw(10)<<"二分次数"<<setw(16)<<"c"<<setw(16)<<"fc"<<'\n'; while(1) { if(fa*fb>0) { cout<<"不能用二分法求解!"; break; } c=(a+b)/2,fc=(*f)(c); if(fabs(fc)<delta) break; else if(fa*fc<0) { b=c,fb=fc; } else{ a=c,fa=fc; } if(b-a<eps) break; cout<<setw(10)<<n++<<setw(16)<<c<<setw(16)<<fc<<'\n'; } return c; } float f(float x) { return x*x*x+x*x-3*x-3; } void main() { float a=1,b=2; float x; x=Bisection(a,b,f); cout<<"方程的根为"<<x<<'\n'; }
)。
，并指出
有效位数。 1） 程序：
#include<iostream.h> #include<math.h> #include<iomanip.h> void main() { double y1,y2,b=1/3.0; int n,N; cout<<"N="; cin>>N; n=N; y2=(1.5-1/N-1/(N+1))/2; cout<<"y["<<N<<"]="<<setw(15)<<y2; while(1) { y1=y2-1/(pow(n,2)-1); cout<<"y["<<n-1<<"]="<<setw(15)<<y1; if(n<=3) break; y2=y1; n--; if(n%3==0) cout<<'\n'; } cout<<endl<<"其精确值为"<<b<<endl; cout<<"误差大约为"<<fabs(y1-b)<<endl; }
实验二 方程求根
实习题 1：求方程 f(x)=x + x − 3x − 3 = 0在 1.5 附近的根。 程序：
#include<iostream.h> #include<math.h> #include<iomanip.h> #define eps 5e-6 #define delta 1e-6
论文报告 PDF 文件使用 "pdfFactory Pro" 试用版本创建
运行结果：
论文报告 PDF 文件使用 "pdfFactory Pro" 试用版本创建
4
计算方法
实习题 2：设S = ∑ 1） 2） 3）
编制按从大到小的顺序计算S 的程序； 编制按从小到大的顺序计算S 的程序； 按 2 种顺序分别计算S ，S ，S
，已知其精确值为 ( − −
3）
N=1000，顺序 运行结果：
论文报告 PDF 文件使用 "pdfFactory Pro" 试用版本创建
6
计算方法
4）
N=1000，倒序 运行结果：
N=10000，顺序 运行结果：
论文报告 PDF 文件使用 "pdfFactory Pro" 试用版本创建
n
≈ 1.644834，分析其误
运行结果：
论文报告 PDF 文件使用 "pdfFactory Pro" 试用版本创建
3
计算方法
程序二（倒序） ：
#include<iostream.h> #include<math.h> #include<iomanip.h> void main() { double y1,y2=1.644834; int n=10000; cout<<"y[10000]="<<setw(8)<<y2; while(1) { y1=y2-1/pow(n,2); cout<<"y["<<n-1<<"]="<<setw(8)<<y2; if(n<=2) break; y2=y1; n--; if(n%3==0) cout<<'\n'; } }
计算方法
东 南 大 学
计算方法与实习
上机实验报告
学号：02607113 姓名：周金波 指导老师：李翠平 实验时间：2009.12.26
论文报告 PDF 文件使用 "pdfFactory Pro" 试用版本创建
计算方法
目录
1 舍入误差与数值稳定性 ................................................................................................... 3 1.1 舍入误差与数值稳定性 ............................................................................................ 3 1.2 实习题 ...................................................................................................................... 3 2 方程求根 .......................................................................................................................... 9 2.1 二分法 ...................................................................................................................... 9 2.2 牛顿迭代法 ............................................................................................................ 11 3 线性方程组数值解法 ..................................................................................................... 17 3.1 列主元高斯消去法 ................................................................................................. 17 3.2 追赶法 .................................................................................................................... 19 3.3 迭代法 ............................................................................................................................... 21 3.3.1 雅可比迭代法 .......................................................................................................... 21 3.3.2 高斯-赛德尔迭代法 ....................................................................................... 24 4 插值法 ............................................................................................................................ 26 4.1 拉格朗日插值多项式 ............................................................................................. 26 4.2 牛顿插值多项式 .................................................................................................... 28 5 曲线拟合 ........................................................................................................................ 29 5.1 最小二乘法............................................................................................................. 29 6 数值积分 ........................................................................................................................ 33 6.1 复化梯形公式与复化辛卜生公式的自适应算法 ................................................... 33 6.1.1 复化辛卜生公式 ............................................................................................ 33 6.1.2 自适应梯形公式 ............................................................................................ 34