北航數(shù)值分析大作業(yè)一(存儲(chǔ)版)
【正文】 [tk+2][k]。 } //先將矩陣Doolittle分解 LU(array)。 //求yk1 for(i=0。i++) { for(j=max2(0,(i2))。j++) { u[i]=array[ij+2][j]*u[j]。 /*直接返回距離原點(diǎn)P最接近的A的特征值*/}//主函數(shù)main(){ int i。 d=mifa(u,MatrixC,0)。 printf(λ501=%.12e\n,d501)。 chuzhi(u)。 for(i=0。i=500。以冪法為例（反冪法原理相同），選取初始迭代向量ui=ei(i=0,1，…500)；即u[j]=0,j≠i；u[j]=1,j=i。 getchar()。 assignment(MatrixC)。 LU(MatrixC)。i++) { a=d1+(i*(d501d1))/40。 d1=a。//將矩陣A轉(zhuǎn)存為MatrixCassignment(MatrixC)。 c=b。i) { for(j=i+1。 } for(i=1。i++) { a=a+u[i]*u[i]。i=500。j++) { for(t=max3(0,k2,j2)。k=500。i++) { a[i]=(*rand()/RAND_MAX)。i++) { b+=y[i]*u[i]。i=500。 /*選用第一種迭代格式*///求ηk1for(i=0。 j++。j=500。 } i++。 return(t)。二． 源程序includeincludeinclude include includeincludeincludedefine E /*定義全局變量相對(duì)誤差限*/int max2(int a,int b) /*求兩個(gè)整型數(shù)最大值的子程序*/{ if(ab) return a?！稊?shù)值分析B》大作業(yè)一SY1103120 朱舜杰一． 算法設(shè)計(jì)方案：將帶狀線性矩陣A[501][501]轉(zhuǎn)存為一個(gè)矩陣MatrixC[5][501] .由于C語(yǔ)言中數(shù)組角標(biāo)都是從0開始的，所以在數(shù)組MatrixC[5][501]中檢索A的帶內(nèi)元素aij的方法是：A的帶內(nèi)元素aij=C中的元素cij+2,j
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1