數(shù)值分析 C語言 代碼

1、信息與軟件工程學(xué)院標(biāo) 準(zhǔn) 實(shí) 驗(yàn) 報(bào) 告（實(shí)驗(yàn)）課程名稱 計(jì)算數(shù)學(xué)基礎(chǔ) 電子科技大學(xué)教務(wù)處制表電 子 科 技 大 學(xué)實(shí) 驗(yàn) 報(bào) 告學(xué)生姓名：何安迪 學(xué) 號(hào)：2013220701020 指導(dǎo)教師：顧小豐實(shí)驗(yàn)地點(diǎn)：實(shí)驗(yàn)中心306 實(shí)驗(yàn)時(shí)間：2015. 11.28一、實(shí)驗(yàn)室名稱：實(shí)驗(yàn)中心306二、實(shí)驗(yàn)項(xiàng)目名稱： 計(jì)算數(shù)學(xué)（數(shù)值分析）的核心算法實(shí)現(xiàn)三、實(shí)驗(yàn)原理通過計(jì)算數(shù)學(xué)基礎(chǔ)的課堂學(xué)習(xí)，參照數(shù)值分析簡(jiǎn)明教程的內(nèi)容，熟練掌握插值方法、數(shù)值積分、方程求根的迭代法、線性方程組的迭代法、線性方程組的直接法等解題技巧。使用算法對(duì)例題進(jìn)行求解，并且通過在Visual Studio中編寫C程序?qū)崿F(xiàn)解題的算法。四、實(shí)

2、驗(yàn)?zāi)康幕谡n程的學(xué)習(xí)，使用Visual Studio編寫C程序?qū)崿F(xiàn)計(jì)算數(shù)學(xué)（數(shù)值分析）的核心算法，熟練掌握課堂上學(xué)習(xí)的各種算法，充分理解課本上的知識(shí)，并將之用于實(shí)例求解。五、實(shí)驗(yàn)內(nèi)容1、最小二乘法：在-1, 1區(qū)間上取 N = 21個(gè)等距結(jié)點(diǎn)，計(jì)算出以相應(yīng)節(jié)點(diǎn)上 ex 的值做為數(shù)據(jù)樣本，作出 m = 3，5，7，9次的最小二乘多項(xiàng)式。并計(jì)算出不同階最小二乘多項(xiàng)式給出的最小偏差六、實(shí)驗(yàn)環(huán)境OS:Microsoft Windows 7 64bit 環(huán)境： Microsoft Visual Studio 2012七、實(shí)驗(yàn)步驟1. 實(shí)驗(yàn)一：最小二乘法，同樣用到上面構(gòu)造的高斯函數(shù)，根據(jù)最小二乘法的分劃和

3、多項(xiàng)式擬合公式構(gòu)建實(shí)驗(yàn)二的函數(shù)void problem2()，代碼如下。void problem2() / 最小二乘法double x21 = 0 , y21 = 0 ;double a44 = 0 , result_a4 = 0 ;double b66 = 0 , result_b6 = 0 ;double c88 = 0 , result_c8 = 0 ;double d1010 = 0 , result_d10 = 0 ;double right10 = 0 ;double *result4 = result_a, result_b, result_c, result_d ;for (i

4、nt i = 0; i < 21; i+)xi = -1 + i * 2 / 20.0;yi = exp(xi);for (int i = 0; i < 10; i+)for (int j = 0; j < 10; j+)for (int k = 0; k < 21; k+)dij += pow(xk, i + j);if (i < 4 && j < 4)aij = dij; bij = dij; cij = dij;else if (i < 6 && j < 6)bij = dij; cij = dij;else

5、 if (i < 8 && j < 8)cij = dij;for (int i = 0; i < 10; i+)for (int j = 0; j < 21; j+)righti += pow(xj, i)*yj;Gauss(*a, right, result_a, 4);Gauss(*b, right, result_b, 6);Gauss(*c, right, result_c, 8);Gauss(*d, right, result_d, 10);for (int i = 0; i < 4; i+)int k = i * 2 + 3;prin

6、tf("在-1,1內(nèi)ex的%d階最小二乘多項(xiàng)式為：n", k);for (int j = 0; j < k; j+)if (j = 0) printf(" %10.8f*x%d", resultij, j);elseif (j % 4 = 0)printf("n");printf("+%10.8f*x%d", resultij, j);printf("nn");八、實(shí)驗(yàn)數(shù)據(jù)及結(jié)果分析1. 主界面，供選擇模擬哪個(gè)實(shí)驗(yàn)的算法。圖1 主界面2. 錄入2，選擇實(shí)驗(yàn)一，計(jì)算得3、5、6、9階的最小二乘多項(xiàng)式結(jié)果如下，根據(jù)驗(yàn)證結(jié)果正確。圖3 實(shí)驗(yàn)2結(jié)果九、實(shí)驗(yàn)結(jié)論在計(jì)算數(shù)學(xué)基礎(chǔ)的這個(gè)小實(shí)驗(yàn)中，我重新系統(tǒng)地復(fù)習(xí)了數(shù)值分析中的最小二乘法。成功地用C/C+模擬了上述算法，對(duì)本學(xué)期的課程所學(xué)的內(nèi)容有了進(jìn)一步的鞏固和提