exercise 6.9

c0=298792586   % m/s

 

d=.02;a=.02286;b=.01016      % cavity dimensions [m]

sigma_Al=3.816*10^7    % [S/m]

 

e0=8.854*10^-12     % [F/m] free space permittivity

mu0=4*pi*10^-7       % [H/m] free space permeability

er=1;mur=1

 

p_101=1

 

p_102=2

 

 

etha0=120*pi;etha=etha0/er^.5

 

 

fc=@(n,m,p)...

c0/(2*pi*er*mur)*((n*pi/a)^2+(m*pi/b)^2+(p*pi/d)^2)^.5

 

 

% cut-off

 

fc_101=fc(1,0,1)

 

fc_102=fc(1,0,2)

 

 

 

Rs= @(f,s) (2*pi*f*mu0/(2*s))^.5    % surface resistivity

 

 

 

Rs_101=Rs(fc_101,sigma_Al)

 

Rs_102=Rs(fc_102,sigma_Al)

 

 

 

lambda_101=c0/(er^.5*fc_101)

 

lambda_102=c0/(er^.5*fc_102)

 

 

 

k_101=2*pi*fc_101*(mu0*e0*er*mur)^.5   % wave number

 

k_102=2*pi*fc_102*(mu0*e0*er*mur)^.5

 

 

Qc_101=...

(kc_101*a*d)^3*b*etha/(2*pi^2*Rs_101)*1/...

(2*p_101^2*a^3*b+2*b*d^3+p_101^2*a^3*d+a*d^3)

 

Qc_102=...

(kc_102*a*d)^3*b*etha/(2*pi^2*Rs_102)*1/...

(2*p_102^2*a^3*b+2*b*d^3+p_102^2*a^3*d+a*d^3)