24小时热门版块排行榜    

查看: 1339  |  回复: 5

lixy1217

木虫 (著名写手)


[交流] 谱方法中多项式求根的问题

基于多项式(Legendre、Hermite等)的谱方法求解偏微分方程,难免要涉及到对多项式求根的问题

我用32位计算机C++编程计算,采用通常的方法来求Hermite多项式,超过70次以后,多项式的拟合就已经失真了。而用二分法对Hermite多项式求根,超过40次就解不出了,如果用QR算法那效果更差。舍入误差在多项式中体现得淋漓尽致。

但是我想,求解一个偏微分方程,100个以上的基函数也不算多吧?虽然不要求能够像Fourier基函数那么好的稳定性。可总应该有什么好的算法能够有效模拟100次以上的带权高斯多项式并对其求根,在此请高人指点指点。如果说是要换计算机那就免了。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cool_smile

木虫 (著名写手)


★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
lixy1217: 金币+2, 欢迎交流 2013-12-31 19:37:49
方法一、可以利用特征值方法求正交多项式的零点。譬如:N次Hermite多项式的零点恰好就是一个三对角(对称)矩阵的特征值。
方法二、利用迭代法求根,关键是初值的选取。初值怎么选取可以参考文献 Pan V. (1997, 39(2), SIAM Review)
2楼2013-12-31 19:09:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lixy1217

木虫 (著名写手)


引用回帖:
2楼: Originally posted by cool_smile at 2013-12-31 19:09:43
方法一、可以利用特征值方法求正交多项式的零点。譬如:N次Hermite多项式的零点恰好就是一个三对角(对称)矩阵的特征值。
方法二、利用迭代法求根,关键是初值的选取。初值怎么选取可以参考文献 Pan V. (1997, 39 ...

不知道你试着算过没有。其实求多项式的根本身不是太难,可是次数一旦高起来,舍入误差的影响将会非常大。一个70次的Hermite多项式,它的系数在数量级上差别能达到10的40次方。而且x的70次方,就算是一点小小的扰动,带来的差别也是巨大的。用32位计算机操作,不说求根,如何来把这个多项式表示出来,都是个问题。
3楼2013-12-31 19:47:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cool_smile

木虫 (著名写手)


★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
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
不知道你试着算过没有。其实求多项式的根本身不是太难,可是次数一旦高起来,舍入误差的影响将会非常大。一个70次的Hermite多项式,它的系数在数量级上差别能达到10的40次方。而且x的70次方,就算是一点小小的扰动 ...

此言差矣,在实际的数值计算过程中,一般都不需要把多项式直接表示出来(否则就会导致你所说的舍入误差积累),如果是正交多项式,一可以利用递推公式去计算n次正交多项式的函数值或者导数值。

我可以提供matlab子程序给你,用于计算Hermite求积节点和权, 当多项式次数很高时(n>100),你可以检验下,计算结果非常精确。代码见附件。

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱: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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lixy1217

木虫 (著名写手)


引用回帖:
4楼: Originally posted by cool_smile at 2013-12-31 21:22:11
此言差矣,在实际的数值计算过程中,一般都不需要把多项式直接表示出来(否则就会导致你所说的舍入误差积累),如果是正交多项式,一可以利用递推公式去计算n次正交多项式的函数值或者导数值。

我可以提供matla ...

听你这么一说,好像有点意思~~看来直观上先得到系数在研究多项式的方式在这里是不适用了。
不禁追问一下,三对角对称矩阵表示Hermite多项式的根,这是个什么三对角矩阵?迭代法求根,是常用的牛顿迭代吗?如果这么说来,更简单的二分法是否也会适用,因为我只需要知道实根?
5楼2013-12-31 22:35:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lixy1217

木虫 (著名写手)


引用回帖:
4楼: Originally posted by cool_smile at 2013-12-31 21:22:11
此言差矣,在实际的数值计算过程中,一般都不需要把多项式直接表示出来(否则就会导致你所说的舍入误差积累),如果是正交多项式,一可以利用递推公式去计算n次正交多项式的函数值或者导数值。

我可以提供matla ...

果然不应该求多项式系数,而应该使用Hermite基函数(不是Hermite多项式)迭代公式直接逐点表示。用最原始野蛮的二分法,可以求到750阶以内的Hermite基函数的根。但是对于更高阶的多项式,由于根的范围突破了38.5这个坎,意味着exp(-0.5*x*x)会低于C程序数据范围的下限,从而直接等于0,也就无法求根了。
不知道这时候有什么好的处理办法(呵呵,我这个人有点贪心)
6楼2014-01-01 20:04:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lixy1217 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +5 yunziaaaaa 2026-03-01 6/300 2026-03-01 23:57 by ccp273206157
[考研] 一志愿郑大材料学硕298分,求调剂 +5 wsl111 2026-03-01 5/250 2026-03-01 23:45 by 暮雨星晴
[考研] 材料学硕318求调剂 +5 February_Feb 2026-03-01 5/250 2026-03-01 23:31 by L135790
[考研] 0856调剂 +5 刘梦微 2026-02-28 5/250 2026-03-01 22:30 by wang_dand
[硕博家园] 博士自荐 +7 科研狗111 2026-02-26 11/550 2026-03-01 22:24 by 哲平L
[考研] 26考研报考西工大材料308分求调剂 +3 weizhong123 2026-03-01 3/150 2026-03-01 21:42 by 公瑾逍遥
[考研] 272求调剂 +6 田智友 2026-02-28 6/300 2026-03-01 21:40 by 公瑾逍遥
[考研] 0856求调剂285 +10 吕仔龙 2026-02-28 10/500 2026-03-01 21:37 by 公瑾逍遥
[考研] 一志愿中南大学理学化学 +4 15779376950 2026-03-01 5/250 2026-03-01 19:00 by Fff-1
[考博] 26申博 +4 想申博! 2026-02-26 6/300 2026-03-01 17:32 by 想申博!
[考研] 328求调剂 +3 aaadim 2026-03-01 5/250 2026-03-01 17:29 by njzyff
[考研] 281求调剂 +4 2026计算机_诚心 2026-03-01 7/350 2026-03-01 17:20 by 2026计算机_诚心
[考研] 321求调剂一志愿东北林业大学材料与化工英二数二 +4 虫虫虫虫虫7 2026-03-01 7/350 2026-03-01 16:52 by caszguilin
[考研] 0856材料求调剂 +4 麻辣鱿鱼 2026-02-28 4/200 2026-03-01 16:51 by caszguilin
[考研] 285求调剂 +8 满头大汗的学生 2026-02-28 8/400 2026-03-01 16:47 by caszguilin
[考研] 311求调剂 +6 亭亭亭01 2026-03-01 6/300 2026-03-01 15:41 by 324616
[考研] 求调剂 +6 repeatt?t 2026-02-28 6/300 2026-03-01 14:37 by Sakura绘
[考研] 课题组接收材料类调剂研究生 +3 gaoxiaoniuma 2026-02-28 4/200 2026-03-01 14:30 by jjj三跨
[硕博家园] 2025届双非化工硕士毕业,申博 +3 更多的是 2026-02-27 4/200 2026-03-01 10:04 by ztg729
[基金申请] 面上模板改不了页边距吧? +5 ieewxg 2026-02-25 6/300 2026-03-01 00:10 by addressing
信息提示
请填处理意见