北航数值分析实习题目第一题(4)

北航数值分析实习题目第一题

sumAY=sumAY+C[i-j+2][j]*y[j];

u[i]=sumAY;

sumAY=0;

}

//求幂法中的βk

β=βk; //将β(k-1)放在β中

βk=0;

for(i=0;i<=500;i++) //求βk

βk=βk+y[i]*u[i];

if(k>=2)

ε=fabs(βk-β)/fabs(βk);

k++;

}

return(βk);

}

//带状矩阵的LU分解

void daizhuangLU(double A[5][501])

{

int i,j,k,m,t;

double sumukj=0,sumlik=0;

for(k=0;k<=500;k++)

{

for(j=k;j<=min(k+2,500);j++) //求ukj并存在A[k-j+2][j]中 {

for(t=max3(0,k-2,j-2);t<=k-1;t++)

sumukj=sumukj+A[k-t+2][t]*A[t-j+2][j];

A[k-j+2][j]=A[k-j+2][j]-sumukj;

sumukj=0;

}

if(k<500)

for(i=k+1;i<=min(k+2,500);i++) //求lik并存在A[i-k+2][k]中 {

for(m=max3(0,i-2,k-2);m<=k-1;m++)

sumlik=sumlik+A[i-m+2][m]*A[m-k+2][k];

A[i-k+2][k]=(A[i-k+2][k]-sumlik)/A[2][k];

sumlik=0;

}

}

}

//反幂法函数：用反幂法求矩阵的模最小的特征值

double fanmifa(double M[5][501])

{

搜索“diyifanwen.net”或“第一范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读，第一范文网，提供最新外语学习北航数值分析实习题目第一题(4)全文阅读和word下载服务。