function?get_peak（）;
directory=‘E:\MATLAB701\work\yl\data\‘;
dirname=strcat（directory‘*.txt‘）;%.txt‘）;
namelist=dir（dirname）;
a=length（namelist）;
answer=zeros（8a-1）;%記錄結果
%answer（:1）=[‘名稱‘‘峰高‘‘半高寬‘‘基線‘‘峰位‘‘半高位1‘‘半高位2‘‘半高位縱坐標‘];
for?k=1:a;
????if?strcmp（namelist（k）.name‘base.txt‘）;
????????filename=strcat（directorynamelist（k）.name）;
????????fid=fopen（filename‘rt‘）;?
????????filecontent=textscan（fid‘%n%n‘）;
????????m=length（filecontent{1}）;%行數
????????data_base=zeros（m2）;
????????for?n=1:2;
????????????data_base（:n）=filecontent{n};
????????end;
????????fclose（fid）;
????end;
end;
kkk=0;
for?k=1:a;
????if?strcmp（namelist（k）.name‘base.txt‘）;
????else
????????kkk=kkk+1;
????????filename=strcat（directorynamelist（k）.name）;
????????fid=fopen（filename‘rt‘）;
????????filecontent=textscan（fid‘%n%n‘）;
????????m=length（filecontent{1}）;%行數
????????data=zeros（m2）;
????????for?n=1:2;
????????????data（:n）=filecontent{n};
????????end;
????????data（:2）=data（:2）./data_base（:2）;
????????data（:2）=smooth（data（:2）10）;%平滑曲線！！！！！！！！！！！！
????????data_write=data‘;
????????directory_temp=‘E:\MATLAB701\work\yl\output\‘;
????????outputname=strcat（directory_tempnamelist（k）.name）;
????????fid2=fopen（outputname‘wt‘）;
????????fprintf（fid2‘%6f\t%6f\n‘data_write）;
????????fclose（fid2）;
????????x2=data（11）:0.2:data（m1）;
????????y2=interp1（data（:1）data（:2）x2‘spline‘）;
????????%baseline=（data（m-42）+data（m-32）+data（m-22）+data（m-12）+data（m2））/5;%基線
????????baseline=（data（12）+data（22）+data（32）+data（m-22）+data（m-12）+data（m2））/6;%基線
????????%baseline=0;
????????max_y2=y2（1）;
????????max_x2=1;
????????for?ii=2:length（x2）;
????????????if?y2（ii）>max_y2;
????????????????max_y2=y2（ii）;
????????????????max_x2=ii;
????????????end;
????????end;
????????height_peak=max_y2-baseline;%峰高
????????peakline=baseline+0.9*height_peak;
????????num_temp=0;
????????peak_pos=0;%峰位的序號
????????for?kk=1:length（x2）;
????????????if?y2（kk）>peakline;
????????????????num_temp=num_temp+1;
????????????????peak_pos=peak_pos+kk;
????????????end;
????????end;
????????peak_pos=fix（peak_pos/num_temp）;
????????height_peak2=mean（y2（peak_pos-3:peak_pos+3））-baseline;%修正后的峰值
????????half_y=baseline+height_peak/2;
????????half_x1=x2（1）;
????????half_x2=x2（length（x2））;
????????for?jj=2:max_x2;
????????????if?y2（jj-1）<=half_y&&y2（jj）>=half_y;
????????????????half_x1=x2（jj）;
????????????end;
????????end;
????????for?jj=max_x2:length（x2）;
????????????if?y2（jj-1）>=half_y&&y2（jj）<=half_y;
????????????????half_x2=x2（jj）;
????????????end;
????????end;
????????%????half_x1=interp1（y2（1:max_x2）x2（1:max_x2）half_y）;
????????%????half_x2=interp1（y2（max_x2:length（x2））x2（max_x2:length（x2））half_y）;
????????half_width_peak=half_x2-half_x1;
????????%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
????????data2=data;
????????data2（:2）=data2（:2）-baseline;
????????data2（:2）=data2（:2）/height_peak2;
????????data_write2=data2

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????4623??2011-12-14?16:31??get_peak.m