当前位置：首页 > 北航数值分析报告计算实习报告材料一
实用文案
int max_iterations;//最大值迭代次数
int abs_min_iterations;//求绝对值最小的迭代次数 };
int main() {
FinalValue main_num= {0,0,0,0}; double temp;//两值交换中间变量 int temp1;
double u[501]={0};//,y[501];//为u0赋初值 double Norm_u=0;//范数 double C[5][501] = {0}; InitC(C);//初始化C int i=0,*iterations; Initu(u);
iterations = &main_num.min_iterations; main_num.min = Get_Fabs_Eigenvalue(C[2],u,iterations);//将求得的绝对值最大特征值放到min变量中
main_num.abs_max = main_num.min; A_sub_minI(C[2],main_num.min); for(i=0;i<501;i++) u[i]=i+1;
iterations = &main_num.max_iterations; main_num.max = Get_Fabs_Eigenvalue(C[2],u,iterations);//将求得的绝对值最大特征值放到min变量中
main_num.max += main_num.min; if(main_num.min>main_num.max) {
temp = main_num.min;
main_num.min = main_num.max; main_num.max = temp;
temp1 = main_num.min_iterations;
main_num.min_iterations = main_num.max_iterations; main_num.max_iterations = temp1; }
printf(\最小特征值为：%.12e,迭代次数为：%d\\n\ printf(\最大特征值为：%.12e,迭代次数为：%d\\n\ /**********************************/ /*以下利用反幂法求解模最小的特征值*/ /**********************************/ InitC(C);//初始化C Initu(u);
标准文档
实用文案
DoolittleC(C,501,2,2); main_num.detA = detA(C);
iterations = &main_num.abs_min_iterations;
main_num.abs_min = Get_min_Eigenvalue(C,u,501,2,2,iterations); main_num.cond2 = fabs(main_num.abs_max/main_num.abs_min);
printf(\绝对值最小特征值为：%.12e,迭代次数为：%d\\n\ printf(\的行列式为：%.12e；\
printf(\的条件数cond(A)2为：%.12e\\n\ /**********************************************/ /*以下利用反幂法求与数列Uk中元素最相近的特征值*/ /**********************************************/ double Uk[39];//保存Uk的值并保存与Uk[i]最接近的λ double B[39]; int diedai;
iterations = &diedai; for( i=1;i<=39;i++) {
Uk[i-1] = main_num.min + i*(main_num.max - main_num.min)/40; InitC(C); Initu(u);
for(int j =0;j<501;j++) C[2][j] -= Uk[i-1]; DoolittleC(C,501,2,2);
B[i-1] = Get_min_Eigenvalue(C,u,501,2,2,iterations); B[i-1] +=Uk[i-1];
printf(\μ%-2d=%-3.12e,与μ%-2d最相近的特征=%-3.12e\\n\ }
return 0; }
double detA(double C[5][501]) {
int i =0; double e =1;
for(i =0;i<501;i++) e*=C[2][i]; return e; }
void InitMatrix(double *p) {
标准文档
值λ实用文案
for(int i=1;i<=501;i++) {
*p=(1.64-0.024*i)*sin(0.2*i)-0.64*exp(0.1/i); p++; } }
void Initu(double *p) {
for(int i=1;i<=501;i++) {
if(i<=500) *p=0; else *p=1; p++; } // *p=1; }
void A_sub_minI(double *a,double min) {
for(int i=1;i<=501;i++) {
*a=*a-min; a++; } }
double NeiJi(double *a,double *b) {
double e=0.0;
for(int i=0;i<501;i++) {
e=e+(*a)*(*b); a++; b++; }
return e; }
void get_y(double *y,double *u) { double Norm_u=sqrt(NeiJi(u,u)); for(int i=0;i<501;i++) {
标准文档
实用文案
*y = *u/Norm_u; y++; u++; } }
void get_u(double *u,double *y,double *a) {
u[0]=a[0]*y[0]+b*y[1]+c*y[2];
u[1]=b*y[0]+a[1]*y[1]+b*y[2]+c*y[3]; for(int i=2;i<499;i++)
u[i]=c*y[i-2]+b*y[i-1]+a[i]*y[i]+b*y[i+1]+c*y[i+2]; u[499]=c*y[497]+b*y[498]+a[499]*y[499]+b*y[500]; u[500]=c*y[498]+b*y[499]+a[500]*y[500]; }
void InitC(double C[5][501]) {
int i;
for(i = 2;i < 501;i++) {
C[0][i] = -0.064; C[4][i-2] = -0.064; }
for(i = 1;i < 501;i++) {
C[1][i] = 0.16; C[3][i-1] = 0.16; }
for(i = 1;i <= 501;i++) {
C[2][i-1] = (1.64-0.024*i)*sin(0.2*i)-0.64*exp(0.1/i); } }
double Get_Fabs_Eigenvalue(double *a,double *u,int *iterations) {
double y[501],B_k0=0,B_k1=0; double wucha; int i=0; while(1) {
++i;
标准文档
本文作者：...共分享92篇相关文档
