function main2010715.m %% Data load load ('/Users/koji/Documents/MATLAB/gw_analysis/raw_data/data_samp10cut1.dat'); data = (data_samp10cut1 - mean(data_samp10cut1)); samp = 10; cal = 4/(4*pi*0.050/1064e-9)*1/(0.097*2); l1 = length(data); t = (1:l1)/samp; %% Rejection broken-data data = data(601:2.6e+5); [spe,fr]=fft_data2(data,2^10,samp,1); [fr,amp_s,pha_s,amp_l,pha_l]=filter_design20090119(fr); %% loglog(fr,abs(spe').*cal*(1+amp_l/28.9)); % Select 0.2 Hz A=abs(spe(21,:))*cal.*(1+amp_l(21)/28.9); % Data quality cut by eye B=A(111:230); % Data quality cut large 3% data out=B<0.6e-8; C=B(out).^2; % plot histgrum [y,x]=hist(C,2^4); m = mean(C) s = std(C) plot(x,y,x,exp(-1.8e+17*x+3.8))