Line | |
---|
1 | clear all; close all; clc
|
---|
2 | %inductive grating
|
---|
3 | %phys. const
|
---|
4 | c=299792458;%speed of light, m/s
|
---|
5 | sigma=1.45e+6*1e-7*c*c;%stainless steel conductivity, [s]
|
---|
6 | j=0;
|
---|
7 | nu_filter=[20 23 24 25 26 27 30 40 50].*1e+9;
|
---|
8 | nu=(0.001:0.01:1000)*1e+11;
|
---|
9 | for f=nu_filter
|
---|
10 | j=j+1;
|
---|
11 | lamb0=c/f;
|
---|
12 | ratio=0.02;%a/g ratio
|
---|
13 | omega0=(1-0.27*(ratio));
|
---|
14 | g=lamb0*omega0;
|
---|
15 | a=ratio*g;%wire width
|
---|
16 |
|
---|
17 |
|
---|
18 |
|
---|
19 | disp(['g=' num2str(g*1e+3) 'mm a=' num2str(a*1e+3) 'mm nu=' num2str((1-0.27*ratio)*c/g*1e-9) ' GHz'])
|
---|
20 | disp(['g=' num2str(round(g*1e+6)*1e-3) 'mm a=' num2str(round(a*1e+6)*1e-3) 'mm nu=' num2str((1-0.27*ratio)*c/g*1e-9) ' GHz'])
|
---|
21 | disp(['Plate=' num2str(round(g*1e+6)*1e-3-2*round(a*1e+6)*1e-3) ' mm and ' num2str(round(1000/(round(g*1e+6)*1e-3))/10) ' squares on line'])
|
---|
22 | disp([' '])
|
---|
23 | omega=nu*g/c;
|
---|
24 | lamb=g./omega;
|
---|
25 | eta=1/(1-2*ratio);
|
---|
26 | R=eta./2.*sqrt(c./lamb./sigma);
|
---|
27 | Om=lamb0./lamb-lamb./lamb0;
|
---|
28 | Z0=(2*log(csc(a*pi/(2*g))));
|
---|
29 | t2=(R.^2+(Z0./Om).^2)./((1+R).^2+(Z0./Om).^2);
|
---|
30 |
|
---|
31 |
|
---|
32 | semilogy(nu/1e+9,t2,'color',rand([1 3]),'linew',2)
|
---|
33 |
|
---|
34 | hold on
|
---|
35 | end
|
---|
36 | set(gca,'fontsize',16)
|
---|
37 | xlabel('Frequency, [GHz]')
|
---|
38 | ylabel('Transmissivity')
|
---|
39 | legend('20 GHz','23 GHz','24 GHz','25 GHz','26 GHz','27 GHz','30 GHz','40 GHz','50 GHz')
|
---|
40 | grid on
|
---|
41 | xlim([10 100])
|
---|
Note: See
TracBrowser
for help on using the repository browser.