24小时热门版块排行榜

返回列表

当前主题已经存档。

cxd15505

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 569.5
散金: 4
红花: 3
帖子: 120
在线: 55.4小时
虫号: 593982
注册: 2008-09-06

[交流] 【求助】用matlab 模拟高斯拟合已有2人参与

先在想把光谱数据，用matlab用高斯拟合重新模拟求出最大峰值，但是看了help还是不知道如何使用。求教。
举个简单的例子x=[1,2,3,4,5,6,7]
y=[1,3,5,10,6,4,2]
请问用matlab用高斯拟合如何实现，请写明编写代码，谢谢

回复此楼

» 猜你喜欢

论文投稿，期刊推荐已经有6人回复
寻求一种能扛住强氧化性腐蚀性的容器密封件已经有3人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有4人回复
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗已经有6人回复
孩子确诊有中度注意力缺陷已经有14人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有5人回复
2025冷门绝学什么时候出结果已经有3人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复

1楼2010-04-19 09:42:23

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

anyuezhiji

银虫 (正式写手)

星空行者

仿真EPI: 5
应助: 0 (幼儿园)
贵宾: 0.208
金币: 467.3
红花: 3
帖子: 416
在线: 89.3小时
虫号: 943880
注册: 2010-01-17
性别: GG
专业: 凝聚态物性 II ：电子结构

★ ★
adu886886(金币+2):谢谢交流 2010-04-20 08:19
cxd15505(金币+5): 2010-04-20 08:54

很简单的

引用回帖:

function gussfit
%正态分布拟合
%%PSL@CSU
%%QQ:547423688
%%Email:anyuezhiji@qq.com
%%Edit @ 2010.4.19
x=[1,2,3,4,5,6,7];
y=[1,3,5,10,6,4,2];
%这样
% R=[];
% for i=1:length(x)
%    R=[R repmat(x(i),1,y(i))];
% end
% [muhat,sigmahat,muci,sigmaci] = normfit(R)
%或
alpha=0.05;
censoring=[];
[muhat,sigmahat,muci,sigmaci] = normfit(x,alpha,censoring,y);
disp(['在',num2str(100*(1-alpha)),'%的置信度下：']);
disp(['峰值mu为  : ',num2str(muhat,'%10.5f')]);
disp(['置信区间为:  [',num2str(muci','%10.5f'),']']);
disp(['sigma为  : ',num2str(sigmahat,'%10.5f')]);
disp(['置信区间为:  [',num2str(sigmaci','%10.5f'),']']);
[mu1,sigma1]= normpdffit(x,y/sum(y),muhat,sigmahat);
disp(['继续调用自定义函数拟合：']);
disp(['峰值mu为  : ',num2str(mu1,'%10.5f')]);
disp(['sigma为  : ',num2str(sigma1,'%10.5f')]);
figure;
%绘制直方图
bar(x,y/sum(y));
hold on;
%绘制对应点
plot(x,y/sum(y),'ro');
%定义绘图区间
t=min(x)

max(x)-min(x))/50:max(x);
%拟合分布曲线数据
nyf=normpdf(t,muhat,sigmahat);
plot(t,nyf,'r-');
mynyf=normpdf(t,mu1,sigma1);
plot(t,mynyf,'k-.','LineWidth',2);
legend('hist','hist value','normfit','mynormfit');
function [muhat,sigmahat]= normpdffit(x,y,mu0,sigma0)
%正态分布pdf参数估计
p=[mu0 sigma0];
opt=optimset('TolX',0.001,'Display','off');
[x,fval]=fminsearch(@mynormfit,p,opt,x,y);
muhat=x(1);
sigmahat=x(2);
function Variance= mynormfit(p,x,y)
%计算对于任何一组参数p(1)，p(2)，给出当前数据下的方差来。
Variance=var(y-normpdf(x,p(1),p(2)));

引用回帖:

在95%的置信度下：
峰值mu为  : 4.19355
置信区间为:  [3.65430 4.73280]
sigma为  : 1.47013
置信区间为:  [1.17480 1.96509]
继续调用自定义函数拟合：
峰值mu为  : 4.17697
sigma为  : 1.45444