數(shù)值代數(shù)上機(jī)報(bào)告(文件)
【正文】 ][1]0) b[i][j][1]=a[i][j][1]。k++) { for( max=0,i=1。j++) { if(ijamp。 d=j。 if(u=0) x=2*b[c][d][k]。 f=*x/t*(1/g)。j=n。amp。j==c) r[c][c][k]=g。 if(i==camp。amp。i!=camp。 for(i=1。j++) { a[i][c][k+1]=b[i][c][k]*g+b[i][d][k]*f。 if(i!=camp。i!=damp。 v=b[c][d][k]*(g*gf*f)。 a[d][d][k+1]=b[c][c][k]*f*f2*b[c][d][k]*f*g+b[d][d][k]*g*g。j=n。 for(i=1。j++) { if(i!=j) w=w+b[i][j][k+1]*b[i][j][k+1]。 else continue。三、 源程序：includeincludefloat solution(float a,float b,float c,float d)。 scanf(%f%f%f%f,amp。d)。 f1=3*a*x0*x0+2*b*x0+c。 //限制條件10的負(fù)六次方 return x。 x0=x。}float solution(float a,float b,float c,float d) //定義求解函數(shù){ float x,x0=1,f,f1。b,amp。 printf(\n函數(shù)f(x)=a*x*x*x+b*x*x+c*x+d\n\n)。逐步線性化是構(gòu)造迭代函數(shù)最常用的方法，Newton迭代法就是用這種方法構(gòu)造而成的。 } printf(x=%lf,w)。i++) for(j=1。 else b[i][j][k+1]=a[i][j][k+1]。i=n。 a[d][c][k+1]=h+v。j!=c) a[i][j][k+1]=b[i][j][k]。j!=damp。 a[i][d][k+1]=b[i][c][k]*f+b[i][d][k]*g。i++) for(j=1。j!=d) r[i][j][k]=1。 if(i==jamp。j==d) r[c][d][k]=f。amp。 if(i==camp。amp。i=n。 t=sqrt(x*x+y*y)。 if(u=0) y=u。b[i][j][k]max+) { max=b[i][j][k]。i++) for(j=1。 } for(k=1。i++) for(j=1。a[i][j][1])。i=n。 scanf(%d,amp。二、 算法：1) 輸入系數(shù)矩陣 A；2) 取，形成一個(gè)相似矩陣序列，k=1,2,…3) 選取，使得元素=0，這時(shí)應(yīng)滿足 4) 令，因此可以得到： ，三、 源程序：include includedefine P void main(){ double a[20][20][20],b[20][20][20],r[20][20][20],l[20][20][20]。i=n。 printf(最接近的特征值為：\n)。i++) { y0[i]=y1[i]。 if(s==n)flag=0。 /*判斷：若|Z[k]Z[k1]|=e，則進(jìn)行下一輪迭代 ；否則得解。i++) if(fabs(y1[i])fabs(y1[i1]))m=y1[i]。k++) q=q+r[i][k]*y1[k]。 } for(i=n。i++) { for(q=0,k=1。p++) q=q+l[i][p]*r[p][k]。 } for(i=k+1。i++) { for(q=0,p=1。 while(flag==1) { /*用Doolittle解方程 LR*y2=z1中的y2，其中B=LR */ for(k=1。 for(i=1。 y1[i]=(y[i]q)/r[i][i]。i=1。ki。 l[i][k]=(R[i][k]q)/r[k][k]。i=n。pk。k=n。 /*第一步迭代，用Doolittle計(jì)算R*y1=y0中的y1*/ for(i=1。i=n。j=n。 L[i][k]=(B[i][k]q)/R[k][k]。i=n。pk。k=n。j=n。lm)。j++) scanf(%lf,amp。 for(i=1。 printf(請(qǐng)輸入矩陣的階數(shù)：n=)。三、 源程序：includeincludedefine e define N 100void main(){ int i,j,k,n,p,s,kk=1,flag=1。 //return 0。 for(i=0。jn。 for(j=0。 } xmax=。j++) { nx[j]=0。 for(i=0。 for(i=0。i++) for(j=0。 return 1。n)。
研究報(bào)告相關(guān)推薦
鄂ICP備17016276號(hào)-1