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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 求用最小二乘拟合时间序列
北京石油化工学院2026年研究生招生接收调剂公告
41/1返回列表
查看: 2176  |  回复: 3
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

梧桐雨9008

新虫 (初入文坛)

[求助] 求用最小二乘拟合时间序列

现有一组时间序列，想对其进行最小二乘拟合，确定其参数，参数模型是：
y=a+b*x+c*sin((pi/14)*x+d)+e*sin((pi/947)*x+f)+g*sin((pi/2016)*x+h)+i*sin((pi/9624)*x+j)+k*sin((pi/20381)*x+l),现在想拟合出a,b,c,d,e,f,g,h,i,j,k,l 这十二个参数，求帮忙，试了好几种方法都算不出来好结果。
数据：
x     y
1        5
2        11
3        16
4        23
5        36
6        58
7        29
8        20
9        10
10        8
11        3
12        0
13        0
14        2
15        11
16        27
17        47
18        63
19        60
20        39
21        28
22        26
23        22
24        11
25        21
26        40
27        78
28        122
29        103
30        73
31        47
32        35
33        11
34        5
35        16
36        34
37        70
38        81
39        111
40        101
41        73
42        40
43        20
44        16
45        5
46        11
47        22
48        40
49        60
50        80.9
51        83.4
52        47.7
53        47.8
54        30.7
55        12.2
56        9.6
57        10.2
58        32.4
59        47.6
60        54
61        62.9
62        85.9
63        61.2
64        45.1
65        36.4
66        20.9
67        11.4
68        37.8
69        69.8
70        106.1
71        100.8
72        81.6
73        66.5
74        34.8
75        30.6
76        7
77        19.8
78        92.5
79        154.4
80        125.9
81        84.8
82        68.1
83        38.5
84        22.8
85        10.2
86        24.1
87        82.9
88        132
89        130.9
90        118.1
91        89.9
92        66.6
93        60
94        46.9
95        41
96        21.3
97        16
98        6.4
99        4.1
100        6.8
101        14.5
102        34
103        45
104        43.1
105        47.5
106        42.2
107        28.1
108        10.1
109        8.1
110        2.5
111        0
112        1.4
113        5
114        12.2
115        13.9
116        35.4
117        45.8
118        41.1
119        30.1
120        23.9
121        15.6
122        6.6
123        4
124        1.8
125        8.5
126        16.6
127        36.3
128        49.6
129        64.2
130        67
131        70.9
132        47.8
133        27.5
134        8.5
135        13.2
136        56.9
137        121.5
138        138.3
139        103.2
140        85.7
141        64.6
142        36.7
143        24.2
144        10.7
145        15
146        40.1
147        61.5
148        98.5
149        124.7
150        96.3
151        66.6
152        64.5
153        54.1
154        39
155        20.6
156        6.7
157        4.3
158        22.7
159        54.8
160        93.8
161        95.8
162        77.2
163        59.1
164        44
165        47
166        30.5
167        16.3
168        7.3
169        37.6
170        74
171        139
172        111.2
173        101.6
174        66.2
175        44.7
176        17
177        11.3
178        12.4
179        3.4
180        6
181        32.3
182        54.3
183        59.7
184        63.7
185        63.5
186        52.2
187        25.4
188        13.1
189        6.8
190        6.3
191        7.1
192        35.6
193        73
194        85.1
195        78
196        64
197        41.8
198        26.2
199        26.7
200        12.1
201        9.5
202        2.7
203        5
204        24.4
205        42
206        63.5
207        53.8
208        62
209        48.5
210        43.9
211        18.6
212        5.7
213        3.6
214        1.4
215        9.6
216        47.4
217        57.1
218        103.9
219        80.6
220        63.6
221        37.6
222        26.1
223        14.2
224        5.8
225        16.7
226        44.3
227        63.9
228        69
229        77.8
230        64.9
231        35.7
232        21.2
233        11.1
234        5.7
235        8.7
236        36.1
237        79.7
238        114.4
239        109.6
240        88.8
241        67.8
242        47.5
243        30.6
244        16.3
245        9.6
246        33.2
247        92.6
248        151.6
249        136.3
250        134.7
251        83.9
252        69.4
253        31.5
254        13.9
255        4.4
256        38
257        141.7
258        190.2
259        184.8
260        159
261        112.3
262        53.9
263        37.6
264        27.9
265        10.2
266        15.1
267        47
268        93.8
269        105.9
270        105.5
271        104.5
272        66.6
273        68.9
274        38
275        34.5
276        15.5
277        12.6
278        27.5
279        92.5
280        155.4
281        154.6
282        140.4
283        115.9
284        66.6
285        45.9
286        17.9
287        13.4
288        29.4
289        100.2
290        157.6
291        142.6
292        145.7
293        94.3
294        54.6
295        29.9
296        17.5
297        8.6
298        21.5
299        64.3
300        93.3
301        119.6
302        111
303        104
304        63.7
305        40.4
306        29.8
307        15.2
308        7.5
309        2.9
310        3.1
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-10-13 20:57:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sg47

铜虫 (小有名气)

【答案】应助回帖

★
感谢参与，应助指数 +1
dbb627: 金币+1, 感谢参与讨论 2013-10-14 09:42:54
用matlab不就行了，最小化目标函数
赞一下 回复此楼
2楼2013-10-13 23:08:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shenying1215

铜虫 (小有名气)
用R软件，季节模型
赞一下(1人) 回复此楼
3楼2013-11-24 21:31:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主
镜像延拓至512长度，做FFT,参数就出来了，不必死硬拟合。
赞一下(1人) 回复此楼
4楼2013-11-25 13:21:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 梧桐雨9008 的主题更新
41/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料专硕322 +4 哈哈哈吼吼吼哈 2026-04-05 4/200 2026-04-05 23:24 by 来看流星雨10
[考研] 一志愿北京化工085600 310分求调剂 +14 0856材料与化工3 2026-04-04 16/800 2026-04-05 22:57 by Hdyxbekcb
[考研] 化学调剂 +15 艾志恒 2026-04-03 16/800 2026-04-05 22:36 by syh9288
[考研] 377求调剂 +6 by.ovo 2026-04-05 6/300 2026-04-05 22:18 by dongzh2009
[考研] 285求调剂 +11 哦呦呼o 2026-04-04 11/550 2026-04-05 17:59 by 猪会飞
[考研] 318求调剂 +11 ykyhsa 2026-04-05 13/650 2026-04-05 12:44 by aidndnjck
[考研] 材料调剂 +11 一样YWY 2026-04-02 13/650 2026-04-04 23:10 by 无际的草原
[考研] 325求调剂 +4 春风不借意 2026-04-04 4/200 2026-04-04 14:46 by 湘农储能材料
[考研] 一志愿北京科技大学材料工程085601，求调剂 +17 cdyw 2026-04-02 18/900 2026-04-04 11:14 by w_xuqing
[考研] 305求调剂 +3 77Qi 2026-04-03 3/150 2026-04-03 23:01 by qzxyhcsy
[考研] 化工求调剂 +11 荔香芝士椰奶 2026-04-03 11/550 2026-04-03 22:06 by 啵啵啵0119
[考研] 310求调剂 +18 争取九点睡 2026-03-30 18/900 2026-04-03 18:35 by ls刘帅
[考研] 求调剂 +8 akdhjs 2026-04-03 8/400 2026-04-03 18:17 by 戴维ING
[考研] 考研调剂 +8 不爱喝饮料 2026-04-03 8/400 2026-04-03 16:40 by Mistake-J
[考研] 303求调剂 +3 一色清羽 2026-04-02 4/200 2026-04-03 10:22 by 蓝云思雨
[考研] 调剂 +3 好好读书。 2026-04-01 6/300 2026-04-02 15:49 by liumengping
[考研] 275学硕081000服从调剂到其他专业，保不住本专业了 +7 一只小小水牛 2026-04-02 8/400 2026-04-02 14:23 by alice-2022
[考研] 286分调剂 +20 Faune 2026-03-30 22/1100 2026-04-02 13:24 by clyblh
[考研] 085410 一志愿211 22408分数359求调剂 +3 123456789qw 2026-03-31 4/200 2026-04-02 00:06 by 义文wang
[考研] 江苏科技大学招材料研究生 +4 Su032713. 2026-04-01 5/250 2026-04-01 22:03 by cccchenso
信息提示
关闭
请填处理意见
关闭