作者teddy0819 (在一步就是天堂地狱)
看板MATLAB
标题[问题] import raw data然後用Gaussian找peak
时间Thu Jan 23 12:26:12 2014
向大家请教一些很基本的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
1F：→ teddy0819:以解决 01/24 02:44