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

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

return(z);

}

//将矩阵A压缩存储在矩阵C中

void YaSuoJZ(double C[5][501],double a[501],double b,double c)

{

int i;

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

{

if(i>=2) C[0][i]=c;

else C[0][i]=0;

if(i>=1) C[1][i]=b;

else C[1][i]=0;

if(i<=499) C[3][i]=b;

else C[3][i]=0;

if(i<=498) C[4][i]=c;

else C[4][i]=0;

C[2][i]=a[i];

}

}

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

double mifa(double C[5][501])

{

double u[501];

double y[501]={0},η=0;

double β,βk=0;

double ε=1; // ε为精度

double sumu=0,sumAY=0;

int i,j,k=1; //k为循环次数

for (i=0;i<=500;i++) //取任一非零向量u0

u[i] = 1.0;

while(ε>=1e-12)

{

for(i=0;i<=500;i++) //求u(k-1)的2范数η

sumu=sumu+u[i]*u[i];

η=sqrt(sumu);

sumu=0;

for(i=0;i<=500;i++) //求y(k-1)

y[i]=u[i]/η;

for(i=0;i<=500;i++) //求u(k)的各分量u[i]

{

for(j=max2(0,i-2);j<=min(i+2,500);j++)

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