392 lines
17 KiB
Matlab
392 lines
17 KiB
Matlab
|
|
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.5];
|
|
params.awg_vpp = [0.35];
|
|
params.eq_mode = [2];
|
|
params.i_atten = [0:4:40];
|
|
|
|
wh = DataStorage(params);
|
|
|
|
wh.addStorage("ber");
|
|
wh.addStorage("pd_in");
|
|
wh.addStorage("m");
|
|
wh.addStorage("sir");
|
|
wh.addStorage("s_pow");
|
|
wh.addStorage("i_pow");
|
|
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 = 2;
|
|
|
|
M = 4;
|
|
pn_key = 2;
|
|
usemrds = 0;
|
|
fsym = 68e9;
|
|
fdac = 256e9;
|
|
awg_vpp = 0.35;
|
|
fadc = 160e9;
|
|
rrcalpha = 0.05;
|
|
v_bias = 2.25;
|
|
pd_in_set = 6;
|
|
|
|
looptotal = prod(wh.dim);
|
|
iterationTimes = zeros(looptotal, 1); % Preallocate for speed
|
|
|
|
disp(['Start Measurement of ',num2str(looptotal),' loops...'])
|
|
hWaitbar = waitbar(0, 'Starting measurement...', 'Name', 'Processing Progress');
|
|
loopcnt = 0;
|
|
|
|
estimatedTimeRemaining = 0;
|
|
estimatedTotalTime = 0;
|
|
for eq_mode = wh.parameter.eq_mode.values
|
|
for i_atten = wh.parameter.i_atten.values
|
|
for v_bias = wh.parameter.vbias.values
|
|
for awg_vpp = wh.parameter.awg_vpp.values
|
|
|
|
iterationStartTime = tic;
|
|
|
|
loop_name = ['_iatten_',num2str(i_atten)];
|
|
|
|
loopcnt = loopcnt+1;
|
|
progressFraction = loopcnt / looptotal;
|
|
waitbar(progressFraction, hWaitbar, ...
|
|
sprintf('Progress: %d/%d\nEstimated time remaining: %.2f hours\nEstimated time remaining: %.2f hours', ...
|
|
loopcnt, looptotal, estimatedTimeRemaining/60/60, estimatedTotalTime/60/60));
|
|
|
|
|
|
|
|
switch eq_mode
|
|
case 1
|
|
ffe_only = 1;
|
|
postfilter_approach = 0;
|
|
db_channel_approach = 0;
|
|
db_coding_approach = 0;
|
|
|
|
db_precode = db_coding_approach || db_channel_approach;
|
|
|
|
case 2
|
|
ffe_only = 0;
|
|
postfilter_approach = 1;
|
|
db_channel_approach = 0;
|
|
db_coding_approach = 0;
|
|
|
|
db_precode = db_coding_approach || db_channel_approach;
|
|
|
|
case 3
|
|
ffe_only = 0;
|
|
postfilter_approach = 0;
|
|
db_channel_approach = 1;
|
|
db_coding_approach = 0;
|
|
|
|
db_precode = db_coding_approach || db_channel_approach;
|
|
case 4
|
|
ffe_only = 0;
|
|
postfilter_approach = 0;
|
|
db_channel_approach = 0;
|
|
db_coding_approach = 1;
|
|
|
|
db_precode = db_coding_approach || db_channel_approach;
|
|
end
|
|
|
|
|
|
|
|
%%%%% SET Voltages %%%%%%
|
|
dcs = DC_supply("active",[1,1],"voltage",[v_bias, 9]);
|
|
dcs.set("voltage",[v_bias, 9]);
|
|
|
|
%%%%% SET Attenuator %%%%%%
|
|
voa = OptAtten("active",[1,2,1,1],"value",[0,pd_in_set,0,i_atten],"wavelength",[1310,1310,1310,1310]);
|
|
voa.set('active',[1,2,1,1],'value',[0,pd_in_set,0,i_atten]);
|
|
% voa.readvals();
|
|
|
|
%%%%% 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]);
|
|
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);
|
|
|
|
[Digi_sig,Symbols,Bits] = PAMsource("fsym",fsym,"M",M,"order",19,"useprbs",1,...
|
|
"fs_out",fdac,"applyclipping",0,"clipfactor",1.7,...
|
|
"applypulseform",0,"pulseformer",Pform,"randkey",pn_key,...
|
|
"db_precode",db_precode,...
|
|
"mrds_code",usemrds,"mrds_blocklength",512,"db_encode",db_coding_approach).process();
|
|
|
|
%%%%% Precompensation Routine %%%%%%
|
|
if precomp_mode == 1 % measure channel
|
|
freqresp = ChannelFreqResp("Nacq",1024,"Navg",64,"Ncp",63,'f_ref',Digi_sig.fs);
|
|
Digi_sig = freqresp.buildOFDM();
|
|
elseif precomp_mode == 2 % apply precomp
|
|
freqresp = ChannelFreqResp("Nacq",1024,"Navg",64,"Ncp",63,'f_ref',Digi_sig.fs);
|
|
Digi_sig = freqresp.precomp(Digi_sig,'maxampdb',precomp_amp_max,'loadPath',precomp_path,'fileName',precomp_fn);
|
|
end
|
|
|
|
%%%%% 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");
|
|
end
|
|
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
|
|
|
%%%%% AWG --> Scope %%%%%%
|
|
[~,Scpe_sig,~,~] = A2S.process("signal2",Digi_sig);
|
|
|
|
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
|
Scpe_sig.spectrum("displayname","Scope PSD","fignum",20);
|
|
|
|
% Scpe_sig.plot("displayname","Scope raw signal","fignum",25);
|
|
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
|
|
|
|
|
%%%%%% Sample to 2x fsym %%%%%%
|
|
Scpe_sig = Scpe_sig.resample("fs_in",160e9,"fs_out",2*fsym);
|
|
|
|
%%%%% Precompensation Routine %%%%%%
|
|
if precomp_mode == 1
|
|
freqresp.estimate(Scpe_sig,"save",true,"savePath",precomp_path,"fileName",precomp_fn);
|
|
freqresp.plot();
|
|
end
|
|
|
|
%%%%%% Sync Rx signal with reference %%%%%%
|
|
[Scpe_sig,S] = Scpe_sig.tsynch("reference",Symbols,"fs_ref",fsym);
|
|
|
|
%%%%%% SNR CHEAT - Avg. the measured signal occurences %%%%%%
|
|
average_signals = 0;
|
|
if average_signals
|
|
scope_mean = zeros(size(S{1}.signal));
|
|
for n=1:numel(S)
|
|
scope_mean = scope_mean + S{n}.signal;
|
|
end
|
|
scope_mean = scope_mean ./ n;
|
|
Scpe_sig.signal = scope_mean;
|
|
end
|
|
|
|
%%%%% Plot and Save Routine 2 %%%%%%%%%%%%%%%%%%%%%%%%%
|
|
save([folderpath,experiment_name,'rx_signal',loop_name],"S");
|
|
% 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 = 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 %%%%%%%%%%%%%%%%%%%%%%%%%%%
|
|
|
|
[EQ_sig] = Eq.process(Scpe_sig,Symbols);
|
|
|
|
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);
|
|
|
|
disp(['FFE: ',sprintf('%.1E',ber),'| SIR: ',num2str(sir),' dB | PD_in: ',num2str(pd_in),' dBm']);
|
|
|
|
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));
|
|
hold on
|
|
histogram(received(lvl,:),1000,"EdgeAlpha",0);
|
|
end
|
|
end
|
|
|
|
elseif postfilter_approach %%%%%%%%%%%%%%%%%%%%%%%%%%%
|
|
|
|
[EQ_sig, Noi] = Eq.process(Scpe_sig_,Symbols);
|
|
|
|
% EQ_sig.plot("displayname",'After VNLE','fignum',90,'clear',1);
|
|
|
|
% 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);
|
|
[~,~,ber_vnle,~] = calc_ber(Rx_bits.signal,Bits.signal,"skip_front",100,"skip_end",150,"returnErrorLocation",1);
|
|
|
|
EQ_sig.normalize('mode','rms').spectrum('displayname','EQ Out','fignum',1234,'normalizeTo0dB',1,'normalizeToNyquist',1,'color',cols(nc,:));
|
|
|
|
Noi.normalize('mode','rms').spectrum('displayname','Noise PSD optimal','fignum',1234,'normalizeTo0dB',1,'normalizeToNyquist',1,'color',cols(nc+1,:));
|
|
|
|
for nc = 1:3
|
|
|
|
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
|
|
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
|
|
|
|
disp(['FFE: ',sprintf('%.1E',ber_vnle),' -> PF -> MLSE: ',sprintf('%.1E',ber(nc)),' dB | PD_in: ',num2str(pd_in),' dBm']);
|
|
|
|
elseif db_channel_approach %%%%%%%%%%%%%%%%%%%%%%%%%%%
|
|
|
|
[EQ_sig, Noi] = Eq.process(Scpe_sig,Duobinary().encode(Symbols));
|
|
|
|
EQ_sig.plot("fignum",50,"displayname",'After EQ','clear',1);
|
|
|
|
EQ_sig = MLSE("DIR",[1,1],"duobinary_output",1,"M",M,"trellis_states",PAMmapper(M,0).levels).process(EQ_sig);
|
|
EQ_sig = Duobinary().decode(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([' DB Precode -> Channel -> FFE -> Decode/ Mod ',sprintf('%.1E',ber),' | PD_in: ',num2str(pd_in),' dBm']);
|
|
|
|
elseif db_coding_approach %%%%%%%%%%%%%%%%%%%%%%%%%%%
|
|
|
|
[EQ_sig, Noi] = Eq.process(Scpe_sig,Symbols);
|
|
EQ_sig = MLSE("DIR",[1,1],"duobinary_output",1,"M",M,"trellis_states",PAMmapper(M,0).levels).process(EQ_sig);
|
|
EQ_sig = Duobinary().decode(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);
|
|
|
|
% EQ_sig.plot("fignum",50,"displayname",'After EQ');
|
|
|
|
disp([' DB Precode -> DB Code -> Channel -> FFE -> Decode/ Mod ',sprintf('%.1E',ber),' | PD_in: ',num2str(pd_in),' dBm']);
|
|
|
|
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);
|
|
wh.addValueToStorage(sir,'sir',v_bias,awg_vpp,eq_mode,i_atten);
|
|
wh.addValueToStorage(s_pow,'s_pow',v_bias,awg_vpp,eq_mode,i_atten);
|
|
wh.addValueToStorage(i_pow,'i_pow',v_bias,awg_vpp,eq_mode,i_atten);
|
|
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
|
|
end
|
|
|
|
close(hWaitbar);
|
|
|
|
wh.save([folderpath,experiment_name,'wh']);
|
|
|
|
|
|
cols = linspecer(8);
|
|
|
|
i_atten_vals = wh.parameter.i_atten.values;
|
|
v_bias = wh.parameter.vbias.values(1);
|
|
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
|
|
plot(sirs, bers, "LineWidth", 0.5, "LineStyle", "-", "Marker", ".", "MarkerSize", 15, "DisplayName", experiment_name);
|
|
|
|
|
|
% Continue with the rest of your plot settings
|
|
yline(3.8e-3, 'DisplayName', 'HD-FEC', 'LineStyle', '--', 'HandleVisibility', 'off');
|
|
xlabel('Signal to Interference Ratio in dB');
|
|
ylabel('Bit Error Rate (BER)');
|
|
title('Bit Error Rate vs. SIR');
|
|
set(gca, 'yscale', 'log');
|
|
set(gca, 'Box', 'on');
|
|
grid on;
|
|
grid minor;
|
|
legend('Interpreter', 'none');
|
|
|
|
|
|
autoArrangeFigures(3,3,2)
|
|
|
|
disp("measurement done") |