24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1406  |  回复: 5
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

saji

铁虫 (初入文坛)

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

各位前辈:
小弟使用1stopt拟合一等式,原始形式如下:

得到的拟合结果正常,图形如下:


但是,使用同样的样本,进行变式的拟合时,却得到怪异的图形。。。不知问题在哪里?
变式如下:

得到的函数图形如下:


下面是我使用变式的拟合代码:
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

1stOpt拟合代码的问题求助
未命名.JPG


1stOpt拟合代码的问题求助-1
error2.JPG
回复此楼

» 猜你喜欢

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

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

ybkooo

至尊木虫 (著名写手)

^^

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
abcd
5楼2015-06-12 11:54:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
公式能写出显示就用显示形式,隐式方程拟合很麻烦的,写成:y=f(x)这种形式
3楼2015-06-12 11:35:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ybkooo

至尊木虫 (著名写手)

^^

【答案】应助回帖

感谢参与,应助指数 +1
最好写成显式方程
abcd
4楼2015-06-12 11:35:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

saji

铁虫 (初入文坛)

引用回帖:
3楼: Originally posted by dingd at 2015-06-12 11:35:06
公式能写出显示就用显示形式,隐式方程拟合很麻烦的,写成:y=f(x)这种形式

显式的拟合不稳定,在数学版发帖后,版主建议用隐式拟合,具体原因不明。
6楼2015-06-12 13:26:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见