[問題] import raw data然後用Gaussian找peak

[teddy0819 (在一步就是天堂地獄)]

向大家請教一些很基本的Matlab問題

目前我需要一個script可以直接輸入實驗取到的raw data file

需要用Gaussian取到這個data的Peak值，half height width

然後作圖(raw data vs Gaussian)

所以我參考了很多有的沒有

結果得到很雜亂的東西

--

function [ sigma, mu, A ] = GaussFit( filename )
%UNTITLED9 Summary of this function goes here
%   Detailed explanation goes here
data=dlmread(filename,'',2,1);
X=data(:,1);
Y=data(;,2);

xlim_(1) = min(xlim_(1),min(x));
xlim_(2) = max(xlim_(2),max(x));

ylog=log(Y);
xlog=X;
p=polyfit(xlog,ylog,2);
A2=p(1);
A1=p(2);
A0=p(3);
sigma=sqrt(-1/(2*A2));
mu=A1*sigma^2;
A=exp(A0+mu^2/(2*sigma^2));

end

---

因為是亂湊一通的 所以我也不太會自己去debug

麻煩大家有經驗的幫我看一下是出了什麼問題

(我知道是少了作圖的部份...)

--
※ 發信站: 批踢踢實業坊(ptt.cc) 
◆ From: 98.251.48.87
→ teddy0819:以解決