lixy1217

木蟲 (著名寫手)

應(yīng)助: 4 (幼兒園)
金幣: 4618.6
帖子: 1381
在線: 234.2小時(shí)
蟲號(hào): 1125891

[交流] 譜方法中多項(xiàng)式求根的問題

基于多項(xiàng)式（Legendre、Hermite等）的譜方法求解偏微分方程，難免要涉及到對(duì)多項(xiàng)式求根的問題

我用32位計(jì)算機(jī)C++編程計(jì)算，采用通常的方法來求Hermite多項(xiàng)式，超過70次以后，多項(xiàng)式的擬合就已經(jīng)失真了。而用二分法對(duì)Hermite多項(xiàng)式求根，超過40次就解不出了，如果用QR算法那效果更差。舍入誤差在多項(xiàng)式中體現(xiàn)得淋漓盡致。

但是我想，求解一個(gè)偏微分方程，100個(gè)以上的基函數(shù)也不算多吧？雖然不要求能夠像Fourier基函數(shù)那么好的穩(wěn)定性�？煽倯�(yīng)該有什么好的算法能夠有效模擬100次以上的帶權(quán)高斯多項(xiàng)式并對(duì)其求根，在此請(qǐng)高人指點(diǎn)指點(diǎn)。如果說是要換計(jì)算機(jī)那就免了。

回復(fù)此樓

» 猜你喜歡

306求調(diào)劑已經(jīng)有8人回復(fù)
招08考數(shù)學(xué) 已經(jīng)有12人回復(fù)
接收2026碩士調(diào)劑(學(xué)碩+專碩) 已經(jīng)有4人回復(fù)
北科281學(xué)碩材料求調(diào)劑已經(jīng)有9人回復(fù)
306求0703調(diào)劑一志愿華中師范已經(jīng)有6人回復(fù)
276求調(diào)劑。有半年電池和半年高分子實(shí)習(xí)經(jīng)歷已經(jīng)有8人回復(fù)
材料與化工考研調(diào)劑已經(jīng)有3人回復(fù)
081700 調(diào)劑 267分已經(jīng)有3人回復(fù)
一志愿重慶大學(xué)085700資源與環(huán)境，總分308求調(diào)劑已經(jīng)有3人回復(fù)
336化工調(diào)劑已經(jīng)有4人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

求助個(gè)多項(xiàng)式二次型的問題已經(jīng)有8人回復(fù)
請(qǐng)問計(jì)算中的譜方法和擬譜方法的本質(zhì)區(qū)別是什么？已經(jīng)有5人回復(fù)
matlaB求解多項(xiàng)式反函數(shù)遇阻已經(jīng)有7人回復(fù)

» 搶金幣啦！回帖就可以得到:

查看全部散金貼

1樓 2013-12-31 12:02:24

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

cool_smile

木蟲 (著名寫手)

應(yīng)助: 23 (小學(xué)生)
貴賓: 0.01
金幣: 4916.4
帖子: 1033
在線: 966.9小時(shí)
蟲號(hào): 503053

★ ★ ★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖
lixy1217: 金幣+2, 歡迎交流 2013-12-31 19:37:49

方法一、可以利用特征值方法求正交多項(xiàng)式的零點(diǎn)。譬如：N次Hermite多項(xiàng)式的零點(diǎn)恰好就是一個(gè)三對(duì)角（對(duì)稱）矩陣的特征值。
方法二、利用迭代法求根，關(guān)鍵是初值的選取。初值怎么選取可以參考文獻(xiàn) Pan V. (1997, 39(2), SIAM Review)

贊一下

回復(fù)此樓

2樓2013-12-31 19:09:43

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

lixy1217

木蟲 (著名寫手)

應(yīng)助: 4 (幼兒園)
金幣: 4618.6
帖子: 1381
在線: 234.2小時(shí)
蟲號(hào): 1125891

引用回帖:

2樓: Originally posted by cool_smile at 2013-12-31 19:09:43
方法一、可以利用特征值方法求正交多項(xiàng)式的零點(diǎn)。譬如：N次Hermite多項(xiàng)式的零點(diǎn)恰好就是一個(gè)三對(duì)角（對(duì)稱）矩陣的特征值。
方法二、利用迭代法求根，關(guān)鍵是初值的選取。初值怎么選取可以參考文獻(xiàn) Pan V. (1997, 39 ...

不知道你試著算過沒有。其實(shí)求多項(xiàng)式的根本身不是太難，可是次數(shù)一旦高起來，舍入誤差的影響將會(huì)非常大。一個(gè)70次的Hermite多項(xiàng)式，它的系數(shù)在數(shù)量級(jí)上差別能達(dá)到10的40次方。而且x的70次方，就算是一點(diǎn)小小的擾動(dòng)，帶來的差別也是巨大的。用32位計(jì)算機(jī)操作，不說求根，如何來把這個(gè)多項(xiàng)式表示出來，都是個(gè)問題。

贊一下

回復(fù)此樓

3樓2013-12-31 19:47:42

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

cool_smile

木蟲 (著名寫手)

應(yīng)助: 23 (小學(xué)生)
貴賓: 0.01
金幣: 4916.4
帖子: 1033
在線: 966.9小時(shí)
蟲號(hào): 503053

★ ★ ★ ★ ★ ★ ★ ★ ★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖
feixiaolin: 金幣+3 2013-12-31 22:37:10
lixy1217: 金幣+5 2014-01-01 19:58:05

引用回帖:

3樓: Originally posted by lixy1217 at 2013-12-31 19:47:42
不知道你試著算過沒有。其實(shí)求多項(xiàng)式的根本身不是太難，可是次數(shù)一旦高起來，舍入誤差的影響將會(huì)非常大。一個(gè)70次的Hermite多項(xiàng)式，它的系數(shù)在數(shù)量級(jí)上差別能達(dá)到10的40次方。而且x的70次方，就算是一點(diǎn)小小的擾動(dòng) ...

此言差矣，在實(shí)際的數(shù)值計(jì)算過程中，一般都不需要把多項(xiàng)式直接表示出來（否則就會(huì)導(dǎo)致你所說的舍入誤差積累），如果是正交多項(xiàng)式，一可以利用遞推公式去計(jì)算n次正交多項(xiàng)式的函數(shù)值或者導(dǎo)數(shù)值。

我可以提供matlab子程序給你，用于計(jì)算Hermite求積節(jié)點(diǎn)和權(quán), 當(dāng)多項(xiàng)式次數(shù)很高時(shí)（n>100），你可以檢驗(yàn)下，計(jì)算結(jié)果非常精確。代碼見附件。

贊一下

回復(fù)此樓

» 本帖附件資源列表

歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對(duì)版權(quán)有異議，請(qǐng)聯(lián)系郵箱：xiaomuchong@tal.com
附件 1 : hegs.m

2013-12-31 21:20:49, 692 bytes

附件 2 : hefun.m

2013-12-31 21:21:13, 1.45 K

附件 3 : hepolyn.m

2013-12-31 21:21:32, 1000 bytes

4樓2013-12-31 21:22:11

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

lixy1217

木蟲 (著名寫手)

應(yīng)助: 4 (幼兒園)
金幣: 4618.6
帖子: 1381
在線: 234.2小時(shí)
蟲號(hào): 1125891

引用回帖:

4樓: Originally posted by cool_smile at 2013-12-31 21:22:11
此言差矣，在實(shí)際的數(shù)值計(jì)算過程中，一般都不需要把多項(xiàng)式直接表示出來（否則就會(huì)導(dǎo)致你所說的舍入誤差積累），如果是正交多項(xiàng)式，一可以利用遞推公式去計(jì)算n次正交多項(xiàng)式的函數(shù)值或者導(dǎo)數(shù)值。

我可以提供matla ...

聽你這么一說，好像有點(diǎn)意思~~看來直觀上先得到系數(shù)在研究多項(xiàng)式的方式在這里是不適用了。
不禁追問一下，三對(duì)角對(duì)稱矩陣表示Hermite多項(xiàng)式的根，這是個(gè)什么三對(duì)角矩陣？迭代法求根，是常用的牛頓迭代嗎？如果這么說來，更簡單的二分法是否也會(huì)適用，因?yàn)槲抑恍枰缹?shí)根？

贊一下

回復(fù)此樓

5樓2013-12-31 22:35:07

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

lixy1217

木蟲 (著名寫手)

應(yīng)助: 4 (幼兒園)
金幣: 4618.6
帖子: 1381
在線: 234.2小時(shí)
蟲號(hào): 1125891

引用回帖:

果然不應(yīng)該求多項(xiàng)式系數(shù)，而應(yīng)該使用Hermite基函數(shù)（不是Hermite多項(xiàng)式）迭代公式直接逐點(diǎn)表示。用最原始野蠻的二分法，可以求到750階以內(nèi)的Hermite基函數(shù)的根。但是對(duì)于更高階的多項(xiàng)式，由于根的范圍突破了38.5這個(gè)坎，意味著exp(-0.5*x*x)會(huì)低于C程序數(shù)據(jù)范圍的下限，從而直接等于0，也就無法求根了。
不知道這時(shí)候有什么好的處理辦法（呵呵，我這個(gè)人有點(diǎn)貪心）

贊一下

回復(fù)此樓

6樓2014-01-01 20:04:20

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 北科281學(xué)碩材料求調(diào)劑 +8	tcxiaoxx 2026-03-20	9/450	2026-03-23 12:16 by tcxiaoxx
[考研] 求調(diào)劑材料學(xué)碩080500，總分289分 5+3	@taotao 2026-03-19	21/1050	2026-03-23 10:17 by 冠c哥
[考研] 求調(diào)劑 +5	Zhangbod 2026-03-21	7/350	2026-03-22 13:13 by Zhangbod
[考研] 280求調(diào)劑 +11	咕嚕曉曉 2026-03-18	12/600	2026-03-21 22:40 by ACS Nano——
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +3	邱gl 2026-03-20	7/350	2026-03-21 19:05 by 15709483992
[考研] 一志愿南大，0703化學(xué)，分?jǐn)?shù)336，求調(diào)劑 +3	收到VS 2026-03-21	3/150	2026-03-21 18:42 by 學(xué)員8dgXkO
[考研] 工科0856求調(diào)劑 +3	沐析汀汀 2026-03-21	3/150	2026-03-21 18:30 by 學(xué)員8dgXkO
[考研] 278求調(diào)劑 +9	煙火先于春 2026-03-17	9/450	2026-03-21 17:47 by 學(xué)員8dgXkO
[考研] 297求調(diào)劑 +11	戲精丹丹丹 2026-03-17	12/600	2026-03-21 17:47 by ColorlessPI
[考研] 296求調(diào)劑 +4	www_q 2026-03-20	4/200	2026-03-21 17:26 by 學(xué)員8dgXkO
[考研] 299求調(diào)劑 +4	某某某某位 2026-03-21	4/200	2026-03-21 16:30 by barlinike
[考研] 085601調(diào)劑 358分 +3	zzzzggh 2026-03-20	4/200	2026-03-21 10:21 by luoyongfeng
[考研] 070300化學(xué)319求調(diào)劑 +7	錦鯉0909 2026-03-17	7/350	2026-03-21 03:46 by JourneyLucky
[考研] 一志愿西南交大，求調(diào)劑 +5	材化逐夢(mèng)人 2026-03-18	5/250	2026-03-21 00:26 by JourneyLucky
[考研] 321求調(diào)劑 +9	何潤采123 2026-03-18	11/550	2026-03-20 23:19 by JourneyLucky
[考研] 290求調(diào)劑 +7	^O^乜 2026-03-19	7/350	2026-03-20 21:43 by JourneyLucky
[考研] 295復(fù)試調(diào)劑 +8	簡木ChuFront 2026-03-19	8/400	2026-03-20 20:44 by zhukairuo
[考研] 085601求調(diào)劑 +4	Du.11 2026-03-16	4/200	2026-03-17 17:08 by ruiyingmiao
[考研] 一志愿南京大學(xué)，080500材料科學(xué)與工程，調(diào)劑 +4	Jy? 2026-03-16	4/200	2026-03-17 11:02 by gaoqiong
[考研] 283求調(diào)劑 +3	聽風(fēng)就是雨； 2026-03-16	3/150	2026-03-17 07:41 by 熱情沙漠