% fcc_ieeec951_levels_5g_bands_rev_3.m
% USA
E(1) = 61.4
Z = 120*pi
P(1) = E(1)^2/Z
% CANADA
E(2) = 40.7
P(2) = E(2)^2/Z
% Russia, China, Eastern Europe, Switzerland
E(3) = 6.14
P(3) = E(3)^2/Z
% BIO Initiative Report 2007 Salzburg June 2000
E(4) = 0.614
P(4) = E(4)^2/Z
% European Parliament May 27, 2011
E(5) = 0.2
P(5) = E(5)^2/Z
% Average Indoor Urban Exposure, Toronto, CANADA
Eindoor = [.3 .4 .5 .6 .7 .8 .9 1.0 1.1 1.2 1.3 1.4]
Pindoor = Eindoor.^2/Z

freq_mhz = [0.3
0.4
0.5
0.6
0.7
0.8
0.9
1
1.1
1.2
1.3
1.4
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
30
40
50
60
70
80
90
100
110
120
130
140
150
160
170
180
190
200
210
220
230
240
250
260
270
280
290
300
300
400
500
600
700
800
900
1000
1100
1200
1300
1400
1500
2000
3000
4000
5000
6000
7000
8000
9000
10000
11000
12000
13000
14000
15000
16000
17000
18000
19000
20000
21000
22000
23000
24000
25000
26000
27000
28000
29000
30000
31000
32000
33000
34000
35000
36000
37000
38000
39000
40000
41000
42000
43000
44000
45000
46000
47000
48000
49000
50000
51000
52000
53000
54000
55000
56000
57000
58000
59000
60000
61000
62000
63000
64000
65000
66000
67000
68000
69000
70000
71000
72000
73000
74000
75000
76000
77000
78000
79000
80000
81000
82000
83000
84000
85000
86000
87000
88000
89000
90000
91000
92000
93000
94000
95000
96000
97000
98000
99000
100000];

S1 = [100
100
100
100
100
100
100
100
100
100
100
100
180/freq_mhz(13)^2
180/freq_mhz(14)^2
180/freq_mhz(15)^2
180/freq_mhz(16)^2
180/freq_mhz(17)^2
180/freq_mhz(18)^2
180/freq_mhz(19)^2
180/freq_mhz(20)^2
180/freq_mhz(21)^2
180/freq_mhz(22)^2
180/freq_mhz(23)^2
180/freq_mhz(24)^2
180/freq_mhz(25)^2
180/freq_mhz(26)^2
180/freq_mhz(27)^2
180/freq_mhz(28)^2
180/freq_mhz(29)^2
180/freq_mhz(30)^2
180/freq_mhz(31)^2
180/freq_mhz(32)^2
180/freq_mhz(33)^2
180/freq_mhz(34)^2
180/freq_mhz(35)^2
180/freq_mhz(36)^2
180/freq_mhz(37)^2
180/freq_mhz(38)^2
180/freq_mhz(39)^2
180/freq_mhz(40)^2
180/freq_mhz(41)^2
.2
.2
.2
.2
.2
.2
.2
.2
.2
.2
.2
.2
.2
.2
.2
.2
.2
.2
.2
.2
.2
.2
.2
.2
.2
.2
.2
.2
freq_mhz(70)/1500
freq_mhz(71)/1500
freq_mhz(72)/1500
freq_mhz(73)/1500
freq_mhz(74)/1500
freq_mhz(75)/1500
freq_mhz(76)/1500
freq_mhz(77)/1500
freq_mhz(78)/1500
freq_mhz(79)/1500
freq_mhz(80)/1500
freq_mhz(81)/1500
freq_mhz(82)/1500
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1];
% S is in unite of mW/cm^2

S2 = S1*E(2)^2/E(1)^2
S3 = S1*E(3)^2/E(1)^2
S4 = S1*E(4)^2/E(1)^2
S5 = S1*E(5)^2/E(1)^2
Sindoor = S1*Eindoor.^2/E(1)^2

c = 299792458
freq_1 = 700E6
Pt_max_d = -5.7% Transmit Power, dBm
Pt_max = 10^((Pt_max_d - 30)/10)% Transmit Power, Watts
% Variable Antenna Gains
%Gain_d = [1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16]
Gain_d = 3
Gain = 10.^(Gain_d/10)
d = 1% Minimum Distance from Transmitter
Pwr_Den_1 = .1*Pt_max*Gain./(4*pi*d.^2)
Pwr_Den_1_d = 10*log10(Pwr_Den_1)

freq_2 = 900E6
Pt_max_d = -4.7% Transmit Power, dBm
Pt_max = 10^((Pt_max_d - 30)/10)% Transmit Power, Watts
Pwr_Den_2 = .1*Pt_max*Gain./(4*pi*d.^2)
Pwr_Den_2_d = 10*log10(Pwr_Den_2)

freq_3 = 1800E6
Pt_max_d = -3% Transmit Power, dBm
Pt_max = 10^((Pt_max_d - 30)/10)% Transmit Power, Watts
Pwr_Den_3 = .1*Pt_max*Gain./(4*pi*d.^2)
Pwr_Den_3_d = 10*log10(Pwr_Den_3)

freq_4 = 2100E6
Pt_max_d = -3% Transmit Power, dBm
Pt_max = 10^((Pt_max_d - 30)/10)% Transmit Power, Watts
Pwr_Den_4 = .1*Pt_max*Gain./(4*pi*d.^2)
Pwr_Den_4_d = 10*log10(Pwr_Den_4)

freq_5 = 2600E6
Pt_max_d = -3% Transmit Power, dBm
Pt_max = 10^((Pt_max_d - 30)/10)% Transmit Power, Watts
Pwr_Den_5 = .1*Pt_max*Gain./(4*pi*d.^2)
Pwr_Den_5_d = 10*log10(Pwr_Den_5)

figure(1)
loglog(freq_mhz,S1,'-r','LineWidth',2.0)
hold on;
loglog(freq_mhz,S2,'-g','LineWidth',2.0)
hold on;
loglog(freq_mhz,S3,'-b','LineWidth',2.0)
hold on;
loglog(freq_mhz,S4,'-k','LineWidth',2.0)
hold on;
loglog(freq_mhz,S5,'-c','LineWidth',2.0)
hold on;
loglog(freq_mhz,Sindoor,'-m','LineWidth',2.0)
hold on;
loglog(freq_1/1E6,Pwr_Den_1,'--s','LineWidth',2.0)
hold off;
xlabel('Frequency (MHz)')
ylabel('Power Density (mW/cm^2)')
title('GP Mult Stds,700 MHz:1m distance, Gain 3 dBi, -5.7 dBm, Power Density mW/cm^2')

figure(2)
loglog(freq_mhz,S1,'-r','LineWidth',2.0)
hold on;
loglog(freq_mhz,S2,'-g','LineWidth',2.0)
hold on;
loglog(freq_mhz,S3,'-b','LineWidth',2.0)
hold on;
loglog(freq_mhz,S4,'-k','LineWidth',2.0)
hold on;
loglog(freq_mhz,S5,'-c','LineWidth',2.0)
hold on;
loglog(freq_mhz,Sindoor,'-m','LineWidth',2.0)
hold on;
loglog(freq_2/1E6,Pwr_Den_2,'--s','LineWidth',2.0)
hold off;
xlabel('Frequency (MHz)')
ylabel('Power Density (mW/cm^2)')
title('GP Mult Stds,900 MHz:1m distance, Gain 3 dBi, -4.7 dBm, Power Density mW/cm^2')

figure(3)
loglog(freq_mhz,S1,'-r','LineWidth',2.0)
hold on;
loglog(freq_mhz,S2,'-g','LineWidth',2.0)
hold on;
loglog(freq_mhz,S3,'-b','LineWidth',2.0)
hold on;
loglog(freq_mhz,S4,'-k','LineWidth',2.0)
hold on;
loglog(freq_mhz,S5,'-c','LineWidth',2.0)
hold on;
loglog(freq_mhz,Sindoor,'-m','LineWidth',2.0)
hold on;
loglog(freq_3/1E6,Pwr_Den_3,'--s','LineWidth',2.0)
hold off;
xlabel('Frequency (MHz)')
ylabel('Power Density (mW/cm^2)')
title('GP Mult Stds, 1800 MHz:1m distance, Gain 3 dBi, -3 dBm, Power Density mW/cm^2')

figure(4)
loglog(freq_mhz,S1,'-r','LineWidth',2.0)
hold on;
loglog(freq_mhz,S2,'-g','LineWidth',2.0)
hold on;
loglog(freq_mhz,S3,'-b','LineWidth',2.0)
hold on;
loglog(freq_mhz,S4,'-k','LineWidth',2.0)
hold on;
loglog(freq_mhz,S5,'-c','LineWidth',2.0)
hold on;
loglog(freq_mhz,Sindoor,'-m','LineWidth',2.0)
hold on;
loglog(freq_4/1E6,Pwr_Den_4,'--s','LineWidth',2.0)
hold off;
xlabel('Frequency (MHz)')
ylabel('Power Density (mW/cm^2)')
title('GP Mult Stds, 2100 MHz:1m distance, Gain 3 dBi, -3 dBm, Power Density mW/cm^2')

figure(5)
loglog(freq_mhz,S1,'-r','LineWidth',2.0)
hold on;
loglog(freq_mhz,S2,'-g','LineWidth',2.0)
hold on;
loglog(freq_mhz,S3,'-b','LineWidth',2.0)
hold on;
loglog(freq_mhz,S4,'-k','LineWidth',2.0)
hold on;
loglog(freq_mhz,S5,'-c','LineWidth',2.0)
hold on;
loglog(freq_mhz,Sindoor,'-m','LineWidth',2.0)
hold on;
loglog(freq_5/1E6,Pwr_Den_5,'--s','LineWidth',2.0)
hold off;
xlabel('Frequency (MHz)')
ylabel('Power Density (mW/cm^2)')
title('GP Mult Stds, 2600 MHz:1m distance, Gain 3 dBi, -3 dBm, Power Density mW/cm^2')