質問編集履歴
1
データの追加
test
CHANGED
|
File without changes
|
test
CHANGED
|
@@ -1,13 +1,14 @@
|
|
|
1
1
|
最小二乗法で一次関数をデータにフィッティングしたいです。
|
|
2
|
-
具体的には
|
|
2
|
+
具体的には下においてます、date.txtのデータです。
|
|
3
3
|
これをその一つ上のfitting.cppでフィッティングしたのですが、フィッティングできていません。
|
|
4
4
|
fitting.cppは
|
|
5
5
|
https://www.mk-mode.com/blog/2014/03/02/cpp-least-squares-method/
|
|
6
6
|
こちらのサイトを参考にさせていただいております。
|
|
7
7
|
実際に、他のデータで試してみるとフィッティングはうまくいきます。
|
|
8
|
-
しかし以下のdate.txtのデータ
|
|
8
|
+
しかし以下のdate.txtのデータではできません。
|
|
9
|
-
|
|
9
|
+
一次の係数として出力される値は約90です。
|
|
10
10
|
何が原因か、わかる方いらっしゃいましたらどうかよろしくお願いします。
|
|
11
|
+
[補足]date2.txtでも一次の係数として出力される値は約19で、うまくいきませんでした。
|
|
11
12
|
```fitting.txt
|
|
12
13
|
#include <iostream>
|
|
13
14
|
#include <stdio.h>
|
|
@@ -245,4 +246,186 @@
|
|
|
245
246
|
-125.825 -883.022
|
|
246
247
|
-125.699 -893.022
|
|
247
248
|
```
|
|
248
|
-
|
|
249
|
+
|
|
250
|
+
```date2.txt
|
|
251
|
+
-124.491 876.8
|
|
252
|
+
-124.626 866.801
|
|
253
|
+
-124.758 856.801
|
|
254
|
+
-124.889 846.802
|
|
255
|
+
-125.019 836.803
|
|
256
|
+
-125.146 826.804
|
|
257
|
+
-125.273 816.805
|
|
258
|
+
-125.399 806.806
|
|
259
|
+
-125.523 796.806
|
|
260
|
+
-125.645 786.807
|
|
261
|
+
-125.766 776.808
|
|
262
|
+
-125.884 766.809
|
|
263
|
+
-126.001 756.809
|
|
264
|
+
-126.116 746.81
|
|
265
|
+
-126.232 736.811
|
|
266
|
+
-126.345 726.811
|
|
267
|
+
-126.457 716.812
|
|
268
|
+
-126.567 706.812
|
|
269
|
+
-126.675 696.813
|
|
270
|
+
-126.781 686.814
|
|
271
|
+
-126.886 676.814
|
|
272
|
+
-126.991 666.815
|
|
273
|
+
-127.094 656.815
|
|
274
|
+
-127.195 646.816
|
|
275
|
+
-127.294 636.816
|
|
276
|
+
-127.392 626.817
|
|
277
|
+
-127.488 616.817
|
|
278
|
+
-127.582 606.818
|
|
279
|
+
-127.677 596.818
|
|
280
|
+
-127.769 586.818
|
|
281
|
+
-127.86 576.819
|
|
282
|
+
-127.949 566.819
|
|
283
|
+
-128.036 556.82
|
|
284
|
+
-128.122 546.82
|
|
285
|
+
-128.207 536.82
|
|
286
|
+
-128.291 526.821
|
|
287
|
+
-128.373 516.821
|
|
288
|
+
-128.453 506.821
|
|
289
|
+
-128.532 496.822
|
|
290
|
+
-128.609 486.822
|
|
291
|
+
-128.684 476.822
|
|
292
|
+
-128.759 466.823
|
|
293
|
+
-128.832 456.823
|
|
294
|
+
-128.904 446.823
|
|
295
|
+
-128.973 436.823
|
|
296
|
+
-129.041 426.824
|
|
297
|
+
-129.108 416.824
|
|
298
|
+
-129.172 406.824
|
|
299
|
+
-129.237 396.824
|
|
300
|
+
-129.3 386.824
|
|
301
|
+
-129.361 376.825
|
|
302
|
+
-129.42 366.825
|
|
303
|
+
-129.478 356.825
|
|
304
|
+
-129.534 346.825
|
|
305
|
+
-129.589 336.825
|
|
306
|
+
-129.644 326.825
|
|
307
|
+
-129.696 316.826
|
|
308
|
+
-129.747 306.826
|
|
309
|
+
-129.795 296.826
|
|
310
|
+
-129.843 286.826
|
|
311
|
+
-129.888 276.826
|
|
312
|
+
-129.933 266.826
|
|
313
|
+
-129.977 256.826
|
|
314
|
+
-130.019 246.826
|
|
315
|
+
-130.059 236.826
|
|
316
|
+
-130.097 226.826
|
|
317
|
+
-130.134 216.826
|
|
318
|
+
-130.171 206.827
|
|
319
|
+
-130.206 196.827
|
|
320
|
+
-130.239 186.827
|
|
321
|
+
-130.27 176.827
|
|
322
|
+
-130.3 166.827
|
|
323
|
+
-130.328 156.827
|
|
324
|
+
-130.354 146.827
|
|
325
|
+
-130.38 136.827
|
|
326
|
+
-130.405 126.827
|
|
327
|
+
-130.427 116.827
|
|
328
|
+
-130.448 106.827
|
|
329
|
+
-130.467 96.827
|
|
330
|
+
-130.485 86.827
|
|
331
|
+
-130.501 76.827
|
|
332
|
+
-130.516 66.827
|
|
333
|
+
-130.53 56.827
|
|
334
|
+
-130.542 46.827
|
|
335
|
+
-130.553 36.827
|
|
336
|
+
-130.562 26.827
|
|
337
|
+
-130.569 16.827
|
|
338
|
+
-130.576 6.82704
|
|
339
|
+
-130.581 -3.17296
|
|
340
|
+
-130.584 -13.173
|
|
341
|
+
-130.586 -23.173
|
|
342
|
+
-130.586 -33.173
|
|
343
|
+
-130.584 -43.173
|
|
344
|
+
-130.581 -53.173
|
|
345
|
+
-130.577 -63.173
|
|
346
|
+
-130.572 -73.173
|
|
347
|
+
-130.565 -83.173
|
|
348
|
+
-130.556 -93.1729
|
|
349
|
+
-130.546 -103.173
|
|
350
|
+
-130.534 -113.173
|
|
351
|
+
-130.52 -123.173
|
|
352
|
+
-130.504 -133.173
|
|
353
|
+
-130.488 -143.173
|
|
354
|
+
-130.471 -153.173
|
|
355
|
+
-130.452 -163.173
|
|
356
|
+
-130.431 -173.173
|
|
357
|
+
-130.408 -183.173
|
|
358
|
+
-130.384 -193.173
|
|
359
|
+
-130.357 -203.173
|
|
360
|
+
-130.331 -213.173
|
|
361
|
+
-130.303 -223.173
|
|
362
|
+
-130.274 -233.173
|
|
363
|
+
-130.242 -243.173
|
|
364
|
+
-130.209 -253.173
|
|
365
|
+
-130.174 -263.172
|
|
366
|
+
-130.138 -273.172
|
|
367
|
+
-130.099 -283.172
|
|
368
|
+
-130.061 -293.172
|
|
369
|
+
-130.021 -303.172
|
|
370
|
+
-129.979 -313.172
|
|
371
|
+
-129.935 -323.172
|
|
372
|
+
-129.89 -333.172
|
|
373
|
+
-129.843 -343.172
|
|
374
|
+
-129.794 -353.172
|
|
375
|
+
-129.743 -363.172
|
|
376
|
+
-129.691 -373.171
|
|
377
|
+
-129.638 -383.171
|
|
378
|
+
-129.584 -393.171
|
|
379
|
+
-129.528 -403.171
|
|
380
|
+
-129.471 -413.171
|
|
381
|
+
-129.411 -423.171
|
|
382
|
+
-129.35 -433.17
|
|
383
|
+
-129.288 -443.17
|
|
384
|
+
-129.223 -453.17
|
|
385
|
+
-129.157 -463.17
|
|
386
|
+
-129.089 -473.17
|
|
387
|
+
-129.021 -483.169
|
|
388
|
+
-128.951 -493.169
|
|
389
|
+
-128.879 -503.169
|
|
390
|
+
-128.806 -513.169
|
|
391
|
+
-128.731 -523.168
|
|
392
|
+
-128.654 -533.168
|
|
393
|
+
-128.575 -543.168
|
|
394
|
+
-128.495 -553.167
|
|
395
|
+
-128.413 -563.167
|
|
396
|
+
-128.329 -573.167
|
|
397
|
+
-128.244 -583.166
|
|
398
|
+
-128.158 -593.166
|
|
399
|
+
-128.071 -603.166
|
|
400
|
+
-127.982 -613.165
|
|
401
|
+
-127.891 -623.165
|
|
402
|
+
-127.799 -633.164
|
|
403
|
+
-127.705 -643.164
|
|
404
|
+
-127.609 -653.163
|
|
405
|
+
-127.511 -663.163
|
|
406
|
+
-127.411 -673.162
|
|
407
|
+
-127.31 -683.162
|
|
408
|
+
-127.207 -693.161
|
|
409
|
+
-127.102 -703.161
|
|
410
|
+
-126.996 -713.16
|
|
411
|
+
-126.888 -723.16
|
|
412
|
+
-126.778 -733.159
|
|
413
|
+
-126.666 -743.158
|
|
414
|
+
-126.553 -753.158
|
|
415
|
+
-126.439 -763.157
|
|
416
|
+
-126.324 -773.157
|
|
417
|
+
-126.207 -783.156
|
|
418
|
+
-126.088 -793.155
|
|
419
|
+
-125.968 -803.154
|
|
420
|
+
-125.846 -813.154
|
|
421
|
+
-125.722 -823.153
|
|
422
|
+
-125.596 -833.152
|
|
423
|
+
-125.469 -843.151
|
|
424
|
+
-125.34 -853.15
|
|
425
|
+
-125.209 -863.15
|
|
426
|
+
-125.076 -873.149
|
|
427
|
+
-124.942 -883.148
|
|
428
|
+
-124.806 -893.147
|
|
429
|
+
|
|
430
|
+
```
|
|
431
|
+
|