[問題] import raw data然後用Gaussian找peak
向大家請教一些很基本的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
→
01/24 02:44, , 1F
01/24 02:44, 1F
MATLAB 近期熱門文章
PTT數位生活區 即時熱門文章
