北航數(shù)值分析計(jì)算實(shí)習(xí)報(bào)告一-資料下載頁(yè)
2025-07-22 01:55本頁(yè)面
【正文】 wucha // wucha=(fabs(B_k1B_k0))/(fabs(B_k1))。//get wucha if(wucha=Precision) break。 else if(i10000) { printf(迭代次數(shù)超長(zhǎng)，請(qǐng)更改初始向量\n)。 break。 } else B_k0 = B_k1。 } *iterations = i。 return (1/B_k1)。}void Doolittle_getx(double C[5][501],double y[501],double u[501],int n,int s,int r){ int i,t。 double e。 for(i = 2。i = n。i++) { e = 0。 for(t = max(1,ir)。t = i1。t++) { e+=C[it+s+11][t1]*y[t1]。 } y[i1] = e。 } u[n1] = y[n1]/C[s+11][n1]。 for(i = n1。i=1。i) { e = 0。 for(t = i+1。t=min(i+s,n)。t++) e += C[it+s+11][t1]*u[t1]。 u[i1] = (y[i1] e)/C[s+11][i1]。 }}void copy(double b[501],double y[501]){ for(int i=0。i501。i++) { b[i] = y[i]。 }}void DoolittleC(double C[5][501],int n,int s,int r){ int k,j,i,t。 double e。 for(k = 1。k = n。k++) { for(j = k。j = min(k+s,n)。j++) { e = 0。 for(t = max(max(1,kr),js)。t = k1。t++) { e = e + C[kt+s+11][t1]*C[tj+s+11][j1]。 } C[kj+s+11][j1] = C[kj+s+11][j1] e。 } if(k == n) break。 for(i = k+1。i = min(k+r,n)。i++) { e = 0。 for(t = max(max(1,ir),ks)。t = k1。t++) { e = e + C[it+s+11][t1]*C[tk+s+11][k1]。 } C[ik+s+11][k1] = (C[ik+s+11][k1] e)/C[s+11][k1]。 } } }int min(int a,int b){ if(ab) return a。 else return b。}int max(int a,int b){ if(a=b) return a。 else return b。} 附錄：部分實(shí)驗(yàn)程序截圖 另，其余都為0 ，其余為0 ，其余為01其余為0
高考資料相關(guān)推薦
鄂ICP備17016276號(hào)-1