应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » 其它 » 【致H兄】数列题的万能多项式拟合求法
北京石油化工学院2026年研究生招生接收调剂公告
292/3返回列表上一页123下一页
查看: 2279  |  回复: 28
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
本帖产生 2 个 程序强帖 ，点击这里进行查看

wangww2011

木虫 (著名写手)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
xzhdty(金币+2): 中秋快乐，欢迎常来 2011-09-09 23:09:51
引用回帖:
6楼: Originally posted by sudo at 2011-09-09 19:42:55:
这个拟合出来的式子也可以无限长啊

哦，莫非你觉得不是整数不爽？简单改改就行：

把n从0开始计数，设给定的数列长度为N-1，同理使用上面的方法得到一个多项式F(n)之后，再令

f(n) = F(n mod N)

不就 ...

不是这个意思，举个例子(这台电脑上的MATLAB被卸载了,就用Mathematica来说明吧)，譬如数列：
CODE:
-1,1,-1,1,-1,1,-1,1,...

有 a(n)=(-1)^n
CODE:
FindSequenceFunction[{-1, 1, -1, 1, -1, 1, -1, 1}, n]

结果为
CODE:
(-1)^n

有了这个通项公式后，可以轻易的知道a(50)=1,a(51)=-1,...
但是如果用多项式拟合，
CODE:
InterpolatingPolynomial[{-1, 1, -1, 1, -1, 1, -1, 1}, x] // Expand

结果为
CODE:
-255 + (64864 x)/105 - (2856 x^2)/5 + (12056 x^3)/45 - 70 x^4 + (464 x^5)/45 - (4 x^6)/5 + (8 x^7)/315

如果另
CODE:
y(x)=-255 + (64864 x)/105 - (2856 x^2)/5 + (12056 x^3)/45 - 70 x^4 + (464 x^5)/45 - (4 x^6)/5 + (8 x^7)/315

可以知道y(1)到y(8)都是符合给定序列的，但是
y(50)=10158083585
y(51)=11832465919
就不符合了，所以多项式拟合不靠谱
希望我表达清楚了

[ Last edited by wangww2011 on 2011-9-9 at 20:15 ]
赞一下(2人) 回复此楼
11楼2011-09-09 20:14:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sudo

木虫 (正式写手)
★
余泽成(金币+1): 鼓励交流！ 2011-09-09 23:41:17
引用回帖:
11楼: Originally posted by wangww2011 at 2011-09-09 20:14:46:
不是这个意思，举个例子(这台电脑上的MATLAB被卸载了,就用Mathematica来说明吧)，譬如数列：
CODE:
-1,1,-1,1,-1,1,-1,1,...

有 a(n)=(-1)^n
[code]
FindSequenceFunction[{-1, 1, -1, 1, - ...

-1,1,-1,1,-1,1,-1,1,...

何以见得是(-1)^n？

或者，你能定义“最符合人类直觉”这种事情么？
赞一下(1人) 回复此楼
12楼2011-09-09 20:33:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangww2011

木虫 (著名写手)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
余泽成(金币+2): 鼓励讨论！ 2011-09-09 23:40:58
引用回帖:
12楼: Originally posted by sudo at 2011-09-09 20:33:32:
-1,1,-1,1,-1,1,-1,1,...

何以见得是(-1)^n？

或者，你能定义“最符合人类直觉”这种事情么？

这是个很open的话题，瞎扯几句吧
呵呵，我开始觉得是思维定势，可是为什么会有如此的思维定势？受训练的结果吗？
我觉得这个源之于人类的本能之一----控制欲，人们希望它是可控的----需要有规律才能被控制----可预测性(具有周期性的最简单，好预测)。
其实无论是(-1)^n(这里可以用sin或者其他任意函数扩充,但要至少满足周期性)或者是多项式拟合都具有可预测性，只不过预测的结果不同罢了。对于-1,1,-1,1,-1,1,-1,1,...这个数列，这里又多了一项本能的追求---对称性，如果画出图就会发现，多项式不满足周期性或者对称性。

而对于其他更普通的数列，不会严格满足周期性，但是人类还是会本能的希望其某些方面尽量满足，譬如震荡但振幅增加----我们会希望振幅增加着但是也要一直震荡着-----如果是多项式拟合的话，当x很大时，总是最高幂项起到关键作用，震荡会消失------这可能是人们本能上不想看到的吧。

[ Last edited by wangww2011 on 2011-9-9 at 22:07 ]
赞一下(2人) 回复此楼
13楼2011-09-09 20:38:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sudo

木虫 (正式写手)
★
余泽成(金币+1): 2011-09-10 20:26:39
引用回帖:
13楼: Originally posted by wangww2011 at 2011-09-09 20:38:11:
这是个很open的话题，瞎扯几句吧
呵呵，我开始觉得是思维定势，可是为什么会有如此的思维定势？受训练的结果吗？
我觉得这个源之于人类的本能之一----控制欲，人们希望它是可控的----需要有规律才能被控制- ...

=,=咳，取模的话，周期性完美出现...然后对称性也有了....

其实我才是那个刻板的教条主义者

[ Last edited by sudo on 2011-9-9 at 22:47 ]
赞一下(1人) 回复此楼
14楼2011-09-09 22:46:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangww2011

木虫 (著名写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
余泽成(金币+1): 鼓励交流！ 2011-09-10 20:26:56
引用回帖:
14楼: Originally posted by sudo at 2011-09-09 22:46:22:
=,=咳，取模的话，周期性完美出现...然后对称性也有了....

其实我才是那个刻板的教条主义者
[ Last edited by sudo on 2011-9-9 at 22:47 ]

呵呵 你这个取模的话会出现不连续性，
赞一下(2人) 回复此楼
15楼2011-09-10 10:01:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sudo

木虫 (正式写手)
★
余泽成(金币+1): 鼓励讨论！ 2011-09-10 20:28:35
引用回帖:
15楼: Originally posted by wangww2011 at 2011-09-10 10:01:06:
呵呵 你这个取模的话会出现不连续性，

确实，把定义域只能是整数域了

不过实话讲，对于数列题，定义域能在R内能连续的“通项公式”其实也不多吧

比如

(-1)^n = exp(iπ)^n = exp(inπ)

并不能在R上连续啊，只能在C上连续~

PS：我又想了下...莫非你所指的“连续性”不是实变函数论中的定义？

[ Last edited by sudo on 2011-9-10 at 10:11 ]
赞一下(1人) 回复此楼
16楼2011-09-10 10:10:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangww2011

木虫 (著名写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
余泽成(金币+1): 鼓励交流！ 2011-09-10 20:28:50
引用回帖:
16楼: Originally posted by sudo at 2011-09-10 10:10:03:
确实，把定义域只能是整数域了

不过实话讲，对于数列题，定义域能在R内能连续的“通项公式”其实也不多吧

比如

(-1)^n = exp(iπ)^n = exp(inπ)

并不能在R上连续啊，只能在C上连续~

...

我说的就是实数域，像你前面提到的，不需要用(-1)^n，可以用三角函数，其实也就是Re((-1)^n)

差点搞混了，譬如复函数f(x)=exp(I*pi*x)中定义域可以是实数域啊，这个复函数在R上是连续的啊
赞一下(2人) 回复此楼
17楼2011-09-10 10:50:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sudo

木虫 (正式写手)
★
余泽成(金币+1): 鼓励交流！ 2011-09-10 20:29:01
引用回帖:
17楼: Originally posted by wangww2011 at 2011-09-10 10:50:08:
我说的就是实数域，像你前面提到的，不需要用(-1)^n，可以用三角函数，其实也就是Re((-1)^n)

差点搞混了，譬如复函数f(x)=exp(I*pi*x)中定义域可以是实数域啊，这个复函数在R上是连续的啊

啊说错了...确实exp(i*pi*x)是定义域R内连续的...

其实我当时想表达的是另一个意思，数列题好像没见过有出现值域在复数域的东西...

就是一般直接讨论

f(x)=(-1)^x

的定义域的话，都是建立在f(x)定义在R上的前提再讨论x的吧？

[ Last edited by sudo on 2011-9-10 at 10:58 ]
赞一下(1人) 回复此楼
18楼2011-09-10 10:56:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangww2011

木虫 (著名写手)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
xzhdty(金币+2): 中秋快乐，欢迎常来 2011-09-10 19:45:09
引用回帖:
18楼: Originally posted by sudo at 2011-09-10 10:56:27:
啊说错了...确实exp(i*pi*x)是定义域R内连续的...

其实我当时想表达的是另一个意思，数列题好像没见过有出现值域在复数域的东西...

就是一般直接讨论

f(x)=(-1)^x

的定义域的话，都是建立在f ...

呵呵  好像也见过带虚数的题 两个一组的计算 那么多出题的人 什么怪题出不出来呀 哎
赞一下(2人) 回复此楼
19楼2011-09-10 11:02:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sudo

木虫 (正式写手)
引用回帖:
19楼: Originally posted by wangww2011 at 2011-09-10 11:02:39:
呵呵  好像也见过带虚数的题 两个一组的计算 那么多出题的人 什么怪题出不出来呀 哎

=,=居然连复数域的题都有...搞死人...

话说这种题意义几乎为零，不知道为什么考公务员要考
赞一下 回复此楼
20楼2011-09-10 11:08:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 sudo 的主题更新
292/3返回列表上一页123下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 336材料与化工085600求调剂 +7 水星记infp 2026-04-05 7/350 2026-04-06 00:57 by fmesaito
[考研] 材料调剂 +14 壹贰贰亿 2026-04-04 14/700 2026-04-05 23:31 by 来看流星雨10
[考研] 280求调剂 +4 李rien 2026-04-04 4/200 2026-04-05 18:44 by imissbao
[考研] 284求调剂 +7 徐同学_001 2026-04-04 13/650 2026-04-05 17:19 by yulian1987
[考研] 一志愿北交大材料工程总分358求调剂 +6 cs0106 2026-04-05 6/300 2026-04-05 16:34 by imissbao
[考研] 调剂求助 +10 想换手机不想解� 2026-04-02 13/650 2026-04-05 09:41 by sam3303
[考研] 413求调剂 +4 柯某某 2026-03-31 4/200 2026-04-04 22:18 by 学员6BFVa3
[考研] 一志愿211，化学学硕，310分，本科重点双非，求调剂 +11 努力奋斗112 2026-04-04 11/550 2026-04-04 20:51 by 蓝云思雨
[考研] 085601，一志愿厦大334复试被刷求调剂 +13 曾仰之 2026-04-03 15/750 2026-04-04 20:13 by dongzh2009
[论文投稿] 求文献 5+3 ys879651$ 2026-04-02 3/150 2026-04-04 17:22 by bobvan
[考研] 293分求调剂，外语为俄语 +6 加一一九 2026-03-31 6/300 2026-04-04 14:57 by 聪明的大松鼠
[考研] 材料科学与工程考研 +10 拯救皮特托先生 2026-04-02 10/500 2026-04-03 23:57 by userper
[考研] 一志愿华中农业071010，总分320求调剂 +7 困困困困坤坤 2026-04-02 7/350 2026-04-03 17:26 by Yuena_Wang
[考研] 初试成绩337找调剂 +3 ??? ?. ? 2026-04-03 3/150 2026-04-03 11:43 by 土木硕士招生
[考研] 279求调剂 +6 qazplm0852 2026-04-02 6/300 2026-04-03 10:03 by 蓝云思雨
[考研] 调剂 +3 好好读书。 2026-04-01 6/300 2026-04-02 15:49 by liumengping
[考研] 能源动力 调剂 +3 不破不立0 2026-04-02 3/150 2026-04-02 12:46 by ffffjjjj
[考研] 一志愿9初试366 本双非求调剂 +4 运气来得若有似� 2026-04-02 4/200 2026-04-02 09:56 by guanxin1001
[考研] 考研调剂 +12 Amber00 2026-03-31 12/600 2026-04-02 09:04 by sanrepian
[考研] 求调剂，一志愿南京师范大学计算机专硕，初试373，六级通过， +3 计算机追梦人 2026-04-01 3/150 2026-04-02 07:57 by fxue1114
信息提示
关闭
请填处理意见
关闭