【正文】 =M。j++) B[i]+=BF[i][j]*Y[j]。jM+1。 數(shù)值計算課程設計 29 //求 F的轉(zhuǎn)置 for(i=0。i++) cinY[i]。 cinn。 for(j=0。 for(j=0。 for(i=0。go39。 coutendl。 else cout t*(x x[i])^3。 i n。 三次樣條插值算法 24 }while(ch == 39。 i++) { a[i] = h[i1] / (h[i] + h[i1])。 break。 fxym[0] = 6*((y[1] y[0]) / (x[1] x[0]) f0) / h[0]。 float f0, f1。 coutPlease put in Yi39。 char ch。 for(i = 1。 float b = (y[x2] y[x1]) / (x[x2] x[x1])。S239。 y3=*(x3 3).^3 + *(x3 2).^3 *(x3 3) + *(x3 2)。 cinab。 龍貝格求積分算法 18 } p = fabs(q y[m1])。 } p = (y[0] + h*p)/。 ep = epsilon + 。 double s, q。iN。iN。i++) { for(j=0。 牛頓法解非線性方程組 14 coutendl。i) for(j=N21。iN。i0。 } } for (i=0。 for (i=0。 else aug[i][j]=0。 for (i=0。i++) {for(j=0。 y=xx[1]。 for( i=0。 return 0。 for (i=0。 coutendl。 void inv_jacobian(float yy[N][N],float inv[N][N])。}數(shù)值計算課程設計 9 1122( , ) ( , )()( , ) ( , )k k k kkk k k kf p q f p qxyJPf p q f p qxy??????? ????12( , )() ( , )kkk kkf p qFP f p q??? ????1kkP P P? ? ?? 牛頓法解非線性方程組 、 算法說明 設 kP 已知。 //system(pause)。 for(j=i+1。k++) { l[k][i]=a[k][i]/a[i][i]。 a[i][n]=a[c][n]。 for(i=0。 int main() { int i,j。 |*(head+(m+1)*k+i)|max K=k+1。 數(shù)值計算課程設計 5 高斯列主元法解線性方程組 、 算法說明 首先將線性方程組做成增光矩陣，對增廣矩陣進行行變換。 for(i=0。 cinab。} return(result)。} else { double x1,y1。 //f 為函數(shù)的入口地址， x0、 y0為初值， xn為所求點， step 為計算次數(shù) double Runge_Kuta( double (*f)(double x, double y), double x0, double y0, double xn, int step ) { double k1,k2,k3,k4,result。該算法是構(gòu)建在數(shù)學支持的基礎之上的。 k3=f(x0+h/2, y0+h*k2/2)。 k3=f(x1+h/2, y1+h*k2/2)。 cout請輸入初值 x0,y0:。 cinstep。 } double f(double x, double y) { double r。 、 高斯列主元算法流程圖 Y N max i!=i temp=*(head*(m+1)*i+k)。 const N=20。 coutendl。 /*找列最大元素 */ 高斯列主元法解線性方程組 8 for(n=0。 coutendl。 for(i=m2。 for(i=0。i++) for(j=0。 const int N2=2*N。 for (i=0。 //jis uan jacobian ju zhen de ni juzhen invjacobian inv_jacobian(jacobian,invjacobian)。iN。 //非線性方程組 yy[0]=x*x2*xy+。 } void ffjacobian(float xx[N],float yy[N][N]) { float x,y。 cout雅克比矩陣： endl。 } void inv_jacobian(float yy[N][N],float inv[N][N]) {float aug[N][N2],L。 for(j=N。j++) coutaug[i][j] 。 for(j=i。 coutendl。j=0。 } coutendl。i++) {for(j=0。 cout雅克比矩陣的逆： endl。 } void newdim(float x0[N], float inv[N][N],float y0[N],float x1[N]) { int i,j。 x1[i]=x0[i]sum。 、 龍貝格求積分算法代碼 includeiostream includecmath using namespace std。 y[0] = h*(f(aa) + f(bb))/。 in。 q = (s*p y[k1])/(s )。 } return (q)。 x2=1::2。) gtext(39。 float x[max], y[max], h[max]。 i n。 i) fxym[i] = fxym[i] B[i]*fxym[i+1]。 i++) { coutPlease put in Xi39。 i++) //求步長 h[i] = x[i+1] x[i]。 cinf0f1。 c[0] = a[n] = 0。 i++) fxym[i] = 6 * f(i1, i, i+1)。 printout(n)。 return 0。 else coutt*(x x[i+1])^3。 t = (y[i+1] fxym[i+1]*h[i]*h[i]/6)/h[i]。 x=3::3。) 、 M 次多項式曲線 擬合算法代碼 includeiostream includecmath define MAX 20 using namespace std。 } for(i=0。kn。 double X[MAX]={0},Y[MAX]={0},F[MAX][MAX]={0},B[MAX]={0}。 cout\n請輸入 n個點的 Y坐標序列 :\n。i++) for(k=1。 } } //計算其轉(zhuǎn)置的 BF 與 F 的乘 for(i=0。iM+1。jn。 for(i=M。 includecmath double f(double x) { return x*x2*x1 } double df(double x) { return 2*x2。k++) { p1=p0f(p0)/df(p0)。 cout迭代次數(shù)為 :kendl。 圖 81 程序調(diào)試圖 不動點法解非線性方程 34 cout不動點法解非線性方程 f(x)=1x^2/2endl。 err=fabs(P[k1]P[k2])。 cout求解不動點近似值的序列 :endl。 int main() { double Sn=0。 for(i=0。i++) { Sn=Sn+func(X,i,x,n)*y[i]。i++) { if(i==k) continue。 double err,norm,A[n][n],B[n],P[n],X[n]。 for (i=0。kmax1。 } norm=0。i++) P[i]=X[i]。 我今天認真的進行課程設計，學會腳踏實地邁開這一步，就是為明天能穩(wěn)健地在社會大潮中奔跑打下堅實的基礎 。 參考文獻 40 參 考 文 獻 [1]. John H. Mathews. Peking: Numerical Methods Using MATLAB Fourth Edition. [M].電子工業(yè)出版社 . [2]. 譚忠強 . C++程序設計 .[M].清華大學出版社 . 2021. [3]. C++程序設計題解與上機指導 . [M] . 清華大學出版社 .2021. [4]. John 等著 .數(shù)值方法 .（ MATLAB 版）（第四版） .[M]. 電子工業(yè)出版社 .2