function?idemo1
%?Self-contained?demonstration?function?for?comparing?the?iPeak?and
%?Peakfit?functions?applied?to?a?test?signal?consisting?of?four?peaks.
%??T.?C.?O‘Haver?July?2012
increment=1;
x=[1:increment:500];
%?For?each?simulated?peak?specify?the?amplitude?position?and?width
%?and?the?overall?random?noise?level.
amp=[100?200?300?400];??%?Amplitudes?of?the?peaks??（Change?if?desired）
pos=[100?200?300?400];???%?Positions?of?the?peaks?（Change?if?desired）
wid=[50?50?50?50];???%?Widths?of?the?peaks?（Change?if?desired）
Noise=2;?%?Amount?of?random?noise?added?to?the?signal.?（Change?if?desired）?

%?A?=?matrix?containing?one?of?the?unit-amplidude?peak?in?each?of?its?rows
A?=?zeros（length（pos）length（x））;
ActualPeaks=[0?0?0?0?0];
p=1;
for?k=1:length（pos）
??????%?Create?a?series?of?peaks?of?different?x-positions
??????A（k:）=exp（-（（x-pos（k））./（0.6005615.*wid（k）））.^2）;?%?Gaussian?peaks
??????%?alternatively?A（k:）=ones（size（x））./（1+（（x-pos（k））./（0.5.*wid（k）））.^2）;??%?Lorentzian?peaks
??????%?Assembles?actual?parameters?into?ActualPeaks?matrix:?each?row?=?1
??????%?peak;?columns?are?Peak?#?Position?Height?Width?Area
??????ActualPeaks（p:）?=?[p?pos（k）?amp（k）?wid（k）?1.0646.*amp（k）.*wid（k）];?
??????p=p+1;
end
z=amp*A;??%?Multiplies?each?row?by?the?corresponding?amplitude?and?adds?them?up
y=z+Noise.*randn（size（z））;??%?Optionally?adds?random?noise
%?y=y+gaussian（x04000）;?%?Optionally?adds?a?broad?background?signal
demodata=[x‘?y‘];?%?Assembles?x?and?y?vectors?into?data?matrix
?
%?Compare?MeasuredPeaks?to?ActualPeaks
disp（‘-----------------------------------------------------------------‘）
disp（‘Demonstration?script?for?comparing?the?iPeak?and?Peakfit?functions.‘）
disp（‘Detection?and?measurement?of?four?peaks?with?heights?100?200?‘）
disp（‘300?and?400.?Data?are?contained?in?the?2?x?n?matrix?demodata‘）
disp（‘?‘）
disp（[‘Random?noise?in?signal?=?‘?num2str（Noise）]）
disp（‘?‘）
disp（‘?????????Peak?#????Position??????Height???????Width????????Area‘）
ActualPeaks
disp（‘?‘）
disp（‘Using?iPeak:‘）
PeakD=6;
disp（[‘>>?MeasuredPeaks=ipeak（demodata‘?num2str（PeakD）?‘）‘]）
disp（‘?‘）
disp（‘iPeak?finds?all?the?peaks?in?a?signal?at?once?returning?a‘）
disp（‘table?of?estimated?peak?positions?heights?widths?and?areas.‘）
disp（‘?‘）
%?Now?call?iPeak?with?specified?value?of?the?peak?detection?factor?（PeakD）.
%?The?numerical?argument?（PeakD）?is?an?estimate?of?the?ratio?of?the?typical
%?peak?width?to?the?length?of?the?entire?data?record?（PeakD）.?Small?values
%?detect?fewer?peaks;?larger?values?detect?more?peaks.
tic;
MeasuredPeaks=ipeak（demodataPeakD）;
iPeakTime=toc;
NumPeaks=max（MeasuredPeaks（:1））;
disp（[‘Number?of?peaks?detected?=?‘?num2str（NumPeaks）]）?
disp（‘?????????Peak?#????Position??????Height???????Width????????Area‘）
MeasuredPeaks
if?size（ActualPeaks）==size（MeasuredPeaks）
????AverageiPeakErrors=mean（abs（100.*（ActualPeaks-MeasuredPeaks）./ActualPeaks））
end
disp（[‘Elapsed?time?is?‘?num2str（iPeakTime）?‘

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件??????15310??2013-07-21?16:40??findpeaks\idemos\hs_err_pid820.log

?????文件??????89415??2012-07-31?07:07??findpeaks\idemos\idemo1.m

?????文件??????91682??2012-07-30?06:25??findpeaks\idemos\idemo2.m

?????文件?????107015??2013-07-21?16:18??findpeaks\idemos\ipf.m

?????文件??????37350??2013-06-10?10:34??findpeaks\idemos.zip

?????文件????1397560??2013-06-10?10:12??findpeaks\PeakFind2.wmv

????..A.SH.??????5632??2013-06-25?22:17??findpeaks\Thumbs.db

?????目錄??????????0??2013-07-24?01:24??findpeaks\idemos

?????目錄??????????0??2013-06-25?22:17??findpeaks

-----------?---------??----------?-----??----

??????????????1743964????????????????????9