數(shù)值分析上機(jī)實(shí)習(xí)報(bào)告(參考版)
2025-01-21 21:49本頁面
【正文】 。（5）問題討論矩陣A是嚴(yán)格對(duì)角占優(yōu)陣，而且還是實(shí)數(shù)對(duì)稱陣，滿足Householder變換的條件，并且得到的B就是與A相似的三對(duì)角陣。x=mgauss2(A,b)。,]。,。,。,。endformat longA=[,。x(n)=b(n)/A(n,n)。 A(k+1:n,k)=zeros(nk,1)。 A(k+1:n,k+1:n)=A(k+1:n,k+1:n)m*A(k,k+1:n)。 b([k p],:)=b([p k],:)。 p=p+k1。endn=length(b)。j++) { printf(%f,A[i][j])。 for(j=0。i9。 } } } printf(\nHousehold變換矩陣B為:\n)。j9。i9。i++) { q[i]=y[i]k*u[i]。 for(i=0。i++) { k=k+u[i]*y[i]。 for(i=0。 } y[i]=y[i]/l。j9。i++) { y[i]=0。 } for(i=0。 else if(i==r+1) u[i]=A[r+1][r]+sign*s。i9。 else if(A[r+1][r]0) sign=1。 l=s*s+fabs(A[r+1][r])*s。i++) { s=s+A[i][r]*A[i][r]。 for(i=r+1。r7。j++) { A0[i][j]=A[i][j]。i++) { for(j=0。 for(i=0。 double A[9][9]= { {,}, {,}, {,}, {,}, {,}, {,}, {,}, {,}, {,} }。s1 , 0, … ,0)T,則d1=H1a1 , 這里為了增加計(jì)算的靈活性，避免同號(hào)數(shù)相減，選取符號(hào)使2a21s10,這樣只需要把上式右端第二項(xiàng)取為2sign(a21)a21 s1即可，從而歸納起來算法步驟為：（2）程序主體include include void main( ){ int i,j,m,r,sign。s1 , 0 , 0 , … , 0）T,這里,前面的正負(fù)表示鏡像可取為兩個(gè)方向相反的向量，這增加計(jì)算的靈活性。 第五題 第 20 頁 A=b=( )T (1).用列主元素消去法求解Bx=b. （1）理論依據(jù)用Househloder變換，把A化為三對(duì)角陣。end，精度高，可根據(jù)有變化的情況與需要的精度自動(dòng)修改步長，誤差小且程序簡單，存儲(chǔ)量少。Z(2)=0*x+0*y1+0*y2+1*y3。endfunction Z=df1(x,y1,y2,y3)Z=zeros(3,1)。 K(:,4)=feval(df1,x,y1,y2,y3)。y2=S(2,1)。S(:,1)=Y+h.*K(:,3)。 K(:,3)=feval(df1,x,y1,y2,y3)。y2=S(2,1)。 S(:,1)=Y+*h.*K(:,2)。y3=S(3,1)。 y1=S(1,1)。 x=x+*h。y3=a。y1=a。K=zeros(3,4)。Y=zeros(3,1)。 } if(i==200) printf(\ny1()=%f y2()=%f y3()=%f\n\n,y1,y2,y3)。 } if(i==170) { printf(\ny1()=%f y2()=%f y3()=%f,y1,y2,y3)。 } if(i==90) { printf(\ny1()=%f y2()=%f y3()=%f,y1,y2,y3)。 if(i==50) { printf(\ny1()=%f y2()=%f y3()=%f,y1,y2,y3)。 k4=h*(10001000*y2100*(y3+k3))。 k2=h*(10001000*y2100*(y3+*k1))。 y2+=(k1+2*k2+2*k3+k4)/6。 y1+=(k1+2*k2+2*k3+k4)/6。i=200。 double y1=,y2=,y3=。 double h=。 第四題,打印, Ru
教學(xué)課件相關(guān)推薦
鄂ICP備17016276號(hào)-1