設萬維讀者為首頁 萬維讀者網 -- 全球華人的精神家園 廣告服務 聯繫我們 關於萬維
 
首  頁 新  聞 視  頻 博  客 論  壇 分類廣告 購  物
搜索>> 發表日誌 控制面板 個人相冊 給我留言
幫助 退出
     
  木樁的博客
  我來自金星
我的名片
木樁
註冊日期: 2014-11-04
訪問總量: 739,355 次
點擊查看我的個人資料
Calendar
我的公告欄
最新發布
· 公益貼 ---- 談談川普的關稅
· 2025 回國見聞 - 上海淮海路的今
· 2025 回國見聞 - 西安和杭州
· 2025 重返故鄉上海
· 閒談應用數學 --- 計算流體力學
· 文化長廊
· 一隻小鳥的詠嘆
友好鏈接
分類目錄
【科學與學術 討論- 3】
· 公益貼 ---- 談談川普的關稅
· 閒談應用數學 --- 計算流體力學
【科學與學術討論 - 2】
· 飛機安全嗎?--- 分析幾次重大空
· 談談世界頂尖科學家的排名榜
· 公益貼 --- 再談疫苗
· 公益帖 --- 新冠疫苗的最新發展
· 公益貼 --- 談談疫苗
· 科學計算 (Scientific Computing
· 從三峽大壩想到水利工程
· 檢測新冠病毒的3 種方法 --- 分
· 波音737- MAX8 出事原因在哪裡?
· 介紹一個美國頂級應用數學系
【科學與學術討論 -1】
· 格林威治 -- 國際標準時間
· 美國的福利保險費是怎樣算出來的
· 美國的401(K)養老金
· 關於量子保密通訊
· 揭露中國學者在期刊上造假作弊
· 中國和美國教育學術界的若干比較
· 計算數學家在幹些什麼?
· 學術雜誌的“影響因子”
· 有感中國人的高校情結
· 首爾-2014世界數學家大會
【回國感想 - I】
· 2025 重返故鄉上海
· 2015 回國見聞 - 爽爽貴陽
· 2015 回國見聞 - 貴州織金洞
· 2015 回國見聞-黃果樹大瀑布
· 2015 回國見聞 -雲南石林
· 回國見聞 -上海受騙上當記
· 回國見聞 -國內的幾個有趣現象
【回國感想- II】
· 2025 回國見聞 - 上海淮海路的今
· 2025 回國見聞 - 西安和杭州
【茗香茶語】
· 天鵝的夢囈
· 母親節,櫻花又開了
· 禮輕情意重
【回憶錄】
· 童年的回憶
· 西南聯大逸聞軼事
· 魂縈上海灘(5)-文革中的那些鄰居
· 魂縈上海灘(4)--我的母親
· 魂縈上海灘(3)-媽媽給我的禮物
· 魂縈上海灘(2)–綠房子裡的故事
· 魂縈上海灘(1)-萬航渡路
· 童年舊事 -我的小提琴老師
· 父親留下的傳家寶
【藝術-1】
· 藝術座談 (2) - 藝術的兩重性
· 藝術座談 (1) --- 怎樣理解音樂
· 藝術也能變戲法
· 我的藝術創作 --- 雕塑, 油畫,
· 話說故宮
· 漫談藝術——西方畫 vs 中國畫
· 裝修設計 - 讓電腦為人類服務
· 家庭住宅的裝修和設計
· 淺談人體雕塑
【藝術-2】
【鋼琴演奏】
· 愛吧!能愛就愛吧!—- 李斯特的
· 聆聽古典音樂,欣賞現代建築
· 肖邦夜曲,情愁悽美,孤獨浪漫
· 介紹鋼琴夜曲的祖宗 John Field
· 舒曼的鋼琴曲--《異國陌人》
· 音樂和數學 - 十二平均律
· 2015 聖誕之夜-《少女的祈禱》
· 鋼琴獨奏 -《 夢幻》
【鋼琴演奏 -2】
【世界各地】
· 倫敦白金漢宮皇家衛隊換崗儀式
· 希臘羅德島太陽神巨像失蹤之謎
· 巴西--寫在2016 夏季奧運會之前
· 西班牙的小偷和鬥牛士
· 法國的香農舍 -女人的城堡
· 印度的硅谷
· 莫斯科之行
· 巴黎的墓地
【歌劇】
· 庚子年 --- 媽媽,生日快樂!
· 二首上海經典老歌給萬維朋友拜年
· 祝萬維朋友 平安夜好!聖誕快樂
· 歌劇中的詠嘆調《我親愛的爸爸》
· 歌劇欣賞 (2)-我聽歌劇的體會
· 歌劇欣賞(1)- 什麼是歌劇
【詩詞歌賦】
· 文化長廊
· 七言絕句三首 -《春天的思念》
· 七言絕句二首 - 《訪友》
· 楓葉又紅了 ---- 七言 . 思君
· 名家誦讀 --- 唐詩二首
【芸芸眾生】
· 一隻小鳥的詠嘆
· 為生活在戰爭中的人們祈禱 !
· 深秋的祭念
· Memorial Day - 告慰全世界新冠
· 20世紀最大的藝術偽造案 -米格倫
· 龐氏騙局
· 從數學家納什聯想到 。。。
· 我的露露和咪咪
【核桃樹 --- 萬維博客科學沙龍】
· 藝術也能變戲法
· 紀念切爾諾貝利核事故三十周年:
· 瘋子數學家 -- 約翰.納什
【高山流水論壇年慶】
· 【金牛賀春】貝多芬鋼琴奏鳴曲《
· 【金豬賀歲】上海懷舊老歌《月圓
· 【金狗旺旺】— 獻給愛麗絲
· 【高山年慶】-- 向朋友們致敬!
存檔目錄
02/01/2026 - 02/28/2026
10/01/2025 - 10/31/2025
09/01/2025 - 09/30/2025
08/01/2025 - 08/31/2025
01/01/2025 - 01/31/2025
12/01/2024 - 12/31/2024
04/01/2024 - 04/30/2024
04/01/2022 - 04/30/2022
03/01/2022 - 03/31/2022
01/01/2022 - 01/31/2022
12/01/2021 - 12/31/2021
04/01/2021 - 04/30/2021
03/01/2021 - 03/31/2021
02/01/2021 - 02/28/2021
12/01/2020 - 12/31/2020
11/01/2020 - 11/30/2020
10/01/2020 - 10/31/2020
09/01/2020 - 09/30/2020
08/01/2020 - 08/31/2020
06/01/2020 - 06/30/2020
05/01/2020 - 05/31/2020
04/01/2020 - 04/30/2020
03/01/2020 - 03/31/2020
02/01/2020 - 02/29/2020
01/01/2020 - 01/31/2020
12/01/2019 - 12/31/2019
11/01/2019 - 11/30/2019
10/01/2019 - 10/31/2019
07/01/2019 - 07/31/2019
03/01/2019 - 03/31/2019
02/01/2019 - 02/28/2019
12/01/2018 - 12/31/2018
11/01/2018 - 11/30/2018
10/01/2018 - 10/31/2018
08/01/2018 - 08/31/2018
07/01/2018 - 07/31/2018
02/01/2018 - 02/28/2018
12/01/2017 - 12/31/2017
11/01/2017 - 11/30/2017
10/01/2017 - 10/31/2017
09/01/2017 - 09/30/2017
08/01/2017 - 08/31/2017
05/01/2017 - 05/31/2017
04/01/2017 - 04/30/2017
03/01/2017 - 03/31/2017
02/01/2017 - 02/28/2017
12/01/2016 - 12/31/2016
11/01/2016 - 11/30/2016
10/01/2016 - 10/31/2016
09/01/2016 - 09/30/2016
07/01/2016 - 07/31/2016
06/01/2016 - 06/30/2016
05/01/2016 - 05/31/2016
04/01/2016 - 04/30/2016
03/01/2016 - 03/31/2016
02/01/2016 - 02/29/2016
01/01/2016 - 01/31/2016
12/01/2015 - 12/31/2015
11/01/2015 - 11/30/2015
10/01/2015 - 10/31/2015
09/01/2015 - 09/30/2015
08/01/2015 - 08/31/2015
07/01/2015 - 07/31/2015
06/01/2015 - 06/30/2015
05/01/2015 - 05/31/2015
04/01/2015 - 04/30/2015
03/01/2015 - 03/31/2015
02/01/2015 - 02/28/2015
01/01/2015 - 01/31/2015
12/01/2014 - 12/31/2014
發表評論
作者:
用戶名: 密碼: 您還不是博客/論壇用戶?現在就註冊!
     
評論:
科學計算 (Scientific Computing)
   

常常碰到網友問我,什麼是計算數學?什麼是科學計算?這篇文章企圖闡明數學在“科學計算”里的作用。

首先,什麼是科學計算?在這方面有着不同的定義。在我看來,這應該是用數學作為工具,設計有效的算法,然後在計算機上實現。的確,科學計算和計算機科學有着緊密的聯繫,但不是一回事。科學計算着重於數學的角色,強調數學的作用。通過數學,保證設計方法的有效性,有足夠的精確度,足夠的穩定性,而且提高效率。

還是通過幾個有名的例子來說明問題吧。


例一:快速傅立葉變換 (Fast Fourier transform, FFT)

我想,這個方法大家都不陌生吧?在很多科學技術領域裡,應用的那個部分里,都需要做傅立葉變換,事實上,當今世界上的計算機在運算的時侯,很多都在作快速傅立葉變換。這裡,你就可以看到數學對它的貢獻了。其實,這僅僅是初等數學,而且還是比較簡單的數學。

一維空間:從  2π 之間,均勻地分布個點。

頻譜空間:傅立葉級數展開。 

傅立葉級數的係數和點空間之間是有關係的。這樣說吧,如果你知道了這些點值,就能算出這些係數,反之,知道了這些係數,也能算出那些點值。 

計算機在點值和係數的空間來回計算,這是要有代價的。如果已知係數要算一個點值,這是 N 數量級的工作量;把所有N個點值全算出來  則是 N2 數量級的工作量。反之,如果已知點值,要把所有N個係數全算出來,也  N2  數量級的工作量。當 N 很大時,這是相當大的計算量,這不,如果 N = 100000,再平方一下,這數字太可怕了,而且,計算機的大量運算會產生誤差積累,導致算出的結果不夠精確。

能不能把 N降下來?能快速傅立葉變換( FFT) 的思想其實很簡單:仔細觀察公式,發現有很多重複計算。在算點值時,打個比方(純粹胡說):第一個點,第四個點,第七個點。。。有重複計算。FFT 把重複計算的部分,一次性算好,然後在第一個點,第四個點,第七個點。。。同時使用。就這麼點簡單的思想,你只要把它做到極致,一個 N2  數量級的工作量就降到了N log N 數量級了,相比於 Nlog N 是一個很小的數,幾乎等同於常數。這,是一個革命性的創舉。

可以這麼說,如果沒有得益於快速傅立葉變換( FFT)的應用,各種信號處理,電子工程,絕對不會變得像今天這麼廣泛。FFT 在數學上是一個很簡單的東西,這件事情之所以做得非常飄亮,就是在數學上並沒有做新的事情,也沒有改變原來的公式,思路也是簡單的,只是把傅立葉公式算得很快,多快好省幹革命。

FFT的第一篇文章,1965年發表在Mathematics of Computation , 這是美國數學學會的頂級雜誌,也是非常古老的計算數學雜誌。這篇FFT 論文的影響力是巨大的。

 

例二:快速多極法 fast multi-pole method 

這種方法是用來快速計算 N 個粒子之間兩兩相互作用的效果。 比如,兩兩粒子間的引力作用。對一個固定的粒子來說,它的運動取決於它和其他所有粒子兩兩相互作用的效果之和。把這些相互作用都算清楚,需要 N 數量級的工作量。那麼,要算好所有 N 個粒子的運動規律,則需要 N2  數量級的工作量。

這次運氣不好,和 FFT 不同,我們似乎找不到明顯的重複運算 (除了兩兩相互作用只需算一次,這只能減少一半運算量,沒法在數量級上減少)。但這難不倒數學家。數學家們發現,一個粒子和距離它較遠的一堆粒子的相互作用,不用兩兩算清楚,只需要算這個粒子和這一堆粒子(看成一個大粒子)的相互作用,也就八九不離十了。當然,誤差是多少,一堆粒子到底可以多大,如何定義 “較遠”,這些都要分析清楚。這次沒有辦法像 FFT 一樣得到和原來公式一樣的數學效果,但是,如果給定能容忍的誤差,比如千分之一,則上述的辦法只需要 N log N 數量級的工作量。實際上,如果把算法和分析再做得精細些,工作量能降到 N 數量級。這種方法叫做 “快速多極法”。 

快速多極法的發明,歸功於 耶魯大學的 V. Rokhlin L. Greengard1985 年,他們發表在 Journal of Computational Physics 的一篇文章中,首創性地提出和分析了快速多極法。和FFT 類似,這個方法也在科學和工程中得到廣泛的應用。因為這個工作,這兩位數學家得到了美國數學學會 2001 年的 一個重量級獎項 (the Leroy P. Steele Prize)。這項獎通常都是授予純粹數學家的,應用數學家極少得到該獎,這也說明了快速多極法在數學和應用中的雙重重要性。


例三:多重網格法 multigrid method

怎樣解大型的線性方程組

                         Ax = b 

A N x N矩陣,是向量。這個問題看上去十分簡單,用中學學到的Cramer’s Rule 就能得到解。不過,用Cramer’s Rule 的計算量是 N 的階乘(N!)的數量級。當矩陣很大的時候,比如, N=1000 N! 是個非常龐大的數,不信你在計算器上摁一摁,會發現計算器都無法顯示出這麼龐大的結果。不管今天的計算機有多快,都無法在短時間內實現這樣的計算。所以,Cramer’s Rule 在解大型的線性方程組時,只有理論上的意義,沒有實際意義,紙上談兵也。

用高斯消除法(Gaussian elimination)求解,可以將工作量從 N! 數量級降到N3數量級。這已經大大提高了計算效率。對一般的線性方程組,這經常是實際計算中唯一可行的方法。不過, N3  還是實在太大了,數學家希望再進一步降低工作量。對於一些實際應用中的矩陣,有些特殊的性質,比如對稱正定矩陣,何不利用一下?應用數學知識可以設計 N 數量級的算法,如果能如願以賞,的確是美事一樁,故而發明了“多重網格法”。要知道,不是所有的矩陣都能採用多重網格法,矩陣需要滿足特殊的性質。因為哪怕A是對角矩陣,求解也需要 N 數量級的計算,再也沒有油水可榨了,N數量級已經達到了極限。

多重網格法的主要貢獻者是以色列數學家Achi Brandt , 為此,2005 年,他得到了美國工業與應用數學學會(SIAM)和美國計算機學會 (ACM) 共同授予的 “計算科學和工程獎”  (Prize in Computational Science and Engineering ) 這是個非常了不起的大獎。


向科學家們致以崇高的敬禮!


 


 
關於本站 | 廣告服務 | 聯繫我們 | 招聘信息 | 網站導航 | 隱私保護
Copyright (C) 1998-2026. Creaders.NET. All Rights Reserved.