多層誘電体板による反射・透過特性
モジュールのインポートや初期設定など
import numpy as np
import matplotlib.pyplot as plt
np.set_printoptions(precision=4)
CCC = 299.79
qq = 180.0/np.pi
import scienceplots
#Warning : As of version 2.0.0, you need to add import scienceplots before setting the style (plt.style.use('science')).
plt.style.use(['science', 'notebook'])
#plt.rcParams['font.family'] = 'Times New Roman' # font familyの設定
plt.rcParams['font.family'] = 'serif' # font familyの設定
plt.rcParams['mathtext.fontset'] = 'cm' # math fontの設定
plt.rcParams["font.size"] = 10
ユーザ関数
多層誘電体板の反射・透過係数を求める関数は，
def func_slab_tr(nt, ep, td, t, sk, th, ic, ip):
zk00 = ep[0]*(1.0-1j*td[0])
zk00 = sk*np.sqrt(zk00)
zsin00 = np.sin(th)
zk1 = ep[nt+1]*(1.0-1j*td[nt+1])
zk1 = sk*np.sqrt(zk1)
zsin1 = zk00/zk1*zsin00
zcos1 = np.sqrt(1.0-zsin1**2)
za1, zf1 = 1.0, 1.0
zb1, zg1 = 0.0, 0.0
ttl = 0.0
zzzg = 1.0
for i in range(nt,-1,-1):
zk2 = zk1
zsin2, zcos2 = zsin1, zcos1
za2, zb2, zf2, zg2 = za1, zb1, zf1, zg1
zk1 = ep[i]*(1.0-1j*td[i])
zk1 = sk*np.sqrt(zk1)
zsin1 = zk2/zk1*zsin2
zcos1 = np.sqrt(1.0-zsin1**2)
zcos21 = zcos2/zcos1
zy2 = zcos21*(zk2/zk1)
zw2 = zcos21*(zk1/zk2)
zgm1 = 1j*t[i]*zk1*zcos1
za1 = np.exp( zgm1)*(za2*(1.0+zy2)+zb2*(1.0-zy2))/2.0
zb1 = np.exp(-zgm1)*(za2*(1.0-zy2)+zb2*(1.0+zy2))/2.0
zf1 = np.exp( zgm1)*(zf2*(1.0+zw2)+zg2*(1.0-zw2))/2.0*(zk2/zk1)
zg1 = np.exp(-zgm1)*(zf2*(1.0-zw2)+zg2*(1.0+zw2))/2.0*(zk2/zk1)
ttl = ttl+t[i]
zzzg = zzzg*np.exp(-1j*zk1*t[i]*zcos1)
ztn = 1.0/za1# transmission coefficient for perpendicular polarization
ztp = 1.0/zf1# transmission coefficient for parallel polarization
zrn = zb1/za1# reflection coefficient for perpendicular polarization
zrp = zg1/zf1# reflection coefficient for parallel polarization
if ip==1:
zrp = -zrp
if ic==0:
zzz = 1j*sk*ttl*np.cos(th)
zzz = np.exp(zzz)
ztn = ztn*zzz
ztp = ztp*zzz
zzzg = zzzg*zzz
cth = np.cos(th)
n2sth2 = np.sqrt(ep[nt+1]-np.sin(th)**2);
cth12 = n2sth2/cth
ztn = ztn*np.sqrt(cth12)
ztp = ztp*np.sqrt(cth12)
return zrn,zrp,ztn,ztp,zzzg
電磁界の複素係数から絶対値，デジベル(dB)，偏角(deg)を計算するため，
def dbdeg(z):
q = 180.0/np.pi
amp = np.abs(z)
db = 20.0*np.log10(amp)
rad = np.angle(z)
# rad = np.unwrap(rad, discont=0.95*np.pi)
# rad = np.unwrap(rad)
deg = rad*q
return amp,db,deg
多層誘電体の寸法，比誘電率，誘電正接をフォーマット付で出力するため，
def output_dielectric(nt,t,ep,td):
print(f"{'No.':>4s} {'Thickness[mm]':>20s} {'Dielectric constant':>20s} {'Loss tangent':>20s}")
for i in range(nt+2):
print(f"{i:4d} {t[i]:20.3f} {ep[i]:20.3f} {td[i]:20.3f}")
多層誘電体による反射・透過係数の計算値をフォーマット付で出力するため，
def output_tr(aa,cha,ndel,ztn,ztp,zrn,zrp,zzzg):
atn,dtn,ptn = dbdeg(ztn)
atp,dtp,ptp = dbdeg(ztp)
arn,drn,prn = dbdeg(zrn)
arp,drp,prp = dbdeg(zrp)
atg,dtg,ptg = dbdeg(zzzg)
print(f"{'':>7s} {'Transmission coefficients':>33s} {'Reflection coefficients':>33s}")
print(f"{'':>8s} {'Perpendicular':>16s} {'Parallel':>16s} {'Perpendicular':>16s} {'Parallel':>16s}")
print(f"{cha:>8s} {'Amp[dB]':>8s} {'Ph[deg]':>8s} {'Amp[dB]':>8s} {'Ph[deg]':>8s} {'Amp[dB]':>8s} {'Ph[deg]':>8s} {'Amp[dB]':>8s} {'Ph[deg]':>8s}")
#Y1 = np.array([all,atn,ptn,atp,ptp,arn,prn,arp,prp,ptg])
#YY1 = Y1.T
for i in range(0,len(aa),ndel):
# print(('{:8.3f} '*len(aa)).format(*YY1[i]))
print(f"{aa[i]:8.2f} {dtn[i]:8.3f} {ptn[i]:8.2f} {dtp[i]:8.3f} {ptp[i]:8.2f} {drn[i]:8.3f} {prn[i]:8.2f} {drp[i]:8.3f} {prp[i]:8.2f} {ptg[i]:8.2f}")
return dtn,ptn,dtp,ptp,drn,prn,drp,prp,dtg,ptg
反射・透過係数の入射角特性を図示するため，
def figure4a(tthdeg,drn1,drp1,prn1,prp1,dtn1,dtp1,ptn1,ptp1,fname):
ymin = -40.0
dymin = np.abs(ymin)/15.0
fig, ax = plt.subplots(nrows=2, ncols=2, figsize=(11, 7)) # グラフ領域の作成，nrows=縦に並べる数，ncols=横に並べる数
ax[0][0].minorticks_on()
ax[0][0].set_xlim(0, 90)
ax[0][0].set_xticks(np.arange(0.0, 100.0, step=10))
ax[0][0].set_ylim(ymin,0.0)
ax[0][0].set_ylabel('Amplitude [dB]')
ax[0][0].plot(tthdeg,drn1, color='blue')
ax[0][0].plot(tthdeg,drp1, color='tab:green')
ax[1][0].set_xlim(0, 90)
ax[1][0].set_ylim(-180, 180)
ax[1][0].set_xticks(np.arange(0.0, 100.0, step=10))
ax[1][0].set_yticks(np.arange(-180.0, 270.0, step=90))
ax[1][0].set_ylabel('Phase [deg]')
ax[1][0].plot(tthdeg,prn1,'.', color='tab:blue')
ax[1][0].plot(tthdeg,prp1,'.', color='tab:green')
ax[0][1].minorticks_on()
ax[0][1].set_xlim(0, 90)
ax[0][1].set_xticks(np.arange(0.0, 100.0, step=10))
ax[0][1].set_ylim(ymin,0.0)
ax[0][1].set_ylabel('Amplitude [dB]')
ax[0][1].plot(tthdeg,dtn1, color='tab:orange')
ax[0][1].plot(tthdeg,dtp1, color='tab:red')
ax[1][1].set_xlim(0, 90)
ax[1][1].set_ylim(-180, 180)
ax[1][1].set_xticks(np.arange(0.0, 100.0, step=10))
ax[1][1].set_yticks(np.arange(-180.0, 270.0, step=90))
ax[1][1].set_ylabel('Phase [deg]')
ax[1][1].plot(tthdeg,ptn1,'.', color='tab:orange')
ax[1][1].plot(tthdeg,ptp1,'.', color='tab:red')
[[ax[i][j].set_xlabel('Incident angle [deg]') for i in range(2)] for j in range(2)]
[ax[i][0].legend([r'$R_\perp$ (TE, perpendicular-pol.)','$R_\parallel$ (TM, parallel-pol.)'], \
ncol=1, loc='lower right', fancybox=False, frameon=True, fontsize=10) \
for i in range(2)]
[ax[i][1].legend([r'$T_\perp$ (TE, perpendicular-pol.)','$T_\parallel$ (TM, parallel-pol.)'], \
ncol=1, loc='lower right', fancybox=False, frameon=True, fontsize=10) \
for i in range(2)]
fig.tight_layout()
fig.savefig(fname+'.pdf')
plt.show()
反射・透過係数の周波数特性を図示するため，
def figure4b(tthdeg,ffreq,drn2,drp2,prn2,prp2,dtn2,dtp2,ptn2,ptp2,fname):
ymin = -40.0
dymin = np.abs(ymin)/15.0
fig, ax = plt.subplots(2,2,figsize=(11, 7)) # グラフ領域の作成，nrows=縦に並べる数，ncols=横に並べる数
ax[0,0].minorticks_on()
ax[0,0].set_xlim(0, 90)
ax[0,0].set_xticks(np.arange(0.0, 100.0, step=10))
ax[0,0].set_ylim(ymin,0.0)
ax[0,0].set_xlabel('Incident angle [deg]')
ax[0,0].set_ylabel('Amplitude [dB]')
ax[0,0].plot(tthdeg,drn1, tthdeg,drp1, tthdeg,dtn1, tthdeg,dtp1)
ax[0,0].legend([r'$R_\perp$ (TE, perpendicular-pol.)','$R_\parallel$ (TM, parallel-pol.)','$T_\perp$ (TE, perpendicular-pol.)','$T_\parallel$ (TM, parallel-pol.)'], ncol=1, loc='lower right', fancybox=False, frameon=True, fontsize=10)
ax[1,0].set_xlim(0, 90)
ax[1,0].set_ylim(-180, 180)
ax[1,0].set_xticks(np.arange(0.0, 100.0, step=10))
ax[1,0].set_yticks(np.arange(-180.0, 270.0, step=90))
ax[1,0].set_xlabel('Incident angle [deg]')
ax[1,0].set_ylabel('Phase [deg]')
ax[1,0].plot(tthdeg,prn1,'.', tthdeg,prp1,'.', tthdeg,ptn1,'.', tthdeg,ptp1,'.')
ax[1,0].legend([r'$R_\perp$ (TE, Perpendicular-pol.)','$R_\parallel$ (TM, parallel-pol.)','$T_\perp$ (TE, Perpendicular-pol.)','$T_\parallel$ (TM, parallel-pol.)'], ncol=1, loc='lower right', fancybox=False, frameon=True, fontsize=10)
ax[0,1].set_xlim(ffreq[0], ffreq[-1])
ax[0,1].set_ylim(ymin,0.0)
ax[0,1].set_xlabel('Frequency [GHz]')
ax[0,1].set_ylabel('Amplitude [dB]')
ax[0,1].plot(ffreq,drn2, ffreq,drp2, ffreq,dtn2, ffreq,dtp2)
ax[0,1].legend([r'$R_\perp$ (TE, perpendicular-pol.)','$R_\parallel$ (TM, parallel-pol.)','$T_\perp$ (TE, perpendicular-pol.)','$T_\parallel$ (TM, parallel-pol.)'], ncol=1, loc='lower right', fancybox=False, frameon=True, fontsize=10)
ax[1,1].set_xlim(ffreq[0], ffreq[-1])
ax[1,1].set_ylim(-180, 180)
ax[1,1].set_yticks(np.arange(-180.0, 270.0, step=90))
ax[1,1].set_xlabel('Frequency [GHz]')
ax[1,1].set_ylabel('Phase [deg]')
ax[1,1].plot(ffreq,prn2,'.', ffreq,prp2,'.', ffreq,ptn2,'.', ffreq,ptp2,'.',)
ax[1,1].legend([r'$R_\perp$ (TE, perpendicular-pol. )','$R_\parallel$ (TM, parallel-pol.)','$T_\perp$ (TE, perpendicular-pol.)','$T_\parallel$ (TM, parallel-pol.)'], ncol=1, loc='lower right', fancybox=False, frameon=True, fontsize=10)
fig.tight_layout()
fig.savefig(fname+'.pdf')
plt.show()
計算条件
周波数については，
start_f, stop_f, nfreq = 5.0, 30.0, 25*20+1# 周波数の設定
ffreq = np.linspace(start_f, stop_f, nfreq)
wwl = CCC/ffreq# 波長[mm]
ssk = 2*np.pi/wwl# 波数
入射角については，
tthdeg = np.linspace(0.0, 90.0, nth+1)# incident angle
tth = tthdeg/qq
中心周波数ほかの設定は，
thdeg = 0.0
freq = 10.0# 中心周波数
nu_th = np.int64(np.round(nth)*thdeg/90)
nu_freq = np.int64(np.round((nfreq-1)*(freq-start_f)/(stop_f-start_f)))
ic = 0# =0:規格化しない, =1:誘電体スラブがない場合の位相で規格化
ip = 1# =1:式のとおり, =-1:reflection for parallel-pol. wave
誘電体境界面による反射・透過
誘電体境界面による反射・透過係数を求めるための計算条件は，
nt = 0# 誘電体スラブの層数
ep_end = 9.0# 終端の媒質の比誘電率
ep = np.array([1.0, ep_end]) # dielectric constant
td = np.array([0.0, 0.0]) # loss tangent
t = np.array([0.0, 0.0]) # thickness [mm], 境界面での反射・透過を求める場合
output_dielectric(nt,t,ep,td)
No. Thickness[mm] Dielectric constant Loss tangent
0 0.000 1.000 0.000
1 0.000 9.000 0.000
誘電体境界面による反射・透過係数を求め，計算値を出力および図示すると，
zrn,zrp,ztn,ztp,zzzg = func_slab_tr(nt, ep, td, t, ssk[nu_freq], tth, ic, ip)
print(f'Frequency={ffreq[nu_freq]:.1f}[GHz], wavelength={wwl[nu_freq]:.1f}[mm]')
dtn1,ptn1,dtp1,ptp1,drn1,prn1,drp1,prp1,dtg1,ptg1 = output_tr(tthdeg,'th[deg]',10,ztn,ztp,zrn,zrp,zzzg)
figure4a(tthdeg,drn1,drp1,prn1,prp1,dtn1,dtp1,ptn1,ptp1,'tr_nt0_angle')
Frequency=10.0[GHz], wavelength=30.0[mm]
Transmission coefficients Reflection coefficients
Perpendicular Parallel Perpendicular Parallel
th[deg] Amp[dB] Ph[deg] Amp[dB] Ph[deg] Amp[dB] Ph[deg] Amp[dB] Ph[deg]
0.00 -1.249 0.00 -1.249 0.00 -6.021 180.00 -6.021 -180.00 0.00
5.00 -1.257 0.00 -1.242 0.00 -5.999 180.00 -6.043 -180.00 0.00
10.00 -1.279 0.00 -1.220 0.00 -5.933 180.00 -6.110 -180.00 0.00
15.00 -1.317 0.00 -1.183 0.00 -5.823 180.00 -6.226 -180.00 0.00
20.00 -1.373 0.00 -1.131 0.00 -5.670 180.00 -6.396 -180.00 0.00
25.00 -1.448 0.00 -1.064 0.00 -5.475 180.00 -6.630 -180.00 0.00
30.00 -1.545 0.00 -0.982 0.00 -5.239 180.00 -6.940 -180.00 0.00
35.00 -1.668 0.00 -0.884 0.00 -4.963 180.00 -7.348 -180.00 0.00
40.00 -1.823 0.00 -0.772 0.00 -4.649 180.00 -7.884 -180.00 0.00
45.00 -2.018 0.00 -0.645 0.00 -4.299 180.00 -8.598 -180.00 0.00
50.00 -2.262 0.00 -0.508 0.00 -3.915 180.00 -9.574 -180.00 0.00
55.00 -2.571 0.00 -0.362 0.00 -3.499 180.00 -10.970 -180.00 0.00
60.00 -2.966 0.00 -0.217 0.00 -3.055 180.00 -13.120 -180.00 0.00
65.00 -3.480 0.00 -0.088 0.00 -2.586 180.00 -16.977 -180.00 0.00
70.00 -4.171 0.00 -0.006 0.00 -2.096 180.00 -28.255 -180.00 0.00
75.00 -5.140 0.00 -0.043 0.00 -1.587 180.00 -20.103 -0.00 0.00
80.00 -6.623 0.00 -0.379 0.00 -1.066 180.00 -10.780 -0.00 0.00
85.00 -9.357 0.00 -1.677 0.00 -0.535 180.00 -4.945 -0.00 0.00
90.00 nan nan nan nan nan nan nan nan 0.00
整合層による反射・透過
整合層による反射・透過係数を求めるための計算条件は，
nt = 1# 誘電体スラブの層数
ep_end = 9.0# 終端の媒質の比誘電率
epl, tdl, tl = [3.0], [0.0], [17.3084/4.0]# 1/4波長整合層
ep = np.array([1.0] + epl + [ep_end])
td = np.array([0.0] + tdl + [0.0])
t = np.array([0.0] + tl + [0.0])
output_dielectric(nt,t,ep,td)
No. Thickness[mm] Dielectric constant Loss tangent
0 0.000 1.000 0.000
1 4.327 3.000 0.000
2 0.000 9.000 0.000
整合層による反射・透過係数の入射角特性を求め，計算値を出力および図示すると，
zrn,zrp,ztn,ztp,zzzg = func_slab_tr(nt, ep, td, t, ssk[nu_freq], tth, ic, ip)
print(f'Frequency={ffreq[nu_freq]:.1f}[GHz], wavelength={wwl[nu_freq]:.1f}[mm]')
dtn1,ptn1,dtp1,ptp1,drn1,prn1,drp1,prp1,dtg1,ptg1 = output_tr(tthdeg,'th[deg]',10,ztn,ztp,zrn,zrp,zzzg)
figure4a(tthdeg,drn1,drp1,prn1,prp1,dtn1,dtp1,ptn1,ptp1,'tr_nt1_matching_angle')
Frequency=10.0[GHz], wavelength=30.0[mm]
Transmission coefficients Reflection coefficients
Perpendicular Parallel Perpendicular Parallel
th[deg] Amp[dB] Ph[deg] Amp[dB] Ph[deg] Amp[dB] Ph[deg] Amp[dB] Ph[deg]
0.00 -0.000 -38.04 -0.000 -38.04 -121.381 -90.00 -121.381 -90.00 -38.04
5.00 -0.000 -38.10 -0.000 -38.10 -56.885 126.55 -56.923 53.53 -38.12
10.00 -0.000 -38.30 -0.000 -38.31 -44.777 127.15 -44.928 53.15 -38.37
15.00 -0.001 -38.63 -0.001 -38.65 -37.627 128.17 -37.962 52.44 -38.80
20.00 -0.002 -39.10 -0.002 -39.15 -32.481 129.60 -33.064 51.37 -39.40
25.00 -0.006 -39.70 -0.005 -39.83 -28.415 131.41 -29.298 49.87 -40.19
30.00 -0.014 -40.43 -0.010 -40.71 -25.018 133.59 -26.236 47.85 -41.17
35.00 -0.027 -41.30 -0.019 -41.82 -22.069 136.12 -23.635 45.23 -42.36
40.00 -0.050 -42.30 -0.032 -43.19 -19.439 138.99 -21.333 41.99 -43.77
45.00 -0.087 -43.45 -0.052 -44.86 -17.042 142.17 -19.212 38.12 -45.42
50.00 -0.146 -44.72 -0.084 -46.85 -14.821 145.63 -17.179 33.73 -47.32
55.00 -0.238 -46.14 -0.134 -49.21 -12.736 149.35 -15.164 28.98 -49.49
60.00 -0.381 -47.68 -0.217 -51.96 -10.755 153.29 -13.120 24.12 -51.96
65.00 -0.606 -49.35 -0.357 -55.15 -8.854 157.44 -11.028 19.35 -54.74
70.00 -0.963 -51.15 -0.601 -58.80 -7.016 161.75 -8.884 14.86 -57.83
75.00 -1.551 -53.05 -1.046 -62.96 -5.224 166.19 -6.696 10.71 -61.26
80.00 -2.599 -55.06 -1.915 -67.65 -3.465 170.74 -4.479 6.92 -65.01
85.00 -4.839 -57.16 -3.942 -72.93 -1.727 175.35 -2.243 3.39 -69.10
90.00 nan nan nan nan nan nan nan nan -73.48
zrn,zrp,ztn,ztp,zzzg = func_slab_tr(nt, ep, td, t, ssk, tth[nu_th], ic, ip)
print(f'Incident angle={tthdeg[nu_th]:.1f}[deg]')
dtn2,ptn2,dtp2,ptp2,drn2,prn2,drp2,prp2,dtg2,ptg2 = output_tr(ffreq,'f[GHz]',20,ztn,ztp,zrn,zrp,zzzg)
figure4b(tthdeg,ffreq,drn2,drp2,prn2,prp2,dtn2,dtp2,ptn2,ptp2,'tr_nt1_matching_angle_freq')
Incident angle=0.0[deg]
Transmission coefficients Reflection coefficients
Perpendicular Parallel Perpendicular Parallel
f[GHz] Amp[dB] Ph[deg] Amp[dB] Ph[deg] Amp[dB] Ph[deg] Amp[dB] Ph[deg]
5.00 -0.669 -14.91 -0.669 -14.91 -8.451 139.11 -8.451 139.11 -19.02
6.00 -0.473 -18.83 -0.473 -18.83 -9.860 129.99 -9.860 129.99 -22.82
7.00 -0.289 -23.16 -0.289 -23.16 -11.919 120.47 -11.919 120.47 -26.63
8.00 -0.136 -27.87 -0.136 -27.87 -15.108 110.57 -15.108 110.57 -30.43
9.00 -0.035 -32.87 -0.035 -32.87 -20.920 100.36 -20.920 100.36 -34.23
10.00 -0.000 -38.04 -0.000 -38.04 -121.381 -90.00 -121.381 -90.00 -38.04
11.00 -0.035 -43.21 -0.035 -43.21 -20.920 -100.36 -20.920 -100.36 -41.84
12.00 -0.136 -48.21 -0.136 -48.21 -15.108 -110.57 -15.108 -110.57 -45.65
13.00 -0.289 -52.92 -0.289 -52.92 -11.919 -120.47 -11.919 -120.47 -49.45
14.00 -0.473 -57.25 -0.473 -57.25 -9.860 -129.99 -9.860 -129.99 -53.25
15.00 -0.669 -61.16 -0.669 -61.16 -8.451 -139.11 -8.451 -139.11 -57.06
16.00 -0.857 -64.68 -0.857 -64.68 -7.469 -147.82 -7.469 -147.82 -60.86
17.00 -1.020 -67.86 -1.020 -67.86 -6.793 -156.19 -6.793 -156.19 -64.67
18.00 -1.144 -70.75 -1.144 -70.75 -6.352 -164.28 -6.352 -164.28 -68.47
19.00 -1.223 -73.46 -1.223 -73.46 -6.102 -172.19 -6.102 -172.19 -72.27
20.00 -1.249 -76.08 -1.249 -76.08 -6.021 180.00 -6.021 180.00 -76.08
21.00 -1.223 -78.69 -1.223 -78.69 -6.102 172.19 -6.102 172.19 -79.88
22.00 -1.144 -81.40 -1.144 -81.40 -6.352 164.28 -6.352 164.28 -83.68
23.00 -1.020 -84.30 -1.020 -84.30 -6.793 156.19 -6.793 156.19 -87.49
24.00 -0.857 -87.47 -0.857 -87.47 -7.469 147.82 -7.469 147.82 -91.29
25.00 -0.669 -90.99 -0.669 -90.99 -8.451 139.11 -8.451 139.11 -95.10
26.00 -0.473 -94.91 -0.473 -94.91 -9.860 129.99 -9.860 129.99 -98.90
27.00 -0.289 -99.23 -0.289 -99.23 -11.919 120.47 -11.919 120.47 -102.70
28.00 -0.136 -103.94 -0.136 -103.94 -15.108 110.57 -15.108 110.57 -106.51
29.00 -0.035 -108.95 -0.035 -108.95 -20.920 100.36 -20.920 100.36 -110.31
30.00 -0.000 -114.12 -0.000 -114.12 -111.838 -90.00 -111.838 -90.00 -114.12
半波長誘電体板
半波長誘電体板による反射・透過係数を求めるための計算条件は，
nt = 1# 誘電体スラブの層数
ep_end = 1.0# 終端の媒質の比誘電率
epl, tdl, tl = [2.25], [0.0], [19.986/2.0]# 半波長レドーム
ep = np.array([1.0] + epl + [ep_end])
td = np.array([0.0] + tdl + [0.0])
t = np.array([0.0] + tl + [0.0])
output_dielectric(nt,t,ep,td)
No. Thickness[mm] Dielectric constant Loss tangent
0 0.000 1.000 0.000
1 9.993 2.250 0.000
2 0.000 1.000 0.000
半波長誘電体板による反射・透過係数の垂直入射における周波数特性を求め，計算値を出力および図示すると，
zrn,zrp,ztn,ztp,zzzg = func_slab_tr(nt, ep, td, t, ssk, tth[nu_th], ic, ip)
print(f'Incident angle={tthdeg[nu_th]:.1f}[deg]')
dtn2,ptn2,dtp2,ptp2,drn2,prn2,drp2,prp2,dtg2,ptg2 = output_tr(ffreq,'f[GHz]',20,ztn,ztp,zrn,zrp,zzzg)
figure4b(tthdeg,ffreq,drn2,drp2,prn2,prp2,dtn2,dtp2,ptn2,ptp2,'tr_nt1_half_angle_freq')
Incident angle=0.0[deg]
Transmission coefficients Reflection coefficients
Perpendicular Parallel Perpendicular Parallel
f[GHz] Amp[dB] Ph[deg] Amp[dB] Ph[deg] Amp[dB] Ph[deg] Amp[dB] Ph[deg]
5.00 -0.695 -30.00 -0.695 -30.00 -8.299 -180.00 -8.299 -180.00 -30.00
6.00 -0.633 -34.70 -0.633 -34.70 -8.674 163.30 -8.674 163.30 -36.00
7.00 -0.467 -39.85 -0.467 -39.85 -9.912 146.15 -9.912 146.15 -42.00
8.00 -0.253 -45.79 -0.253 -45.79 -12.473 128.21 -12.473 128.21 -48.00
9.00 -0.071 -52.61 -0.071 -52.61 -17.876 109.39 -17.876 109.39 -54.00
10.00 0.000 -60.00 0.000 -60.00 -322.453 42.59 -323.628 129.21 -60.00
11.00 -0.071 -67.39 -0.071 -67.39 -17.876 -109.39 -17.876 -109.39 -66.00
12.00 -0.253 -74.21 -0.253 -74.21 -12.473 -128.21 -12.473 -128.21 -72.00
13.00 -0.467 -80.15 -0.467 -80.15 -9.912 -146.15 -9.912 -146.15 -78.00
14.00 -0.633 -85.30 -0.633 -85.30 -8.674 -163.30 -8.674 -163.30 -84.00
15.00 -0.695 -90.00 -0.695 -90.00 -8.299 -180.00 -8.299 -180.00 -90.00
16.00 -0.633 -94.70 -0.633 -94.70 -8.674 163.30 -8.674 163.30 -96.00
17.00 -0.467 -99.85 -0.467 -99.85 -9.912 146.15 -9.912 146.15 -102.00
18.00 -0.253 -105.79 -0.253 -105.79 -12.473 128.21 -12.473 128.21 -108.00
19.00 -0.071 -112.61 -0.071 -112.61 -17.876 109.39 -17.876 109.39 -114.00
20.00 0.000 -120.00 0.000 -120.00 -318.698 61.46 -319.155 112.19 -120.00
21.00 -0.071 -127.39 -0.071 -127.39 -17.876 -109.39 -17.876 -109.39 -126.00
22.00 -0.253 -134.21 -0.253 -134.21 -12.473 -128.21 -12.473 -128.21 -132.00
23.00 -0.467 -140.15 -0.467 -140.15 -9.912 -146.15 -9.912 -146.15 -138.00
24.00 -0.633 -145.30 -0.633 -145.30 -8.674 -163.30 -8.674 -163.30 -144.00
25.00 -0.695 -150.00 -0.695 -150.00 -8.299 180.00 -8.299 180.00 -150.00
26.00 -0.633 -154.70 -0.633 -154.70 -8.674 163.30 -8.674 163.30 -156.00
27.00 -0.467 -159.85 -0.467 -159.85 -9.912 146.15 -9.912 146.15 -162.00
28.00 -0.253 -165.79 -0.253 -165.79 -12.473 128.21 -12.473 128.21 -168.00
29.00 -0.071 -172.61 -0.071 -172.61 -17.876 109.39 -17.876 109.39 -174.00
30.00 0.000 -180.00 0.000 -180.00 -315.765 70.07 -315.992 105.21 -180.00
薄皮サンドイッチ型
薄皮サンドイッチ型に対する反射・透過係数を求める計算条件は，
nt = 3# 誘電体スラブの層数
ep_end = 1.0# 終端の媒質の比誘電率
epl, tdl, tl = [4.0, 1.15, 4.0], [0.015, 0.002, 0.015], [14.9895*0.03, 14.9895*0.47, 14.9895*0.03]# 薄皮サンドイッチ型（対称構造）A.W.Rudge,"Handbook of Antenna Design," Vol.2, p.471.
ep = np.array([1.0] + epl + [ep_end])
td = np.array([0.0] + tdl + [0.0])
t = np.array([0.0] + tl + [0.0])
output_dielectric(nt,t,ep,td)
No. Thickness[mm] Dielectric constant Loss tangent
0 0.000 1.000 0.000
1 0.450 4.000 0.015
2 7.045 1.150 0.002
3 0.450 4.000 0.015
4 0.000 1.000 0.000
薄皮サンドイッチ型に対する反射・透過係数の垂直入射における周波数特性を求め，計算値を出力および図示すると，
zrn,zrp,ztn,ztp,zzzg = func_slab_tr(nt, ep, td, t, ssk, tth[nu_th], ic, ip)
print(f'Incident angle={tthdeg[nu_th]:.1f}[deg]')
dtn2,ptn2,dtp2,ptp2,drn2,prn2,drp2,prp2,dtg2,ptg2 = output_tr(ffreq,'f[GHz]',20,ztn,ztp,zrn,zrp,zzzg)
figure4b(tthdeg,ffreq,drn2,drp2,prn2,prp2,dtn2,dtp2,ptn2,ptp2,'tr_nt3_thin_skin_angle_freq')
Incident angle=0.0[deg]
Transmission coefficients Reflection coefficients
Perpendicular Parallel Perpendicular Parallel
f[GHz] Amp[dB] Ph[deg] Amp[dB] Ph[deg] Amp[dB] Ph[deg] Amp[dB] Ph[deg]
5.00 -0.109 -10.60 -0.109 -10.60 -17.363 -149.25 -17.363 -149.25 -8.46
6.00 -0.113 -12.59 -0.113 -12.59 -17.436 -160.70 -17.436 -160.70 -10.15
7.00 -0.103 -14.62 -0.103 -14.62 -18.377 -172.15 -18.377 -172.15 -11.85
8.00 -0.084 -16.78 -0.084 -16.78 -20.559 176.41 -20.559 176.41 -13.54
9.00 -0.065 -19.09 -0.065 -19.09 -25.310 165.27 -25.310 165.27 -15.23
10.00 -0.060 -21.58 -0.060 -21.58 -49.046 177.25 -49.046 177.25 -16.92
11.00 -0.084 -24.20 -0.084 -24.20 -24.853 -41.98 -24.853 -41.98 -18.62
12.00 -0.146 -26.87 -0.146 -26.87 -18.166 -53.33 -18.166 -53.33 -20.31
13.00 -0.251 -29.48 -0.251 -29.48 -14.348 -65.18 -14.348 -65.18 -22.00
14.00 -0.396 -31.89 -0.396 -31.89 -11.786 -76.98 -11.786 -76.98 -23.69
15.00 -0.565 -34.00 -0.565 -34.00 -9.980 -88.52 -9.980 -88.52 -25.39
16.00 -0.741 -35.75 -0.741 -35.75 -8.698 -99.72 -8.698 -99.72 -27.08
17.00 -0.901 -37.14 -0.901 -37.14 -7.813 -110.57 -7.813 -110.57 -28.77
18.00 -1.023 -38.21 -1.023 -38.21 -7.250 -121.10 -7.250 -121.10 -30.46
19.00 -1.091 -39.08 -1.091 -39.08 -6.969 -131.42 -6.969 -131.42 -32.16
20.00 -1.092 -39.88 -1.092 -39.88 -6.958 -141.66 -6.958 -141.66 -33.85
21.00 -1.022 -40.77 -1.022 -40.77 -7.242 -151.98 -7.242 -151.98 -35.54
22.00 -0.885 -41.96 -0.885 -41.96 -7.886 -162.54 -7.886 -162.54 -37.23
23.00 -0.697 -43.62 -0.697 -43.62 -9.034 -173.51 -9.034 -173.51 -38.93
24.00 -0.485 -45.93 -0.485 -45.93 -10.981 175.03 -10.981 175.03 -40.62
25.00 -0.294 -48.99 -0.294 -48.99 -14.457 163.27 -14.457 163.27 -42.31
26.00 -0.176 -52.76 -0.176 -52.76 -22.407 153.52 -22.407 153.52 -44.00
27.00 -0.182 -57.00 -0.182 -57.00 -25.470 -54.20 -25.470 -54.20 -45.69
28.00 -0.343 -61.31 -0.343 -61.31 -14.653 -62.08 -14.653 -62.08 -47.39
29.00 -0.653 -65.21 -0.653 -65.21 -10.182 -74.48 -10.182 -74.48 -49.08
30.00 -1.070 -68.34 -1.070 -68.34 -7.578 -86.70 -7.578 -86.70 -50.77
1/4波長表皮サンドイッチ型
1/4波長表皮サンドイッチ型に対する反射・透過係数を求めるための計算条件は，
nt = 3# 誘電体スラブの層数
ep_end = 1.0# 終端の媒質の比誘電率
epl, tdl, tl = [3.0, 9.0, 3.0], [0.0, 0.0, 0.0], [17.3804*0.25, 9.993*0.5, 17.3804*0.25]# 1/4 波長表皮サンドイッチ型（対称構造）
ep = np.array([1.0] + epl + [ep_end])
td = np.array([0.0] + tdl + [0.0])
t = np.array([0.0] + tl + [0.0])
output_dielectric(nt,t,ep,td)
No. Thickness[mm] Dielectric constant Loss tangent
0 0.000 1.000 0.000
1 4.345 3.000 0.000
2 4.997 9.000 0.000
3 4.345 3.000 0.000
4 0.000 1.000 0.000
1/4波長表皮サンドイッチ型に対する反射・透過係数の垂直入射における周波数特性を求め，計算値を出力および図示すると，
zrn,zrp,ztn,ztp,zzzg = func_slab_tr(nt, ep, td, t, ssk, tth[nu_th], ic, ip)
print(f'Incident angle={tthdeg[nu_th]:.1f}[deg]')
dtn2,ptn2,dtp2,ptp2,drn2,prn2,drp2,prp2,dtg2,ptg2 = output_tr(ffreq,'f[GHz]',20,ztn,ztp,zrn,zrp,zzzg)
figure4b(tthdeg,ffreq,drn2,drp2,prn2,prp2,dtn2,dtp2,ptn2,ptp2,'tr_nt3_quater_skin_angle_freq')
Incident angle=0.0[deg]
Transmission coefficients Reflection coefficients
Perpendicular Parallel Perpendicular Parallel
f[GHz] Amp[dB] Ph[deg] Amp[dB] Ph[deg] Amp[dB] Ph[deg] Amp[dB] Ph[deg]
5.00 -1.580 -82.12 -1.580 -82.12 -5.158 105.70 -5.158 105.70 -98.20
6.00 -0.293 -105.51 -0.293 -105.51 -11.856 65.87 -11.856 65.87 -117.84
7.00 -0.013 -131.36 -0.013 -131.36 -25.307 -156.40 -25.307 -156.40 -137.48
8.00 -0.167 -153.66 -0.167 -153.66 -14.228 164.86 -14.228 164.86 -157.11
9.00 -0.102 -174.48 -0.102 -174.48 -16.342 127.60 -16.342 127.60 -176.75
10.00 -0.000 163.49 -0.000 163.49 -42.445 -90.86 -42.445 -90.86 163.61
11.00 -0.118 141.55 -0.118 141.55 -15.708 -129.25 -15.708 -129.25 143.97
12.00 -0.174 120.80 -0.174 120.80 -14.069 -166.42 -14.069 -166.42 124.33
13.00 -0.010 98.38 -0.010 98.38 -26.365 154.71 -26.365 154.71 104.69
14.00 -0.324 72.41 -0.324 72.41 -11.427 -67.69 -11.427 -67.69 85.05
15.00 -1.650 49.24 -1.650 49.24 -5.003 -107.30 -5.003 -107.30 65.41
16.00 -3.058 34.11 -3.058 34.11 -2.963 -138.86 -2.963 -138.86 45.77
17.00 -3.599 24.24 -3.599 24.24 -2.492 -165.17 -2.492 -165.17 26.13
18.00 -2.923 14.32 -2.923 14.32 -3.100 168.48 -3.100 168.48 6.49
19.00 -1.185 -2.67 -1.185 -2.67 -6.219 135.05 -6.219 135.05 -13.15
20.00 -0.001 -33.02 -0.001 -33.02 -36.426 -91.73 -36.426 -91.73 -32.79
21.00 -1.275 -62.74 -1.275 -62.74 -5.945 -137.89 -5.945 -137.89 -52.43
22.00 -2.956 -79.16 -2.956 -79.16 -3.066 -170.75 -3.066 -170.75 -72.07
23.00 -3.546 -89.03 -3.546 -89.03 -2.534 162.96 -2.534 162.96 -91.71
24.00 -2.933 -99.24 -2.933 -99.24 -3.089 136.31 -3.089 136.31 -111.34
25.00 -1.510 -115.02 -1.510 -115.02 -5.321 104.09 -5.321 104.09 -130.98
26.00 -0.263 -138.62 -0.263 -138.62 -12.304 64.05 -12.304 64.05 -150.62
27.00 -0.016 -164.33 -0.016 -164.33 -24.422 -158.09 -24.422 -158.09 -170.26
28.00 -0.160 173.48 -0.160 173.48 -14.408 163.29 -14.408 163.29 170.10
29.00 -0.087 152.58 -0.087 152.58 -17.044 125.95 -17.044 125.95 150.46
30.00 -0.002 130.47 -0.002 130.47 -32.907 -92.59 -32.907 -92.59 130.82
2項整合変成器
10 GHz を中心周波数とする2項整合変成器の格段のインピーダース特性（「伝送線路論」参照）より，
wl = CCC/freq
nt = 5# 誘電体スラブの層数
zL = 3.0# 終端のインピーダンス
znl = [1.0349, 1.2287, 1.7321, 2.4415, 2.8988]# 反射の最大値-30dBとする各段の特性インピーダンス
zn = np.array([1.0] + znl + [zL])
td = np.zeros(nt+2)# 誘電正接
ep = zn**2# 比誘電率
t = wl/zn/4# 1/4波長
t[0] = 0.0
t[-1] = 0.0
output_dielectric(nt,t,ep,td)
No. Thickness[mm] Dielectric constant Loss tangent
0 0.000 1.000 0.000
1 7.242 1.071 0.000
2 6.100 1.510 0.000
3 4.327 3.000 0.000
4 3.070 5.961 0.000
5 2.585 8.403 0.000
6 0.000 9.000 0.000
2項整合変成器に対する反射・透過係数の垂直入射における周波数特性を求め，計算値を出力および図示すると，
zrn,zrp,ztn,ztp,zzzg = func_slab_tr(nt, ep, td, t, ssk, tth[nu_th], ic, ip)
print(f'Incident angle={tthdeg[nu_th]:.1f}[deg]')
dtn2,ptn2,dtp2,ptp2,drn2,prn2,drp2,prp2,dtg2,ptg2 = output_tr(ffreq,'f[GHz]',20,ztn,ztp,zrn,zrp,zzzg)
figure4b(tthdeg,ffreq,drn2,drp2,prn2,prp2,dtn2,dtp2,ptn2,ptp2,'tr_nt2_binominal_angle_freq')
Incident angle=0.0[deg]
Transmission coefficients Reflection coefficients
Perpendicular Parallel Perpendicular Parallel
f[GHz] Amp[dB] Ph[deg] Amp[dB] Ph[deg] Amp[dB] Ph[deg] Amp[dB] Ph[deg]
5.00 -0.046 -82.12 -0.046 -82.12 -19.739 -42.15 -19.739 -42.15 -84.96
6.00 -0.008 -100.04 -0.008 -100.04 -27.600 -88.03 -27.600 -88.03 -101.95
7.00 -0.001 -117.68 -0.001 -117.68 -38.500 -133.49 -38.500 -133.49 -118.94
8.00 -0.000 -135.15 -0.000 -135.15 -54.112 -177.56 -54.112 -177.56 -135.93
9.00 -0.000 -152.55 -0.000 -152.55 -75.505 156.30 -75.505 156.30 -152.93
10.00 -0.000 -169.92 -0.000 -169.92 -84.328 180.00 -84.328 180.00 -169.92
11.00 -0.000 172.71 -0.000 172.71 -75.505 -156.30 -75.505 -156.30 173.09
12.00 -0.000 155.32 -0.000 155.32 -54.112 177.56 -54.112 177.56 156.10
13.00 -0.001 137.84 -0.001 137.84 -38.500 133.49 -38.500 133.49 139.11
14.00 -0.008 120.20 -0.008 120.20 -27.600 88.03 -27.600 88.03 122.12
15.00 -0.046 102.29 -0.046 102.29 -19.739 42.15 -19.739 42.15 105.12
16.00 -0.173 84.16 -0.173 84.16 -14.081 -3.97 -14.081 -3.97 88.13
17.00 -0.437 66.33 -0.437 66.33 -10.188 -49.81 -10.188 -49.81 71.14
18.00 -0.801 49.60 -0.801 49.60 -7.738 -94.55 -7.738 -94.55 54.15
19.00 -1.122 34.36 -1.122 34.36 -6.428 -137.79 -6.428 -137.79 37.16
20.00 -1.249 20.17 -1.249 20.17 -6.021 -180.00 -6.021 -180.00 20.17
21.00 -1.122 5.97 -1.122 5.97 -6.428 137.79 -6.428 137.79 3.17
22.00 -0.801 -9.27 -0.801 -9.27 -7.738 94.55 -7.738 94.55 -13.82
23.00 -0.437 -26.00 -0.437 -26.00 -10.188 49.81 -10.188 49.81 -30.81
24.00 -0.173 -43.83 -0.173 -43.83 -14.081 3.97 -14.081 3.97 -47.80
25.00 -0.046 -61.96 -0.046 -61.96 -19.739 -42.15 -19.739 -42.15 -64.79
26.00 -0.008 -79.87 -0.008 -79.87 -27.600 -88.03 -27.600 -88.03 -81.78
27.00 -0.001 -97.51 -0.001 -97.51 -38.500 -133.49 -38.500 -133.49 -98.78
28.00 -0.000 -114.98 -0.000 -114.98 -54.112 -177.56 -54.112 -177.56 -115.77
29.00 -0.000 -132.38 -0.000 -132.38 -75.505 156.30 -75.505 156.30 -132.76
30.00 -0.000 -149.75 -0.000 -149.75 -84.328 -180.00 -84.328 -180.00 -149.75
チェビシェフ変成器
10 GHz を中心周波数とし，反射係数が $-30$ dB等リプルとなるチェビシェフ変成器の格段のインピーダース特性（「伝送線路論」参照）より，
nt = 5# 誘電体スラブの層数
#zL = 2.0# 終端のインピーダンス
#znl = [1.0806, 1.2146, 1.4142, 1.6466, 1.8508]# 反射の最大値-30dBとする各段の特性インピーダンス
zL = 3.0# 終端のインピーダンス
znl = [1.1067, 1.3367, 1.7321, 2.2444, 2.7108]# 反射の最大値-30dBとする各段の特性インピーダンス
zn = np.array([1.0] + znl + [zL])
td = np.zeros(nt+2)# 誘電正接
ep = zn**2# 比誘電率
t = wl/zn/4# 1/4波長
t[0] = 0.0
t[-1] = 0.0
output_dielectric(nt,t,ep,td)
No. Thickness[mm] Dielectric constant Loss tangent
0 0.000 1.000 0.000
1 6.772 1.225 0.000
2 5.607 1.787 0.000
3 4.327 3.000 0.000
4 3.339 5.037 0.000
5 2.765 7.348 0.000
6 0.000 9.000 0.000
チェビシェフ変成器に対する反射・透過係数の垂直入射における周波数特性を求め，計算値を出力および図示すると，
zrn,zrp,ztn,ztp,zzzg = func_slab_tr(nt, ep, td, t, ssk, tth[nu_th], ic, ip)
print(f'Incident angle={tthdeg[nu_th]:.1f}[deg]')
dtn2,ptn2,dtp2,ptp2,drn2,prn2,drp2,prp2,dtg2,ptg2 = output_tr(ffreq,'f[GHz]',20,ztn,ztp,zrn,zrp,zzzg)
figure4b(tthdeg,ffreq,drn2,drp2,prn2,prp2,dtn2,dtp2,ptn2,ptp2,'tr_nt2_Chebyshev_angle_freq')
Incident angle=0.0[deg]
Transmission coefficients Reflection coefficients
Perpendicular Parallel Perpendicular Parallel
f[GHz] Amp[dB] Ph[deg] Amp[dB] Ph[deg] Amp[dB] Ph[deg] Amp[dB] Ph[deg]
5.00 -0.002 -86.23 -0.002 -86.23 -34.128 136.78 -34.128 136.78 -88.04
6.00 -0.003 -104.37 -0.003 -104.37 -31.319 91.27 -31.319 91.27 -105.65
7.00 -0.000 -122.39 -0.000 -122.39 -52.168 -133.99 -52.168 -133.99 -123.26
8.00 -0.004 -140.33 -0.004 -140.33 -30.879 -179.44 -30.879 -179.44 -140.87
9.00 -0.003 -158.20 -0.003 -158.20 -31.510 135.31 -31.510 135.31 -158.48
10.00 -0.000 -176.09 -0.000 -176.09 -97.928 180.00 -97.928 180.00 -176.09
11.00 -0.003 166.02 -0.003 166.02 -31.510 -135.31 -31.510 -135.31 166.30
12.00 -0.004 148.15 -0.004 148.15 -30.879 179.44 -30.879 179.44 148.70
13.00 -0.000 130.21 -0.000 130.21 -52.168 133.99 -52.168 133.99 131.09
14.00 -0.003 112.19 -0.003 112.19 -31.319 -91.27 -31.319 -91.27 113.48
15.00 -0.002 94.06 -0.002 94.06 -34.128 -136.78 -34.128 -136.78 95.87
16.00 -0.014 75.44 -0.014 75.44 -24.927 -2.83 -24.927 -2.83 78.26
17.00 -0.174 56.30 -0.174 56.30 -14.059 -49.36 -14.059 -49.36 60.65
18.00 -0.580 37.95 -0.580 37.95 -9.031 -95.10 -9.031 -95.10 43.04
19.00 -1.046 21.91 -1.046 21.91 -6.695 -138.53 -6.695 -138.53 25.44
20.00 -1.249 7.83 -1.249 7.83 -6.021 -180.00 -6.021 -180.00 7.83
21.00 -1.046 -6.25 -1.046 -6.25 -6.695 138.53 -6.695 138.53 -9.78
22.00 -0.580 -22.29 -0.580 -22.29 -9.031 95.10 -9.031 95.10 -27.39
23.00 -0.174 -40.65 -0.174 -40.65 -14.059 49.36 -14.059 49.36 -45.00
24.00 -0.014 -59.79 -0.014 -59.79 -24.927 2.83 -24.927 2.83 -62.61
25.00 -0.002 -78.41 -0.002 -78.41 -34.128 136.78 -34.128 136.78 -80.22
26.00 -0.003 -96.54 -0.003 -96.54 -31.319 91.27 -31.319 91.27 -97.83
27.00 -0.000 -114.56 -0.000 -114.56 -52.168 -133.99 -52.168 -133.99 -115.43
28.00 -0.004 -132.50 -0.004 -132.50 -30.879 -179.44 -30.879 -179.44 -133.04
29.00 -0.003 -150.37 -0.003 -150.37 -31.510 135.31 -31.510 135.31 -150.65
30.00 -0.000 -168.26 -0.000 -168.26 -97.928 -180.00 -97.928 -180.00 -168.26
前のページに戻る
「目次」のページに戻る