Begin with High Speed Measurement

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
+                    % SET Ref clock out (which is conn. to scope) to 10 MHz
-                    obj.writeNcheck(v, ':OUTPut:FREQuency ''M1.RefClkOut'',100000000');
+                    obj.writeNcheck(v, ':OUTPut:FREQuency ''M1.RefClkOut'',10000000');
 
                     % :OUTP:GLOB 'M1.System', OFF
                     obj.writeNcheck(v, ':OUTP:GLOB ''M1.System'', ON');

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
+                    case scope_model.UXR1104B
-                        v=visadev('TCPIP0::134.245.243.223::inst0::INSTR');
+                        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
 

Datatypes/scope_model.m

@@ -3,7 +3,7 @@ classdef scope_model < int32
     enumeration
         DSAZ634A (1)
         UXR1102A (2)
-        UXR1104A (3)
+        UXR1104B (3)
     end
 
 end

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();

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);

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)
+