silas/imdd_silas
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

setupt_simulation changed

Browse Source
This commit is contained in:
Silas Oettinghaus
2023-10-09 15:26:03 +02:00
parent cc4a66adae
commit 9ca9db0722
1 changed files with 48 additions and 46 deletions
Download Patch File Download Diff File Expand all files Collapse all files

94
setup_simulation.m
View File

@@ -9,11 +9,11 @@ for lp = 1
%% A) Set Simulation Variables
sir = -18;
delay = 10; %mpi delay in meter
delay = 50; %mpi delay in meter
fiblen = 0; %main link in km
laser_linewidth =5e6;
laser_linewidth =1e6;
O = 18; %order of prbs
N = 2^(O-1); %length of prbs
@@ -72,8 +72,8 @@ for lp = 1
eq = EQ("K",2,"plottrain",0,"plotfinal",0,...
"training_length",4096,"training_loops",2,...
"Ne",[50,5,3],"Nb",[3,2,2],...
"DCmu",0.05,"DDmu",[0.0004 0.0005 0.0006 0.0007 ],"DFEmu",0.005,"FFEmu",0.005,...
"Ne",[50,5,3],"Nb",[3,0,0],...
"DCmu",0.05,"DDmu",[0.0004 0.0005 0.0006 0.0007 ],"DFEmu",0.005,"FFEmu",0.00,...
"dd_loops",2,"epsilon",[10 100 1000 ],"M",2,...
"thres",[0.005 0.004 0.0005 ],"l1act",0,"delay",0,"rho",0.0005,"ideal_dfe",0,"DB_aim",0);
@@ -128,51 +128,53 @@ for lp = 1
% Delay the reflected signal
[R,n] = R.delay("delay_meter",delay);
% % Add together
% col = cbrewer2('qual','Paired',8);
%
% xax = (1:X.length);% / fsimu * (physconst("LightSpeed")/1.4677);
% xax_mtr = (1:X.length) / fsimu * (physconst("LightSpeed")/1.4677);
% thresh = 0.7;
% phaseX = phase(X.signal);
% phaseR = phase(R.signal);
% phasediff = wrapToPi(phaseX-phaseR);
% [pos_high] = find(abs(phasediff)>thresh);
% [pos_low] = find(abs(phasediff)<=thresh);
%
% figure()
% subplot(3,1,1)
% scatter(xax_mtr,wrapToPi(phaseX),4,'.','MarkerEdgeColor',col(4,:),'DisplayName','Phase of original Signal')
% hold on
% scatter(xax_mtr,wrapToPi(phaseR),4,'.','MarkerEdgeColor',col(7,:),'DisplayName','Phase of delayed Signal')
% hold on
% xline(delay,'LineWidth',4,'LineStyle','--','HandleVisibility','off')
% lg = legend;
% lg.Location = "southwest";
%
% subplot(3,1,2)
% scatter(xax_mtr(pos_low),phasediff(pos_low),4,'.','MarkerEdgeColor',col(6,:),'DisplayName','Phase Difference')
% hold on
% scatter(xax_mtr(pos_high),phasediff(pos_high),4,'.','MarkerEdgeColor',col(4,:),'HandleVisibility','off')
% xline(delay,'LineWidth',4,'LineStyle','--','HandleVisibility','off')
% yline(thresh,'LineWidth',2,'LineStyle','-','HandleVisibility','off')
% yline(-thresh,'LineWidth',2,'LineStyle','-','HandleVisibility','off')
% lg = legend;
% lg.Location = "southwest";
% Add together
col = cbrewer2('qual','Paired',8);
X = X+R;
xax = (1:X.length);% / fsimu * (physconst("LightSpeed")/1.4677);
xax_sec = (1:X.length) / fsimu .* 1e6;%* (physconst("LightSpeed")/1.4677);
xax_sec = xax_sec(n:end).';
thresh = pi/4;
phaseX = phase(X.signal(n:end));
phaseR = phase(R.signal(n:end));
phasediff = wrapToPi(phaseX-phaseR);
[pos_high] = find(abs(phasediff)>thresh);
[pos_low] = find(abs(phasediff)<=thresh);
% subplot(3,1,3)
% scatter(xax_mtr(pos_low),abs(X.signal(pos_low).^2),4,'.','MarkerEdgeColor',col(6,:),'DisplayName','Constructive Interference');
% hold on
% scatter(xax_mtr(pos_high),abs(X.signal(pos_high).^2),4,'.','MarkerEdgeColor',col(4,:),'DisplayName','Destructive Interference');
% xline(delay,'LineWidth',4,'LineStyle','--','HandleVisibility','off')
%
% lg = legend;
% lg.Location = "southwest";
figure()
subplot(3,1,1)
scatter(xax_sec.',wrapToPi(phaseX),4,'.','MarkerEdgeColor',col(2,:),'DisplayName','Phase of original Signal')
hold on
scatter(xax_sec,wrapToPi(phaseR),4,'.','MarkerEdgeColor',col(1,:),'DisplayName','Phase of delayed Signal')
hold on
% xline(delay/(physconst("LightSpeed")/1.4677),'LineWidth',4,'LineStyle','--','HandleVisibility','off')
lg = legend;
lg.Location = "southwest";
subplot(3,1,2)
scatter(xax_sec(pos_low),phasediff(pos_low),4,'.','MarkerEdgeColor',col(5,:),'DisplayName','Phase Difference')
hold on
scatter(xax_sec(pos_high),phasediff(pos_high),4,'.','MarkerEdgeColor',col(2,:),'HandleVisibility','off')
% xline(delay/(physconst("LightSpeed")/1.4677),'LineWidth',4,'LineStyle','--','HandleVisibility','off')
yline(thresh,'LineWidth',2,'LineStyle','-','HandleVisibility','off')
yline(-thresh,'LineWidth',2,'LineStyle','-','HandleVisibility','off')
lg = legend;
lg.Location = "southwest";
X.signal = X.signal(n:end);
R.signal = R.signal(n:end);
X = X+R;
subplot(3,1,3)
scatter(xax_sec(pos_low),abs(X.signal(pos_low).^2),4,'.','MarkerEdgeColor',col(5,:),'DisplayName','Constructive Interference');
hold on
scatter(xax_sec(pos_high),abs(X.signal(pos_high).^2),4,'.','MarkerEdgeColor',col(2,:),'DisplayName','Destructive Interference');
% xline(delay/(physconst("LightSpeed")/1.4677),'LineWidth',4,'LineStyle','--','HandleVisibility','off')
lg = legend;
lg.Location = "southwest";
disp(['SIR ',num2str(10*log10(X.power/R.power))]);
@@ -214,7 +216,7 @@ for lp = 1
Eq_in.signal = Eq_in.signal - 1/wl .* movsum( Eq_in.signal,[wl/2,wl/2]);
% Equalize Signal
[Eq_out] = eq2.process(Eq_in,digimod_out);
[Eq_out] = eq.process(Eq_in,digimod_out);
@@ -379,5 +381,5 @@ for lp = 1
end
end