From bafc7f12b778c71f1356e4ea912635e2ca80f245 Mon Sep 17 00:00:00 2001 From: Silas Labor Zizou Date: Thu, 24 Oct 2024 19:41:59 +0200 Subject: [PATCH] Commit before first big baudrate sweep --- Classes/00_signals/Signal.m | 14 +- Classes/01_transmit/ChannelFreqResp.m | 12 +- Classes/01_transmit/PAMsource.m | 43 +- Classes/05_Lab/Awg2Scope.m | 28 +- Classes/05_Lab/DC_supply.m | 126 ++-- Classes/05_Lab/OptAtten.m | 5 +- Classes/05_Lab/ScopeKeysight.m | 37 +- Functions/showCurrentMeasurement.m | 4 +- projects/400G_FTN_setups/imdd_mpi_dsp.m | 70 ++- .../bias_optimization.m | 536 +++++++++++------- projects/Lab_2024/bias_sweep_evaluation.m | 31 +- 11 files changed, 588 insertions(+), 318 deletions(-) diff --git a/Classes/00_signals/Signal.m b/Classes/00_signals/Signal.m index a2f70ae..3ecb6e4 100644 --- a/Classes/00_signals/Signal.m +++ b/Classes/00_signals/Signal.m @@ -510,7 +510,7 @@ classdef Signal end %% - function [obj,S] = tsynch(obj,options) + function [obj,S,isFlipped] = tsynch(obj,options) % time sync and cut arguments obj Signal @@ -525,14 +525,16 @@ classdef Signal q = obj.fs/options.fs_ref; b = options.reference.resample("fs_in",options.fs_ref,"fs_out",obj.fs).normalize("mode","oneone").signal; + max_occurences = floor(length(a)/length(b)); + %estimate delay between signals [co,lags] = xcorr(a,b); - [~,pos] = max(co); + [~,pos] = max(abs(co)); D = lags(pos); %estimate start pos of signal maxpeaknum = floor(length(a)/length(b)); - [pks,pkpos] = findpeaks(co./max(co),'MinPeakDistance',length(b)/2,'MinPeakHeight',0.2,'NPeaks',maxpeaknum); + [pks,pkpos] = findpeaks(abs(co./max(co)),'MinPeakDistance',length(b)/2,'MinPeakHeight',0.2,'NPeaks',maxpeaknum); shifts = lags(pkpos); %Cut occurences of ref signal from signal (only positive shifts) @@ -543,6 +545,12 @@ classdef Signal S{end+1,1} = sig; end + % + isFlipped=0; + if all(sign(co(pkpos))) + isFlipped = 1; + end + %return/keep the sinal with the highest correlation (only within positive shifts) [~,idx]=max(pks(shifts>0)); obj.signal = S{idx}.signal; diff --git a/Classes/01_transmit/ChannelFreqResp.m b/Classes/01_transmit/ChannelFreqResp.m index e77cf43..586e580 100644 --- a/Classes/01_transmit/ChannelFreqResp.m +++ b/Classes/01_transmit/ChannelFreqResp.m @@ -202,19 +202,19 @@ classdef ChannelFreqResp < handle function plot(obj) figure(55); - clf; + %clf; Havg = obj.H; %1) - subplot(2,1,1);hold all;box on;title('Magnitude Freq. Response'); + subplot(2,1,1);hold on;box on;title('Magnitude Freq. Response'); plot(obj.faxis/1e9, 20*log10(abs(obj.H_all)),'linewidth',0.1,'LineStyle','-','Color','#808080') ; xlim([0.2 .5*max(obj.faxis)*1e-9]); plot(obj.faxis/1e9, 20*log10(abs(Havg)),'LineWidth',2); grid on; %2) - subplot(2,1,2); hold all; box on; title('Phase Freq. Response'); + subplot(2,1,2); hold on; box on; title('Phase Freq. Response'); plot(obj.faxis/1e9, angle(obj.H_all),'linewidth',0.1,'LineStyle','-','Color','#808080') ; plot(obj.faxis/1e9, unwrap(angle(Havg)),'LineWidth',2) ; xlim([0.2 .5*max(obj.faxis)*1e-9]); @@ -228,7 +228,7 @@ classdef ChannelFreqResp < handle Havg = Havg./mean(Havg(2:10)); %3) - subplot(2,1,1); hold all; box on; title('Inverse Magnitude Freq. Response'); + subplot(2,1,1); hold on; box on; title('Inverse Magnitude Freq. Response'); plot(obj.faxis/1e9, 20*log10(abs(1./Havg)),"LineWidth",2,"Color",[0.3467 0.5360 0.6907]) ; xlim([0.2 .5*max(obj.faxis)*1e-9]); grid on; ylim([-1 15]); @@ -236,12 +236,12 @@ classdef ChannelFreqResp < handle yline(3,'LineWidth',2,'LineStyle','--'); %4) - subplot(2,1,2); hold all; box on; title('Inverse Phase Freq. Response'); + subplot(2,1,2); hold on; box on; title('Inverse Phase Freq. Response'); plot(obj.faxis/1e9, unwrap(angle(1./Havg)),"LineWidth",2,"Color",[0.3467 0.5360 0.6907]) ; xlim([0.2 .5*max(obj.faxis)*1e-9]); grid on; %%% plot for publication - figure(98989);hold all;box on;title('Magnitude Freq. Response'); + figure(98989);hold on;box on;title('Magnitude Freq. Response'); xlim([0.2 .5*max(obj.faxis)*1e-9]); ylim([-20, 2]); plot(obj.faxis/1e9, 20*log10(abs(Havg)),'LineWidth',2); diff --git a/Classes/01_transmit/PAMsource.m b/Classes/01_transmit/PAMsource.m index 099b6be..4e96911 100644 --- a/Classes/01_transmit/PAMsource.m +++ b/Classes/01_transmit/PAMsource.m @@ -77,9 +77,42 @@ classdef PAMsource bitpattern=[]; if obj.useprbs - for i = 1:log2(obj.M) - [bitpattern(:,i),seed] = prbs(O,N,seed); + % for i = 1:log2(obj.M) + % [bitpattern(:,i),seed] = prbs(O,N,seed); + % end + + %%%%% MOVE-IT PRMS %%%% + + state = struct(); + + para = struct(); + + if obj.M == 6 + para.bl = 2^(obj.order-2); + para.dimension = 5; + else + para.bl = 2^(obj.order-1); + para.dimension = log2(obj.M); %2.5bits/sym -> 2 bit/sym end + + para.rand = 0; + + para.order = floor(obj.order / log2(obj.M)); + para.skip =0; + para.bruijn = 0; + para.reset_prms = 0; + para.method = 1; + + data_in = []; + global loop; + loop = 0; + [data_out,state_] = prms(data_in, state, para); + loop = 1; + [data_out,state_out] = prms(data_in, state_, para); + bitpattern = data_out'; + + %%%%% END MOVE-IT %%%%%%% + else s = RandStream('twister','Seed',obj.randkey); for i = 1:log2(obj.M) @@ -88,7 +121,7 @@ classdef PAMsource end if obj.M == 6 - bitpattern = reshape(bitpattern,[],1); + bitpattern = reshape(bitpattern',[],1); bitpattern = bitpattern(1:end-mod(length(bitpattern),5)); end @@ -108,7 +141,7 @@ classdef PAMsource end % figure(12);hold on;histogram(symbols.signal,'Normalization','probability'); - + if obj.mrds_code symbols = MRDS_coding("blocklength",obj.mrds_blocklength).encode(symbols); end @@ -128,7 +161,7 @@ classdef PAMsource %%%%% Re-sample to f DAC %%%%%% digi_sig = digi_sig.resample("fs_in",digi_sig.fs,"fs_out",obj.fs_out,"n",10,"beta",5); - % digi_sig.spectrum("fignum",111,"displayname","after pulseforming"); + % digi_sig.spectrum("fignum",111,"displayname","after pulseforming"); %%%%% Hard clip digital signal to PAM range before DAC %%%%%% if obj.applyclipping diff --git a/Classes/05_Lab/Awg2Scope.m b/Classes/05_Lab/Awg2Scope.m index 765546c..e097bb9 100644 --- a/Classes/05_Lab/Awg2Scope.m +++ b/Classes/05_Lab/Awg2Scope.m @@ -7,11 +7,12 @@ classdef Awg2Scope Scope mapping; + waitUntilClick end methods (Access=public) - function obj = Awg2Scope(Awg,Scope,mapping) + function obj = Awg2Scope(Awg,Scope,mapping,options) %Simple class to call the Awg and Scope and map the signals %accordingly in the correct formats with correct l % ogbook @@ -21,16 +22,19 @@ classdef Awg2Scope Awg Scope mapping + + options.waitUntilClick = 0; end obj.Awg = Awg; obj.Scope = Scope; obj.mapping = mapping; % AWG CH [1,2,3,4] -> Scope CH [0,0,0,1] + obj.waitUntilClick = options.waitUntilClick; end - function [S1,S2,S3,S4] = process(obj,channels) + function [S1,S2,S3,S4] = process(obj,channels, options) arguments obj @@ -41,15 +45,33 @@ classdef Awg2Scope channels.signal4 Informationsignal = Informationsignal([]) % add new optional arguments here + options.waitUntilClick = obj.waitUntilClick; end + + + %%% UPLOAD TO AWG %%% + [S1,S2,S3,S4]=obj.Awg.upload("signal1",channels.signal1,... "signal2",channels.signal2,... "signal3",channels.signal3,... "signal4",channels.signal4... ); + + %%% UPLOAD TO AWG %%% - scpe_sig_cell = obj.Scope.read(); + + + %%% READ FROM SCOPE %%% + + scpe_sig_cell = obj.Scope.read("waitUntilClick",options.waitUntilClick); + + %%% READ FROM SCOPE %%% + + + + + % Map Scope measurement to output signal % mapping index is the AWG chanel and mapping number is the diff --git a/Classes/05_Lab/DC_supply.m b/Classes/05_Lab/DC_supply.m index 5c17d68..99bef76 100644 --- a/Classes/05_Lab/DC_supply.m +++ b/Classes/05_Lab/DC_supply.m @@ -51,90 +51,90 @@ classdef DC_supply < handle success = [0,0]; - try - % %connect to device - % if exist('v','var') - % %v = visadev("GPIB1::19::INSTR"); - % v = obj.connectDevice(); - % else - % - % end + % %connect to device + % if exist('v','var') + % %v = visadev("GPIB1::19::INSTR"); + % v = obj.connectDevice(); + % else + % + % end - v = visadev("GPIB1::19::INSTR"); + v = visadev("GPIB1::19::INSTR"); - debug = 0; - if debug - disp(['Connected to Instrument: ',char(v.Vendor),' ',char(v.Model),' SerNo:',char(v.SerialNumber)]); - end + debug = 0; + if debug + disp(['Connected to Instrument: ',char(v.Vendor),' ',char(v.Model),' SerNo:',char(v.SerialNumber)]); + end - cmd = 'INST:SEL?'; + cmd = 'INST:SEL?'; + writeline(v, cmd); + prev_selected_channel = readline(v); + + for ch = 1:2 + + % choose channel + cmd = ['INST:SEL OUT',num2str(ch)]; writeline(v, cmd); - prev_selected_channel = readline(v); - for ch = 1:2 + % get current voltage level + cmd = 'VOLT?'; + writeline(v, cmd); + act_volt = str2num(readline(v)); - % choose channel - cmd = ['INST:SEL OUT',num2str(ch)]; - writeline(v, cmd); + % desired voltage (round to two digits after comma) + des_volt = round(options.voltage(ch),2); + + cnt = 0; + while abs(act_volt-des_volt) > 0 + + selected_channel = ['OUT',num2str(ch)]; % get current voltage level cmd = 'VOLT?'; writeline(v, cmd); act_volt = str2num(readline(v)); - % desired voltage (round to two digits after comma) - des_volt = round(options.voltage(ch),2); + % difference + diff_volt = act_volt-des_volt; - cnt = 0; - while abs(act_volt-des_volt) > 0 - - selected_channel = ['OUT',num2str(ch)]; - - % get current voltage level - cmd = 'VOLT?'; - writeline(v, cmd); - act_volt = str2num(readline(v)); - - % difference - diff_volt = act_volt-des_volt; - - % set new voltage - increment_voltage = -0.01* sign(diff_volt) ; - cmd = ['VOLT ',num2str(act_volt+increment_voltage)]; - writeline(v, cmd); - - cnt = cnt+1; - if mod(cnt,100) == 0 - wait(1); - end - - % get current voltage level - cmd = 'VOLT?'; - writeline(v, cmd); - act_volt = str2num(readline(v)); + % set new voltage + increment_voltage = -0.01* sign(diff_volt) ; + cmd = ['VOLT ',num2str(act_volt+increment_voltage)]; + writeline(v, cmd); + cnt = cnt+1; + if mod(cnt,100) == 0 + pause(1); end - % check if voltage is set - if act_volt ~= des_volt - hMsgBox = msgbox('An error occurred in dc supply module. Check if voltage is set correctly.'); - uiwait(hMsgBox); - else - success(ch) = 1; - end + % get current voltage level + cmd = 'VOLT?'; + writeline(v, cmd); + act_volt = str2num(readline(v)); end - % choose channel - cmd = ['INST:SEL ',char(strtrim(prev_selected_channel))]; - writeline(v, cmd); + % check if voltage is set + if act_volt ~= des_volt + hMsgBox = msgbox('An error occurred in dc supply module. Check if voltage is set correctly.'); + uiwait(hMsgBox); + else + success(ch) = 1; + end - %disconnect - delete(v); - - catch - end + + % choose channel + cmd = ['INST:SEL ',char(strtrim(prev_selected_channel))]; + writeline(v, cmd); + + %disconnect + delete(v); + + + + + end function [voltage,current] = readVals(obj) diff --git a/Classes/05_Lab/OptAtten.m b/Classes/05_Lab/OptAtten.m index 84e98b5..93deccc 100644 --- a/Classes/05_Lab/OptAtten.m +++ b/Classes/05_Lab/OptAtten.m @@ -193,8 +193,9 @@ classdef OptAtten < handle differ = abs(state.outputpower(i)-options.value(i)); %create msgbox when necessary if differ >= 2 - hMsgBox = msgbox(['Output Power at attenuator slot ' num2str(2) ' differs by 2dB or more!']); - uiwait(hMsgBox); + warning(['Output Power at attenuator slot ' num2str(2) ' differs by ', num2str(differ), 'or more!']) + % hMsgBox = msgbox(['Output Power at attenuator slot ' num2str(2) ' differs by 2dB or more!']); + % uiwait(hMsgBox); end end end diff --git a/Classes/05_Lab/ScopeKeysight.m b/Classes/05_Lab/ScopeKeysight.m index 4565e69..3635a30 100644 --- a/Classes/05_Lab/ScopeKeysight.m +++ b/Classes/05_Lab/ScopeKeysight.m @@ -12,6 +12,8 @@ classdef ScopeKeysight interpolate recordLen IPaddress + + waitUntilClick end methods (Access=public) @@ -30,6 +32,7 @@ classdef ScopeKeysight options.interpolate logical = 0; options.recordLen double = 1000000; options.IPaddress + options.waitUntilClick = 0; end % @@ -55,6 +58,7 @@ classdef ScopeKeysight arguments obj options.channel logical = obj.channel + options.waitUntilClick = obj.waitUntilClick; end %%%%%%%%%%%%%%%%%%%%%%%% @@ -119,14 +123,43 @@ classdef ScopeKeysight obj.writeNcheck(v,sprintf(':CHANnel%u:DISPlay ON',n));% display captured data trace end + set(v,'Timeout',150); if obj.autoScale obj.writeNcheck(v,':AUTOscale'); for n = 1:4 range = str2double(obj.writeReceiveCheck(v,sprintf(':CHANnel%u:RANGe?',n))); - obj.writeNcheck(v,sprintf(':CHANnel%u:RANGe %.3f',n,range/1.2)); + + %use this to finetune autoscaling - VERY helpful + % higher value leads to higher scaling + fintunefactor = 1.4; % within [1,...,2] + obj.writeNcheck(v,sprintf(':CHANnel%u:RANGe %.3f',n,range/fintunefactor)); end else - obj.writeNcheck(v,':SINGLE'); + %obj.writeNcheck(v,':SINGLE'); + end + + % After Autoscale, let the user adjust the scope scaling... + % "options.waitUntilClick" was a shit name but now this is it :-) + if options.waitUntilClick + % Create a dialog box with the desired text + d = dialog('Name', 'Scale Scope then press continue', 'Position', [300, 300, 300, 180]); + + % Add a text label with the instruction + uicontrol('Parent', d, ... + 'Style', 'text', ... + 'Position', [20, 100, 260, 40], ... + 'String', 'Please adjust scope scaling now, then press continue', ... + 'HorizontalAlignment', 'center'); + + % Add a button to close the dialog and resume execution + uicontrol('Parent', d, ... + 'Style', 'pushbutton', ... + 'Position', [100, 40, 100, 40], ... + 'String', 'Continue', ... + 'Callback', 'uiresume(gcbf); delete(gcbf)'); + + % Pause execution until the dialog box is closed + uiwait(d); end if obj.extRef diff --git a/Functions/showCurrentMeasurement.m b/Functions/showCurrentMeasurement.m index b9710b8..7205f6f 100644 --- a/Functions/showCurrentMeasurement.m +++ b/Functions/showCurrentMeasurement.m @@ -39,9 +39,9 @@ function showCurrentMeasurement(varargin) % Convert values to strings for display for i = 1:length(values) - varNameLower = lower(names{i}); % Convert variable name to lowercase for case-insensitive comparison + varNameLower = names{i}; % Convert variable name to lowercase for case-insensitive comparison if isnumeric(values{i}) - if strcmp(varNameLower, 'ber') + if any(strcmpi(varNameLower, {'ber','mlse','ffe','db'})) % Format 'ber' values in exponential notation with two decimal places values{i} = sprintf('%.2e', values{i}); else diff --git a/projects/400G_FTN_setups/imdd_mpi_dsp.m b/projects/400G_FTN_setups/imdd_mpi_dsp.m index d9464c2..761d473 100644 --- a/projects/400G_FTN_setups/imdd_mpi_dsp.m +++ b/projects/400G_FTN_setups/imdd_mpi_dsp.m @@ -2,31 +2,31 @@ %% Parameter to simulate and save params = struct; -params.M = [4]; -params.datarate = [448]; -params.rop = [0]; +params.M = [6]; +params.datarate = [250]; +params.rop = [-10]; params.sir = 40;%15:1:40; -params.random_key_laser_phase = 10:20; +params.random_key_laser_phase = 1; -precomp_mode = 0; %0=do nothing ; 1= measure; 2=precomp active -postfilter = 0; % noise whiten. approach -> Postfilter + MLSE +precomp_mode = 1; %0=do nothing ; 1= measure; 2=precomp active +postfilter = 1; % noise whiten. approach -> Postfilter + MLSE -db_precode = 1; +db_precode = 0; db_encode = 0; -db_channelapproach = 1; +db_channelapproach = 0; -laser_linewidth = 50e5; +laser_linewidth = 0e5; random_key_sequence = 15; random_key_laser_phase = 66; -sir = 20; +sir = 60; if ismac precomp_path = "/Users/silasoettinghaus/Documents/MATLAB/imdd_simulation/projects/standard_system"; else - precomp_path = "C:\Users\Silas\Documents\MATLAB\imdd_simulation\projects\standard_system\"; + precomp_path = "C:\Users\sioe\Documents\MATLAB\imdd_simulation\projects\standard_system\"; end -precomp_fn = "400G_simulative_setup"; +precomp_fn = "400G_simulative_setup_meas"; usemrds = 0; @@ -78,9 +78,9 @@ for M = wh.parameter.M.values % Digi_sig.normalize("mode","rms").spectrum("displayname",'Tx Signal','fignum',10); if precomp_mode == 1 %measure - freqresp = ChannelFreqResp("Nacq",1024,"Navg",64,"Ncp",63,'f_ref',Digi_sig.fs); - Digi_sig = freqresp.buildOFDM(); - Digi_sig_I = freqresp.buildOFDM(); + precomp_est = ChannelFreqResp("Nacq",1024,"Navg",64,"Ncp",63,'f_ref',Digi_sig.fs); + Digi_sig = precomp_est.buildOFDM(); + Digi_sig_I = precomp_est.buildOFDM(); elseif precomp_mode == 2 %apply Digi_sig = ChannelFreqResp("Nacq",1024,"Navg",64,"Ncp",63,'f_ref',Digi_sig.fs).precomp(Digi_sig,'maxampdb',3,'loadPath',precomp_path,'fileName',precomp_fn); Digi_sig_I = ChannelFreqResp("Nacq",1024,"Navg",64,"Ncp",63,'f_ref',Digi_sig_I.fs).precomp(Digi_sig_I,'maxampdb',3,'loadPath',precomp_path,'fileName',precomp_fn); @@ -117,7 +117,7 @@ for M = wh.parameter.M.values % MAIN SIGNAL %%%%% MODULATE E/O CONVERSION %%%%%% - vbias_rel = 0.5; + vbias_rel = 0.7; u_pi = 2.9; vbias = -vbias_rel*u_pi; @@ -180,10 +180,12 @@ for M = wh.parameter.M.values "adcresolution",10,"quantbuffer",0.1,'block_dc',1,'lpf_active',1,'H_lpf',Lp_scpe).process(Rx_sig); if precomp_mode == 1 - freqresp.estimate(Scpe_sig,"save",true,"savePath",precomp_path,"fileName",precomp_fn); - freqresp.plot(); + precomp_est.estimate(Scpe_sig,"save",true,"savePath",precomp_path,"fileName",precomp_fn); + precomp_est.plot(); end + + % Scpe_sig_normalized = Scpe_sig.normalize("mode","rms"); % Scpe_sig.normalize("mode","rms").spectrum("displayname",'After Scope','fignum',10); @@ -193,6 +195,8 @@ for M = wh.parameter.M.values %%%%%% Sync Rx signal with reference %%%%%% [Scpe_sig,S] = Scpe_sig.tsynch("reference",Symbols,"fs_ref",fsym); + + Scpe_sig.eye(fsym,M,"fignum",50,"displayname",'Simulated after Scope'); %%%%% EQUALIZE %%%%%% Eq = FFE("epochs_tr",5,"epochs_dd",5,"len_tr",4096*2,"mu_dd",1e-4,"mu_tr",0,"order",25,"sps",2,"decide",0); @@ -237,8 +241,34 @@ for M = wh.parameter.M.values % 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(i,j),~] = calc_ber(Rx_bits.signal,Bits.signal,"skip_front",100,"skip_end",150,"returnErrorLocation",1); + [~,~,ber_ffe(i,j),~] = calc_ber(Rx_bits.signal,Bits.signal,"skip_front",100,"skip_end",150,"returnErrorLocation",1); + + if 1 + figure(55); + clf + title(sprintf('PAM %d ; BER: %1.2e',M, ber_ffe(i,j))); + constellation = unique(Symbols.signal); + received = NaN(numel(constellation),length(Symbols)); + for lvl = 1:numel(constellation) + %Separate the equalized signal into the + %respective levels based on the actually + %transmitted level! + received(lvl,Symbols.signal==constellation(lvl)) = EQ_sig.signal(Symbols.signal==constellation(lvl)); + intermediate = received(lvl,:); + cnt(lvl) = numel(intermediate(~isnan(intermediate))); + hold on + histogram(received(lvl,:),1000,"EdgeAlpha",0,'DisplayName',['Lvl ',num2str(lvl),' | ',num2str(cnt(lvl)),' entries']); + end + legend + end + + + + + + + cols = linspecer(8); 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,:)); @@ -256,7 +286,7 @@ for M = wh.parameter.M.values % EQ_sig.spectrum("displayname","Signal Spectrum after MLSE","fignum",1234); if 1 - cols = linspecer(12); + diff --git a/projects/HighSpeedExperiment_2024/bias_optimization.m b/projects/HighSpeedExperiment_2024/bias_optimization.m index fecc926..287c768 100644 --- a/projects/HighSpeedExperiment_2024/bias_optimization.m +++ b/projects/HighSpeedExperiment_2024/bias_optimization.m @@ -1,6 +1,9 @@ folderpath = 'C:\Users\sioe\Documents\High_Speed_Measurement_2024\bias_testing\'; -experiment_name = 'PAM4_b2b_'; +experiment_name = 'PAM6_b2b_rop_referenz_Thormax'; +currentTime = datetime('now', 'Format', 'yyyyMMdd_HHmmss'); +timeStr = char(currentTime); +experiment_name = [experiment_name, timeStr]; ffe_only = 0; postfilter_approach = 1; @@ -12,12 +15,17 @@ db_precode = db_coding_approach || db_channel_approach; %%% SIR Sweep for MPI Experiment %%% params = struct; -params.vbias = [3]; +params.vbias = [2.67]; % PAM6=2.3V %PAM8=2.68V params.awg_vpp = [2.7]; +params.precomp_amp_max = [5]; +params.rop_atten = [0]; +params.M = [8]; wh = DataStorage(params); -wh.addStorage("ber"); +wh.addStorage("ber_ffe"); +wh.addStorage("ber_mlse"); +wh.addStorage("ber_db"); wh.addStorage("pd_in"); wh.addStorage("rop"); wh.addStorage("m"); @@ -26,233 +34,343 @@ wh.addStorage("signals"); precomp_path = "C:\Users\sioe\Documents\High_Speed_Measurement_2024\precomp\"; precomp_fn = "lab_high_speed"; precomp_mode = 2; %0=do nothing ; 1= measure; 2=precomp active + precomp_amp_max = 4; +random_key = 2; +pd_in_set = 7; -M = 4; -pn_key = 2; -usemrds = 0; -fsym = 170e9; -fdac = 256e9; -awg_vpp = 0.35; -fadc = 256e9; -rrcalpha = 0.05; -v_bias = 2.25; -pd_in_set = 6; -rop_atten = 0; +looptotal = prod(wh.dim); + +disp(['Start Measurement of ',num2str(looptotal),' loops...']) +iterationTimes = zeros(looptatal, 1); % Preallocate for speed +if ~exist('hWaitbar', 'var') || ~isvalid(hWaitbar) + hWaitbar = waitbar(0, sprintf('Starting %d measurements',looptotal), 'Name', 'Processing Progress'); +else + waitbar(0, hWaitbar, sprintf('Starting %d measurements',looptotal)); +end -looptatal = prod(wh.dim); -disp(['Start Measurement of ',num2str(looptatal),' loops...']) -hWaitbar = waitbar(0, 'Starting measurement...', 'Name', 'Processing Progress'); loopcnt = 0; +estimatedTimeRemaining = 0; +estimatedTotalTime = 0; for v_bias = wh.parameter.vbias.values - for awg_vpp = wh.parameter.awg_vpp.values + for rop_atten = wh.parameter.rop_atten.values + for precomp_amp_max = wh.parameter.precomp_amp_max.values + for M = wh.parameter.M.values + + for awg_vpp = wh.parameter.awg_vpp.values + + iterationStartTime = tic; + loopcnt = loopcnt+1; - loopcnt = loopcnt+1; - progressFraction = loopcnt / looptatal; - waitbar(progressFraction, hWaitbar, ... - sprintf('Progress: %d/%d', loopcnt, looptatal)); + if M == 4 + fsym = 170e9; + elseif M == 6 + fsym = 144e9; + elseif M == 8 + fsym = 160e9; + end - loop_name = ['_fsym_',num2str(fsym)]; + %%%%% Loop Preps + %fsym = round(targetrate/log2(M)); + loop_name = ['_fsym_',num2str(fsym)]; + + %%%%% SET Voltages %%%%%% + dcs = DC_supply("active",[1,0],"voltage",[v_bias, 0]); + dcs.set("voltage",[v_bias, 0]); + + %%%%% SET Attenuator %%%%%% + voa = OptAtten("active",[1,2,1,1],"value",[rop_atten,pd_in_set,0,0],"wavelength",[1310,1310,1310,1310]); + voa.set('active',[1,2,1,1],'value',[rop_atten,pd_in_set,0,0]); + % voa.readvals(); + + %%%%% Construct AWG and Scope Modules %%%%%% + fdac = 256e9; + fadc = 256e9; + SCP = ScopeKeysight("model","UXR1104B",'autoscale',1,"fadc","GSa_256","channel",[0,1,0,0],"recordLen",2000000,"removeDC",1); + AWG = AwgKeysight("model","M8199B","fdac",fdac,"scaletodac",[1,1],"skews",[0,0],"voltages",[0,awg_vpp]); + A2S = Awg2Scope(AWG,SCP,[0,2,0,0],"waitUntilClick",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",random_key,... + "db_precode",db_precode,"db_encode",db_coding_approach,... + "mrds_code",usemrds,"mrds_blocklength",512).process(); + + %%%%% Precompensation Routine %%%%%% + if precomp_mode == 1 % measure channel + precomp_est = ChannelFreqResp("Nacq",1024,"Navg",64,"Ncp",63,'f_ref',fdac); + Digi_sig = precomp_est.buildOFDM(); + elseif precomp_mode == 2 % apply precomp + precomp_est = ChannelFreqResp("Nacq",1024,"Navg",64,"Ncp",63,'f_ref',Digi_sig.fs); + Digi_sig = precomp_est.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); + + + %%%%% Plot and Save Routine 1 %%%%%%%%%%%%%%%%%%%%%%%%% + %Digi_sig.spectrum("displayname","Normal Tx","fignum",10); + + save([folderpath,[experiment_name,'_bits'],loop_name],"Bits"); + save([folderpath,[experiment_name,'_symbols'],loop_name],"Symbols"); + %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% + + %%%%% AWG --> Scope %%%%%% + [~,Scpe_sig_raw,~,D] = A2S.process("signal2",Digi_sig,"waitUntilClick",0); + + % Scpe_sig_raw.spectrum("displayname","Scope PSD","fignum",20); + % Scpe_sig_raw.plot("displayname","Scope raw signal","fignum",25,"clear",1); + %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% + + %%%%%% Sample to 2x fsym %%%%%% + Scpe_sig_resampled = Scpe_sig_raw.resample("fs_in",fadc,"fs_out",2*fsym); + + %%%%% Precompensation Routine %%%%%% + if precomp_mode == 1 + precomp_est.estimate(Scpe_sig_resampled,"save",true,"savePath",precomp_path,"fileName",precomp_fn); + precomp_est.plot(); + end + + voa.readvals(); + rop = voa.power_state(1); + pd_in = voa.power_state(2); + disp(['ROP: ',num2str(rop),' dBm || PD in: ',num2str(pd_in), ' dBm']); + + %%%%%% Sync Rx signal with reference (S is a cell array with all occurences) %%%%%% + [Scpe_sig_syncd,S,isFlipped] = Scpe_sig_resampled.tsynch("reference",Symbols,"fs_ref",fsym); + + %%%%%% SNR CHEAT - Avges the measured signal occurences found after correlation in "tsynch" %%%%%% + 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_syncd.signal = scope_mean; + end + + %%%%% Plot and Save Routines: SAVE RECEIVED SIGNALS %%%%%%%%%%%%%%%%%%%%%%%%% + save([folderpath,experiment_name,'_rx_signal',loop_name],"S"); + + Scpe_sig_syncd.eye(fsym,M,"fignum",40,"displayname",' after Scope'); + + %%%%% 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.05,"DDmu",[0.0004 0.0004 0.0004 0.0004 ],"DFEmu",0.005,"FFEmu",0,"plotfinal",0,"ideal_dfe",1); + + % set to minus one not zero not avoid confusion if BER is acutally zero + ber_ffe = -1; + ber_mlse = -1; + ber_db = -1; + + if ffe_only %%%%%%%%%%%%%%%%%%%%%%%%%%% + + [EQ_sig] = Eq.process(Scpe_sig_syncd,Symbols); + + EQ_sig.plot("fignum",50,"displayname",'After EQ'); + + Rx_bits = PAMmapper(M,0).demap(EQ_sig); + + [~,errors_bm,ber_ffe,errors] = calc_ber(Rx_bits.signal,Bits.signal,"skip_front",100,"skip_end",150,"returnErrorLocation",1); + + disp(['FFE: ',sprintf('%.1E',ber_ffe),'| ROP: ',num2str(rop),' dB | PD_in: ',num2str(pd_in),' dBm']); + - %%%%% SET Voltages %%%%%% - dcs = DC_supply("active",[1,0],"voltage",[v_bias, 0]); - dcs.set("voltage",[v_bias, 0]); + if 1 - %%%%% SET Attenuator %%%%%% - voa = OptAtten("active",[1,2,1,1],"value",[rop_atten,pd_in_set,0,0],"wavelength",[1310,1310,1310,1310]); - voa.set('active',[1,2,1,1],'value',[rop_atten,pd_in_set,0,0]); - % voa.readvals(); + EQ_sig.plot("fignum",50,"displayname",'After EQ','clear',1); - %%%%% Construct AWG and Scope Modules %%%%%% - SCP = ScopeKeysight("model","UXR1104B",'autoscale',0,"fadc","GSa_256","channel",[0,1,0,0],"recordLen",2000000,"removeDC",1); - AWG = AwgKeysight("model","M8199B","fdac",fdac,"scaletodac",[1,1],"skews",[0,0],"voltages",[0,awg_vpp]); - A2S = Awg2Scope(AWG,SCP,[0,2,0,0]); + + figure(56); + clf + title(sprintf('PAM %d ; BER: %1.2e',M, ber_ffe)); + constellation = unique(Symbols.signal); + received = NaN(numel(constellation),length(Symbols)); + for lvl = 1:numel(constellation) + %Separate the equalized signal into the + %respective levels based on the actually + %transmitted level! + received(lvl,Symbols.signal==constellation(lvl)) = EQ_sig.signal(Symbols.signal==constellation(lvl)); + intermediate = received(lvl,:); + cnt(lvl) = numel(intermediate(~isnan(intermediate))); + hold on + histogram(received(lvl,:),1000,"EdgeAlpha",0,'DisplayName',['Lvl ',num2str(lvl),' | ',num2str(cnt(lvl)),' entries']); + end + legend + end + + elseif postfilter_approach %%%%%%%%%%%%%%%%%%%%%%%%%%% + + [EQ_sig] = Eq.process(Scpe_sig_syncd,Symbols); + + EQ_sig.plot("fignum",50,"displayname",'After EQ','clear',1); + + Noi = EQ_sig-Symbols; + + Rx_bits = PAMmapper(M,0).demap(EQ_sig); + + [~,num_errors,ber_ffe,pos_errors] = calc_ber(Rx_bits.signal,Bits.signal,"skip_front",100,"skip_end",150,"returnErrorLocation",1); - %%%%% 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",18,"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 - precomp_est = ChannelFreqResp("Nacq",1024,"Navg",64,"Ncp",63,'f_ref',fdac); - Digi_sig = precomp_est.buildOFDM(); - elseif precomp_mode == 2 % apply precomp - precomp_est = ChannelFreqResp("Nacq",1024,"Navg",64,"Ncp",63,'f_ref',Digi_sig.fs); - Digi_sig = precomp_est.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); - - - %%%%% Plot and Save Routine 1 %%%%%%%%%%%%%%%%%%%%%%%%% - % Digi_sig.spectrum("displayname","Normal Tx","fignum",10); - - % save([folderpath,[experiment_name,'bits'],loop_name],"Bits"); - % save([folderpath,[experiment_name,'symbols'],loop_name],"Symbols"); - %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% - - %%%%% AWG --> Scope %%%%%% - [~,Scpe_sig,~,D] = 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",fadc,"fs_out",2*fsym); - - %%%%% Precompensation Routine %%%%%% - if precomp_mode == 1 - precomp_est.estimate(Scpe_sig,"save",true,"savePath",precomp_path,"fileName",precomp_fn); - precomp_est.plot(); - end - - voa.readvals(); - rop = voa.power_state(1); - pd_in = voa.power_state(2); - disp(['ROP: ',num2str(rop),' dBm || PD in: ',num2str(pd_in), ' dBm']); - - %%%%%% 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(); - rop = voa.power_state(1); - pd_in = voa.power_state(2); - - %%%%% 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.05,"DDmu",[0.0004 0.0004 0.0004 0.0004 ],"DFEmu",0.005,"FFEmu",0,"plotfinal",0,"ideal_dfe",1); - - 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),'| ROP: ',num2str(rop),' dB | PD_in: ',num2str(pd_in),' dBm']); - - if 1 - 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); + + nc = 2; + burg_coeff = arburg(Noi.signal,nc); + + EQ_sig = EQ_sig.filter(burg_coeff,1); + + if 0 + 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 + + 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); + + [~,num_errors,ber_mlse,pos_errors] = calc_ber(Rx_bits.signal,Bits.signal,"skip_front",100,"skip_end",150,"returnErrorLocation",1); + + disp(['FFE: ',sprintf('%.1E',ber_ffe),' -> PF -> MLSE: ',sprintf('%.1E',ber_mlse),' dB | PD_in: ',num2str(pd_in),' dBm']); + + elseif db_channel_approach %%%%%%%%%%%%%%%%%%%%%%%%%%% + + [EQ_sig, Noi] = Eq.process(Scpe_sig_syncd,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); + [~,num_errors,ber_db,pos_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_db),' | PD_in: ',num2str(pd_in),' dBm']); + + elseif db_coding_approach %%%%%%%%%%%%%%%%%%%%%%%%%%% + + [EQ_sig, Noi] = Eq.process(Scpe_sig_syncd,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_db,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_db),' | PD_in: ',num2str(pd_in),' dBm']); + + end + + + %%%%% Store measurement into measurement "warehouse" %%%%%% + wh.addValueToStorage(ber_ffe,'ber_ffe',v_bias,awg_vpp,precomp_amp_max,rop_atten,M); + wh.addValueToStorage(ber_mlse,'ber_mlse',v_bias,awg_vpp,precomp_amp_max,rop_atten,M); + wh.addValueToStorage(ber_db,'ber_db',v_bias,awg_vpp,precomp_amp_max,rop_atten,M); + + wh.addValueToStorage(rop,'rop',v_bias,awg_vpp,precomp_amp_max,rop_atten,M); + wh.addValueToStorage(pd_in,'pd_in',v_bias,awg_vpp,precomp_amp_max,rop_atten,M); + % wh.addValueToStorage(Rx_bits,'signals',v_bias,awg_vpp,precomp_amp_max,rop_atten); + wh.addValueToStorage(M,'m',v_bias,awg_vpp,precomp_amp_max,rop_atten,M); + %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% + + %%%%% Plot stuff into Table (feel free to add own values in -> 'name',value <- notation. Must be closed when table is changed)%%%%%%%%%%%%%%%%%%%%%% + showCurrentMeasurement('FFE', ber_ffe,'MLSE',ber_mlse, 'Fsym',fsym.*1e-9, 'ROP', rop, 'PD in', pd_in, 'PAM',M, 'Vbias', v_bias, 'AWG Vpp', awg_vpp, 'Precomp MaxAmp',precomp_amp_max); + + %%%%% Arrange Figures %%%%%%%%%%%%%%%%%%%%%% + autoArrangeFigures(3,3,2); + + iterationTimes(loopcnt) = toc(iterationStartTime); + averageTimePerIteration = mean(iterationTimes(1:loopcnt)); + estimatedTotalTime = averageTimePerIteration * looptatal; + estimatedTimeRemaining = estimatedTotalTime - sum(iterationTimes(1:loopcnt)); + progressFraction = loopcnt / looptotal; + waitbar(progressFraction, hWaitbar, ... + sprintf('Loop: %d of %d \n Runtime: %.1f min | %.1f sec per Loop |Time to go: %.1f min ', ... + loopcnt, looptotal, sum(iterationTimes(1:loopcnt))/60, averageTimePerIteration, estimatedTimeRemaining/60 )); + + wh.save([folderpath,experiment_name,'_wh']); + end end - - elseif postfilter_approach %%%%%%%%%%%%%%%%%%%%%%%%%%% - - [EQ_sig] = Eq.process(Scpe_sig,Symbols); - - EQ_sig.plot("fignum",50,"displayname",'After EQ','clear',1); - - Noi = EQ_sig-Symbols; - - 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); - - nc = 3; - burg_coeff = arburg(Noi.signal,nc); - - EQ_sig = EQ_sig.filter(burg_coeff,1); - - if 0 - 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 - - 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']); - - 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,'ber',v_bias,awg_vpp); - wh.addValueToStorage(rop,'rop',v_bias,awg_vpp); - wh.addValueToStorage(pd_in,'pd_in',v_bias,awg_vpp); - wh.addValueToStorage(Rx_bits,'signals',v_bias,awg_vpp); - wh.addValueToStorage(M,'m',v_bias,awg_vpp); - - showCurrentMeasurement('BER', ber, 'ROP', rop, 'PD in', pd_in, 'PAM',M, 'Vbias', v_bias, 'AWG Vpp', awg_vpp, 'Precomp MaxAmp',precomp_amp_max); - autoArrangeFigures(3,3,2); - - end end + close(hWaitbar); wh.save([folderpath,experiment_name,'_wh']); autoArrangeFigures(3,3,2) + + + + +%%% ROP PLOT + +if 0 + + rop_vals = wh.parameter.rop_atten.values; + + ber_ffe = wh.getStoValue('ber_ffe',v_bias,awg_vpp,precomp_amp_max,rop_vals,M); + ber_mlse = wh.getStoValue('ber_mlse',v_bias,awg_vpp,precomp_amp_max,rop_vals,M); + rop_measured = wh.getStoValue('rop',v_bias,awg_vpp,precomp_amp_max,rop_vals,M); + pd_in_measured = wh.getStoValue('pd_in',v_bias,awg_vpp,precomp_amp_max,rop_vals,M); + + legendname = ['Thormax: ',num2str(fsym.*1e-9),' GBd | PAM',num2str(M),' | PD: ',num2str(pd_in_set),'| Vbias: ', num2str(v_bias),'V | ']; + + figure(230); + hold on; % Retain the plot so new points can be added without complete redraw + + % Plot the data and get the line handle + ffeLine = plot(rop_measured, ber_ffe, "LineWidth", 0.5, "LineStyle", "-", "Marker", ".", "MarkerSize", 15, "DisplayName", [legendname,' + FFE']); + mlseLine = plot(rop_measured, ber_mlse, "LineWidth", 0.5, "LineStyle", "-", "Marker", ".", "MarkerSize", 15, "DisplayName", [legendname,' +MLSE']); + + % Store pd_in_measured in the ZData property + ffeLine.ZData = pd_in_measured; + % Customize the data tips + % Set labels for existing data tip rows + ffeLine.DataTipTemplate.DataTipRows(1).Label = 'ROP'; + ffeLine.DataTipTemplate.DataTipRows(2).Label = 'FFE'; + ffeLine.DataTipTemplate.DataTipRows(2).Format = '%.2e'; % Format BER as "3e-4" + % Add a new data tip row for PDin + pdinRow = dataTipTextRow('PDin', 'ZData'); + ffeLine.DataTipTemplate.DataTipRows(3) = pdinRow; + + % Store pd_in_measured in the ZData property + mlseLine.ZData = pd_in_measured; + % Customize the data tips + % Set labels for existing data tip rows + mlseLine.DataTipTemplate.DataTipRows(1).Label = 'ROP'; + mlseLine.DataTipTemplate.DataTipRows(2).Label = 'MLSE'; + mlseLine.DataTipTemplate.DataTipRows(2).Format = '%.2e'; % Format BER as "3e-4" + % Add a new data tip row for PDin + pdinRow = dataTipTextRow('PDin', 'ZData'); + mlseLine.DataTipTemplate.DataTipRows(3) = pdinRow; + + % Continue with the rest of your plot settings + yline(3.8e-3, 'DisplayName', 'HD-FEC', 'LineStyle', '--', 'HandleVisibility', 'off'); + xlabel('Received Optical Power (dBm)'); + ylabel('Bit Error Rate (BER)'); + title('Bit Error Rate vs. ROP'); + set(gca, 'yscale', 'log'); + set(gca, 'Box', 'on'); + grid on; + grid minor; + legend('Interpreter', 'none'); + +end \ No newline at end of file diff --git a/projects/Lab_2024/bias_sweep_evaluation.m b/projects/Lab_2024/bias_sweep_evaluation.m index 46ebf62..b465e05 100644 --- a/projects/Lab_2024/bias_sweep_evaluation.m +++ b/projects/Lab_2024/bias_sweep_evaluation.m @@ -1,6 +1,8 @@ -filename = "C:\Users\Silas\Nextcloud4\Dokumente\02_Ablage_Office\Lab_Data_24\bias_sweep\PAM4_10km_ffe__wh.mat"; + +filename = "C:\Users\sioe\Documents\High_Speed_Measurement_2024\bias_testing\PAM4_b2b_bias_sweep_20241023_130342_wh.mat"; + a = load(filename); wh2 = a.obj; @@ -9,13 +11,19 @@ m = wh2.getStoValue('m',wh2.parameter.vbias.values(1),wh2.parameter.awg_vpp.valu v_bias_vals = wh2.parameter.vbias.values; awg_vpp_vals = wh2.parameter.awg_vpp.values; -bers = []; +ber_ffe= []; +ber_mlse= []; +rop_measured= []; +pd_in_measured= []; + rop_measured = []; cnt = 0; for awg_vpp_cur = awg_vpp_vals cnt = cnt+1; - bers(cnt,:) = wh2.getStoValue('ber',v_bias_vals,awg_vpp_cur); + ber_ffe(cnt,:) = wh2.getStoValue('ber_ffe',v_bias_vals,awg_vpp_cur)'; + ber_mlse(cnt,:) = wh2.getStoValue('ber_mlse',v_bias_vals,awg_vpp_cur); rop_measured(cnt,:) = wh2.getStoValue('rop',v_bias_vals,awg_vpp_cur); + pd_in_measured(cnt,:) = wh2.getStoValue('pd_in',v_bias_vals,awg_vpp_cur); end [bestber,bestindex] = min(bers,[],'all'); @@ -25,6 +33,23 @@ bestvbias=v_bias_vals(v_bias_pos); disp(['Best Vpp: ',num2str(bestvbias),' V; Best Vpp AWG: ',num2str(bestawgvpp),' V' ]) + +figure(100) +hold on +plot(v_bias_vals,ber_ffe,'DisplayName',['FFE only']); +plot(v_bias_vals,ber_mlse,'DisplayName',['MLSE']); +yline(3.8e-3, 'DisplayName', 'HD-FEC', 'LineStyle', '--', 'HandleVisibility', 'off'); +xlabel('V bias'); +ylabel('Bit Error Rate (BER)'); +title('Bit Error Rate vs. V bias'); +set(gca, 'yscale', 'log'); +set(gca, 'Box', 'on'); +grid on; +grid minor; +legend('Interpreter', 'none'); + + + figure(); sgtitle(['PAM ', num2str(m)]) subplot1 = subplot(1,2,1);