Fortran Coder

标题: 高斯勒让德求积公式的误差估计 [打印本页]

作者: 岸边的鱼 时间: 2014-7-14 20:27
标题: 高斯勒让德求积公式的误差估计
今天早上2AC的那谁，坚决要我搞定自动变步长的高斯勒让德积分公式的误差估计。2AC的，自己教过数值分析自己不清楚啊。
于是我问百度阿姨，百度阿姨也不愿意说。Google大叔也上不去。问数值分析，数值分析只给出了在标准积分区间上的误差
表达式，于是我无奈的选择使用复合辛普森公式控制了误差，晚上2AC的又召见我，居然让他分析出来了已知份数的高斯勒让德
数值积分公式的误差估计。现在整理为下文，供群友们查阅，以防止再有类似的2AC残害群友。

作者: fcode 时间: 2014-7-14 22:11
楼主研究勒让德很久了啊

作者: 岸边的鱼 时间: 2014-7-15 09:59
本帖最后由岸边的鱼于 2014-7-15 10:12 编辑

fcode 发表于 2014-7-14 22:11
楼主研究勒让德很久了啊

被那个2AC逼得，我就不懂了，我一个物理海洋学的屌丝，研究数值算法有毛用。2AC的以前说使用5个高斯点的公式计算的效果不好，点太少，我被逼的找了任意个高斯点的计算程序。现在我说：“群里有人说，太多高斯点的话计算截断误差会比较大”，2AC的说“恩，那是肯定的，我们就用5点的吧”。早2AC去了！SAC的，成天嫌吧我公式输入的不好看（说什么虚数单位i在公式编辑器中要使用文本格式，试问这个群里谁知道？而且改起来又麻烦，关键这就是平时看，我使用直接输入的怎么了？要是需要投稿，你说改成啥格式我都没话说啊），现在又不投稿，我写好看了，你吃啊

作者: fcode 时间: 2014-7-15 10:23

岸边的鱼发表于 2014-7-15 09:59
被那个2AC逼得，我就不懂了，我一个物理海洋学的屌丝，研究数值算法有毛用。2AC的以前说使用5个高斯点的公 ...

脑残的人很多，不必太在意。按照它的意思做一份，再按照自己的意思做一份既可。

作者: mangix2010 时间: 2014-7-17 15:16
mark，小鱼好厉害~~

作者: 岸边的鱼 时间: 2014-7-19 09:16

mangix2010 发表于 2014-7-17 15:16
mark，小鱼好厉害~~

小鱼，厉害个屁，成天被虐，一会又要去被虐了，那天小鱼不被虐，才是真的厉害

作者: vvt 时间: 2014-7-19 09:21
程序员注定被虐，不是被老板，就是被电脑。永无翻身

作者: 王培杰 时间: 2014-7-23 17:41
我有1000阶2000阶5000阶10000阶的高斯勒让德零点和权系数资料，楼主要不要试试？如果函数在x轴附近震荡的话，还是高斯积分点多一些比较好。可以看看我的附件里面的计算效果。

gsld.rar

8.91 KB, 下载次数: 18

作者: 岸边的鱼 时间: 2014-7-29 09:21

王培杰发表于 2014-7-23 17:41
我有1000阶2000阶5000阶10000阶的高斯勒让德零点和权系数资料，楼主要不要试试？如果函数在x轴附近震荡的话 ...

谢谢，我先看看哈，最近几天回家了，刚回来

作者: 岸边的鱼 时间: 2014-7-29 10:11
本帖最后由岸边的鱼于 2014-7-29 10:13 编辑

王培杰发表于 2014-7-23 17:41
我有1000阶2000阶5000阶10000阶的高斯勒让德零点和权系数资料，楼主要不要试试？如果函数在x轴附近震荡的话 ...

感谢句号大神的友情支援，我刚看了下，惊出一身冷汗，一直以为只要能计算出来结果，那么应该就是差不多的，看了你的统计才知道，原来点少了计算出来了有时候也是错的，我老板现在又不让使用这个方法计算了，点数不好控制，选多了计算量势必很大，选少了计算结果没有保障，现在我们该有自动变步长高斯勒让德计算方法了，计算的结果我发上来你看下，其中m的数量是划分的份数，你应该一看就明白的，效果和真实值比起来还是可以的。
（1）9.99999999999999999998601381929170
M=32
（2）1.57047801697321581899140776826504
M=8192
（3）1.25000000000000000000073298194594
M=4096

作者: 王培杰 时间: 2014-7-29 23:52

岸边的鱼发表于 2014-7-29 10:11
感谢句号大神的友情支援，我刚看了下，惊出一身冷汗，一直以为只要能计算出来结果，那么应该就是差不多的 ...

附件里是高斯勒让德零点和对应权系数，共享出来给需要的人用。。。

GSLD 10000 .dat (859.38 KB, 下载次数: 5)

GSLD 5000 .dat (419.92 KB, 下载次数: 6)

GSLD 2000 .dat (167.97 KB, 下载次数: 1)

GSLD 1000 .dat (99.61 KB, 下载次数: 1)

GSLD 500 .dat (49.8 KB, 下载次数: 1)

GSLD 200 .dat (19.92 KB, 下载次数: 1)

GSLD 100 .dat (9.96 KB, 下载次数: 1)

GSLD 50 .dat (4.98 KB, 下载次数: 1)

GSLD 20 .dat (1.99 KB, 下载次数: 2)

GSLD 10 .dat (1020 Bytes, 下载次数: 2)

GSLD 5 .dat (510 Bytes, 下载次数: 1)

欢迎光临 Fortran Coder (http://bbs.fcode.cn/)