| 查看: 1615 | 回复: 5 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
谱方法中多项式求根的问题
|
|||
|
基于多项式（Legendre、Hermite等）的谱方法求解偏微分方程，难免要涉及到对多项式求根的问题 我用32位计算机C++编程计算，采用通常的方法来求Hermite多项式，超过70次以后，多项式的拟合就已经失真了。而用二分法对Hermite多项式求根，超过40次就解不出了，如果用QR算法那效果更差。舍入误差在多项式中体现得淋漓尽致。 但是我想，求解一个偏微分方程，100个以上的基函数也不算多吧？虽然不要求能够像Fourier基函数那么好的稳定性。可总应该有什么好的算法能够有效模拟100次以上的带权高斯多项式并对其求根，在此请高人指点指点。如果说是要换计算机那就免了。 |
回复此楼» 猜你喜欢
基金系统什么内容也没有
已经有6人回复
-
有多少人是今天查系统知道结果的？
已经有17人回复
-
国自然面上复盘~欢迎讨论
已经有12人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有9人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识
已经有6人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
面上合作单位盖章
已经有5人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有7人回复
-
申请删除本帖
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助个多项式二次型的问题
已经有8人回复
- 请问计算中的谱方法和拟谱方法的本质区别是什么？
已经有5人回复
- matlaB求解多项式反函数遇阻
已经有7人回复
6楼2014-01-01 20:04:20
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
