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

24小时热门版块排行榜    

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

sudo

木虫 (正式写手)

[交流] 【致H兄】数列题的万能多项式拟合求法 已有2人参与

中秋快到了无心工作，就娱乐一下吧~

不知道是怎么回事的虫友请看这个帖子：
http://muchong.com/bbs/viewthread.php?tid=3573558&fpage=0&view=&highlight=&page=1

分析也在那个帖子里面啦，就不多说了，放MATLAB代码：
CODE:
function av = myPolyFit(y)
%A * av = y
    n = length(y);
    y = y(:); %make `y` a column vector
   
    %construct `A`
    A = (1:n)' * ones(1,n);
    for j = n-2 :-1: 1
        A(:,j) = A(:,j) .* A(:,j+1);
    end
    A(:,n) = ones(n,1);
   
    av = A\y;

其中y是序列，av是所找到的多项式的系数。

就以H兄的题目为例子吧
22   26   （）   23   24    21  （）    18

括号里面的数据可以乱填，比如这里就都填20吧，那么就有这段MATLAB程序：
CODE:
y=[22 26 20 23 24 21 20 18];
av=myPolyFit(y);
n=1:0.001:9;
val=polyval(av, n);
poly2str(av, 'n')
plot(n, val)

输出的多项式为：
CODE:
   0.0047619 n^7 - 0.18889 n^6 + 3.025 n^5 - 25.1806 n^4 + 116.3583 n^3
   - 293.6306 n^2 + 365.6119 n - 144

输出的图为：


[ Last edited by sudo on 2011-9-9 at 16:13 ]
回复此楼

» 本帖已获得的红花（最新10朵）

huycwork

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-09-09 16:11:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huycwork

金虫 (著名写手)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
xzhdty(金币+2): 欢迎常来程序语言看看 2011-09-09 23:08:48
余泽成: 你们都是高手，惺惺相惜啊，哈哈！ 2011-09-09 23:38:20
侬这动静闹得这般震撼，风口浪尖儿的
一般情况下如果等差数列，比如1 2 3 4 5 6。
俺表示淡定，说，这个用一个项就可以得出答案了！
Matlab表示淡定，想，sudo同学给了俺一个段儿呢！
计算机表示淡定，合计着，啥时把Matlab调过来，请sudo吃个月饼，一顿饭就解决了！
广大虫友表示淡定，这有啥，这方法不是那儿有嘛！
出考题的老师表示不够淡定，啥？六次方！
赞一下(3人) 回复此楼
漩涡的中心有一块空地，空空的。
2楼2011-09-09 17:51:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangww2011

木虫 (著名写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
xzhdty(金币+1): 中秋快乐！ 2011-09-09 23:09:05
引用回帖:
1楼: Originally posted by sudo at 2011-09-09 16:11:12:
中秋快到了无心工作，就娱乐一下吧~

不知道是怎么回事的虫友请看这个帖子：
http://muchong.com/bbs/viewthread.php?tid=3573558&fpage=0&view=&highlight=&page=1

分析也在那 ...

中秋快乐！

你这个纯多项式拟合不靠谱，因为这是局部的拟合，而实际上那个数列可以无限长的啊
赞一下(2人) 回复此楼
3楼2011-09-09 18:03:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huycwork

金虫 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖
送鲜花一朵
余泽成(程序强帖+1): 鼓励交流！ 2011-09-09 23:39:00
依我看，待定系数法就差不离了，能解决90%以上此类题，分两类，一类是指数式：
f2n = ab^n
另外一类是多项式，这个最高项由数据的多数决定，上面的式子占去两个参数，剩下的就归它：
f1n = an^3+bn^2+cn+d
将两类混合在一起建立方程，就是：
fn = ab^n+cn^3+dn^2+en+f
然后代入自然数列出求解矩阵，剩下的就是解方程组的工夫了。
赞一下(2人) 回复此楼
漩涡的中心有一块空地，空空的。
4楼2011-09-09 18:39:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sudo

木虫 (正式写手)
★ ★
余泽成(金币+2): 鼓励交流！ 2011-09-09 23:39:13
引用回帖:
2楼: Originally posted by huycwork at 2011-09-09 17:51:44:
侬这动静闹得这般震撼，风口浪尖儿的
一般情况下如果等差数列，比如1 2 3 4 5 6。
俺表示淡定，说，这个用一个项就可以得出答案了！
Matlab表示淡定，想，sudo同学给了俺一个段儿呢！
计算机表示淡定，合 ...

原本就符合多项式规律的数列，算出来就是对的啦~

等差数列的话（如1 2 3 4 5 6），最后算出来就是a(n)=n，如果用符号运算的话

MATLAB会因为精度问题保留高次项，不过那些系数都相当小就是~
CODE:
   1.8575e-018 n^5 - 2.8578e-017 n^4 + 1.6489e-016 n^3 - 4.2998e-016 n^2
   + 1 n - 2.0817e-016

其实只有一个'n'
赞一下(1人) 回复此楼
5楼2011-09-09 19:40:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sudo

木虫 (正式写手)
★ ★
余泽成(金币+2): 鼓励交流！ 2011-09-09 23:39:22
引用回帖:
3楼: Originally posted by wangww2011 at 2011-09-09 18:03:00:
中秋快乐！

你这个纯多项式拟合不靠谱，因为这是局部的拟合，而实际上那个数列可以无限长的啊

这个拟合出来的式子也可以无限长啊

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

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

f(n) = F(n mod N)

不就能满足那个整数的要求了嘛~
赞一下(1人) 回复此楼
6楼2011-09-09 19:42:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sudo

木虫 (正式写手)
★ ★
余泽成(金币+2): 鼓励交流！ 2011-09-09 23:39:39
引用回帖:
4楼: Originally posted by huycwork at 2011-09-09 18:39:45:
依我看，待定系数法就差不离了，能解决90%以上此类题，分两类，一类是指数式：
f2n = ab^n
另外一类是多项式，这个最高项由数据的多数决定，上面的式子占去两个参数，剩下的就归它：
f1n = an^3+bn^2+cn+d
将 ...

我弄这个主要是为了表示对这种题目的无语！

如果面试的时候被考官问到类似的题目，可以先详细地问问他“有什么具体的限定条件没有”

如果他不加上诸如“那个通项公式的定义域如果换成实数域，那么它在实数域里面连续可导”这些BT的条件的话，这个可以算通用解法之一啦，哈哈哈
赞一下(1人) 回复此楼
7楼2011-09-09 19:47:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sudo

木虫 (正式写手)
哦，差点忘了，祝各位中秋快乐~~~
赞一下 回复此楼
8楼2011-09-09 20:06:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huycwork

金虫 (著名写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
余泽成(金币+1): 鼓励交流！ 2011-09-09 23:40:01
引用回帖:
7楼: Originally posted by sudo at 2011-09-09 19:47:19:
我弄这个主要是为了表示对这种题目的无语！

如果面试的时候被考官问到类似的题目，可以先详细地问问他“有什么具体的限定条件没有”

如果他不加上诸如“那个通项公式的定义域如果换成实数域，那么它在实 ...

其实我觉得现在好多面试笔试也是扯蛋的，有些算法或者什么观察规律的，反应过来了就很容易，没反应过来就麻烦了，但是现在那些面试官又好这个，很多人也专门研究这些题目，比如那个什么1块蛋糕切三刀五人分或者。唉，这不是没办法的事情嘛，一次没进去就得等明年了。。。。
不过值得庆幸的是还好没有比如鸽巢原理之类的数学应用的先例，要不然就玩完了~
赞一下(2人) 回复此楼
漩涡的中心有一块空地，空空的。
9楼2011-09-09 20:09:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huycwork

金虫 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖
引用回帖:
8楼: Originally posted by sudo at 2011-09-09 20:06:32:
哦，差点忘了，祝各位中秋快乐~~~

大家都提前预祝哈~大家同乐，独乐乐，不如与众~
赞一下(1人) 回复此楼
漩涡的中心有一块空地，空空的。
10楼2011-09-09 20:13:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 sudo 的主题更新
291/3返回列表123下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿郑州大学085600求调剂 +17 吃的不少 2026-04-05 20/1000 2026-04-06 00:32 by T可可西里T
[考研] 机械专硕274求调剂，不挑专业学校 +6 泛泛2333 2026-04-05 7/350 2026-04-05 23:20 by chyhaha
[考研] 081200-11408-276学硕求调剂 +4 崔wj 2026-04-05 4/200 2026-04-05 20:29 by 啵啵啵0119
[考研] 求调剂求调剂 +8 121. 2026-04-02 8/400 2026-04-05 20:15 by lys0704
[考研] 301求调剂 +3 XYPLR 2026-04-05 4/200 2026-04-05 19:07 by XYPLR
[考研] 0817化学工程与技术求调剂，一志愿中海洋319 +13 lv945 2026-04-04 13/650 2026-04-05 18:14 by 猪会飞
[考研] 283分求调剂 +9 试试看呗 2026-04-04 9/450 2026-04-05 10:27 by 果冻大王
[考研] 一志愿北京2，材料与化工308求调剂 +10 熊二想上岸 2026-04-04 10/500 2026-04-05 05:20 by houyaoxu
[考研] 考研调剂 +5 四川王涛 2026-04-04 5/250 2026-04-04 22:18 by 啵啵啵0119
[考研] 368求调剂 +5 今华习 2026-04-03 7/350 2026-04-04 18:47 by imissbao
[考研] 319求调剂 +4 星星不眨眼喽 2026-04-03 4/200 2026-04-04 16:25 by 中飞院空管学院�
[考研] 兽医调剂 +3 wh119216 2026-04-02 3/150 2026-04-03 19:34 by zrongyan
[考研] 315分 085602 求调剂 +15 26考研上岸版26 2026-04-02 15/750 2026-04-03 12:45 by xingguangj
[考研] 326求调剂 +10 崽崽仔 2026-04-02 10/500 2026-04-03 09:08 by 帕尔马拉特
[考研] 362求调剂 +14 西南交材料专硕3 2026-03-31 14/700 2026-04-02 17:50 by yunlongyang
[考研] 354求调剂 +4 lxb598 2026-03-31 5/250 2026-04-02 09:55 by Jaylen.
[考研] 279求调剂 +6 学而思兮知 2026-04-01 6/300 2026-04-02 09:16 by vgtyfty
[考研] 08工科275求调剂，可跨考。 +5 AaAa7420 2026-03-31 5/250 2026-04-01 15:21 by 159357hjz
[考研] 一志愿西交大080500材料学硕349 +6 jqx1258 2026-03-31 7/350 2026-03-31 21:08 by yuq
[考研] 复试调剂 +7 双马尾痞老板2 2026-03-31 7/350 2026-03-31 19:49 by Dyhoer
信息提示
关闭
请填处理意见
关闭