24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

saji

铁虫 (初入文坛)

应助: 1 (幼儿园)
金币: 33.7
帖子: 22
在线: 3.9小时
虫号: 3883640
注册: 2015-05-21
专业: 信号理论与信号处理

[求助] 1stOpt拟合代码的问题求助已有2人参与

各位前辈：
小弟使用1stopt拟合一等式，原始形式如下：
$y_{i}=A+B(t-x_{i})^m+C(t-x_{i})^mcos(wln(t-x_{i})+p)$
得到的拟合结果正常，图形如下：

但是，使用同样的样本，进行变式的拟合时，却得到怪异的图形。。。不知问题在哪里？
变式如下：
$(y_{i}-A)(t-x_{i})^{-m}=B+Ccos(wln(t-x_{i})+p)$
得到的函数图形如下：

下面是我使用变式的拟合代码：
Title "A-share 07-Bubble LPPL Fit";
Parameters a,b,c,t,m,w,p;
Variable y,x;
Function (y-a)*(t-x)^-m = b+c*cos(w*ln(t-x)+p);
Data;
6.941555559 1
6.938227256 2
7.017131521 3
7.030903452 4
7.010569957 5
7.008768258 6
6.997105918 7
6.978060888 8
6.990265709 9
6.98989732 10
7.017165583 11
7.004415456 12
7.004911935 13
6.997323551 14
7.004774906 15
7.025218263 16
7.010842412 17
7.007593374 18
6.985584462 19
6.961858921 20
6.953949696 21
6.946055376 22
6.94074992 23
6.945626043 24
6.925578498 25
6.919188668 26
6.952882552 27
6.944467078 28
6.947375912 29
6.933535105 30
6.919779695 31
6.922008222 32
6.928588747 33
6.928177322 34
6.953035481 35
6.952154866 36
6.978033857 37
6.993849485 38
6.990760977 39
6.987520717 40
6.992967372 41
7.006729647 42
7.01004017 43
7.00541633 44
7.028859903 45
7.037802358 46
7.049718276 47
7.060504692 48
7.076298178 49
7.062975387 50
7.079189449 51
7.070960273 52
7.084087267 53
7.045807064 54
7.047677209 55
7.054965069 56
7.047479829 57
7.06230902 58
7.066872178 59
7.066347509 60
7.051360262 61
7.044368792 62
7.058584449 63
7.077441512 64
7.08074005 65
7.086925208 66
7.067594321 67
7.085323054 68
7.085408456 69
7.08139341 70
7.080204935 71
7.096029086 72
7.104358532 73
7.104374962 74
7.100808215 75
7.10711994 76
7.100536113 77
7.080016399 78
7.056140801 79
7.04923748 80
7.051891986 81
7.032833418 82
7.031541592 83
7.052275448 84
7.052387083 85
7.037862064 86
7.05375166 87
7.057771423 88
7.049781612 89
7.038386094 90
7.031188986 91
7.039836496 92
7.034042496 93
7.035223704 94
7.039939014 95
7.039810207 96
7.022795 97
7.000481213 98
7.00104706 99
6.985518776 100
6.996514045 101
6.993887103 102
7.007412361 103
6.998758926 104
7.003106367 105
7.003656193 106
7.012247718 107
7.010449944 108
6.992373044 109
6.994111858 110
6.992696431 111
6.991644119 112
6.999322098 113
6.998792707 114
7.018400904 115
7.02102682 116
7.001841983 117
7.008275499 118
7.015150324 119
7.016534345 120
7.012856461 121
7.000321698 122
7.002393415 123
7.001925719 124
6.997858289 125
6.983972525 126
6.991907071 127
7.002715398 128
7.00153975 129
7.015242831 130
7.017834942 131
7.019158898 132
7.025882255 133
7.024253045 134
7.027765167 135
7.03152127 136
7.03556961 137
7.030644368 138
7.0345976 139
7.043047246 140
7.053432734 141
7.051238983 142
7.053615113 143
7.064641072 144
7.057086076 145
7.074085486 146
7.087798409 147
7.097897839 148
7.103048999 149
7.107112567 150
7.099245508 151
7.112086177 152
7.107801325 153
7.092464816 154
7.097083896 155
7.117489323 156
7.132158836 157
7.13513624 158
7.135507395 159
7.132547071 160
7.137315003 161
7.160555492 162
7.156251518 163
7.162440132 164
7.146342166 165
7.156688207 166
7.154334067 167
7.159546149 168
7.169477799 169
7.14727048 170
7.144729149 171
7.144834871 172
7.161171164 173
7.157911481 174
7.161503299 175
7.167705864 176
7.167952582 177
7.169373111 178
7.175172908 179
7.159032151 180
7.164950051 181
7.161587092 182
7.138800331 183
7.131205983 184
7.127024133 185
7.127415169 186
7.138608236 187
7.138107982 188
7.150553996 189
7.150065962 190
7.146348468 191
7.161170388 192
7.162703652 193
7.167582482 194
7.17200776 195
7.166037377 196
7.166972288 197
7.169189111 198
7.174392708 199
7.16604819 200
7.168807144 201
7.184986178 202
7.19278233 203
7.200545781 204
7.200228605 205
7.202635135 206
7.214564015 207
7.216877608 208
7.215337033 209
7.194684568 210
7.214900215 211
7.228833926 212
7.233537726 213
7.241869022 214
7.234105027 215
7.256151146 216
7.251193199 217
7.244061073 218
7.25663381 219
7.25610809 220
7.272553242 221
7.311287854 222
7.33378155 223
7.343223751 224
7.337832996 225
7.379527344 226
7.417032505 227
7.386006677 228
7.393352937 229
7.388501623 230
7.414298747 231
7.413179743 232
7.380599869 233
7.372070258 234
7.37238512 235
7.386405171 236
7.407644114 237
7.41293781 238
7.40324389 239
7.429042984 240
7.42021896 241
7.429297672 242
7.426028699 243
7.371203091 244
7.372425963 245
7.346902715 246
7.347813744 247
7.34503877 248
7.333889959 249
7.335620944 250
7.361672081 251
7.369150714 252
7.372954887 253
7.376580714 254
7.375297131 255
7.381321928 256
7.398449009 257
7.402177703 258
7.40201606 259
7.421547299 260
7.421901982 261
7.436783427 262
7.427469482 263
7.4492391 264
7.462484863 265
7.455986508 266
7.458375874 267
7.464971619 268
7.461640392 269
7.412023248 270
7.41770532 271
7.42826605 272
7.429076827 273
7.405592923 274
7.411630001 275
7.41777738 276
7.418147209 277
7.429794396 278
7.430501964 279
7.423668738 280
7.415795627 281
7.385685534 282
7.378142585 283
7.37830376 284
7.37841994 285
7.358926435 286
7.344354647 287
7.365545249 288
7.364272119 289
7.381542367 290
7.381455816 291
7.359300852 292
7.375040843 293
7.38796292 294
7.379839243 295
7.37651939 296
7.3784774 297
7.386071141 298
7.38547655 299
7.392042042 300
7.392052516 301
7.408795986 302
7.40915098 303
7.411674106 304
7.413752063 305
7.4004318 306
7.413077788 307
7.417026495 308
7.421848159 309
7.415289492 310
7.419370991 311
7.423548146 312
7.435946445 313
7.432122204 314
7.432301543 315
7.450686943 316
7.457601788 317
7.458901588 318
7.457293602 319
7.462158074 320
7.453191583 321
7.451454206 322
7.445789567 323
7.453004938 324
7.459893465 325
7.468755251 326
7.487389357 327
7.48709638 328
7.490028998 329
7.483346211 330
7.486980959 331
7.479468455 332
7.475283113 333
7.488397035 334
7.491228767 335
7.49017423 336
7.472724711 337
7.498416143 338
7.499227367 339
7.501438967 340
7.499520132 341
7.500892051 342
7.516429494 343
7.526024268 344
7.523650936 345
7.53174636 346
7.542376228 347
7.5440834 348
7.531990121 349
7.54775587 350
7.540805792 351
7.530353914 352
7.543399928 353
7.561607944 354
7.571240357 355
7.58669856 356
7.609504853 357
7.619504366 358
7.621369082 359
7.631607692 360
7.625991091 361
7.624265896 362
7.620077929 363
7.627588693 364
7.649353995 365
7.650667862 366
7.678626173 367
7.683995127 368
7.676290041 369
7.676357738 370
7.646658502 371
7.687307653 372
7.704789841 373
7.706817571 374
7.718289861 375
7.729255216 376
7.754667637 377
7.768185249 378
7.771996641 379
7.759160976 380
7.759472502 381
7.798014923 382
7.815905769 383
7.838495906 384
7.850724553 385
7.89188184 386
7.906812022 387
7.879039372 388
7.9036698 389
7.940157962 390
7.946467517 391
7.92664231 392
7.889125636 393
7.935480403 394
7.944852737 395
7.929810445 396
7.921892245 397
7.948811545 398
7.983845847 399
7.989270238 400
7.998042678 401
7.95765515 402
7.966433373 403
7.987953396 404
7.982949492 405
7.932482391 406
7.932158257 407
7.891036025 408
7.868077059 409
7.891965186 410
7.90697992 411
7.914885119 412
7.91219827 413
7.939933562 414
7.948693256 415
7.974220683 416
8.004034181 417
8.005858772 418
8.019809814 419
7.927248961 420
7.965918073 421
7.936370621 422
7.948571068 423
7.932113021 424
7.951620949 425
7.971290843 426
7.982079998 427
7.985453723 428
7.991224826 429
7.994562832 430
7.974647772 431
7.990135877 432
7.982921852 433
8.011170017 434
8.017042059 435
8.025313619 436
8.029831784 437
8.030827958 438
8.046488837 439
8.051604124 440
8.062438464 441
8.070136105 442
8.065888208 443
8.087208418 444
8.099037599 445
8.099111427 446
8.107460992 447
8.108799299 448
8.131221839 449
8.144473357 450
8.159152744 451
8.169344609 452
8.165723818 453
8.187700024 454
8.191980923 455
8.192126543 456
8.145844252 457
8.184291692 458
8.219025941 459
8.221621948 460
8.227898887 461
8.238289278 462
8.232138863 463
8.253558841 464
8.281473643 465
8.297315551 466
8.30639833 467
8.299454507 468
8.305580899 469
8.268521041 470
8.290554288 471
8.30605059 472
8.301585673 473
8.311944812 474
8.321270517 475
8.336560368 476
8.331136588 477
8.338012936 478
8.359863363 479
8.374459357 480
8.307234339 481
8.321094119 482
8.294235123 483
8.208056199 484
8.234061019 485
8.236504475 486
8.266370842 487
8.272094122 488
8.292969056 489
8.311923448 490
8.337224065 491
8.322444902 492
8.326727118 493
8.355461717 494
8.359257625 495
8.338382983 496
8.350151816 497
8.316653487 498
8.27921033 499
8.287370132 500
8.313508826 501
8.272367267 502
8.248189716 503
8.252262075 504
8.26866106 505
8.247001271 506
8.193088323 507
8.237835839 508
8.264418955 509
8.256612795 510
8.259904007 511
8.272823705 512
8.272416063 513
8.248507147 514
8.267754688 515
8.276409972 516
8.272044544 517
8.3086557 518
8.346015945 519
8.345296071 520
8.371928316 521
8.37711654 522
8.376863198 523
8.398582839 524
8.405374533 525
8.366501223 526
8.391115097 527
8.425247625 528
8.439903424 529
8.444885905 530
8.447449038 531
8.466761631 532
8.465767046 533
8.480542485 534
8.491420921 535
8.490825191 536
8.469146831 537
8.446035264 538
8.49798081 539
8.508194222 540
8.51320023 541
8.523670965 542
8.538498219 543
8.546774906 544
8.555392036 545
8.538842544 546
8.550218599 547
8.56002756 548
8.579426871 549
8.574337883 550
8.577481945 551
8.592979469 552
8.571146385 553
8.585839576 554
8.539730898 555
8.551135962 556
8.570467003 557
8.577757953 558
8.598107888 559
8.598811519 560
8.593276996 561
8.607045778 562
8.604228135 563
8.609774743 564
8.598935377 565
8.58270374 566
8.595894015 567
8.621967715 568
8.646949593 569
8.651005645 570
8.660681053 571
8.684947492 572
8.683260697 573
8.704516552 574
8.714741071 575

未命名.JPG

error2.JPG

回复此楼

» 猜你喜欢

新西兰Robinson研究所招收全奖PhD 已经有0人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有92人回复
笼目材料中量子自旋液体基态的证据已经有0人回复
数学教学论硕士可以读数学物理博士吗？已经有0人回复
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生已经有4人回复
澳门大学应用物理及材料工程研究院潘晖教授课题组诚招博士后已经有11人回复
求助NH4V4O10晶体的CIF文件已经有0人回复
英国全奖博士招聘-深度学习与量子物理已经有0人回复
间接带隙半导体有效质量求助已经有1人回复
投稿chemical physical letters不送审？已经有2人回复

» 本主题相关价值贴推荐，对您同样有帮助:

1stopt程序运行已经有3人回复
常微分方程参数拟合问题，matlab和1stopt均可已经有3人回复
Table Curve 2D 数据拟合求助！已经有5人回复
1stopt自变量超多时是否有简单语法？AND 自变量多数据量大的话1stopt都会算的特慢么已经有7人回复
求大神跑下1stOpt，一个sharemodel，谢谢！已经有7人回复
求大神帮忙拟合一个非线性方程，求出模型参数已经有15人回复
用1stOpt进行复数曲线拟合时，结果出错。急求教！！！已经有5人回复
1stopt拟合问题求助已经有7人回复
求助，用1stopt解方程组，两次得到的解不唯一？什么原因，大神帮忙已经有6人回复
大家在使用1stOpt进行数据拟合的时候有没有遇到过这种情况？已经有8人回复

1楼 2015-06-12 11:09:18

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

ybkooo

至尊木虫 (著名写手)

应助: 32 (小学生)
金币: 11360
散金: 178
红花: 15
帖子: 2310
在线: 543.5小时
虫号: 2710519
注册: 2013-10-09
性别: GG
专业: 控制理论与方法

Parameters a,b,c,t,m,w,p;
Variable x,y;
Function (y-a)*(t-x)^-m = b+c*cos(w*ln(t-x)+p);
Data;
1 6.9416
2 6.9382
3 7.0171
4 7.0309
5 7.0106

………………
迭代数: 37
计算用时(时:分:秒:微秒): 00:18:40:63
优化算法: 麦夸特法(Levenberg-Marquardt) + 通用全局优化法
计算结束原因: 达到收敛判断标准
均方差(RMSE): 4.65595587544745E-14
残差平方和(SSE): 1.24648069406153E-24
相关系数(R): 1
相关系数之平方(R^2): 1
决定系数(DC): 1
卡方系数(Chi-Square): 8.72434972578136E-26
F统计(F-Statistic): -4.30093764413883E17

参数最佳估算
---------- -------------
a       2.50176168046892
b       -1.22882102188015E-13
c       -1.55820635860655E-13
t       762.35234096126
m       9.63049718398193
w       1.73070665931322
p       5.69063841038279

====== 结果输出 =====

No 实测值y 计算值y
1 6.9416 6.9416000
2 6.9382 6.9382000
3 7.0171 7.0171000
4 7.0309 7.0309000
5 7.0106 7.0106000
6 7.0088 7.0088000
7 6.9971 6.9971000
8 6.9781 6.9781000
9 6.9903 6.9903000
10 6.9899 6.9899000
11 7.0172 7.0172000