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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » 如何用matlab实现分峰拟合？
191/2返回列表12下一页
查看: 6584  |  回复: 18
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

boluomi

银虫 (小有名气)

[求助] 如何用matlab实现分峰拟合？

大家好。遇到一个问题。
光谱中得到两个峰有重叠，我想用数学方法，把它分开，想用两个高斯曲线对两个峰进行拟合。
查了一下matlab里有cftool工具箱，用于曲线拟合，但是似乎只能拟合一条曲线。
那我应该用什么工具或者函数去将它们分开呢？？？

期待高手回复啊~
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-10-10 16:28:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

dingd

铁杆木虫 (职业作家)
★
xzhdty: 金币+1, 感谢参与 2013-10-12 06:59:57
两个高斯方程：
y=y0+(A/(w*sqrt(PI/2)))*exp(-2*((x-xc)/w)^2)+(A1/(w1*sqrt(PI/2)))*exp(-2*((x-xc1)/w1)^2);

1stOpt求解拟合：

均方差(RMSE):30.1128758590207
残差平方和(SSE):663766.83411058
相关系数(R): 0.987353824563334
相关系数之平方(R^2): 0.974867574879842
确定系数(DC): 0.974867574879842
卡方系数(Chi-Square): 473.331101403844
F统计(F-Statistic): 4807.8660437913

参数        最佳估算
----------        -------------
y0        407.883083321768
a        211736.791436449
w        326.456566603881
xc        3373.56868156967
a1        45282.3684328329
w1        161.790412752586
xc1        3181.03375043014
如何用matlab实现分峰拟合？
c1.jpg
赞一下(7人) 回复此楼
11楼2013-10-11 11:19:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与，应助指数 +1
可以自定义两个高斯函数拟合就行了。有数据的话放上来看看。
赞一下(2人) 回复此楼
2楼2013-10-10 17:19:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

virtualzx

木虫 (著名写手)

【答案】应助回帖

感谢参与，应助指数 +1
两个峰的话拟合f=c1*f1+c2*f2就好了吧
赞一下(2人) 回复此楼
6楼2013-10-11 08:13:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

change0618

铁杆木虫 (著名写手)

方丈大师
用Origin软件啊
回复此楼
3楼2013-10-10 17:25:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

boluomi

银虫 (小有名气)
引用回帖:
3楼: Originally posted by change0618 at 2013-10-10 17:25:43
用Origin软件啊

Origin 是可以
文章一直用matlab画图，写好程序以后调用比较方便。
赞一下 回复此楼
4楼2013-10-10 23:04:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

boluomi

银虫 (小有名气)
引用回帖:
2楼: Originally posted by dingd at 2013-10-10 17:19:44
可以自定义两个高斯函数拟合就行了。有数据的话放上来看看。

数据有的，我明天放上来，一列x 一列y
拟合的话 我该用哪个命令？
赞一下 回复此楼
5楼2013-10-10 23:06:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

boluomi

银虫 (小有名气)
引用回帖:
2楼: Originally posted by dingd at 2013-10-10 17:19:44
可以自定义两个高斯函数拟合就行了。有数据的话放上来看看。

你好，是这么两列数据。
3100          708.789
3100.96        694.71
3101.92        668.305
3102.87        705.526
3103.83        736.69
3104.79        709.833
3105.75        677.329
3106.7        701.353
3107.66        745.102
3108.62        716.283
3109.58        736.531
3110.53        750.786
3111.49        659.86
3112.45        717.862
3113.41        681.788
3114.36        718.037
3115.32        684.644
3116.28        710.272
3117.24        672.669
3118.19        723.303
3119.15        752.189
3120.11        758.817
3121.07        733.668
3122.02        731.709
3122.98        740.14
3123.94        907.128
3124.9        689.089
3125.85        753.108
3126.81        721.466
3127.77        722.941
3128.73        747.104
3129.69        769.367
3130.64        765.5
3131.6        784.059
3132.56        768.876
3133.52        762.824
3134.47        744.88
3135.43        724.803
3136.39        783.683
3137.35        781.165
3138.3        806.298
3139.26        894.995
3140.22        741.651
3141.18        774.93
3142.13        808.402
3143.09        820.667
3144.05        817.009
3145.01        746.474
3145.96        773.844
3146.92        762.054
3147.88        726.669
3148.84        768.527
3149.79        856.055
3150.75        812.867
3151.71        807.237
3152.67        820.094
3153.63        837.945
3154.58        837.754
3155.54        834.965
3156.5        831.655
3157.46        791.682
3158.41        855.771
3159.37        789.329
3160.33        804.097
3161.29        829.472
3162.24        856.303
3163.2        851.674
3164.16        829.673
3165.12        876.698
3166.07        862.678
3167.03        839.15
3167.99        800.249
3168.95        874.286
3169.9        835
3170.86        860.34
3171.82        898.497
3172.78        903
3173.73        850.578
3174.69        912.219
3175.65        857.739
3176.61        842.704
3177.56        911.516
3178.52        926.254
3179.48        887.496
3180.44        916.296
3181.4        875.405
3182.35        843.221
3183.31        846.213
3184.27        905.443
3185.23        898.627
3186.18        876.333
3187.14        902.046
3188.1        927.713
3189.06        872.089
3190.01        842.035
3190.97        865.138
3191.93        852.283
3192.89        882.401
3193.84        930.399
3194.8        911.123
3195.76        921.416
3196.72        907.867
3197.67        907.867
3198.63        925.442
3199.59        921.069
3200.55        919.839
3201.5        967.699
3202.46        936.218
3203.42        926.273
3204.38        938.049
3205.34        971.491
3206.29        907.132
3207.25        957.596
3208.21        939.047
3209.17        961.981
3210.12        1009.05
3211.08        949.055
3212.04        936.497
3213        960.569
3213.95        973.64
3214.91        986.265
3215.87        956.627
3216.83        1003.65
3217.78        916.009
3218.74        976.638
3219.7        946.083
3220.66        964.206
3221.61        958.146
3222.57        975.891
3223.53        896.753
3224.49        956.705
3225.44        946.263
3226.4        937.721
3227.36        916.9
3228.32        977.805
3229.27        958.338
3230.23        967.141
3231.19        976.897
3232.15        946.288
3233.11        979.427
3234.06        915.644
3235.02        898.451
3235.98        965.439
3236.94        976.109
3237.89        993.969
3238.85        952.39
3239.81        942.018
3240.77        943.222
3241.72        929.966
3242.68        969.186
3243.64        899.074
3244.6        917.11
3245.55        936.979
3246.51        961.768
3247.47        963.255
3248.43        963.853
3249.38        956.267
3250.34        962.121
3251.3        916.094
3252.26        922.077
3253.21        940.478
3254.17        952.606
3255.13        942.247
3256.09        941.809
3257.05        937.923
3258        938.588
3258.96        968.282
3259.92        975.42
3260.88        953.238
3261.83        946.635
3262.79        969.845
3263.75        947.682
3264.71        924.855
3265.66        947.845
3266.62        944.69
3267.58        956.398
3268.54        958.313
3269.49        990.118
3270.45        981.504
3271.41        990.651
3272.37        968.205
3273.32        967.263
3274.28        963.566
3275.24        936.882
3276.2        904.822
3277.15        914.939
3278.11        971.198
3279.07        966.908
3280.03        970.592
3280.98        930.038
3281.94        935.353
3282.9        908.529
3283.86        909.211
3284.82        924.378
3285.77        962.741
3286.73        966.747
3287.69        996.118
3288.65        963.418
3289.6        940.328
3290.56        959.526
3291.52        1021.25
3292.48        930.478
3293.43        945.468
3294.39        963.49
3295.35        978.06
3296.31        931.783
3297.26        919.043
3298.22        939.896
3299.18        934.485
3300.14        928.746
3301.09        902.453
3302.05        893.227
3303.01        945.764
3303.97        919.9
3304.92        884.78
3305.88        924.451
3306.84        937.369
3307.8        967.786
3308.76        985.965
3309.71        909.839
3310.67        915.817
3311.63        933.451
3312.59        981.602
3313.54        938.72
3314.5        954.212
3315.46        950.619
3316.42        973.048
3317.37        956.87
3318.33        939.521
3319.29        961.001
3320.25        937.975
3321.2        921.51
3322.16        928.797
3323.12        958.44
3324.08        980.693
3325.03        944.999
3325.99        962.885
3326.95        918.368
3327.91        953.907
3328.86        952.715
3329.82        926.838
3330.78        955.141
3331.74        954.141
3332.69        922.3
3333.65        936.756
3334.61        950.497
3335.57        950.934
3336.53        935.433
3337.48        944.388
3338.44        1033.43
3339.4        1062.86
3340.36        942.571
3341.31        915.928
3342.27        900.672
3343.23        934.25
3344.19        1034.13
3345.14        930.494
3346.1        946.672
3347.06        944.012
3348.02        937.845
3348.97        934.173
3349.93        944.629
3350.89        941.77
3351.85        959.849
3352.8        969.823
3353.76        931.532
3354.72        908.798
3355.68        917.32
3356.63        965.02
3357.59        951.885
3358.55        929.668
3359.51        951.239
3360.47        943.042
3361.42        947.303
3362.38        954.67
3363.34        930.663
3364.3        928.161
3365.25        915.537
3366.21        948.427
3367.17        913.928
3368.13        1010.82
3369.08        998.558
3370.04        953.647
3371        930.22
3371.96        922.995
3372.91        910.812
3373.87        923.961
3374.83        944.086
3375.79        963.028
3376.74        950.392
3377.7        927.989
3378.66        922.123
3379.62        942.391
3380.57        905.956
3381.53        934.104
3382.49        940.859
3383.45        944.448
3384.41        904.425
3385.36        932.208
3386.32        895.656
3387.28        906.308
3388.24        903.211
3389.19        958.04
3390.15        956.944
3391.11        901.482
3392.07        958.607
3393.02        940.607
3393.98        964.19
3394.94        922.221
3395.9        907.798
3396.85        928.474
3397.81        930.36
3398.77        933.704
3399.73        913.502
3400.68        912.175
3401.64        943.242
3402.6        961.36
3403.56        955.079
3404.51        917.757
3405.47        944.599
3406.43        909.184
3407.39        850.412
3408.34        938.846
3409.3        923.571
3410.26        951.986
3411.22        918.966
3412.18        974.323
3413.13        917.817
3414.09        896.613
3415.05        902.011
3416.01        955.746
3416.96        910.266
3417.92        910.266
3418.88        910.266
3419.84        870.697
3420.79        883.744
3421.75        916.898
3422.71        893.442
3423.67        881.977
3424.62        837.092
3425.58        872.258
3426.54        903.392
3427.5        898.931
3428.45        886.629
3429.41        859.594
3430.37        916.215
3431.33        921.24
3432.28        869.179
3433.24        911.207
3434.2        868.206
3435.16        905.852
3436.11        900.902
3437.07        891.899
3438.03        904.433
3438.99        867.895
3439.95        862.479
3440.9        857.928
3441.86        915.995
3442.82        874.405
3443.78        921.912
3444.73        857.289
3445.69        870.137
3446.65        959.321
3447.61        880.495
3448.56        852.276
3449.52        849.708
3450.48        884.491
3451.44        864.618
3452.39        874.572
3453.35        841.954
3454.31        889.267
3455.27        929.924
3456.22        925.688
3457.18        866.45
3458.14        804.311
3459.1        910.597
3460.05        898.035
3461.01        891.251
3461.97        858.632
3462.93        839.813
3463.89        862.574
3464.84        808.632
3465.8        839.343
3466.76        835.989
3467.72        807.462
3468.67        825.235
3469.63        859.435
3470.59        861.934
3471.55        839.059
3472.5        855.412
3473.46        825.265
3474.42        842.192
3475.38        868.524
3476.33        790.168
3477.29        811.176
3478.25        799.466
3479.21        822.863
3480.16        819.689
3481.12        813.93
3482.08        855.094
3483.04        858.002
3483.99        884.332
3484.95        837.72
3485.91        855.001
3486.87        843.814
3487.82        829.135
3488.78        738.471
3489.74        816.457
3490.7        821.16
3491.66        780.007
3492.61        841.055
3493.57        818.032
3494.53        814.248
3495.49        799.047
3496.44        811.591
3497.4        782.574
3498.36        826.376
3499.32        796.016
3500.27        775.042
3501.23        781.974
3502.19        792.695
3503.15        810.342
3504.1        682.091
3505.06        796.478
3506.02        795.042
3506.98        771.991
3507.93        717.439
3508.89        797.969
3509.85        799.71
3510.81        779.55
3511.76        738.401
3512.72        761.285
3513.68        743.421
3514.64        822.276
3515.59        798.917
3516.55        789.649
3517.51        732.145
3518.47        791.709
3519.43        757.366
3520.38        782.661
3521.34        742.351
3522.3        764.946
3523.26        752.987
3524.21        744.046
3525.17        761.703
3526.13        742.273
3527.09        738.229
3528.04        783.005
3529        751.86
3529.96        727.953
3530.92        745.169
3531.87        706.289
3532.83        696
3533.79        700.243
3534.75        727.388
3535.7        680.45
3536.66        691.197
3537.62        708.119
3538.58        687.465
3539.53        690.36
3540.49        651.937
3541.45        698.326
3542.41        740.562
3543.37        725.907
3544.32        687.902
3545.28        724.447
3546.24        706.984
3547.2        711.5
3548.15        681.23
3549.11        754.745
3550.07        717.572
3551.03        720.078
3551.98        702.874
3552.94        678.616
3553.9        682.581
3554.86        648.752
3555.81        670.038
3556.77        649.259
3557.73        666.898
3558.69        669.835
3559.64        668.721
3560.6        630.763
3561.56        605.627
3562.52        632.388
3563.47        640.26
3564.43        737.514
3565.39        645.838
3566.35        642.604
3567.31        673.982
3568.26        677.727
3569.22        701
3570.18        622.11
3571.14        657.936
3572.09        612.68
3573.05        595.306
3574.01        642.264
3574.97        646.668
3575.92        646.743
3576.88        637.214
3577.84        701.426
3578.8        616.691
3579.75        678.649
3580.71        624.181
3581.67        607.299
3582.63        638.38
3583.58        667.831
3584.54        654.38
3585.5        655.326
3586.46        562.313
3587.41        606.946
3588.37        609.553
3589.33        610.347
3590.29        618.727
3591.24        627.736
3592.2        623.808
3593.16        592.115
3594.12        629.498
3595.08        616.81
3596.03        679.827
3596.99        517.082
3597.95        543.604
3598.91        574.955
3599.86        618.561
3600.82        641.654
3601.78        577.294
3602.74        559.799
3603.69        600.499
3604.65        608.409
3605.61        568.967
3606.57        578.974
3607.52        600.092
3608.48        600.121
3609.44        596.415
3610.4        601.219
3611.35        593.517
3612.31        537.996
3613.27        555.28
3614.23        552.577
3615.18        617.527
3616.14        568.494
3617.1        627.39
3618.06        618.741
3619.02        588.102
3619.97        620.845
3620.93        585.48
3621.89        539.361
3622.85        592.642
3623.8        594.773
3624.76        576.901
3625.72        561.531
3626.68        525.253
3627.63        562.053
3628.59        575.173
3629.55        580.032
3630.51        582.086
3631.46        528.101
3632.42        572.982
3633.38        585.76
3634.34        554.639
3635.29        587.364
3636.25        559.137
3637.21        584.739
3638.17        569.378
3639.12        560.558
3640.08        583.066
3641.04        490
3642        490.387
3642.95        528.376
3643.91        533.206
3644.87        555.325
3645.83        577.968
3646.79        528.855
3647.74        613.819
3648.7        529.24
3649.66        556.709
3650.62        504.597
3651.57        523.8
3652.53        534.456
3653.49        497.625
3654.45        504.273
3655.4        535.98
3656.36        498.532
3657.32        489.712
3658.28        536.229
3659.23        521.667
3660.19        551.4
3661.15        537.162
3662.11        470.741
3663.06        478.611
3664.02        515.125
3664.98        502.328
3665.94        576.625
3666.89        515.164
3667.85        532.282
3668.81        468.522
3669.77        515.455
3670.73        509.786
3671.68        509.048
3672.64        509.508
3673.6        496.407
3674.56        514.068
3675.51        542.452
3676.47        528.127
3677.43        522.245
3678.39        549.25
3679.34        550.664
3680.3        476.66
3681.26        498.09
3682.22        474.516
3683.17        511.26
3684.13        461.523
3685.09        471.802
3686.05        504.334
3687        514.94
3687.96        478.289
3688.92        441.027
3689.88        503.686
3690.83        634.893
3691.79        578.515
3692.75        450.361
3693.71        451.668
3694.66        486.616
3695.62        435.422
3696.58        474.252
3697.54        519.364
3698.5        454.2
3699.45        464.344
3700.41        509.323
3701.37        483
3702.33        471.863
3703.28        470.959
3704.24        475.101
3705.2        514.618
3706.16        482.239
3707.11        419.304
3708.07        498.671
3709.03        501.512
3709.99        511.495
3710.94        456.052
3711.9        432.678
3712.86        453.414
3713.82        455.601
3714.77        421.664
3715.73        422.22
3716.69        454.883
3717.65        441.939
3718.6        474.271
3719.56        466.96
3720.52        460.774
3721.48        441.027
3722.44        440.266
3723.39        410.192
3724.35        426.313
3725.31        429.377
3726.27        422.021
3727.22        480.295
3728.18        487
3729.14        441.895
3730.1        447.645
3731.05        478.235
3732.01        468.704
3732.97        428.597
3733.93        493.424
3734.88        490.757
3735.84        473.123
3736.8        485.484
3737.76        458.253
3738.71        488.276
3739.67        489
3740.63        464.276
3741.59        461.363
3742.54        427.571
3743.5        410.286
3744.46        412.224
3745.42        445.682
3746.37        411.404
3747.33        449.029
3748.29        462.779
3749.25        454.332
3750.21        445.557
3751.16        421.803
3752.12        424.244
3753.08        405.196
3754.04        411.032
3754.99        446.668
3755.95        463.616
3756.91        422.208
3757.87        409.932
3758.82        420.17
3759.78        476.829
3760.74        443.784
3761.7        441.856
3762.65        426.05
3763.61        413.501
3764.57        450.228
3765.53        389.194
3766.48        405.804
3767.44        433.103
3768.4        458.688
3769.36        441.85
3770.31        452.908
3771.27        432.46
3772.23        458.556
3773.19        452.971
3774.15        436.36
3775.1        455.314
3776.06        467.744
3777.02        450.397
3777.98        425.073
3778.93        386.28
3779.89        458.179
3780.85        423.898
3781.81        406.266
3782.76        405.143
3783.72        443.129
3784.68        402.42
3785.64        455.596
3786.59        429.65
3787.55        456.874
3788.51        444.833
3789.47        401.194
3790.42        415.727
3791.38        430.789
3792.34        405.753
3793.3        403.431
3794.25        434.64
3795.21        402.425
3796.17        396.774
3797.13        395.829
3798.08        459.885
3799.04        422.098
3800        423.345
赞一下 回复此楼
7楼2013-10-11 08:36:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

boluomi

银虫 (小有名气)
引用回帖:
6楼: Originally posted by virtualzx at 2013-10-11 08:13:39
两个峰的话拟合f=c1*f1+c2*f2就好了吧

不知道怎么写啊，抓狂，还请高手指点指点！matlab不精通。

[filename,filepath,filterindex]=uigetfile('*.txt','select a txt file');
datafile=strcat(filepath,filename);
data=textread(datafile);
x=data(:,1);
y=data(:,2);
f1=a1*exp(-((x-b1)/c1)^2);
f2=a2*exp(-((x-b2)/c2)^2);
y=d1*f1+d2*f2;
plot(x,y);
hold on
plot(x,f1);
hold on
plot(x,f2);
赞一下 回复此楼
8楼2013-10-11 09:47:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

virtualzx

木虫 (著名写手)

【答案】应助回帖

引用回帖:
8楼: Originally posted by boluomi at 2013-10-10 20:47:28
不知道怎么写啊，抓狂，还请高手指点指点！matlab不精通。

=uigetfile('*.txt','select a txt file');
datafile=strcat(filepath,filename);
data=textread(datafile);
x=data(:,1);
y=data(:,2);
f1=a1 ...

你的峰两边都没扫完，很难确定基线啊
赞一下 回复此楼
9楼2013-10-11 10:06:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jeremyfrc

木虫 (正式写手)

【答案】应助回帖

感谢参与，应助指数 +1
如果是XPS的峰，我可以帮你拟合一下，但是如果不是，就很难弄了。
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

boluomi
哈哈哈哈，来吧
10楼2013-10-11 10:30:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 boluomi 的主题更新
191/2返回列表12下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭