-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathfig1_23.py
119 lines (82 loc) · 2.2 KB
/
fig1_23.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
# Bottom loss curves
from numpy import *
from matplotlib.pyplot import *
from bdryr import *
print('Bottom loss curves...')
theta = linspace(90,0,181)
thetar = theta*pi/180.0
grazing_angle = 90 - theta
cw = 1500.0
cp = 1550.0
cs = 0.0
alphap = 0.5
alphas = 0.0
rhow = 1000.0
rhob = 2000.0
R = bdryr(rhow,rhob,cw,cp,cs,alphap,alphas,thetar)
B1 = -10*log10( abs(R)*abs(R) )
cp = 1600.0
R = bdryr(rhow,rhob,cw,cp,cs,alphap,alphas,thetar)
B2 = -10*log10( abs(R)*abs(R) )
cp = 1800.0
R = bdryr(rhow,rhob,cw,cp,cs,alphap,alphas,thetar)
B3 = -10*log10( abs(R)*abs(R) )
cp = 1600.0
alphap = 0.0
R = bdryr(rhow,rhob,cw,cp,cs,alphap,alphas,thetar)
B4 = -10*log10( abs(R)*abs(R) )
alphap = 0.5
R = bdryr(rhow,rhob,cw,cp,cs,alphap,alphas,thetar)
B5 = -10*log10( abs(R)*abs(R) )
alphap = 1.0
R = bdryr(rhow,rhob,cw,cp,cs,alphap,alphas,thetar)
B6 = -10*log10( abs(R)*abs(R) )
alphap = 0.5
rhob = 1500.0
R = bdryr(rhow,rhob,cw,cp,cs,alphap,alphas,thetar)
B7 = -10*log10( abs(R)*abs(R) )
rhob = 2000.0
R = bdryr(rhow,rhob,cw,cp,cs,alphap,alphas,thetar)
B8 = -10*log10( abs(R)*abs(R) )
rhob = 2500.0
R = bdryr(rhow,rhob,cw,cp,cs,alphap,alphas,thetar)
B9 = -10*log10( abs(R)*abs(R) )
cp = 1600.0
alphap = 0.0
alphas = 0.0
rhob = 2000.0
cs = 0.0
R = bdryr(rhow,rhob,cw,cp,cs,alphap,alphas,thetar)
B10 = -10*log10( abs(R)*abs(R) )
cs = 200.0
R = bdryr(rhow,rhob,cw,cp,cs,alphap,alphas,thetar)
B11 = -10*log10( abs(R)*abs(R) )
cs = 400.0
R = bdryr(rhow,rhob,cw,cp,cs,alphap,alphas,thetar)
B12 = -10*log10( abs(R)*abs(R) )
cs = 600.0
R = bdryr(rhow,rhob,cw,cp,cs,alphap,alphas,thetar)
B13 = -10*log10( abs(R)*abs(R) )
figure(1)
plot(grazing_angle,B1,grazing_angle,B2,'--',grazing_angle,B3,'.')
grid(True)
xlabel('Grazing angle (degrees)')
ylabel('Loss (dB)')
figure(2)
plot(grazing_angle,B4,grazing_angle,B5,'--',grazing_angle,B6,'.')
grid(True)
xlabel('Grazing angle (degrees)')
ylabel('Loss (dB)')
figure(3)
plot(grazing_angle,B7,grazing_angle,B8,'--',grazing_angle,B9,'.')
grid(True)
xlabel('Grazing angle (degrees)')
ylabel('Loss (dB)')
figure(4)
plot(grazing_angle,B10,grazing_angle,B11,'--',
grazing_angle,B12,'.',grazing_angle,B13,'r')
grid(True)
xlabel('Grazing angle (degrees)')
ylabel('Loss (dB)')
show()
print("done.")