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

Adapt AWG to work with M8199B

Browse Source 
Minor changes in Scope
Measurement Codes for OFC 2025
This commit is contained in:
Silas Labor Zizou
2024-10-20 17:23:35 +02:00
parent bf94e3dc2f
commit 1b554f2d25
8 changed files with 315 additions and 136 deletions
Download Patch File Download Diff File Expand all files Collapse all files

164
projects/Lab_2024/lab_sir_sweep.m
View File

@@ -1,13 +1,13 @@
folderpath = 'C:\Users\sioe\Nextcloud\Dokumente\02_Ablage_Office\Lab_Data_24\sir_sweep_pam4\';
experiment_name = 'PAM4_DB_encoded_10km_';
folderpath = 'C:\Users\sioe\Nextcloud\Dokumente\02_Ablage_Office\Lab_Data_24\sir_sweep_sd40\';
experiment_name = 'PAM4_56_v2';
%%% SIR Sweep for MPI Experiment %%%
params = struct;
params.vbias = [2.45];
params.awg_vpp = [0.25];
params.eq_mode = [4];
params.vbias = [2.5];
params.awg_vpp = [0.35];
params.eq_mode = [2];
params.i_atten = [0:4:40];
wh = DataStorage(params);
@@ -23,13 +23,13 @@ wh.addStorage("signals");
precomp_path = "C:\Users\sioe\Documents\MATLAB\imdd_simulation\projects\standard_system\";
precomp_fn = "lab_mpi_setup_2";
precomp_mode = 2; %0=do nothing ; 1= measure; 2=precomp active
precomp_amp_max = 5;
precomp_amp_max = 2;
M = 4;
pn_key = 2;
usemrds = 0;
fsym = 92e9;
fdac = 92e9;
fsym = 68e9;
fdac = 256e9;
awg_vpp = 0.35;
fadc = 160e9;
rrcalpha = 0.05;
@@ -60,21 +60,6 @@ for eq_mode = wh.parameter.eq_mode.values
sprintf('Progress: %d/%d\nEstimated time remaining: %.2f hours\nEstimated time remaining: %.2f hours', ...
loopcnt, looptotal, estimatedTimeRemaining/60/60, estimatedTotalTime/60/60));
try
z
ber = wh2.getStoValue('ber',v_bias,awg_vpp,eq_mode,precomp_amp_max);
assert(~isempty(ber),'err')
rop = wh2.getStoValue('rop',v_bias,awg_vpp,eq_mode,precomp_amp_max);
assert(~isempty(rop))
pd_in = wh2.getStoValue('pd_in',v_bias,awg_vpp,eq_mode,precomp_amp_max);
assert(~isempty(pd_in))
M = wh2.getStoValue('m',v_bias,awg_vpp,eq_mode,precomp_amp_max);
assert(~isempty(M))
catch
switch eq_mode
@@ -93,6 +78,7 @@ for eq_mode = wh.parameter.eq_mode.values
db_coding_approach = 0;
db_precode = db_coding_approach || db_channel_approach;
case 3
ffe_only = 0;
postfilter_approach = 0;
@@ -122,8 +108,9 @@ for eq_mode = wh.parameter.eq_mode.values
%%%%% Construct AWG and Scope Modules %%%%%%
SCP = ScopeKeysight("model","DSAZ634A",'autoscale',1,"fadc",'GSa_160',"channel",[1,0],"recordLen",2000000,"removeDC",1);
AWG = AwgKeysight("model","M8196A","fdac",fdac,"scaletodac",[1,1,1,1],"skews",[0,0,0,0],"voltages",[0,0,0,awg_vpp]);
A2S = Awg2Scope(AWG,SCP,[0,0,0,1]);
%AWG = AwgKeysight("model","M8196A","fdac",fdac,"scaletodac",[1,1,1,1],"skews",[0,0,0,0],"voltages",[0,0,0,awg_vpp]);
AWG = AwgKeysight("model","M8199B","fdac",fdac,"scaletodac",[1,1],"skews",[0,0],"voltages",[0,awg_vpp]);
A2S = Awg2Scope(AWG,SCP,[0,1,0,0]);
%%%%% Symbol Generation %%%%%%
Pform = Pulseformer("fsym",fsym,"fdac",4*fsym,"pulse","rrc","pulselength",16,"rrcalpha",rrcalpha);
@@ -146,10 +133,12 @@ for eq_mode = wh.parameter.eq_mode.values
%%%%% Resample to DAC rate %%%%%%
Digi_sig = Digi_sig.resample("fs_out",AWG.fdac);
Digi_sig = Filter('filtdegree',5,"f_cutoff",1.1*(fsym/log2(M)),"fs",Digi_sig.fs,"filterType",filtertypes.gaussian).process(Digi_sig);
%%%%% Plot and Save Routine 1 %%%%%%%%%%%%%%%%%%%%%%%%%
Digi_sig.spectrum("displayname","Normal Tx","fignum",10);
if loopcnt == 1
save([folderpath,[experiment_name,'bits'],loop_name],"Bits");
save([folderpath,[experiment_name,'symbols'],loop_name],"Symbols");
@@ -157,7 +146,7 @@ for eq_mode = wh.parameter.eq_mode.values
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%% AWG --> Scope %%%%%%
[~,~,~,Scpe_sig] = A2S.process("signal4",Digi_sig);
[~,Scpe_sig,~,~] = A2S.process("signal2",Digi_sig);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Scpe_sig.spectrum("displayname","Scope PSD","fignum",20);
@@ -194,17 +183,43 @@ for eq_mode = wh.parameter.eq_mode.values
% Scpe_sig.eye(fsym,M,"fignum",40,"displayname",' after Scope');
voa.readvals();
pd_in = voa.power_state(2);
s_pow = voa.power_state(3);
i_pow = voa.power_state(4);
sir = s_pow- i_pow;
%%%%% EQUALIZE %%%%%%
Eq = FFE("epochs_tr",5,"epochs_dd",5,"len_tr",4096*2,"mu_dd",1e-4,"mu_tr",0,"order",50,"sps",2,"decide",0);
Eq = VNLE("epochs_tr",5,"epochs_dd",5,"len_tr",4096*2,"mu_dd",[0.0004 0.0005 0.0006],"mu_tr",0,"order",[50,7,7],"sps",2,"decide",1);
% Eq = VNLE("epochs_tr",5,"epochs_dd",5,"len_tr",4096*2,"mu_dd",[0.0004 0.0005 0.0006],"mu_tr",0,"order",[50,7,7],"sps",2,"decide",1);
Eq = EQ("Ne",[50,7,7],"Nb",[0,0,0],"training_length",4096*2,"training_loops",5,"dd_loops",5,"K",2,"DCmu",0.0,"DDmu",[0.0004 0.0004 0.0004 0.0004 ],"DFEmu",0.005,"FFEmu",0,"plotfinal",0,"ideal_dfe",1);
Eq = EQ("Ne",[50,0,0],"Nb",[0,0,0],"training_length",4096*2,"training_loops",5,"dd_loops",5,"K",2,"DCmu",0.0,"DDmu",[0.0004 0.0004 0.0004 0.0004 ],"DFEmu",0.005,"FFEmu",0,"plotfinal",0,"ideal_dfe",1);
S = Scpe_sig.signal;
% recursion
N1 = 201;
% Initialize the moving sum for the first window
half_window = (N1 - 1) / 2;
moving_sum = sum(S(1:N1));
% Calculate the first element of R1
cic_filtrd(1:half_window+1) = S(1:half_window+1) - (moving_sum / N1);
% Loop over the signal and apply the recursive moving average subtraction
for n = (half_window+2):(length(S)-half_window)
% Update the moving sum by subtracting the oldest value and adding the new one
moving_sum = moving_sum - S(n-half_window-1) + S(n+half_window);
% Calculate the new value of R1
cic_filtrd(n) = S(n) - (moving_sum / N1);
end
cic_filtrd(n+1:length(S)) = S(n+1:length(S)) - (moving_sum / N1);
Scpe_sig_ = Scpe_sig;
Scpe_sig_.signal = cic_filtrd';
if ffe_only %%%%%%%%%%%%%%%%%%%%%%%%%%%
@@ -212,7 +227,6 @@ for eq_mode = wh.parameter.eq_mode.values
EQ_sig.plot("fignum",50,"displayname",'After EQ');
Rx_bits = PAMmapper(M,0).demap(EQ_sig);
[~,errors_bm,ber,errors] = calc_ber(Rx_bits.signal,Bits.signal,"skip_front",100,"skip_end",150,"returnErrorLocation",1);
@@ -231,46 +245,62 @@ for eq_mode = wh.parameter.eq_mode.values
elseif postfilter_approach %%%%%%%%%%%%%%%%%%%%%%%%%%%
[EQ_sig] = Eq.process(Scpe_sig,Symbols);
[EQ_sig, Noi] = Eq.process(Scpe_sig_,Symbols);
EQ_sig.plot("fignum",50,"displayname",'After EQ','clear',1);
% EQ_sig.plot("displayname",'After VNLE','fignum',90,'clear',1);
Noi = EQ_sig-Symbols;
% Quantization is too far from orig. symbols ->
% error psd is quite different
% Sym_ = PAMmapper(M,0).quantize(EQ_sig);
% Noi_ = Sym_-EQ_sig;
% Noi_.normalize('mode','rms').spectrum('displayname','Noise PSD','fignum',1234,'normalizeTo0dB',1,'normalizeToNyquist',1,'color',cols(nc+1,:));
Rx_bits = PAMmapper(M,0).demap(EQ_sig);
[~,errors_bm,ber_ffe_only,errors] = calc_ber(Rx_bits.signal,Bits.signal,"skip_front",100,"skip_end",150,"returnErrorLocation",1);
[~,~,ber_vnle,~] = calc_ber(Rx_bits.signal,Bits.signal,"skip_front",100,"skip_end",150,"returnErrorLocation",1);
nc = 2;
burg_coeff = arburg(Noi.signal,nc);
EQ_sig.normalize('mode','rms').spectrum('displayname','EQ Out','fignum',1234,'normalizeTo0dB',1,'normalizeToNyquist',1,'color',cols(nc,:));
EQ_sig = EQ_sig.filter(burg_coeff,1);
Noi.normalize('mode','rms').spectrum('displayname','Noise PSD optimal','fignum',1234,'normalizeTo0dB',1,'normalizeToNyquist',1,'color',cols(nc+1,:));
if 1
Noi.spectrum('displayname','Noise PSD','fignum',123)
[h,w] = freqz(1,burg_coeff,length(Noi),"whole",Noi.fs);
h = h/max(abs(h));
hold on
w_ = (w - Noi.fs/2);
plot(w_.*1e-9,20*log10(fftshift(h)),'DisplayName',['', num2str(nc), ' coefficients for burg alg.']);
end
for nc = 1:3
if 0
figure(53);
constellation = unique(Symbols.signal);
received = NaN(numel(constellation),length(Symbols));
for lvl = 1:numel(constellation)
received(lvl,Symbols.signal==constellation(lvl)) = EQ_sig.signal(Symbols.signal==constellation(lvl));
burg_coeff = arburg(Noi.signal,nc);
EQ_sig_filt = EQ_sig.filter(burg_coeff,1);
% EQ_sig.spectrum("displayname","Signal Spectrum after Postfilter","fignum",1234);
EQ_sig_mlse = MLSE("DIR",burg_coeff,"duobinary_output",0,"M",M,"trellis_states",PAMmapper(M,0).levels).process(EQ_sig_filt);
% EQ_sig.spectrum("displayname","Signal Spectrum after MLSE","fignum",1234);
if 0
cols = linspecer(12);
EQ_sig_filt.normalize('mode','rms').spectrum('displayname','Noise PSD','fignum',1234,'normalizeTo0dB',1,'normalizeToNyquist',1,'color',cols(nc+2,:));
% [h,w] = freqz(1,burg_coeff,length(Noi),"whole",Noi.fs);
% h = h/max(abs(h));
% hold on
% w_ = (w - Noi.fs/2);
% figure(123)
% plot(w_.*1e-9,20*log10(fftshift(h)),'DisplayName',['', num2str(nc), ' coefficients for burg alg.']);
[h,w] = freqz(1,burg_coeff,length(Noi),"whole");
h = h/max(abs(h));
hold on
histogram(received(lvl,:),1000,"EdgeAlpha",0);
w_ = (w - pi);
plot(w_,20*log10(fftshift(h)),'DisplayName',['', num2str(nc), ' coefficients for burg alg.']);
end
Rx_bits = PAMmapper(M,0).demap(EQ_sig_mlse);
[~,errors_bm,ber(nc),errors] = calc_ber(Rx_bits.signal,Bits.signal,"skip_front",100,"skip_end",150,"returnErrorLocation",1);
% disp(['BER: ',sprintf('%.1E',ber_mlse(i,j)),' - - ROP: ',num2str(patten(i)),'dBm - - PAM-',num2str(M),' - - ',num2str(fsym*1e-9),' GBd']);
end
EQ_sig = MLSE("DIR",burg_coeff,"duobinary_output",0,"M",M,"trellis_states",PAMmapper(M,0).levels).process(EQ_sig);
Rx_bits = PAMmapper(M,0).demap(EQ_sig);
[~,errors_bm,ber,errors] = calc_ber(Rx_bits.signal,Bits.signal,"skip_front",100,"skip_end",150,"returnErrorLocation",1);
disp(['FFE: ',sprintf('%.1E',ber_ffe_only),' -> PF -> MLSE: ',sprintf('%.1E',ber),' dB | PD_in: ',num2str(pd_in),' dBm']);
disp(['FFE: ',sprintf('%.1E',ber_vnle),' -> PF -> MLSE: ',sprintf('%.1E',ber(nc)),' dB | PD_in: ',num2str(pd_in),' dBm']);
elseif db_channel_approach %%%%%%%%%%%%%%%%%%%%%%%%%%%
@@ -301,8 +331,8 @@ for eq_mode = wh.parameter.eq_mode.values
end
end
wh.addValueToStorage(ber_vnle,'vnle',v_bias,awg_vpp,eq_mode,i_atten);
wh.addValueToStorage(ber,'ber',v_bias,awg_vpp,eq_mode,i_atten);
wh.addValueToStorage(pd_in,'pd_in',v_bias,awg_vpp,eq_mode,i_atten);
wh.addValueToStorage(Rx_bits,'signals',v_bias,awg_vpp,eq_mode,i_atten);
@@ -312,13 +342,13 @@ for eq_mode = wh.parameter.eq_mode.values
wh.addValueToStorage(M,'m',v_bias,awg_vpp,eq_mode,i_atten);
showCurrentMeasurement('BER', ber, 'PD in', pd_in, 'PAM',M, 'Vbias', v_bias, 'AWG Vpp', awg_vpp, 'Precomp MaxAmp',precomp_amp_max);
iterationTimes(loopcnt) = toc(iterationStartTime);
averageTimePerIteration = mean(iterationTimes(1:loopcnt));
estimatedTotalTime = averageTimePerIteration * looptotal;
estimatedTimeRemaining = estimatedTotalTime - sum(iterationTimes(1:loopcnt));
%autoArrangeFigures(3,3,2);
end
end
end
@@ -337,18 +367,12 @@ awg_vpp = wh.parameter.awg_vpp.values(1);
eq_mode = wh.parameter.eq_mode.values(1);
bers = wh.getStoValue('ber',v_bias,awg_vpp,eq_mode,i_atten_vals);
sirs = wh.getStoValue('sir',v_bias,awg_vpp,eq_mode,i_atten_vals);
figure(90);
hold on; % Retain the plot so new points can be added without complete redraw
% Plot the data and get the line handle
hLine = plot(i_atten_vals, bers, "LineWidth", 0.5, "LineStyle", "-", "Marker", ".", "MarkerSize", 15, "DisplayName", experiment_name);
% Customize the data tips
% Set labels for existing data tip rows
hLine.DataTipTemplate.DataTipRows(1).Label = 'Fsym';
hLine.DataTipTemplate.DataTipRows(2).Label = 'BER';
hLine.DataTipTemplate.DataTipRows(2).Format = '%.2e'; % Format BER as "3e-4"
plot(sirs, bers, "LineWidth", 0.5, "LineStyle", "-", "Marker", ".", "MarkerSize", 15, "DisplayName", experiment_name);
% Continue with the rest of your plot settings