資源簡(jiǎn)介
第1章 線性代數(shù)方程組的解法
1.1全主元高斯-約當(dāng)(Gauss-Jordan)消去法
1.2LU分解法
1.3追趕法
1.4五對(duì)角線性方程組解法
1.5線性方程組解的迭代改善
1.6范德蒙(Vandermonde)方程組解法
1.7托伯利茲(Toeplitz)方程組解法
1.8奇異值分解
1.9線性方程組的共軛梯度法
1.1對(duì)稱方程組的喬列斯基(Cholesky)分解法
1.11矩陣的QR分解
1.12松弛迭代法
第2章 插值
2.1拉格朗日插值
2.2有理函數(shù)插值
2.3三次樣條插值
2.4有序表的檢索法
2.5插值多項(xiàng)式
2.6二元拉格朗日插值
2.7雙三次樣條插值
第3章 數(shù)值積分
3.1梯形求積法
3.2辛普森(Simpson)求積法
3.3龍貝格(Romberg)求積法
3.4反常積分
3.5高斯(Gauss)求積法
3.6三重積分
第4章 特殊函數(shù)
4.1Γ函數(shù)、貝塔函數(shù)、階乘及二項(xiàng)式系數(shù)
4.2不完全Γ函數(shù)、誤差函數(shù)
4.3不完全貝塔函數(shù)
4.4零階、一階和任意整數(shù)階的第一、二類貝塞爾函數(shù)
4.5零階、一階和任意整數(shù)階的第一、二類變形貝塞爾函數(shù)
4.6分?jǐn)?shù)階第一類貝塞爾函數(shù)和變形貝塞爾函數(shù)
4.7指數(shù)積分和定指數(shù)積分
4.8連帶勒讓德函數(shù)
附錄
第5章 函數(shù)逼近
5.1級(jí)數(shù)求和
5.2多項(xiàng)式和有理函數(shù)
5.3切比雪夫逼近
5.4積分和導(dǎo)數(shù)的切比雪夫逼近
5.5用切比雪夫逼近求函數(shù)的多項(xiàng)式逼近
第6章 隨機(jī)數(shù)
6.1均勻分布隨機(jī)數(shù)
6.2變換方法——指數(shù)分布和正態(tài)分布隨機(jī)數(shù)
6.3舍選法——Γ分布、泊松分布和二項(xiàng)式分布隨機(jī)數(shù)
6.4隨機(jī)位的產(chǎn)生
6.5蒙特卡羅積分法
第7章 排序
7.1直接插入法和Shell方法
7.2堆排序
7.3索引表和等級(jí)表
7.4快速排序
7.5等價(jià)類的確定
附錄
第8章 特征值問(wèn)題
8.1對(duì)稱矩陣的雅可比變換
8.2變實(shí)對(duì)稱矩陣為三對(duì)角對(duì)稱矩陣
8.3三對(duì)角矩陣的特征值和特征向量
8.4變一般矩陣為赫申伯格矩陣
8.5實(shí)赫申伯格矩陣的QR算法
第9章 數(shù)據(jù)擬合
9.1直線擬合
9.2線性最小二乘法
9.3非線性最小二乘法
9.4絕對(duì)值偏差最小的直線擬合
第1章 方程求根和非線性方程組的解法
1.1圖解法
1.2逐步掃描法和二分法
1.3割線法和試位法
1.4布倫特(Brent)方法
1.5牛頓-拉斐森(Newton-Raphson)法
1.6求復(fù)系數(shù)多項(xiàng)式根的拉蓋爾(Laguerre)方法
1.7求實(shí)系數(shù)多項(xiàng)式根的貝爾斯托(Bairstou)方法
1.8非線性方程組的牛頓-拉斐森方法
第11章 函數(shù)的極值和最優(yōu)化
11.1黃金分割搜索法
11.2不用導(dǎo)數(shù)的布倫特(Brent)法
11.3用導(dǎo)數(shù)的布倫特(Brent)法
11.4多元函數(shù)的下山單純形法
11.5多元函數(shù)的包維爾(Powell)法
11.6多元函數(shù)的共軛梯度法
11.7多元函數(shù)的變尺度法
11.8線性規(guī)劃的單純形法
第12章 傅里葉變換譜方法
12.1復(fù)數(shù)據(jù)快速傅里葉變換算法
12.2實(shí)數(shù)據(jù)快速傅里葉變換算法(一)
12.3實(shí)數(shù)據(jù)快速傅里葉變換算法(二)
12.4快速正弦變換和余弦變換
12.5卷積和逆卷積的快速算法
12.6離散相關(guān)和自相關(guān)的快速算法
12.7多維快速傅里葉變換算法
第13章 數(shù)據(jù)的統(tǒng)計(jì)描述
13.1分布的矩——均值、平均差、標(biāo)準(zhǔn)差、方差、斜差和峰態(tài)
13.2中位數(shù)的搜索
13.3均值與方差的顯著性檢驗(yàn)
13.4分布擬合的X2檢驗(yàn)
13.5分布擬合的K-S檢驗(yàn)法
第14章 解常微分方程組
14.1定步長(zhǎng)四階龍格-庫(kù)塔(Runge-Kutta)法
14.2自適應(yīng)變步長(zhǎng)的龍格-庫(kù)塔法
14.3改進(jìn)的中點(diǎn)法
14.4外推法
第15章 兩點(diǎn)邊值問(wèn)題的解法
15.1打靶法(一)
15.2打靶法(二)
15.3松弛法
第16章 偏微分方程的解法
16.1解邊值問(wèn)題的松弛法
16.2交替方向隱式方法(ADI)
代碼片段和文件信息
#include?
#include?
#include?
#include?
#include?
#include?
#include?
int?sgn（double?x）
{
int?temp;
if（x>0）?temp=1;
if（x=0）?temp=0;
if（x<0）?temp=-1;
return?temp;
}
double?bessj1（double?x）
{
double?p1p2p3p4p5q1q2q3q4q5;
double?r1r2r3r4r5r6s1s2s3s4s5s6;
double?bbbcccaaatempaxxxzy;
????r1?=?72362614232.0;???????????r2?=?-7895059235.0;
????r3?=?242396853.1;?????????????r4?=?-2972611.439;
????r5?=?15704.4826;??????????????r6?=?-30.16036606;
????s1?=?144725228442.0;??????????s2?=?2300535178.0;
????s3?=?18583304.74;?????????????s4?=?99447.43394;
????s5?=?376.9991397;?????????????s6?=?1.0;
????p1?=?1.0;?????????????????????p2?=?0.00183105;
????p3?=?-0.00003516396496;???????p4?=?0.000002457520174;
????p5?=?-0.000000240337019;
????q1?=?0.04687499995;???????????q2?=?-0.0002002690873;
????q3?=?0.000008449199096;???????q4?=?-0.00000088228987;
????q5?=?0.000000105787412;
????if?（fabs（x）? {
????????y?=?x*x;
????????aaa?=?r1?+?y?*?（r2?+?y?*?（r3?+?y?*?（r4?+?y?*?（r5?+?y?*?r6））））;
????????bbb?=?s1?+?y?*?（s2?+?y?*?（s3?+?y?*?（s4?+?y?*?（s5?+?y?*?s6））））;
????????temp?=?x?*?aaa?/?bbb;
}
????else
{
????????ax?=?fabs（x）;
????????z?=?8.0?/?ax;
????????y?=?z*z;;
????????xx?=?ax?-?2.356194491;
????????aaa?=?p1?+?y?*?（p2?+?y?*?（p3?+?y?*?（p4?+?y?*?p5）））;
????????bbb?=?q1?+?y?*?（q2?+?y?*?（q3?+?y?*?（q4?+?y?*?q5）））;
????????ccc?=?sqrt（0.636619772?/?ax）;
????????temp?=?ccc?*?（cos（xx）?*?aaa?-?z?*?sin（xx）?*?bbb?*?sgn（x））;
}??
return?temp;
}
double?bessy1（double?x）
{
double?p1p2p3p4p5q1q2q3q4q5;
double?r1r2r3r4r5r6s1s2s3s4s5s6s7;
double?bbbcccaaatempxxzy;
????p1?=?1.0;?????????????????????p2?=?0.00183105;
????p3?=?-0.00003516396496;???????p4?=?0.000002457520174;
????p5?=?-0.000000240337019;
????q1?=?0.04687499995;???????????q2?=?-0.0002002690873;
????q3?=?0.000008449199096;???????q4?=?-0.00000088228987;
????q5?=?0.000000105787412;
????r1?=?-4900604943000.0;????????r2?=?1275274390000.0;
????r3?=?-51534381390.0;??????????r4?=?734926455.1;
????r5?=?-4237922.726;????????????r6?=?8511.937935;
????s1?=?24995805700000.0;????????s2?=?424441966400.0;
????s3?=?3733650367.0;????????????s4?=?22459040.02;
????s5?=?102042.605;??????????????s6?=?354.9632885;
????s7?=?1.0;
????if?（x? {
????????y?=?x*x;
????????aaa?=?r1?+?y?*?（r2?+?y?*?（r3?+?y?*?（r4?+?y?*?（r5?+?y?*?r6））））;
????????bbb?=?s4?+?y?*?（s5?+?y?*?（s6?+?y?*?s7））;
????????bbb?=?s1?+?y?*?（s2?+?y?*?（s3?+?y?*?bbb））;
????????ccc?=?bessj1（x）?*?log（x）?-?1.0?/?x;
????????temp?=?x?*?aaa?/?bbb?+?0.636619772?*?ccc;
}
????else
{
????????z?=?8.0?/?x;
????????y?=?z*z;
????????xx?=?x?-?2.356194491;
????????aaa?=?sqrt（0.636619772?/?x）;
????????bbb?=?p1?+?y?*?（p2?+?y?*?（p3?+?y?*?（p4?+?y?*?p5）））;
????????ccc?=?q1?+?y?*?（q2?+?y?*?（q3?+?y?*?（q4?+?y?*?q5）））;
????????temp?=?aaa?*?（sin（xx）?*?bbb?+?z?*?cos（xx）?*?ccc）;
評(píng)論
共有 條評(píng)論
川公網(wǎng)安備 51152502000135號(hào)