1 | % % % figure(1)
|
---|
2 | % % % plot(SP,SS,'or')
|
---|
3 | % % % hold on
|
---|
4 | % % % plot(aFT)
|
---|
5 | % % % xlim([0 5000])
|
---|
6 | clear all
|
---|
7 | close all
|
---|
8 | Ndet=[1 2 3 4 7 11 14 20 30 40 50 70];
|
---|
9 | s=7
|
---|
10 | load(['result/PrN' num2str(1)],'ChiVk','Corrk','ChiVh','Corrh','RPK','RPH','RSP','Profile');
|
---|
11 | aFT=abs(fft(Profile(s,:)));
|
---|
12 | DetAngl=35+110/(Ndet(s)+1)*(1:Ndet(s));
|
---|
13 |
|
---|
14 | samp1=0.3./(0.05.*(1-cosd(DetAngl)));
|
---|
15 | samp2=0.3./(0.25.*(1-cosd(DetAngl)));
|
---|
16 | samp3=0.3./(1.5.*(1-cosd(DetAngl)));
|
---|
17 | samp=[samp1 samp2 samp3];
|
---|
18 | SFr=sort(samp);
|
---|
19 | SP=unique(round(2100.*SFr./max(25.6459)));%ACHTUNG!!! remove same points
|
---|
20 | SP1=unique(round(2100.*sort(samp1)./max(25.6459)));
|
---|
21 | SP2=unique(round(2100.*sort(samp2)./max(25.6459)));
|
---|
22 | SP3=unique(round(2100.*sort(samp3)./max(25.6459)));
|
---|
23 | SS=aFT(SP);
|
---|
24 | figure(2)
|
---|
25 | hold on
|
---|
26 | plot(SP1,aFT(SP1),'og')
|
---|
27 | plot(SP2,aFT(SP2),'om')
|
---|
28 | plot(SP3,aFT(SP3),'oc')
|
---|
29 | grid on
|
---|
30 | hold on
|
---|
31 | plot(aFT)
|
---|
32 | plot(RSP(s,:),'k')
|
---|
33 | xlim([0 5010])
|
---|
34 | %
|
---|
35 | % s=1
|
---|
36 | % load(['result/PrN' num2str(1)],'ChiVk','Corrk','ChiVh','Corrh','RPK','RPH','RSP','Profile');
|
---|
37 | % close all
|
---|
38 | % figure(3)
|
---|
39 | % plot(Profile(s,:),'k','linew',2)
|
---|
40 | % hold on
|
---|
41 | % %plot(RPK(s,:),'r','linew',2)
|
---|
42 | %
|
---|
43 | % plot(RPH(5,:),'r','linew',2)
|
---|
44 | % plot(RPH(6,:),'g','linew',2)
|
---|
45 | % plot(RPH(10,:),'b','linew',2)
|
---|
46 | % grid on
|
---|
47 | % xlim([32600 32950])
|
---|
48 | % ch=mean(Chi,1);
|
---|
49 | %
|
---|
50 | % bar(ch(1:9))
|
---|
51 | % set(gca,'YGrid','on') % horizontal grid
|
---|
52 | % set(gca,'XTickLabel',{'Lw1', 'Lw5', 'Lw10', 'Lw1X','Lw5X', 'Lw10X', 'Tw10', 'Tw1X', 'Tw5X'},'fontsize',10)
|
---|
53 | % colormap([1,0,0]);
|
---|
54 | %
|
---|
55 | % ylabel('\chi^2')
|
---|
56 | %
|
---|
57 | % title('Comparison of different sample methods, Hilbert')
|
---|
58 | % load SData.mat
|
---|
59 | % [Profile,~]=PR_func(4, 2^(-16));
|
---|
60 | % aFT=abs(fft(Profile));
|
---|
61 | % SFr=sort(Lw1);
|
---|
62 | % SP=unique(round(2100.*SFr./max(Tw10)));%ACHTUNG!!! remove same points
|
---|
63 | % SS=aFT(SP);
|
---|
64 | % figure(1)
|
---|
65 | % plot(SP,SS,'or')
|
---|
66 | % hold on
|
---|
67 | % plot(aFT)
|
---|
68 | % xlim([0 5000]) |
---|