【调剂】北京石油化工学院2024年16个专业接受调剂
查看: 2033  |  回复: 18
【悬赏金币】回答本帖问题,作者力学的坤将赠送您 5 个金币

力学的坤

铜虫 (小有名气)

[求助] 关于matlab进行高斯拟合的一些问题,急!!!在线等已有2人参与

大家好,我用matlab中的cftool工具拟合,选用的是高斯拟合公式。
数据如下:x=[0,0.01932 ,0.03864,0.05682,0.07614,0.09545,0.11477,0.13295,0.15227,0.17159,0.19091,0.20909,0.22841,0.24773,0.26705,0.28523,0.325,0.36136,0.39318,0.41932];
                 y=[0,6.58575,27.79516,51.29901,66.48691,71.92571,68.76235,60.72499,50.39318,40.06137,30.58977,22.83914,16.52167,9.05724,3.0265,0.72391,0.7049,0.68752,0.67231,0.65982];
               cftool
图一是matlab中的高斯公式,我查资料,说高斯分布就是正态分布,请问这个说法正确吗?
图二是cftool工具拟合的结果,那个a1,b1,c1参数那后面的括号是啥意思,是上下误差吗?,那最后的a1,b1,c1取什么呢?

最关键的问题来了(困扰我半天了),那么这个拟合公式得到了,它的均值和方差是多少呢?怎么得到的呢?
我是这么想的,网上说高斯分布就是正态分布,那么我能不能和图三这样这么做,将2个式子进行对比得到均值和方差?可是得到的方差不一样,那个是对的压,快哭了,请大佬清楚拟合这方面的救救我吧,急急急!!!在线等

关于matlab进行高斯拟合的一些问题,急!!!在线等
matlab高斯公式.png


关于matlab进行高斯拟合的一些问题,急!!!在线等-1
matlabcftool拟合结果.png


关于matlab进行高斯拟合的一些问题,急!!!在线等-2
与正态公式对比.jpg
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

hzlhm

至尊木虫 (著名写手)

【答案】应助回帖

题主提出,“高斯分布就是正态分布,那么我能不能和图三这样这么做,将2个式子进行对比得到均值和方差?可是得到的方差不一样,那个是对的?”。因为你的数据基本符合高斯分布,其拟合度为0.95565,并不是1,所以,你算出的方差是会略有差别的。
QQ:2120156492
17楼2019-06-25 13:17:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hzlhm

至尊木虫 (著名写手)

【答案】应助回帖

用matlab的拟合函数nlinfit,可以得到
A =          71.2324517256635
B =         0.112643814630324
C =        0.0803833713846708
决定系数R^2=0.95565
关于matlab进行高斯拟合的一些问题,急!!!在线等-3
123.jpg

QQ:2120156492
18楼2019-06-25 13:25:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

hollya

至尊木虫 (文坛精英)

葛城?{美

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
独孤神宇: 金币+2, 鼓励交流 2019-06-21 08:29:13
1。高斯分布是正态分布。2你的想法没有错,你也看到,标准偏差并没有很大差距。如果想确定它的话,可以在cftool中自定义高斯分布拟合。
でも私たちはきっと1000回もメールをやりとりして、たぶん心は1センチくらいしか近づけませんでした。
2楼2019-06-20 21:25:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

力学的坤

铜虫 (小有名气)

引用回帖:
2楼: Originally posted by hollya at 2019-06-20 21:25:53
1。高斯分布是正态分布。2你的想法没有错,你也看到,标准偏差并没有很大差距。如果想确定它的话,可以在cftool中自定义高斯分布拟合。

嗯嗯,有问题再请教您
3楼2019-06-21 22:58:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

力学的坤

铜虫 (小有名气)

引用回帖:
2楼: Originally posted by hollya at 2019-06-20 21:25:53
1。高斯分布是正态分布。2你的想法没有错,你也看到,标准偏差并没有很大差距。如果想确定它的话,可以在cftool中自定义高斯分布拟合。

您好,我试了试您说的方法,可怎么是一条直线呀,没有拟合?
关于matlab进行高斯拟合的一些问题,急!!!在线等-4
c125c414eb4ff206fe47bed8db06f41.png

4楼2019-06-21 23:14:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hollya

至尊木虫 (文坛精英)

葛城?{美

引用回帖:
4楼: Originally posted by 力学的坤 at 2019-06-21 23:14:59
您好,我试了试您说的方法,可怎么是一条直线呀,没有拟合?

c125c414eb4ff206fe47bed8db06f41.png
...

这样吧,我时间,给你写个脚本,然后你运行一下。
でも私たちはきっと1000回もメールをやりとりして、たぶん心は1センチくらいしか近づけませんでした。
5楼2019-06-21 23:16:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

力学的坤

铜虫 (小有名气)

引用回帖:
5楼: Originally posted by hollya at 2019-06-21 23:16:09
这样吧,我时间,给你写个脚本,然后你运行一下。...

谢谢您了
6楼2019-06-21 23:19:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hollya

至尊木虫 (文坛精英)

葛城?{美

【答案】应助回帖

引用回帖:
6楼: Originally posted by 力学的坤 at 2019-06-21 23:19:36
谢谢您了...

不要谢啦。给花花就开森:)
% Gauss fitting
x=[0,0.01932,0.03864,0.05682,0.07614,0.09545,0.11477,0.13295,...
    0.15227,0.17159,0.19091,0.20909,0.22841,0.24773,0.26705,...
    0.28523,0.325,0.36136,0.39318,0.41932]
y=[0,6.58575,27.79516,51.29901,66.48691,71.92571,68.76235,...
    60.72499,50.39318,40.06137,30.58977,22.83914,16.52167,...
    9.05724,3.0265,0.72391,0.7049,0.68752,0.67231,0.65982]
Fgauss=@(c,x) c(3)/(sqrt(2*pi)*c(1))*exp(-(x-c(2)).^2/(2*c(1)^2))
% c(1): sigma, c(2): mean, c(3): coefficient
c0=[0.05,0.1,10]
c=lsqcurvefit(Fgauss,c0,x,y)
plot(x,y,'x',x,Fgauss(c,x))
legend('experimental','fitting')

» 本帖已获得的红花(最新10朵)

でも私たちはきっと1000回もメールをやりとりして、たぶん心は1センチくらいしか近づけませんでした。
7楼2019-06-22 01:53:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

力学的坤

铜虫 (小有名气)

送红花一朵
引用回帖:
7楼: Originally posted by hollya at 2019-06-22 01:53:21
不要谢啦。给花花就开森:)
% Gauss fitting
x=
y=
Fgauss=@(c,x) c(3)/(sqrt(2*pi)*c(1))*exp(-(x-c(2)).^2/(2*c(1)^2))
% c(1): sigma, c(2): mean, c(3): coefficient
c0=
c=lsqcurvefit(Fgauss,c0,x, ...

花花给您,不过还是要谢谢您了,我先仔细看看程序,有不懂的再请教你,嘻嘻
8楼2019-06-22 08:18:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hollya

至尊木虫 (文坛精英)

葛城?{美

引用回帖:
8楼: Originally posted by 力学的坤 at 2019-06-22 08:18:32
花花给您,不过还是要谢谢您了,我先仔细看看程序,有不懂的再请教你,嘻嘻...

高斯分布PDF下面面积积分必须等于1,所以要加一个C3系数。如果你的数据是归一化了的,则c3=1。

发自小木虫IOS客户端
でも私たちはきっと1000回もメールをやりとりして、たぶん心は1センチくらいしか近づけませんでした。
9楼2019-06-22 08:33:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

力学的坤

铜虫 (小有名气)

引用回帖:
9楼: Originally posted by hollya at 2019-06-22 08:33:50
高斯分布PDF下面面积积分必须等于1,所以要加一个C3系数。如果你的数据是归一化了的,则c3=1。
...

嗯嗯,所以意思是c(3)/(sqrt(2*pi)*c(1))这个就是cftool中那个a1是吗?
10楼2019-06-22 08:41:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 力学的坤 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 普通院校药学硕士,做合成的,感觉找不到工作 +8 pom戴墨镜 2024-04-24 15/750 2024-04-25 20:50 by pom戴墨镜
[考博] 25年博士申请 +5 Changzixuan 2024-04-25 5/250 2024-04-25 17:37 by 那年的旧城
[基金申请] 基金开始函评了吗? +11 wych1103 2024-04-25 11/550 2024-04-25 16:54 by 未岸亡灵
[教师之家] 期末给学生划重点都是什么话术啊 +16 luokereng 2024-04-20 18/900 2024-04-25 15:46 by BusyGer
[硕博家园] 聊天 +11 暮色恋伊人 2024-04-22 12/600 2024-04-25 13:53 by UCTS
[博后之家] 南京大学-广州大学联合招聘博士后 欢迎广大优秀人才!!! +4 黑魔变身啾 2024-04-20 12/600 2024-04-25 11:18 by dodonaomi
[基金申请] 前同事拿我以前未中的国自然标书申请今年的国自然,现在基金委发来调查函 +33 bjdxyxy 2024-04-19 53/2650 2024-04-25 08:05 by 大雁高飞
[考研] 0854-0855调剂 +7 shangannum1 2024-04-21 10/500 2024-04-25 01:31 by 啊廖sh
[考研] 博士生招生!! +5 2458827884 2024-04-19 5/250 2024-04-24 22:28 by 不懂的风清
[基金申请] 企业博后是否能申请CSC博士后项目?谢谢 66+3 龙包子 2024-04-19 6/300 2024-04-24 21:59 by 3115321
[考博] 24年 申博 化学/材料 一作6篇sci +9 wangyp123 2024-04-23 11/550 2024-04-24 19:01 by bangbangbiu
[硕博家园] 考研,求职还是考编? +13 xizj 2024-04-21 21/1050 2024-04-24 17:23 by jxdic
[论文投稿] 期刊推荐 20+4 木颜尘ip 2024-04-22 7/350 2024-04-24 10:06 by bobvan
[考博] 申博成果界定是根据Jcr分区还是中科院分区 +4 我属驴核动力驴 2024-04-22 5/250 2024-04-24 08:47 by 晓目崇
[教师之家] 找工作过程中碰到好多人问南京某新换校长211的政策,规定如下 +15 ZHONGWU_U 2024-04-21 17/850 2024-04-24 05:12 by juta256
[公派出国] 東京大学 应用化学专业柳田研究室 招收修士/博士生 +5 蝶梦123 2024-04-19 5/250 2024-04-23 22:40 by 90472jvvh
[教师之家] 发现督导在听课,需要跟TA招呼示意吗?有同学看手机、课件有疏漏算大问题吗? +11 河西夜郎 2024-04-19 11/550 2024-04-23 15:01 by mylifecg
[高分子] 请问UV灯是365nm的,那么选光引发剂的波长选多少的?要完全一致吗? +4 engledd2004 2024-04-21 4/200 2024-04-22 16:08 by wangcz23
[论文投稿] 研究光催化的,好中的三四区 20+3 sl.0117 2024-04-20 3/150 2024-04-22 09:53 by bobvan
[硕博家园] 材料工程,二本硕士,值得读吗 +3 简风23 2024-04-20 6/300 2024-04-21 12:19 by 简风23
信息提示
请填处理意见