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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 求用最小二乘拟合时间序列
51/1返回列表
查看: 2160  |  回复: 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的回帖

shenying1215

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

sg47

铜虫 (小有名气)

【答案】应助回帖

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

feixiaolin

荣誉版主 (文坛精英)

优秀版主
镜像延拓至512长度，做FFT,参数就出来了，不必死硬拟合。
赞一下(1人) 回复此楼
4楼2013-11-25 13:21:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料专硕306英一数二 +8 z1z2z3879 2026-03-16 10/500 2026-03-18 00:01 by 星空星月
[考研] 302求调剂 +6 呼呼呼。。。。 2026-03-17 6/300 2026-03-17 23:56 by 星空星月
[考博] 26博士申请 +3 1042136743 2026-03-17 3/150 2026-03-17 23:30 by 轻松不少随
[考研] 考研求调剂 +3 橘颂. 2026-03-17 4/200 2026-03-17 21:43 by 有只狸奴
[考研] 302求调剂 +9 负心者当诛 2026-03-11 9/450 2026-03-17 17:13 by ruiyingmiao
[考研] 26考研求调剂 +6 丶宏Sir 2026-03-13 6/300 2026-03-17 16:13 by 醉在风里
[考研] 一志愿苏州大学材料工程（085601）专硕有科研经历三项国奖两个实用型专利一项省级立项 +6 大火山小火山 2026-03-16 8/400 2026-03-17 15:05 by 无懈可击111
[考研] 085600材料与化工 +4 安全上岸！ 2026-03-16 4/200 2026-03-17 14:02 by 勇敢太监王公公
[考研] 321求调剂 +5 大米饭！ 2026-03-15 5/250 2026-03-16 16:33 by houyaoxu
[考研] 0703一志愿211 285分求调剂 +5 ly3471z 2026-03-13 5/250 2026-03-16 16:16 by 哦哦123
[考研] 一志愿华中师范071000，325求调剂 +6 RuitingC 2026-03-12 6/300 2026-03-16 14:50 by 可淡不可忘
[考研] 0703化学调剂 290分有科研经历，论文在投 +7 腻腻gk 2026-03-14 7/350 2026-03-16 10:12 by houyaoxu
[考研] 本科南京大学一志愿川大药学327 +3 麦田耕者 2026-03-14 3/150 2026-03-14 20:04 by 外星文明
[基金申请] 现在如何回避去年的某一个专家，不知道名字 +3 zk200107 2026-03-12 6/300 2026-03-14 17:13 by zk200107
[考研] 一志愿哈工大材料324分求调剂 +5 闫旭东 2026-03-14 5/250 2026-03-14 14:53 by 木瓜膏
[考研] 330求调剂 +3 ?酱给调剂跪了 2026-03-13 3/150 2026-03-14 10:13 by JourneyLucky
[考研] 0703化学调剂 +4 快乐的香蕉 2026-03-11 4/200 2026-03-13 22:41 by JourneyLucky
[考研] 308求调剂 +5 是Lupa啊 2026-03-11 5/250 2026-03-13 22:13 by JourneyLucky
[考研] 求b区学校调剂 +3 周56 2026-03-11 3/150 2026-03-13 16:20 by JourneyLucky
[考研] 工科调剂 +4 Jiang191123！ 2026-03-11 4/200 2026-03-13 15:15 by Miko19
信息提示
关闭
请填处理意见
关闭