数值计算方法上机实习报告

实习报告。

课程名称：数值计算方法上机。

学院：电气信息工程学院。

姓名：李勇。

学号： 3100501063

2011年12月30日。

数值计算方法上机习题。

1.用newton迭代法求解方程在2.0附近的实根，。

源程序：#include<>

#include<>

main()

double x,f,f1;

scanf("%lf",&x);

for(;;x=x-f/f1)

printf("r=%lf",x);

getch ()

return 0;

输入：2.0

运行结果：r=1.8793852. 设。

试用高斯－约旦消去法求方程组的解。

1）源程序：

#include<>

main()

float a[5][5]=,b[5]=,s,t,am;

int i,j,l,k;

clrscr();

for(k=1;k<=4;k++)

if(a[l][k]==0)

printf("qi yi ju zhen.")break; }

if(l!=k)

for(j=k;j<=4;j++)

t=a[k][j];a[k][j]=a[l][j];a[l][j]=t; }

s=b[l]; b[l]=b[k]; b[k]=s;}

for(j=k+1;j<=4;j++)

a[k][j]=a[k][j]/a[k][k];

b[k]=b[k]/a[k][k];

for(i=1;i<=4;i++)

printf("x1=%f",b[1]);

printf("x2=%f",b[2]);

printf("x3=%f",b[3]);

printf("x4=%f",b[4]);

运行结果：x1=19.999809

x2=-11.999882

x3=-4.999953

x4=2.999972

2）源程序：

#include<>

main()

float a[4][4]=,b[4]=,s,t,am;

int i,j,l,k;

clrscr();

for(k=1;k<=3;k++)

if(a[l][k]==0)

printf("qi yi ju zhen.")break; }

if(l!=k)

s=b[l]; b[l]=b[k]; b[k]=s;}

for(j=k+1;j<=3;j++)

a[k][j]=a[k][j]/a[k][k];

b[k]=b[k]/a[k][k];

for(i=1;i<=3;i++)

printf("x1=%f",b[1]);

printf("x2=%f",b[2]);

printf("x3=%f",b[3]);

运行结果：x1=-11.000001

x2=8.000000

x3=-3.000000

3.已知函数的数表。

试用拉格朗日插值法求时的函数值。

源程序：#include<>

void main()

float x[4]=,y[4]=;

float u,v,l;

int i,j; clrscr();

scanf("%f",&u);

v=0;for(j=0;j<=3;j++)

v=v+l*y[j];

printf("u=%fv=%f",u,v);

运行结果：u=0.332400

v=0.326313

4.用复化simpson方法求积分，要求误差不超过。

源程序：#include<>

#include<>

void main()

int i,n=2;

float s;

float f(float);

for(i=0;i<=2;i++)

s=simpson(f,0,1,n);

printf("s(%d)=%f",n,s);

n*=2;getch();

float simpson(float (*f)(float),float a,float b,int n)

int k;

float s,s1,s2=0.0;

float h=(b-a)/n;

s1=f(a+h/2);

for(k=1;k<=n-1;k++)

s1+=f(a+k*h+h/2);

s2+=f(a+k*h);}

s=h/6*(f(a)+4*s1+2*s2+f(b));

return s;

float f(float x)

return 5/(1+x*x);

运行结果：s(2)=3.926961

s(4)=3.926991

s(8)=3.926991

5. 用阿达姆斯**－校正法求解。

的初值问题。

分别取步长和，分析其结果，保留5位有效数字。

源程序：#include<>

#include<>

double f(double x)

return cos(sin(x));

main()

double h,x,y=0,y2;

clrscr();

scanf("%lf",&h);

printf(“h=%f”,h);

for(x=0;x<=1.001;x=x+h)

运行结果：h=0.200000

x=0.000000 y=0.198720

x=0.200000 y=0.389827

x=0.400000 y=0.567151

x=0.600000 y=0.727046

x=0.800000 y=0.868845

x=1.000000 y=0.994710

运行结果：h=0.050000

x=0.000000 y=0.049979

x=0.050000 y=0.099834

x=0.100000 y=0.149441

x=0.150000 y=0.198680

x=0.200000 y=0.247436

x=0.250000 y=0.295600

x=0.300000 y=0.343069

x=0.350000 y=0.389749

x=0.400000 y=0.435555

x=0.450000 y=0.480414

x=0.500000 y=0.524262

x=0.550000 y=0.567047

x=0.600000 y=0.608729

x=0.650000 y=0.649279

x=0.700000 y=0.688682

x=0.750000 y=0.726934

x=0.800000 y=0.764041

x=0.850000 y=0.800023

x=0.900000 y=0.834910

x=0.950000 y=0.868740

x=1.000000 y=0.901564

心得体会。为期四天的数值计算方法上机实习已经接近尾声，尽管时间很短，却让我获益良多。

首先，做什么事都不可能一蹴而就，必要的准备和持之不懈的努力，以及敢于面对失败，愈挫愈勇的精神。编程问题是非常繁琐的，要想正确的编出程序，必须要在充分了解c语言程序的基础上透彻的分析题目的的要求和目的，在进行程序的编程，然而在编程的过程中不可避免的要面临诸多的困难和阻碍，想要一下子就变出正确的程序几乎不可能，面对错误也不要失去耐心，寻根纠结，发现并改正错误，最终取得成功。其次，我们要注意要虚心求教，人并不是万能的，我们总有一些知识漏点，或者一些不能察觉的错误，这时我们就需要向别人求教，当局者迷，旁观者清，别人能过客观的知名自己的错误，让自己茅塞顿开。

数值计算方法上机实习报告

电池计算方法

疲劳分析的数值计算方法及实例部分理论可打印

储量计算方法

螺旋箍筋计算方法

党费新标准计算方法

分辨力的计算方法

姓名五格计算方法

流失率的计算方法

其他用户还读了