北航數(shù)值分析大作業(yè)一-預(yù)覽頁(yè)
【正文】 j=2。 array[4][k]=。) { k=500j。 }//第2行賦值 for(j=0。 array[2][k]=(*j)*sin((double)(*j))*exp((double)())。 /*p為平移量*/ for(。i=500。i=500。i++) { u[i]=0。 } u[i]=u[i]p*y[i]。 } if(fabs((bc)/b)=E) /*達(dá)到精度水平，迭代終止*/ break。 srand((int)time(0))。 /*生成0~10的隨機(jī)數(shù)*/ }} void chuzhi2(double a[],int j) /*令初始迭代向量為ei*/{ int i。 } a[j]=1。k++) { for(j=k。t++) { array[kj+2][j]=array[kt+2][t]*array[tj+2][j]。t=(k1)。 double a,b,c=0。i++) { array[2][i]=p。 b=0。 } a=sqrt(a)。 } //回帶過(guò)程，求解uk for(i=0。i=500。 } } u[500]=u[500]/array[2][500]。j=min2((i+2),500)。i=500。 } return (p+(1/b))。 double MatrixC[5][501]。//用帶原點(diǎn)平移的冪法求解λ1，λ501 chuzhi(u)。 if(d0) { d1=d。 } printf(λ1=%.12e\n,d1)。 printf(λs=%.12e\n,ds)。 assignment(MatrixC)。 }//求A的條件數(shù) d=fabs((d1/ds))。 a=1。 } printf(行列式detA=%.12e\n,a)。 for(i=0。 printf(u%03d,λmax=%+e ,i,d1)。}三. 程序結(jié)果：四．分析初始向量選擇對(duì)計(jì)算結(jié)果的影響 矩陣的初始向量選擇，對(duì)結(jié)果的影響很大，選擇不同的初始向量可能會(huì)得到的特征值。這是因?yàn)橐訟的n個(gè)線性無(wú)關(guān)的特征向量為一組基，將初始向量線性表出時(shí)，按摸最大的特征值λ1對(duì)應(yīng)的特征向量x1的系數(shù)α1如果為0，就無(wú)法求出特征值
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1