folderpath = 'C:\Users\sioe\Documents\High_Speed_Measurement_2024\8km_bitrate_rop_master\'; experiment_name = ''; currentTime = datetime('now', 'Format', 'yyyyMMdd_HHmmss'); timeStr = char(currentTime); experiment_name = [experiment_name, timeStr]; if 1 %%% BITRATE Sweep for MPI Experiment %%% awg_vpp = 2.7; pd_in_set = 8; random_key = 0; params = struct; params.M = [4,6,8]; params.lambda = flip([1293, 1302, 1310, 1318, 1327.4]); %calcWavelengthPlan(16, 400e9 , 1310); params.bitrate = [300:30:480].*1e9; params.duobinary = [0,1]; params.rop_atten = [0]; end if 0 %%% ROP SWEEP awg_vpp = 2.7; pd_in_set = 6; random_key = 0; params = struct; params.M = [8,6,4]; params.lambda = [1310]; %calcWavelengthPlan(16, 400e9 , 1310); params.bitrate = [300:30:480].*1e9; params.duobinary = [1,0]; params.rop_atten = [0:1.5:7.5]; end wh = DataStorage(params); wh.addStorage("ber_ffe"); wh.addStorage("ber_ffe_mlse"); wh.addStorage("ber_vnle"); wh.addStorage("ber_vnle_mlse"); wh.addStorage("ber_db"); wh.addStorage("FFE"); wh.addStorage("VNLE"); wh.addStorage("pd_in"); wh.addStorage("rop"); wh.addStorage("m"); wh.addStorage("dcs"); wh.addStorage("pdfa"); wh.addStorage("exfo"); wh.addStorage("voa"); 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 looptotal = prod(wh.dim); disp(['Start Measurement of ',num2str(looptotal),' loops...']) iterationTimes = zeros(looptotal, 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 loopcnt = 0; estimatedTimeRemaining = 0; estimatedTotalTime = 0; for M = wh.parameter.M.values dcs = DC_supply("active",[1,0],"voltage",[2.3, 0]); if M == 4 v_bias_for_pam = 2.3; dcs.set("voltage",[v_bias_for_pam, 0]); pulsef = 1; elseif M == 6 %pause(7*60); %wait 30 minutes for stable bias v_bias_for_pam=2.3; dcs.set("voltage",[v_bias_for_pam, 0]); pulsef = 0; elseif M == 8 v_bias_for_pam=2.6; dcs.set("voltage",[v_bias_for_pam, 0]); pause(7*60); %wait 30 minutes for stable bias pulsef = 0; end for lambda = wh.parameter.lambda.values exfo = Exfo_laser("serialport_number",'COM8','mainframe_channel',1,'safety_mode',0); pdfa = Thor_PDFA("safety_mode",0); exfo.getLaserInfo; if ~(exfo.cur_wavelength == lambda) % 1) pdfa.disablePDFA; % 2) exfo.setWavelength(lambda); % 3) pdfa.enablePDFA(); % 4) pdfa.setPumpLevel(100); end for bitrate = wh.parameter.bitrate.values fsym = floor( bitrate*1e-9./log2(M) ).*1e9; for db = wh.parameter.duobinary.values if db == 1 ffe_only = 0; postfilter_approach = 0; db_channel_approach = 1; db_coding_approach = 0; db_precode = db_coding_approach || db_channel_approach; if M == 4 pulsef=1; precomp_amp_max = -50; elseif M == 6 pulsef=0; precomp_amp_max = -50; elseif M == 8 pulsef=0; precomp_amp_max = -50; end elseif db == 0 ffe_only = 0; postfilter_approach = 1; db_channel_approach = 0; db_coding_approach = 0; db_precode = db_coding_approach || db_channel_approach; if M == 4 pulsef=1; precomp_amp_max = -38; elseif M == 6 pulsef=0; precomp_amp_max = -34; elseif M == 8 pulsef=0; precomp_amp_max = -34; end end %%%%% Construct AWG and Scope Modules %%%%%% fdac = 256e9; fadc = 256e9; SCP = ScopeKeysight("model","UXR1104B",'autoscale',1,"fadc","GSa_256","channel",[0,1,0,0],"recordLen",4000000,"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 %%%%%% rcalpha = 0.05; Pform = Pulseformer("fsym",fsym,"fdac",4*fsym,"pulse","rrc","pulselength",16,"rrcalpha",rcalpha); Pamsource = PAMsource(... "fsym",fsym,"M",M,"order",19,"useprbs",1,... "fs_out",fdac,... "applyclipping",0,"clipfactor",1.5,... "applypulseform",pulsef,"pulseformer",Pform,... "randkey",random_key,... "db_precode",db_precode,"db_encode",db_coding_approach,... "mrds_code",0,"mrds_blocklength",512); [Digi_sig,Symbols,Bits] = Pamsource.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); % Digi_sig.spectrum("displayname","TX After precomp","fignum",30,"normalizeToNyquist",0,"normalizeTo0dB",1); for rop_atten = wh.parameter.rop_atten.values %%%%% Loop Preps iterationStartTime = tic; loopcnt = loopcnt+1; loop_name = ['_PAM_',num2str(M),'_L_',num2str(lambda),'_R_',num2str(bitrate),'_DB_',num2str(db),'_ROP_',num2str(rop_atten)]; loop_name = strrep(loop_name,'.','_'); %%%%% READ Voltages %%%%%% dcs.readVals(); %%%%% SET Attenuator %%%%%% voa = OptAtten("active",[1,2,1,1],"value",[rop_atten,pd_in_set,0,0],"wavelength",[1310,1310,1310,1310],"speed",[1000,100,1000,1000]); voa.set('active',[1,2,1,1],'value',[rop_atten,pd_in_set,0,0]); voa.readvals(); %%%% SIGNAL USUALLY HERE, NOW ABOVE ROP_ATTEN %%% %%%%% Plot and Save Routine 1 - same for all rops, thus save only once %%%%%%%%%%%%%%%%%%%%%%%%% if rop_atten == 0 save([folderpath,experiment_name,loop_name,'_bits'],"Bits"); save([folderpath,experiment_name,loop_name,'_symbols'],"Symbols"); end %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%% AWG --> Scope %%%%%% [~,Scpe_sig_raw,~,D] = A2S.process("signal2",Digi_sig,"waitUntilClick",0); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % Scpe_sig_raw = Filter('filtdegree',5,"f_cutoff",0.55.*fsym,"fs",fadc,"filterType",filtertypes.gaussian,"active",true).process(Scpe_sig_raw); % % Scpe_sig_raw.plot("displayname","Scope raw signal","fignum",20,"clear",1); % Scpe_sig_raw.spectrum("displayname","Scope PSD","fignum",30,"normalizeTo0dB",1); % Scpe_sig_raw.eye(fsym,M,"displayname",'eye','fignum',200); %%%%%% 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); %%%%%% 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); %%%%% Plot and Save Routines: SAVE RECEIVED SIGNALS %%%%%%%%%%%%%%%%%%%%%%%%% save([folderpath,experiment_name,loop_name,'_rx_signal'],"S"); save([folderpath,experiment_name,loop_name,'_raw_signal'],"Scpe_sig_raw"); %%%%% EQUALIZE %%%%%% % set to minus one not zero not avoid confusion if BER is acutally zero ber_vnle = [-1]; ber_vnle_mlse = [-1]; ber_ffe_mlse =[-1]; ber_ffe = [-1]; ber_db = [-1]; ffe = EQ("Ne",[50,0,0],"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); vnle = 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 postfilter_approach %%%%%%%%%%%%%%%%%%%%%%%%%%% eq_values = min(numel(S),8); Noi = cell(eq_values,1); EQ_vnle= cell(eq_values,1); EQ_ffe= cell(eq_values,1); if 0 parfor s = 1:eq_values if 0 %FFE LINEAR Scpe_sig_syncd = S{s}; [EQ_ffe{s}] = ffe.process(Scpe_sig_syncd,Symbols); Noi{s} = EQ_ffe{s}-Symbols; Rx_bits = PAMmapper(M,0).demap(EQ_ffe{s}); [~,~,ber_ffe(s),~] = calc_ber(Rx_bits.signal,Bits.signal,"skip_front",100,"skip_end",150,"returnErrorLocation",1); end if 0 %FFE + MLSE nc = 2; burg_coeff = arburg(Noi{s}.signal,nc); EQ_ffe{s} = EQ_ffe{s}.filter(burg_coeff,1); EQ_mlse = MLSE("DIR",burg_coeff,"duobinary_output",0,"M",M,"trellis_states",PAMmapper(M,0).levels).process(EQ_ffe{s}); Rx_bits = PAMmapper(M,0).demap(EQ_mlse); [~,~,ber_ffe_mlse(s),~] = calc_ber(Rx_bits.signal,Bits.signal,"skip_front",100,"skip_end",150,"returnErrorLocation",1); end if 1 %VNLE Scpe_sig_syncd = S{s}; [EQ_vnle{s}] = vnle.process(Scpe_sig_syncd,Symbols); Noi{s} = EQ_vnle{s}-Symbols; Rx_bits = PAMmapper(M,0).demap(EQ_vnle{s}); [~,~,ber_vnle(s),~] = calc_ber(Rx_bits.signal,Bits.signal,"skip_front",100,"skip_end",150,"returnErrorLocation",1); end %VNLE + MLSE if 1 nc = 2; burg_coeff = arburg(Noi{s}.signal,nc); EQ_mlse = EQ_vnle{s}.filter(burg_coeff,1); EQ_mlse = MLSE("DIR",burg_coeff,"duobinary_output",0,"M",M,"trellis_states",PAMmapper(M,0).levels).process(EQ_mlse); Rx_bits = PAMmapper(M,0).demap(EQ_mlse); [~,~,ber_vnle_mlse(s),~] = calc_ber(Rx_bits.signal,Bits.signal,"skip_front",100,"skip_end",150,"returnErrorLocation",1); end end disp(['FFE EQ: BEST BER: ',sprintf('%.1E',min(ber_ffe)),' AVG BER: ',sprintf('%.1E',mean(ber_ffe)),' WORST:',sprintf('%.1E',max(ber_ffe)),'. Out of ',num2str(numel(ber_ffe))]); disp(['FFE + MLSE EQ: BEST BER: ',sprintf('%.1E',min(ber_ffe_mlse)),' AVG BER: ',sprintf('%.1E',mean(ber_ffe_mlse)),' WORST:',sprintf('%.1E',max(ber_ffe_mlse)),'. Out of ',num2str(numel(ber_ffe_mlse))]); disp(['VNLE EQ: BEST BER: ',sprintf('%.1E',min(ber_vnle)),' AVG BER: ',sprintf('%.1E',mean(ber_vnle)),' WORST:',sprintf('%.1E',max(ber_vnle)),'. Out of ',num2str(numel(ber_vnle))]); disp(['VNLE+MLSE EQ: BEST BER: ',sprintf('%.1E',min(ber_vnle_mlse)),' AVG BER: ',sprintf('%.1E',mean(ber_vnle_mlse)),' WORST:',sprintf('%.1E',max(ber_vnle_mlse)),'. Out of ',num2str(numel(ber_ffe))]); [~,i] = min(ber_vnle); figure(56); clf title(sprintf('PAM %d ; BER: %1.2e',M, ber_vnle(i))); 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_vnle{i}.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 db_channel_approach %%%%%%%%%%%%%%%%%%%%%%%%%%% ffe = EQ("Ne",[50,0,0],"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 0 parfor s = 1:numel(S) Scpe_sig_syncd = S{s}; [EQ_sig, Noi] = ffe.process(Scpe_sig_syncd,Duobinary().encode(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); [~,num_errors,ber_db(s),pos_errors] = calc_ber(Rx_bits.signal,Bits.signal,"skip_front",100,"skip_end",150,"returnErrorLocation",1); end disp(['DB EQ: BEST BER: ',sprintf('%.1E',min(ber_db)),' AVG BER: ',sprintf('%.1E',mean(ber_db)),' WORST:',sprintf('%.1E',max(ber_db)),'. Out of',num2str(numel(ber_db))]); else disp('Disabled MLSE for DB in all cases, due to time in measurement loop') end elseif db_coding_approach parfor s = 1:numel(S) S{s}=S{s}.normalize("mode","rms"); ffe = 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); [EQ_sig, Noi] = ffe.process(S{s},Symbols); EQ_sig.plot("displayname",'After EQ','fignum',112); 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(s),pos_errors] = calc_ber(Rx_bits.signal,Bits.signal,"skip_front",100,"skip_end",150,"returnErrorLocation",1); end disp(['DB EQ: BEST BER: ',sprintf('%.1E',min(ber_db)),' AVG BER: ',sprintf('%.1E',mean(ber_db)),' WORST:',sprintf('%.1E',max(ber_db)),'. Out of',num2str(numel(ber_db))]); end %%%%% Store measurement into measurement "warehouse" %%%%%% wh.addValueToStorage({ber_ffe},'ber_ffe',M,lambda,bitrate,db,rop_atten); wh.addValueToStorage({ber_ffe_mlse},'ber_ffe_mlse',M,lambda,bitrate,db,rop_atten); wh.addValueToStorage({ber_vnle},'ber_vnle',M,lambda,bitrate,db,rop_atten); wh.addValueToStorage({ber_vnle_mlse},'ber_vnle_mlse',M,lambda,bitrate,db,rop_atten); wh.addValueToStorage({ber_db},'ber_db',M,lambda,bitrate,db,rop_atten); wh.addValueToStorage(rop,'rop',M,lambda,bitrate,db,rop_atten); wh.addValueToStorage(pd_in,'pd_in',M,lambda,bitrate,db,rop_atten); wh.addValueToStorage(M,'m',M,lambda,bitrate,db,rop_atten); wh.addValueToStorage(dcs,'dcs',M,lambda,bitrate,db,rop_atten); wh.addValueToStorage(pdfa,'pdfa',M,lambda,bitrate,db,rop_atten); wh.addValueToStorage(exfo,'exfo',M,lambda,bitrate,db,rop_atten); wh.addValueToStorage(voa,'voa',M,lambda,bitrate,db,rop_atten); wh.addValueToStorage(ffe,'FFE',M,lambda,bitrate,db,rop_atten); wh.addValueToStorage(vnle,'VNLE',M,lambda,bitrate,db,rop_atten); iterationTimes(loopcnt) = toc(iterationStartTime); averageTimePerIteration = mean(iterationTimes(1:loopcnt)); estimatedTotalTime = averageTimePerIteration * looptotal; 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 end end end close(hWaitbar); wh.save([folderpath,experiment_name,'_wh']);