FlintyLemming/PLFM_RADAR
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Add files via upload

Browse Source
This commit is contained in:
NawfalMotii79
2026-03-15 22:03:54 +00:00
committed by GitHub
parent 70e8366dc4
commit e0ff3b1bc5
6 changed files with 310 additions and 0 deletions
Download Patch File Download Diff File Expand all files Collapse all files
5_Simulations/Matlab/16_8.m
+52
View File
@@ -0,0 +1,52 @@
%Size of PCB
pcbThickness = 0.578e-3;
pcbLength = 200e-3;
pcbWidth = 200e-3;
f = 10.5e+9;
c = 3e+8;
lambda = c/f;
Length_patch = 8.76e-3;
Width_patch = 9.545e-3;
dist_patch = lambda/2;
feeder_width = 0.94e-3;
%Specifying Material of PCB
pcbEpsilonR = 3.48;
%Creating dielectic Material
d = dielectric(Name="RO4350B",EpsilonR=pcbEpsilonR,LossTangent=0.003,Thickness=pcbThickness);
GndPlane = antenna.Rectangle('Length',pcbLength,'Width',pcbWidth); %Creating Ground Plane of Antenna
%Creating Different Shapes of antenna
Rec1 = antenna.Rectangle('Length',dist_patch-Length_patch,'Width',feeder_width,'Center',[0,(dist_patch-Length_patch)/2]); %feeder
Rec2 = antenna.Rectangle('Length',Length_patch,'Width',Width_patch,'Center',[0,(dist_patch-Length_patch)+(Length_patch)/2]); %patch
for i=1:1:16
for j=1:1:8
Rec1 = Rec1 + antenna.Rectangle('Length',dist_patch-Length_patch,'Width',feeder_width,'Center',[i-1,j*(dist_patch)+(dist_patch-Length_patch)/2]);
Rec2 = Rec2 + antenna.Rectangle('Length',Length_patch,'Width',Width_patch,'Center',[i-1,j*(dist_patch)+(dist_patch-Length_patch)+(Length_patch)/2]);
end
end
%Adding all different shapes of antenna
AntennaPlane = Rec1 + Rec2;
%%Creating PCB Stack
p = pcbStack;
p.Name = 'Strip-fed slot';
p.BoardShape = GndPlane;
p.BoardThickness = pcbThickness;
p.Layers = {AntennaPlane,d,GndPlane}; %[x Cordinate,y Cordinate,startLayer stopLayer]
p.FeedLocations = [0,0,1,3;0,dist_patch,1,3;0,dist_patch*2,1,3;0,dist_patch*3,1,3;0,dist_patch*4,1,3;0,dist_patch*5,1,3;0,dist_patch*6,1,3;0,dist_patch*7,1,3;0,dist_patch*8,1,3;0,dist_patch*9,1,3;0,dist_patch*10,1,3;0,dist_patch*11,1,3;0,dist_patch*12,1,3;0,dist_patch*13,1,3;0,dist_patch*14,1,3;0,dist_patch*15,1,3];
%Plotting Different patterns and graphs
figure(1);
show(p); %Display Antenna
% figure(2);
% pattern(p,1.943e9); %Display Radiation Pattern at 1.943GHZ
% figure(3);
% impedance(p,1.6e9:2e7:2.2e9); %Display Impedance Graph from 1.6GHz to 2.2GHz
% freq = linspace(1.6e9, 2.2e9, 50); % Creating Frequency Vector
% s = sparameters(p,freq,50); % Calalculate S11 for all frequencys
% figure(4);
% rfplot(s);%Diplay S11 Plot
5_Simulations/Matlab/Antenna16_8.m
+108
View File
@@ -0,0 +1,108 @@
n_array = 1;
m_array = 2;
delta_phy = 15; %phase difference between adjacent patches (in degrees)
f = 10.5e+9;
c = 3e+8;
lambda = c/f;
Width_patch = 8.76e-3;
Length_patch = 9.545e-3;
dist_patch = lambda/2;
feeder_width = dist_patch - Width_patch;
feeder_length = 0.94e-3;
% Define angles (break into smaller chunks)
theta = -90:1:90;
phi = -180:1:180;
chunk_size = 30;
% Define file name for saving
matFileName = 'pattern_data.mat';
% Create a MAT file object for incremental writing
mf = matfile(matFileName, 'Writable', true);
%Size of PCB
pcbThickness = 0.508e-3;
pcbWidth = 0.5*feeder_width + (feeder_width +Width_patch)*(m_array);
pcbLength =dist_patch*n_array;
%Specifying Material of PCB
% pcbMaterial = 'RO4350B';
% pcbEpsilonR = 3.48;
pcbMaterial = 'RO4350B';
pcbEpsilonR = 3.48;
d = dielectric(pcbMaterial);
d.EpsilonR = pcbEpsilonR;
d.Thickness = pcbThickness;
%Creating dielectic Material
GndPlane = antenna.Rectangle('Length',pcbLength,'Width',pcbWidth,'Center',[0,0]); %Creating Ground Plane of Antenna
%Creating Different Shapes of antenna
Rec1 = antenna.Rectangle('Length',feeder_length,'Width',feeder_width,'Center',[-(pcbLength/2)+0.5*dist_patch,-pcbWidth/2+feeder_width/2]);
Rec2 = antenna.Rectangle('Length',Length_patch,'Width',Width_patch,'Center',[-(pcbLength/2)+0.5*dist_patch,-pcbWidth/2 + feeder_width+Width_patch/2]);
for j=1:1:n_array
for i=1:1:m_array
if((i~=1) || (j~=1))
Rec1 = Rec1 + antenna.Rectangle('Length',feeder_length,'Width',feeder_width,'Center',[(j-1)*dist_patch-(pcbLength/2)+0.5*dist_patch,-pcbWidth/2+feeder_width/2+(i-1)*dist_patch]);
Rec2 = Rec2 + antenna.Rectangle('Length',Length_patch,'Width',Width_patch,'Center',[(j-1)*dist_patch-(pcbLength/2)+0.5*dist_patch,-pcbWidth/2+feeder_width+Width_patch/2+ (i-1)*dist_patch]);
end
end
end
AntennaPlane = Rec1 + Rec2;
%%Creating PCB Stack
p = pcbStack;
p.Name = 'Strip-fed slot';
p.BoardShape = GndPlane;
p.BoardThickness = pcbThickness;
p.Layers = {AntennaPlane,d,GndPlane};
%[x Cordinate,y Cordinate,startLayer stopLayer]
p.FeedDiameter = feeder_length/4;
p.FeedLocations = [-(pcbLength/2)+0.5*dist_patch,-pcbWidth/2+feeder_length/1.4,1,3];
%p.FeedLocations = [-(pcbLength/2)+0.5*dist_patch,-pcbWidth/2+feeder_length/1.4,1,3;dist_patch-(pcbLength/2)+0.5*dist_patch,-pcbWidth/2+feeder_length/1.4,1,3;dist_patch*2-(pcbLength/2)+0.5*dist_patch,-pcbWidth/2+feeder_length/1.4,1,3];
%p.FeedLocations = [0,feeder_length/1.4,1,3;dist_patch,feeder_length/1.4,1,3;dist_patch*2,feeder_length/1.4,1,3;dist_patch*3,feeder_length/1.4,1,3;dist_patch*4,feeder_length/1.4,1,3;dist_patch*5,feeder_length/1.4,1,3;dist_patch*6,feeder_length/1.4,1,3;dist_patch*7,feeder_length/1.4,1,3;dist_patch*8,feeder_length/1.4,1,3;dist_patch*9,feeder_length/1.4,1,3;dist_patch*10,feeder_length/1.4,1,3;dist_patch*11,feeder_length/1.4,1,3;dist_patch*12,feeder_length/1.4,1,3;dist_patch*13,feeder_length/1.4,1,3;dist_patch*14,feeder_length/1.4,1,3;dist_patch*15,feeder_length/1.4,1,3];
%p.p.FeedPhase = [0 delta_phy delta_phy*2 delta_phy*3 delta_phy*4 delta_phy*5 delta_phy*6 delta_phy*7 delta_phy*8 delta_phy*9 delta_phy*10 delta_phy*11 delta_phy*12 delta_phy*13 delta_phy*14 delta_phy*15];
figure(1);
show(p); %Display Antenna
%figure(2);
% pattern(p,f,-180:30:180,0);
% pattern3Dfig = figure(2); %array pattern
% pattern(p,f);
% Initialize counters for indexing
index = 1;
% Process and store in chunks
for i = 1:chunk_size:length(theta)
for j = 1:chunk_size:length(phi)
% Simulated computation (replace with actual pattern function)
p_val = pattern(p,f,theta(i),phi(j));
mf.pattern_values(index, 1) = p_val;
% Store computed values directly into the MAT file (avoiding large RAM usage)
mf.theta_values(index, 1) = theta(i);
mf.phi_values(index, 1) = phi(j);
clear p_val theta(i) phi(j);
% Increment index
index = index + 1;
end
end
figure(2);
patternCustom( mf.pattern_values,mf.theta_values ,mf.phi_values);
% Display confirmation message
%M = readtable("pattern_data.mat");
%patternCustom(M(:,3),M(:,2),M(:,1));
5_Simulations/Matlab/Antenna_array.m
+40
View File
@@ -0,0 +1,40 @@
%This Code was Written By Nemeen Shah
%Subscribe to My YouTube Channel: https://Youtube.com/NematicsLab
%Size of PCB
pcbThickness = 1.6e-3; %1.6mm
pcbLength = 152.4e-3; %152.4mm or 6inch
pcbWidth = 101.6e-3; %101.6mm 0r 4inch
%Specifying Material of PCB
pcbMaterial = 'FR4';
pcbEpsilonR = 4.4;
%Creating dielectic Material
d = dielectric(pcbMaterial);
d.EpsilonR = pcbEpsilonR;
d.Thickness = pcbThickness;
AntennaPlane=antenna.Rectangle('Length',0.5e-2,'Width',5e-2,'Center',[0, 0]); %Creating Feed Plane of Antenna
GndPlane = antenna.Rectangle('Length',pcbLength,'Width',pcbWidth); %Creating Ground Plane of Antenna
%Creating Different Shapes of antenna
Rec = antenna.Rectangle('Length',10e-2,'Width',2e-2,'Center',[0,-20e-3]);
Rec1 = antenna.Rectangle('Length',6e-2,'Width',2e-2,'Center',[0,20e-3]);
%%Creating PCB Stack
p = pcbStack;
p.Name = 'Strip-fed slot';
p.BoardShape = GndPlane;
p.BoardThickness = pcbThickness;
p.Layers = {AntennaPlane,d,GndPlane}; %[x Cordinate,y Cordinate,startLayer stopLayer]
p.FeedLocations = [0,(-pcbWidth/2)+6e-2,1,3];
%Adding all different shapes of antenna
AntennaPlane = AntennaPlane + Rec + Rec1;
p.Layers = {AntennaPlane,d,GndPlane};
%Plotting Different patterns and graphs
figure(1);
show(p); %Display Antenna
figure(2);
pattern(p,1.943e9); %Display Radiation Pattern at 1.943GHZ
figure(3);
impedance(p,1.6e9:2e7:2.2e9); %Display Impedance Graph from 1.6GHz to 2.2GHz
freq = linspace(1.6e9, 2.2e9, 50); % Creating Frequency Vector
s = sparameters(p,freq,50); % Calalculate S11 for all frequencys
figure(4);
rfplot(s);%Diplay S11 Plot
5_Simulations/Matlab/DielectricCatalog.csv
+18
View File
@@ -0,0 +1,18 @@
Name,Relative_Permittivity,Loss_Tangent,Frequency,Comments
"RO4350B",4.38,0.003,12000000000,""
"Air",1,0,1000000000,""
"FR4",4.8,0.026,100000000,""
"Teflon",2.1,0.0002,100000000,""
"Foam",1.03,0.00015,50000000,""
"Polystyrene",2.55,0.0001,100000000,""
"Plexiglas",2.59,0.00676,10000000000,""
"Fused quartz",3.78,0.0001,10000000000,""
"E glass",6.22,0.0023,100000000,""
"RO4725JXR",2.55,0.0022,2500000000,""
"RO4730JXR",3,0.0023,2500000000,""
"TMM3",3.45,0.002,10000000000,""
"TMM4",4.7,0.002,10000000000,""
"TMM6",6.3,0.0023,10000000000,""
"TMM10",9.8,0.0022,10000000000,""
"TMM10i",9.9,0.002,10000000000,""
"Taconic RF-35",3.5,0.0018,1900000000,""
1 Name Relative_Permittivity Loss_Tangent Frequency Comments
2 RO4350B 4.38 0.003 12000000000
3 Air 1 0 1000000000
4 FR4 4.8 0.026 100000000
5 Teflon 2.1 0.0002 100000000
6 Foam 1.03 0.00015 50000000
7 Polystyrene 2.55 0.0001 100000000
8 Plexiglas 2.59 0.00676 10000000000
9 Fused quartz 3.78 0.0001 10000000000
10 E glass 6.22 0.0023 100000000
11 RO4725JXR 2.55 0.0022 2500000000
12 RO4730JXR 3 0.0023 2500000000
13 TMM3 3.45 0.002 10000000000
14 TMM4 4.7 0.002 10000000000
15 TMM6 6.3 0.0023 10000000000
16 TMM10 9.8 0.0022 10000000000
17 TMM10i 9.9 0.002 10000000000
18 Taconic RF-35 3.5 0.0018 1900000000
5_Simulations/Matlab/antenna16_antenna8.m
+92
View File
@@ -0,0 +1,92 @@
% Open file to write
csvFileName = 'pattern_data.csv';
fileID = fopen(csvFileName, 'w');
fprintf(fileID, 'Theta,Phi,Pattern\n'); % CSV headers
% Define angles (break into smaller chunks)
theta = -90:1:90;
phi = -180:1:180;
chunk_size = 50;
%Size of PCB
pcbThickness = 0.508e-3;
pcbLength = 100e-3;
pcbWidth = 100e-3;
n_array = 3;
m_array = 2;
f = 10.5e+9;
c = 3e+8;
lambda = c/f;
Length_patch = 8.76e-3;
Width_patch = 9.545e-3;
dist_patch = lambda/2;
feeder_width = 0.94e-3;
delta_phy = 15; %in degrees
%Specifying Material of PCB
pcbMaterial = 'RO4350B';
pcbEpsilonR = 3.48;
%Creating dielectic Material
d = dielectric(pcbMaterial);
d.EpsilonR = pcbEpsilonR;
d.Thickness = pcbThickness;
GndPlane = antenna.Rectangle('Length',pcbLength,'Width',pcbWidth); %Creating Ground Plane of Antenna
%Creating Different Shapes of antenna
%Rec1 = antenna.Rectangle('Length',feeder_width,'Width',dist_patch-Length_patch,'Center',[0,(dist_patch-Length_patch)/2]); %feeder
%Rec2 = antenna.Rectangle('Length',Length_patch,'Width',Width_patch,'Center',[0,(dist_patch-Length_patch)+(Length_patch)/2]); %patch
for i=1:1:n_array
for j=1:1:m_array
Rec1 = Rec1 + antenna.Rectangle('Length',feeder_width,'Width',dist_patch-Length_patch,'Center',[(i-1)*(dist_patch),(j-1)*(dist_patch)+(dist_patch-Length_patch)/2]);
Rec2 = Rec2 + antenna.Rectangle('Length',Length_patch,'Width',Width_patch,'Center',[(i-1)*(dist_patch),(j-1)*(dist_patch)+(dist_patch-Length_patch)+(Length_patch)/2]);
end
end
%Adding all different shapes of antenna
AntennaPlane = Rec1 + Rec2;
%%Creating PCB Stack
p = pcbStack;
p.Name = 'patch slot';
p.BoardShape = GndPlane;
p.BoardThickness = pcbThickness;
p.Layers = {AntennaPlane,d,GndPlane};
%[x Cordinate,y Cordinate,startLayer stopLayer]
% p.FeedLocations = [0,0,1,3;dist_patch,0,1,3;dist_patch*2,0,1,3;dist_patch*3,0,1,3;dist_patch*4,0,1,3;dist_patch*5,0,1,3;dist_patch*6,0,1,3;dist_patch*7,0,1,3;dist_patch*8,0,1,3;dist_patch*9,0,1,3;dist_patch*10,0,1,3;dist_patch*11,0,1,3;dist_patch*12,0,1,3;dist_patch*13,0,1,3;dist_patch*14,0,1,3;dist_patch*15,0,1,3];
% p.FeedPhase = [0,delta_phy, delta_phy*2, delta_phy*3, delta_phy*4, delta_phy*5, delta_phy*6, delta_phy*7, delta_phy*8, delta_phy*9, delta_phy*10, delta_phy*11, delta_phy*12, delta_phy*13, delta_phy*14, delta_phy*15];
p.FeedLocations = [0,0.001,1,3;dist_patch,0.001,1,3;dist_patch*2,0.001,1,3];
%Plotting Different patterns and graphs
figure(1);
show(p); %Display Antenna
%pattern(p,10.5e9); %Display Radiation Pattern at 10.5GHz
% Process and save in chunks
for i = 1:chunk_size:length(theta)
for j = 1:chunk_size:length(phi)
% Simulated computation (replace with actual pattern function)
p_val = rand(); % Example placeholder
% Write to file
fprintf(fileID, '%.2f,%.2f,%.6f\n', theta(i), phi(j), p_val);
end
end
% Close the file
fclose(fileID);
% Display confirmation message
figure(2);
disp('Pattern data saved to CSV.');
% figure(3);
% impedance(p,0.8*f:0.05*f:1.2*f); %Display Impedance Graph from 1.6GHz to 2.2GHz
% freq = linspace(0.8*f, 1.2*f, 50); % Creating Frequency Vector
% s = sparameters(p,freq,50); % Calalculate S11 for all frequencys
% figure(4);
% rfplot(s);%Diplay S11 Plot
5_Simulations/Matlab/pattern_data.mat
BIN
View File
Binary file not shown.