From b7975a96e91fc54653911d2a9d71b2b135d7203b Mon Sep 17 00:00:00 2001 From: Silas Labor Zizou Date: Wed, 23 Oct 2024 09:32:10 +0200 Subject: [PATCH] Begin with High Speed Measurement --- Classes/05_Lab/AwgKeysight.m | 14 +- Classes/05_Lab/ScopeKeysight.m | 14 +- Datatypes/scope_model.m | 2 +- Tests/05_Lab/Scope_minmal_test.m | 6 + Tests/05_Lab/awg_minimal.m | 34 +++ .../bias_optimization.m | 258 ++++++++++++++++++ 6 files changed, 317 insertions(+), 11 deletions(-) create mode 100644 Tests/05_Lab/Scope_minmal_test.m create mode 100644 Tests/05_Lab/awg_minimal.m create mode 100644 projects/HighSpeedExperiment_2024/bias_optimization.m diff --git a/Classes/05_Lab/AwgKeysight.m b/Classes/05_Lab/AwgKeysight.m index 912eb10..1d35328 100644 --- a/Classes/05_Lab/AwgKeysight.m +++ b/Classes/05_Lab/AwgKeysight.m @@ -237,7 +237,7 @@ classdef AwgKeysight obj.writeNcheck(v,sprintf(':FREQuency:RASTer %.15g;', obj.fdac)); case awg_model.M8199B - obj.writeNcheck(v,':ABORt ''M2'''); %war auf M2, das wird wohl der CH sein oder? + obj.writeNcheck(v,':ABORt ''M2'''); obj.writeNcheck(v,sprintf(':FREQ ''M1.ClkGen'', %.15g;', obj.fdac)); end @@ -301,6 +301,14 @@ classdef AwgKeysight end + + else + % If no signal defined, turn off channel + switch obj.model + case awg_model.M8199B + obj.writeNcheck(v, sprintf(':OUTP ''M2.DataOut%d'', OFF',chan)); + end + end end @@ -313,8 +321,8 @@ classdef AwgKeysight case awg_model.M8199B - % SET Ref clock out (which is conn. to scope) to 100 MHz - obj.writeNcheck(v, ':OUTPut:FREQuency ''M1.RefClkOut'',100000000'); + % SET Ref clock out (which is conn. to scope) to 10 MHz + obj.writeNcheck(v, ':OUTPut:FREQuency ''M1.RefClkOut'',10000000'); % :OUTP:GLOB 'M1.System', OFF obj.writeNcheck(v, ':OUTP:GLOB ''M1.System'', ON'); diff --git a/Classes/05_Lab/ScopeKeysight.m b/Classes/05_Lab/ScopeKeysight.m index 5578a32..4565e69 100644 --- a/Classes/05_Lab/ScopeKeysight.m +++ b/Classes/05_Lab/ScopeKeysight.m @@ -67,8 +67,8 @@ classdef ScopeKeysight v=visadev('TCPIP0::keysight.tf.uni-kiel.de::inst0::INSTR');% TCPIP0::127.0.0.1::inst0::INSTR case scope_model.UXR1102A v=visadev('TCPIP0::testscope.tf.uni-kiel.de::inst0::INSTR');% TCPIP0::127.0.0.1::inst0::INSTR - case scope_model.UXR1104A - v=visadev('TCPIP0::134.245.243.223::inst0::INSTR'); + case scope_model.UXR1104B + v=visadev('TCPIP0::134.245.243.239::inst0::INSTR'); end else try @@ -79,7 +79,7 @@ classdef ScopeKeysight end end - debug = 0; + debug = 1; if debug disp(['Connected to Instrument: ',char(v.Vendor),' ',char(v.Model),' SerNo:',char(v.SerialNumber)]); end @@ -87,11 +87,11 @@ classdef ScopeKeysight % Check if Scope is ready to go opdone = 0; acqdone = 0; - procdone = 0; + procdone = 1; while ~opdone || ~acqdone || ~procdone opdone = sscanf(obj.writeReceiveCheck(v,'*opc?'), '%f'); acqdone = sscanf(obj.writeReceiveCheck(v,'ader?'), '%f'); - procdone = sscanf(obj.writeReceiveCheck(v,'pder?'), '%f'); + %procdone = sscanf(obj.writeReceiveCheck(v,'pder?'), '%f'); pause(0.005); end armed = sscanf(obj.writeReceiveCheck(v,'aer?'), '%f'); @@ -111,7 +111,7 @@ classdef ScopeKeysight end else obj.writeNcheck(v,':ACQuire:REDGe 0'); - obj.writeNcheck(v,':ACQuire:SRATe %u',double(obj.fadc).*1e9); + obj.writeNcheck(v,sprintf(':ACQuire:SRATe %u',double(obj.fadc).*1e9)); end % @@ -133,7 +133,7 @@ classdef ScopeKeysight switch obj.model case scope_model.DSAZ634A obj.writeNcheck(v,':TIMebase:REFClock HFR'); - case scope_model.UXR1102A || scope_model.UXR1104A + case scope_model.UXR1102A || scope_model.UXR1104B obj.writeNcheck(v,':TIMebase:REFClock 1'); end diff --git a/Datatypes/scope_model.m b/Datatypes/scope_model.m index d5b7587..4870979 100644 --- a/Datatypes/scope_model.m +++ b/Datatypes/scope_model.m @@ -3,7 +3,7 @@ classdef scope_model < int32 enumeration DSAZ634A (1) UXR1102A (2) - UXR1104A (3) + UXR1104B (3) end end \ No newline at end of file diff --git a/Tests/05_Lab/Scope_minmal_test.m b/Tests/05_Lab/Scope_minmal_test.m new file mode 100644 index 0000000..70318c8 --- /dev/null +++ b/Tests/05_Lab/Scope_minmal_test.m @@ -0,0 +1,6 @@ + + + + SCP = ScopeKeysight("model","UXR1104B",'autoscale',0,"fadc","GSa_256","channel",[1,0],"recordLen",2000000,"removeDC",1); + + scpe_sig_cell = SCP.read(); diff --git a/Tests/05_Lab/awg_minimal.m b/Tests/05_Lab/awg_minimal.m new file mode 100644 index 0000000..cfab03f --- /dev/null +++ b/Tests/05_Lab/awg_minimal.m @@ -0,0 +1,34 @@ + + +M = 4; +usemrds = 0; +fsym = 68e9; +fdac = 256e9; +awg_vpp = 0.1; +rrcalpha = 0.05; + +%%%%% Construct AWG and Scope Modules %%%%%% +SCP = ScopeKeysight("model","UXR1104B",'autoscale',0,"fadc","GSa_256","channel",[0,0,0,1],"recordLen",2000000,"removeDC",1); +AWG = AwgKeysight("model","M8199B","fdac",fdac,"scaletodac",[1,1],"skews",[0,0],"voltages",[awg_vpp,0]); +A2S = Awg2Scope(AWG,SCP,[0,0,0,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",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(); + +%%%%% 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); + + +%%%%% AWG --> Scope %%%%%% +[Scpe_sig,~,~,~] = A2S.process("signal1",Digi_sig); + +Scpe_sig.plot('displayname','Signal from Scope','fignum',10); +Scpe_sig.spectrum("displayname","Scope PSD","fignum",20); \ No newline at end of file diff --git a/projects/HighSpeedExperiment_2024/bias_optimization.m b/projects/HighSpeedExperiment_2024/bias_optimization.m new file mode 100644 index 0000000..fecc926 --- /dev/null +++ b/projects/HighSpeedExperiment_2024/bias_optimization.m @@ -0,0 +1,258 @@ + +folderpath = 'C:\Users\sioe\Documents\High_Speed_Measurement_2024\bias_testing\'; +experiment_name = 'PAM4_b2b_'; + +ffe_only = 0; +postfilter_approach = 1; +db_channel_approach = 0; +db_coding_approach = 0; + +db_precode = db_coding_approach || db_channel_approach; + +%%% SIR Sweep for MPI Experiment %%% +params = struct; + +params.vbias = [3]; +params.awg_vpp = [2.7]; + +wh = DataStorage(params); + +wh.addStorage("ber"); +wh.addStorage("pd_in"); +wh.addStorage("rop"); +wh.addStorage("m"); +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; + +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; + +looptatal = prod(wh.dim); +disp(['Start Measurement of ',num2str(looptatal),' loops...']) +hWaitbar = waitbar(0, 'Starting measurement...', 'Name', 'Processing Progress'); +loopcnt = 0; + +for v_bias = wh.parameter.vbias.values + for awg_vpp = wh.parameter.awg_vpp.values + + loopcnt = loopcnt+1; + progressFraction = loopcnt / looptatal; + waitbar(progressFraction, hWaitbar, ... + sprintf('Progress: %d/%d', loopcnt, looptatal)); + + 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 %%%%%% + 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]); + + %%%%% 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); + 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) +