| 查看: 1604 | 回复: 5 | |||
[交流]
谱方法中多项式求根的问题
|
|
基于多项式（Legendre、Hermite等）的谱方法求解偏微分方程，难免要涉及到对多项式求根的问题 我用32位计算机C++编程计算，采用通常的方法来求Hermite多项式，超过70次以后，多项式的拟合就已经失真了。而用二分法对Hermite多项式求根，超过40次就解不出了，如果用QR算法那效果更差。舍入误差在多项式中体现得淋漓尽致。 但是我想，求解一个偏微分方程，100个以上的基函数也不算多吧？虽然不要求能够像Fourier基函数那么好的稳定性。可总应该有什么好的算法能够有效模拟100次以上的带权高斯多项式并对其求根，在此请高人指点指点。如果说是要换计算机那就免了。 |
回复此楼» 猜你喜欢
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有6人回复
-
我面上完蛋了
已经有10人回复
-
能否退出参与的面上项目解除限项
已经有24人回复
-
filecode，4个jtjc了
已经有17人回复
-
人气不行了
已经有9人回复
-
没有任何消息-是不是就凉了
已经有6人回复
-
2026国自然函评费到账
已经有19人回复
-
建议基金发布提前给出明确的时间点
已经有16人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有21人回复
-
估计是周四
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助个多项式二次型的问题
已经有8人回复
- 请问计算中的谱方法和拟谱方法的本质区别是什么？
已经有5人回复
- matlaB求解多项式反函数遇阻
已经有7人回复
2楼2013-12-31 19:09:43
3楼2013-12-31 19:47:42
★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
feixiaolin: 金币+3 2013-12-31 22:37:10
lixy1217: 金币+5 2014-01-01 19:58:05
小木虫: 金币+0.5, 给个红包，谢谢回帖
feixiaolin: 金币+3 2013-12-31 22:37:10
lixy1217: 金币+5 2014-01-01 19:58:05
|
此言差矣，在实际的数值计算过程中，一般都不需要把多项式直接表示出来（否则就会导致你所说的舍入误差积累），如果是正交多项式，一可以利用递推公式去计算n次正交多项式的函数值或者导数值。 我可以提供matlab子程序给你，用于计算Hermite求积节点和权, 当多项式次数很高时（n>100），你可以检验下，计算结果非常精确。代码见附件。 |
4楼2013-12-31 21:22:11
5楼2013-12-31 22:35:07
6楼2014-01-01 20:04:20
