1 | %%plots |
---|
2 | close all; clear all; clc |
---|
3 | figure(2) |
---|
4 | hold on |
---|
5 | |
---|
6 | load smallFc.mat |
---|
7 | plot(nu_filter*1e-9,F24i,'-.ok','linew',2,'MarkerFaceColor','k') |
---|
8 | load smallF.mat |
---|
9 | plot(nu_filter*1e-9,F24i,'-.ob','linew',2,'MarkerFaceColor','b') |
---|
10 | F0=[350 230 180 104 60 24]/438; |
---|
11 | % plot(nu_filter*1e-9,F0,'-or','linew',2) |
---|
12 | |
---|
13 | F37=[79.2 119.3 85.7 93.2 79.4 88 89.7 81.8]; |
---|
14 | F53=[40.7 38.5 51.2 48.3 47.1 34.6 47.5 53.5]; |
---|
15 | F75=[42.3 23.6 45.6 22.2 30.9 15.2 29.8 25.4]; |
---|
16 | F124=[18.5 8.18 17.4 16.7 6.15 9.84 7.0 7.5]; |
---|
17 | F186=[7.8 11.4 9.5 6.4 4.8 4.1 4.9 4.5 ]; |
---|
18 | F368=[11.03 3.57 5.36 5.7 7.1 2.53 4.5 2.3]; |
---|
19 | Frame=[95.8 92.5 92.7 92.05 92 92.1]; |
---|
20 | Fx=[mean(F37) mean(F53) mean(F75) mean(F124) mean(F186) mean(F368)]./mean(Frame); |
---|
21 | eFx=[std(F37) std(F53) std(F75) std(F124) std(F186) std(F368)]./mean(Frame); |
---|
22 | errorbar(nu_filter*1e-9,Fx,eFx,'^r','linew',2) |
---|
23 | grid on |
---|
24 | set(gca,'fontsize',16) |
---|
25 | xlabel('Frequency, [GHz]') |
---|
26 | ylabel('Transmissivity, [1]') |
---|
27 | legend('Calc. Z_0=2.3', 'Calc. Z_0=1.9','Exp. data') |
---|
28 | grid on |
---|
29 | title('24GHz source') |
---|
30 | |
---|
31 | figure(3) |
---|
32 | hold on |
---|
33 | |
---|
34 | load smallFc.mat |
---|
35 | plot(nu_filter*1e-9,F35i,'-.ok','MarkerFaceColor','k','linew',2) |
---|
36 | load smallF.mat |
---|
37 | plot(nu_filter*1e-9,F35i,'-.ob','MarkerFaceColor','b','linew',2) |
---|
38 | |
---|
39 | % F0=[224 122 88 40 19.8 5.12]/215; |
---|
40 | % plot(nu_filter*1e-9,F0,'-or','linew',2) |
---|
41 | |
---|
42 | F37=[74.9 74.4 81.33 83.06 75.6 83.3 90.74 80.42 79.78 75.7]; |
---|
43 | F53=[63.8 58.6 61.2 57.2 69.7 50.2 94.4 52.01 88.2 54.73]; |
---|
44 | F75=[51.5 35.7 47.13 43.57 53.47 30.16 64.12 41.9 52.5 35.53]; |
---|
45 | F124=[33.1 10.79 11.06 24.49 11.35 17.22 12.3 25 14.7 21.6]; |
---|
46 | F186=[11.38 4.3 9.15 21.64 7.89 6.7 10.3 24.4 13 25.3]; |
---|
47 | F368=[2.6 0.5 0.9 2.12 2.24 2.3 2.8 1.26 7.17 1.68]; |
---|
48 | Frame=[74.5 72.6 73.1 72.8 73.56 74.2 74.1]; |
---|
49 | Fx=[mean(F37) mean(F53) mean(F75) mean(F124) mean(F186) mean(F368)]./mean(Frame); |
---|
50 | eFx=[std(F37) std(F53) std(F75) std(F124) std(F186) std(F368)]./mean(Frame); |
---|
51 | |
---|
52 | errorbar(nu_filter*1e-9,Fx,eFx,'^r','linew',2) |
---|
53 | grid on |
---|
54 | set(gca,'fontsize',16) |
---|
55 | xlabel('Frequency, [GHz]') |
---|
56 | ylabel('Transmissivity, [1]') |
---|
57 | legend('Calc. Z_0=2.3', 'Calc. Z_0=1.9','Exp. data') |
---|
58 | title('35GHz source') |
---|
59 | load bigF.mat |
---|
60 | |
---|
61 | |
---|
62 | |
---|
63 | figure(4) |
---|
64 | plot(nu_filter*1e-9,F24i,'-ok','linew',2) |
---|
65 | nu_f0=[37 53 75 124 186 368]; |
---|
66 | F=[416 368 364 390 420 422 382 352 328]/438; |
---|
67 | F0=[350 230 180 104 60 24]/438; |
---|
68 | |
---|
69 | % Freq=[20 23 24 25 26 27 30 40 50]; |
---|
70 | % Vpp0=[116.8 124.0 120.0 126.4 133.6 124.8 115.2 100 97.6]./[153.6 150.4 148.8 147.2 145.2 145.6 144.8 144.8 143.2]; |
---|
71 | hold on |
---|
72 | % plot(Freq,Vpp0,'-ok') |
---|
73 | % plot(Freq,F,'-og') |
---|
74 | % % plot(nu_f0,F0,'-or') |
---|
75 | % ylim([0 1]) |
---|
76 | F20=[89.1 84.8 83.2 81.2]; |
---|
77 | F23=[72 71.4 70.3 72]; |
---|
78 | F24=[67.7 71 70.5 72.7 ]; |
---|
79 | F25=[74.5 74.6 74.2 74.2 ]; |
---|
80 | F26=[81.7 82 81 82.4]; |
---|
81 | F27=[88 84.2 87.8 88.5]; |
---|
82 | F30=[77.6 77.4 76.6 79.7]; |
---|
83 | F40=[61.3 60.21 60.9 62]; |
---|
84 | F50=[66.2 66 67.3 66]; |
---|
85 | Frame=[94.3 93.8 93 92.3 93.7 93.2 93.1]; |
---|
86 | Fx=[mean(F20) mean(F23) mean(F24) mean(F25) mean(F26) mean(F27) mean(F30) mean(F40) mean(F50)]./mean(Frame); |
---|
87 | eFx=[std(F20) std(F23) std(F24) std(F25) std(F26) std(F27) std(F30) std(F40) std(F50)]./mean(Frame); |
---|
88 | errorbar(nu_filter*1e-9,Fx,eFx,'-.^r','linew',2) |
---|
89 | grid on |
---|
90 | set(gca,'fontsize',16) |
---|
91 | xlabel('Frequency, [GHz]') |
---|
92 | ylabel('Transmissivity, [1]') |
---|
93 | legend('Prediction', 'Experiment','location','southeast') |
---|
94 | title('24GHz source') |
---|
95 | load bigF.mat |
---|
96 | ylim([0 1]) |
---|
97 | |
---|
98 | |
---|
99 | figure(5) |
---|
100 | plot(nu_filter*1e-9,F35i,'-ok','linew',2) |
---|
101 | % Vpp3=[112.8 108.8 111.2 107.2 108.0 112.8 107.2 117.6 118.4]./[122.4 121.6 120.8 120.0 120.0 120.0 120.0 119.2 119.2]; |
---|
102 | % Vpp4=[260 240 240 230 238 212 192 246 252]/226 |
---|
103 | hold on |
---|
104 | % plot(Freq,Vpp4,'-or') |
---|
105 | % plot(Freq,Vpp3,'-ok') |
---|
106 | F20=[94.7 95.4 96.2 94.2 ]; |
---|
107 | F23=[ 90.3 93.3 87 86.6 ]; |
---|
108 | F24=[ 87.6 91.6 88.4 89.5 ]; |
---|
109 | F25=[87.7 88 88.3 87.5]; |
---|
110 | F26=[85.2 83.9 83.1 84.6]; |
---|
111 | F27=[84.9 83 80.5 77.3]; |
---|
112 | F30=[82.6 82 83.3 85.5]; |
---|
113 | F40=[87 89.5 89.2 88]; |
---|
114 | F50=[91.4 94.2 94.8 90.7]; |
---|
115 | Frame=[96.3 93.3 92.7 92 91.8 91.6 91 91.2]; |
---|
116 | Fx=[mean(F20) mean(F23) mean(F24) mean(F25) mean(F26) mean(F27) mean(F30) mean(F40) mean(F50)]./mean(Frame); |
---|
117 | eFx=[std(F20) std(F23) std(F24) std(F25) std(F26) std(F27) std(F30) std(F40) std(F50)]./mean(Frame); |
---|
118 | errorbar(nu_filter*1e-9,Fx,eFx,'-.^r','linew',2) |
---|
119 | grid on |
---|
120 | set(gca,'fontsize',16) |
---|
121 | xlabel('Frequency, [GHz]') |
---|
122 | ylabel('Transmissivity, [1]') |
---|
123 | legend('Prediction', 'Experiment','location','southeast') |
---|
124 | title('35GHz source') |
---|
125 | ylim([0 1.1]) |
---|