北航计算机实习第三次作业

《数值分析 B》计算实习第三题

算法的设计方案

1. 使 用 牛顿迭代 法 ( 简 单迭 代 法 不 收 敛 ) ， 对 原 题 中 给 出 的 ，

，( )的 11*21 组 分别求出原题中方程

组的一组解，于是得到一组和 对应的 。

2.对于已求出的 ，使用分片二次代数插值法对原题中关于 的数表进

行插值得到 。于是产生了 z=f(x,y)的 11*21 个数值解。

3.从 k=1 开始逐渐增大 k 的值，并使用最小二乘法曲面拟合法对 z=f(x,y)进

#define N 4

double e=1e-12;

double gauss(double a[N][N+1],double x[N]) //gauss 求线性方程组

{

int i,j,k;

double p=0;

for(k=0;k<N-1;k++)

{

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

{

}

}

if(a[N-1][N-1]==0)

{

printf("a[%d][%d]=0\n",N-1,N-1);

}

x[N-1]=a[N-1][N]/a[N-1][N-1];

}

x[i]=(a[i][N]-p)/a[i][i];

}

return 0;

}

void A_ni(double A[10][10],int n) //求逆矩阵

{

double a[10][10],b[10][20],c[10][10],t;

int i,j,m;

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

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

for(j=n;j<2*n;j++)

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

b[i][n+i]=1;

for(m=0;m<n;m++)

{

t=b[m][m];

i=m;

while(b[m][m]==0)

for(j=2*n-1;j>=m;j--)

}

m=n-1;

{

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

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

for(j=0;j<n;j++)

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

{

for(j=0;j<n;j++)

}

while(1)

{

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

A[0][0]=-0.5*sin(result[0]);

A[0][1]=1;

A[0][2]=1;

A[0][3]=1;

A[1][0]=1;

A[1][1]=0.5*cos(result[1]);

A[3][0]=1;

A[3][1]=0.5;

A[3][2]=1;

A[3][3]=cos(result[3]);

A[0][4]=-1*(0.5*cos(result[0])+result[1]+result[2]+result[3]-x-2.67);

A[2][4]=-1*(0.5*result[0]+result[1]+cos(result[2])+result[3]-x-3.74);

A[3][4]=-1*(result[0]+0.5*result[1]+result[2]+sin(result[3])-y-0.79);

gauss(A,result);

num1=num2=0;

num1=sqrt(num1);