diff --git a/Classes/01_transmit/ChannelFreqResp.m b/Classes/01_transmit/ChannelFreqResp.m index 592e932..5095539 100644 --- a/Classes/01_transmit/ChannelFreqResp.m +++ b/Classes/01_transmit/ChannelFreqResp.m @@ -192,7 +192,7 @@ classdef ChannelFreqResp < handle % iH(1) is DC ---> iH(end) is High Freq. H_inv = [iH(1) iH fliplr(conj(iH)) conj(iH(1))]; - H_inv = [iH(1) iH iH(end) fliplr(conj(iH)) conj(iH(1))]; + % H_inv = [iH(1) iH iH(end) fliplr(conj(iH)) conj(iH(1))]; obj.H_apply = H_inv; diff --git a/Classes/01_transmit/PAMsource.m b/Classes/01_transmit/PAMsource.m index e72033c..e370af2 100644 --- a/Classes/01_transmit/PAMsource.m +++ b/Classes/01_transmit/PAMsource.m @@ -57,7 +57,7 @@ classdef PAMsource end end - if boolean(obj.applypulseform) && isempty(obj.pulseformer) + if obj.applypulseform && isempty(obj.pulseformer) warning('No Pulseformer given. Proceeding with RRC and alpha 0.05'); options.pulseformer = Pulseformer("fsym",obj.fsym,"fdac",obj.fs_out,"pulse","rrc","pulselength",16,"rrcalpha",0.05); end diff --git a/Classes/05_Lab/AwgKeysight.m b/Classes/05_Lab/AwgKeysight.m new file mode 100644 index 0000000..c37c37f --- /dev/null +++ b/Classes/05_Lab/AwgKeysight.m @@ -0,0 +1,417 @@ +classdef AwgKeysight + + % AWG implementation only for Keysight devices: + % a) M8196A (92 GSa/s 4Ch) + % b) M8199A (128 GSa/s 4Ch) % TODO + % c) M8199A_ILV (256 GSa/s 2Ch) % TODO + % d) M8199B (256 GSa/s 2Ch) % Commands worked in 2023 + + % Summary of SCPI Commands: + % *IDN? — Identify the instrument. + % *RST — Reset the instrument. + % :SOUR:FREQ:RAST — Set the DAC sample frequency. + % :OUTP:STAT ON — Enable output for a channel. + % :VOLT:LEVel[:IMMediate][:AMPLitude] — Set output voltage. + % :TRAC:DATA — Load waveform data. + % :INIT:IMM — Start waveform generation. + + % Construct Class: + % AWG = AwgKeysight("model","M8196A","fdac",92e9,"scaletodac",[1,1,1,1],"skews",[0,0,0,0],"voltages",[0,0,0,0.6]); + + % Upload new Signal (use either normal array of doubles or Silas Informationsignal class) + % AWG.upload("Signal4",Signal); + + + properties(Access=protected) + model awg_model + fdac double + skews double + voltages double + scaletodac logical + + numChannels double + waveformGranularity double + sampleMemorySize double + + numProvidedSignals double + end + + methods (Access=public) + + function obj = AwgKeysight(options) + + + arguments + options.model awg_model + options.fdac double + options.skews double + options.voltages double + options.scaletodac logical + end + + % + fn = fieldnames(options); + for n = 1:numel(fn) + try + obj.(fn{n}) = options.(fn{n}); + end + end + + switch obj.model + case awg_model.M8196A + obj.numChannels = 4; + obj.waveformGranularity = 128; + obj.sampleMemorySize = 512000; + case awg_model.M8199A + obj.numChannels = 4; + obj.waveformGranularity = 512; + obj.sampleMemorySize = 1024000; %to check + case awg_model.M8199A_ILV + obj.numChannels = 2; + obj.waveformGranularity = 512; + obj.sampleMemorySize = 1024000; %to check + case awg_model.M8199B + obj.numChannels = 2; + obj.waveformGranularity = 512; + obj.sampleMemorySize = 1024000; %to check + end + + end + + function upload(obj,channels,options) + + arguments + obj + % leave this as it is! Important for further handling/ parsing + channels.signal1 Informationsignal = Informationsignal([]) + channels.signal2 Informationsignal = Informationsignal([]) + channels.signal3 Informationsignal = Informationsignal([]) + channels.signal4 Informationsignal = Informationsignal([]) + + % add new optional arguments here + options.bla = 1; + end + + % Every signal at this point is an Informationsignal! + % Argument block above casts everything to Informaionsignal... + + % Count the number of nonempty Informationsignal + obj.numProvidedSignals = sum(structfun(@(x) ~isempty(x.signal), channels)); + assert(obj.numProvidedSignals<=obj.numChannels,['Only ',num2str(obj.numChannels),' AWG Channels. User provided ',num2str(obj.numProvidedSignals),' waveforms!']) + + % Additional check for invalid signals beyond available channels + if obj.numChannels == 2 + % Check if Signal3 or Signal4 is nonempty + if ~isempty(channels.signal3.signal) || ~isempty(channels.signal4.signal) + error('Only 2 channels are available. Signals for Channel 3 and 4 should not be defined.'); + end + end + + % Extract the actual signals to cell array + unpackedSignals = {}; + fn = fieldnames(channels); + for s = 1:obj.numChannels + unpackedSignals{s} = channels.(fn{s}).signal; + end + + success = obj.upload_(unpackedSignals); + + end + end + + + + + + methods(Access=private) + + function success = upload_(obj,signal) + + %%%%%%%%%%%%%%%%%%%%%% + %%% PREPARE SIGNALS %% + %%%%%%%%%%%%%%%%%%%%%% + + for s = 1:obj.numChannels + + % assure laying row vector + signal{s} = reshape(signal{s},[1, length(signal{s})]); + + % set to zeros if voltage == 0 + if obj.voltages(s) == 0 + signal{s} = signal{s} .* 0; + end + + %scale to max DAC resolution + if obj.scaletodac(s) + factor = max(abs(signal{s})); + if factor~=0 + signal{s} = signal{s}./factor; + end + end + + % clip warning (if autoscale is off) + if ~isempty(signal{s}) + clip(s) = sum(signal{s}>1); + if clip(s) + warning(['Clipwarning CH ',num2str(s),'. Clipped Samples', num2str(clip(s))]); + end + end + + % The waveform granularity is e.g. 128. This means that the waveform length must be a multiple of this granularity. The minimum waveform length is 128 samples. + if rem(length(signal{s}),obj.waveformGranularity) ~= 0 + signal{s} = [signal{s} signal{s}(:,1:(obj.waveformGranularity-rem(length(signal{s}),obj.waveformGranularity)))]; + warning(['data should have a length multiple of ',num2str(obj.waveformGranularity),'!']); + end + + % Digitally apply skews to the signals + if obj.skews(s) > 0 + delay_sec = obj.skews(s)*1e-12; + signal{s} = delayseq(signal{s},delay_sec,obj.fdac); + end + + % Scale to DAC values + signal{s} = int8(round(127 * signal{s})); + + % Check signal length, truncate if too long for AWG + siglen = min(length(signal{s}),obj.sampleMemorySize); + signal{s} = signal{s}(1:siglen); + if siglen==obj.sampleMemorySize + warning(['Signal ',num2str(s),' was truncated to ',num2str(obj.sampleMemorySize),' to fit into AWG sample memory']) + end + + + end + + %%%%%%%%%%%%%%%%%%%%%% + %%% CONNECT TO AWG %%% + %%%%%%%%%%%%%%%%%%%%%% + + % read/write lock file to our ntserver/scratch + try + obj.check_lock(char(obj.model),1); + obj.lock_device(char(obj.model)); + + catch + warning('Could not reach ntserver to write lock!') + end + + % Open Visa + switch obj.model + case awg_model.M8196A + v = visadev('TCPIP0::zizou.tf.uni-kiel.de::hislip0::INSTR'); + case awg_model.M8196A + v = visadev('TCPIP0::localhost::hislip0::INSTR'); + end + + %check if channel config matches + writeline(v,'*opt?'); + aw=readline(v); + assert(obj.numChannels==sscanf(aw, '%f'),['Wrong Channel config. AWG response: ',num2str(sscanf(aw, '%f')), ' Channels']); + + switch obj.model + case awg_model.M8196A + + obj.writeNcheck(v,':INST:DACM FOUR'); + obj.writeNcheck(v,':ABORt'); + 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,sprintf(':FREQ ''M1.ClkGen'', %.15g;', obj.fdac)); + + end + + %%%%%%%%%%%%%%%%%%%%%% + %%% UPLOAD SIGNALS %%% + %%%%%%%%%%%%%%%%%%%%%% + + for chan = 1:obj.numChannels + + segm_num = 1; %we usually have one segment + segm_len = length(signal{chan});%there are max signl lengths + + if segm_len > 0 + + switch obj.model + case awg_model.M8196A + + % 1) TRAC DEl - delete the previous waveform + obj.writeNcheck(v,sprintf(':TRACe%d:DELete %d', chan, segm_num)); + + % 2) TRAC DEF - define properties of waveform + obj.writeNcheck(v,sprintf(':TRACe%d:DEFine %d,%d', chan, segm_num, segm_len)); + + % 3) TRAC DATA - load the new waveform + offset = 0; + cmd = sprintf(':TRACe%d:DATA %d,%d,', chan, segm_num, offset); + binblockwrite(v, signal{chan}, 'int8', cmd); + aw = obj.writeReceiveCheck(v,'*opc?'); + + % 4) SET VOLTAGE (not possible to set 0v) + obj.writeNcheck(v, sprintf(':VOLTage%d:AMPLitude %g', chan, obj.voltages(chan))); + + % 5) SET OUTPUTS ON/ OFF (generation not started at this point) + obj.writeNcheck(v, sprintf(':OUTPut%d ON', chan)); + + case awg_model.M8199B + + % 1) TRAC DEl - delete the previous waveform + xfprintf(f, sprintf(':TRAC:DEL ''M2.DataOut1'', %d', segm_num), 1); + xfprintf(f, sprintf(':TRAC:DEL ''M2.DataOut2'', %d', segm_num), 1); + + % 2) TRAC DEF - define properties of waveform + xfprintf(f, sprintf(':TRAC:DEF ''M2.DataOut1'', %d,%d', segm_num, len)); + xfprintf(f, sprintf(':TRAC:DEF ''M2.DataOut2'', %d,%d', segm_num, len)); + + % 3) TRAC DATA - load the new waveform + + % Example command: :TRAC:DATA 'M1.DataOut1', 1, 0 , -127,-126,,-1,0,1,,127,-127,,-1,0,1,,127 + % means : :TRAC:DATA 'identifier' , , , | + + cmd = sprintf(':TRAC:DATA ''M2.DataOut1'', %d,%d,', segm_num, offset); + binblockwrite(f, data(1+offset:offset+len), dataFormat, cmd); + + cmd = sprintf(':TRAC:DATA ''M2.DataOut2'', %d,%d,', segm_num, offset); + binblockwrite(f, data(1+offset:offset+len), dataFormat, cmd); + + % 4) SET VOLTAGE (not possible to set 0v) + xfprintf(f, sprintf(':VOLT:AMPL ''M2.DataOut1'', %gv', awg_volt(1))); + xfprintf(f, sprintf(':VOLT:AMPL ''M2.DataOut2'', %gv', awg_volt(2))); + % xfprintf(f, ':VOLT:AMPL ''M2.DataOut2'', 0.1v'); + + % 5) SET OUTPUTS ON/ OFF + xfprintf(f, ':OUTP ''M2.DataOut1'', ON'); %sioe 11.23: Hier habe ich absichtlich auf OFF gestellt, wenn Kanal nicht bentigt! + xfprintf(f, ':OUTP ''M2.DataOut2'', ON'); + + + + end + end + + end + + % Set global output on/ off to start generation + switch obj.model + case awg_model.M8196A + + obj.writeNcheck(v,':INIT:IMMediate'); + + case awg_model.M8199B + + % :OUTP:GLOB 'M1.System', OFF + obj.writeNcheck(v, ':OUTP:GLOB ''M1.System'', ON'); + obj.writeNcheck(v, ':OUTP:GLOB ''M2.System'', ON'); + obj.writeNcheck(v,':INIT:IMMediate ''M2'''); + end + + % operation ceomplete (opc) + obj.writeNcheck(v,'*OPC'); + aw = obj.writeReceiveCheck(v,'*opc?'); + + % close visa connection + delete(v); + + % delete lock file to our ntserver/scratch + try + obj.release_lock(char(obj.model)); + catch + warning('Could not reach ntserver to write lock!') + end + + % at this point the signal upload should be successful :-) + success = 1; + + end + + function writeNcheck(~,visaobj,command) + % 1) send command to intstument + % 2) check if everything is okay + + writeline(visaobj,char(command)); + + writeline(visaobj,':syst:err?'); + err_aw=readline(visaobj); + assert(strfind(err_aw,"No error"),['AWG System Error? after command: ', char(command)]) + + end + + function [aw]=writeReceiveCheck(~,visaobj,command) + % 1) send command to intstument + % 2) receive answer from instrument + % 3) check if everything is okay + + writeline(visaobj,char(command)); + aw = readline(visaobj); + + writeline(visaobj,':syst:err?'); + err_aw=readline(visaobj); + assert(strfind(err_aw,"No error"),['AWG System Error? after command: ', char(command)]) + + end + + function check_lock(~,aDeviceName,lock_timeout) + pause(0.2); + rdy = 0; + for n = 1:(lock_timeout*100) + pause(0.1); + if rem(n,100) == 0 + display(['waiting for device lock! (' num2str(n/10) 's)']); + + button = questdlg(['Lab Measurement is LOCKED! ',... + 'be careful'],... + 'LOCKED','Wait','Unlock','Abort Meas','Wait') ; + + switch button + case {'Abort Meas' ,''} + delete(findobj('name','Simulation Status Window')); + error('Measurement aborted by the user') ; + case 'Unlock' + rdy = 1; + break; + case 'Wait' + %wait another ten secs + end + + end + % file_handle = fopen(['\\ntserver\scratch\labor\lock\' aDeviceName '.lock'],'r'); + file_handle = fopen([filesep, filesep, 'ntserver.tf.uni-kiel.de', filesep, 'scratch', filesep 'Labor', filesep, 'lock', filesep, aDeviceName, '.lock'],'r'); + if file_handle == -1 + rdy = 1; + break; + end + time_flag=datenum(clock); + dev_flag = fread(file_handle,time_flag,'double'); + fclose(file_handle); + + if time_flag > dev_flag+60/(24*60) + rdy = 1; + break; + end + + end + + if ~rdy + + errordlg([aDeviceName ' is locked by someone else!!!']); + errormsg('AWG2SCOPE',[aDeviceName ' is locked by someone else!!!']); + else + end + end + + function lock_device(~,aDeviceName) + file_handle = fopen([filesep, filesep, 'ntserver.tf.uni-kiel.de', filesep, 'scratch', filesep 'Labor', filesep, 'lock', filesep, aDeviceName, '.lock'],'wb'); + + if ~isempty(file_handle) + time_flag=datenum(clock); + fwrite(file_handle,time_flag,'double'); + fclose(file_handle); + end + end + + function release_lock(~,aDeviceName) + pause(0.5); + delete([filesep, filesep, 'ntserver.tf.uni-kiel.de', filesep, 'scratch', filesep 'Labor', filesep, 'lock', filesep, aDeviceName '.lock']); + end + + end +end diff --git a/Classes/05_Lab/DC_supply.m b/Classes/05_Lab/DC_supply.m index 4cb1377..a5a5288 100644 --- a/Classes/05_Lab/DC_supply.m +++ b/Classes/05_Lab/DC_supply.m @@ -1,6 +1,9 @@ classdef DC_supply % Simple control of DC supply via fixed GBIB address. % Agilent E3646A + % No OVP or OCP implementation + % Not tested what happens if this script asks for Voltage in higher V + % Range (8V vs. 22V) properties(Access=public) active @@ -10,10 +13,10 @@ classdef DC_supply methods (Access=public) function obj = DC_supply(options) - + arguments - options.active bool = true + options.active logical = true options.voltage double = [0,0]; end @@ -30,35 +33,40 @@ classdef DC_supply end - function success = set(obj) + function success = set(obj,options) - success = false; + + arguments + obj + options.voltage double = obj.voltage; + end + + success = [0,0]; try %connect to device - v = visadev("GPIB0::19::INSTR"); - % fopen(v); + v = visadev("GPIB1::19::INSTR"); writeline(v, '*IDN?;'); disp(['DC Source: ' readline(v)]); cmd = 'INST:SEL?'; - fprintf(v, cmd); - selected_channel = fscanf(v); + writeline(v, cmd); + prev_selected_channel = readline(v); for ch = 1:2 % choose channel cmd = ['INST:SEL OUT',num2str(ch)]; - fprintf(v, cmd); + writeline(v, cmd); % get current voltage level - cmd = ['VOLT?']; - fprintf(v, cmd); - act_volt = str2num(fscanf(v)); + cmd = 'VOLT?'; + writeline(v, cmd); + act_volt = str2num(readline(v)); % desired voltage (round to two digits after comma) - des_volt = round(obj.voltage(ch),2); + des_volt = round(options.voltage(ch),2); cnt = 0; while abs(act_volt-des_volt) > 0 @@ -66,9 +74,9 @@ classdef DC_supply selected_channel = ['OUT',num2str(ch)]; % get current voltage level - cmd = ['VOLT?']; - fprintf(v, cmd); - act_volt = str2num(fscanf(v)); + cmd = 'VOLT?'; + writeline(v, cmd); + act_volt = str2num(readline(v)); % difference diff_volt = act_volt-des_volt; @@ -76,7 +84,7 @@ classdef DC_supply % set new voltage increment_voltage = -0.01* sign(diff_volt) ; cmd = ['VOLT ',num2str(act_volt+increment_voltage)]; - fprintf(v, cmd); + writeline(v, cmd); cnt = cnt+1; if mod(cnt,100) == 0 @@ -84,9 +92,9 @@ classdef DC_supply end % get current voltage level - cmd = ['VOLT?']; - fprintf(v, cmd); - act_volt = str2num(fscanf(v)); + cmd = 'VOLT?'; + writeline(v, cmd); + act_volt = str2num(readline(v)); end @@ -94,31 +102,21 @@ classdef DC_supply 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 end % choose channel - cmd = ['INST:SEL ',selected_channel]; - fprintf(v, cmd); + cmd = ['INST:SEL ',prev_selected_channel]; + writeline(v, cmd); %disconnect - fclose(v); delete(v); catch - - end - - end - - end - - methods (Access=private) - - - end end diff --git a/Classes/05_Lab/OptAtten.m b/Classes/05_Lab/OptAtten.m new file mode 100644 index 0000000..0356084 --- /dev/null +++ b/Classes/05_Lab/OptAtten.m @@ -0,0 +1,287 @@ +classdef OptAtten + % control optical attenuator Keysight N7764a via LAN + % set parameters of optical attenuator, e.g. attenuation, wavelength, etc. and read out output power from attenuator + + properties(Access=public) + active %'Enable/Disable attenuator slots. Array position indexes detemines attenuator slot: [1->2 3->4 5->6 7->8]. Value determines mode: 0: off; 1: constant attenuation value; 2: constant output power'; + value %'Set attenuation [dB] or constant output power [dBm]. Range attenuation [0, 45], range constant output power [-50, 20]'; + speed %'Enter a desired speed for attenuation changes [dB/s]. Range: [0.1, 1000]'; + wavelength %'Enter the wavelength of the applied light [nm]. Range: [1260, 1640]'; + avgTime %'Enter a desired averaging time [ms]. Higher averaging times will increase stability of the readings but will slow down measurement speed. Range: [2, 1000]'; + dbOffset %'A non-zero attenuation offset shifts the logarithmic attenuation scale. Then, the set attenuation value differs from the internally calibrated attenuation. Range [-200, 200]'; + psetOffset %'A non-zero power offset shifts the logarithnic power scale. Then, the power reading and the set power (Pset) values differ from the internal power calibration. Range [-60, 60]'; + end + + methods (Access=public) + function obj = OptAtten(options) + + arguments + options.active double = [0,0,0,0]; + options.value double = [0,0,0,0]; + options.speed double = [1000 1000 1000 1000]; + options.wavelength double = [1550 1550 1550 1550]; + options.avgTime double = [2 2 2 2]; + options.dbOffset double = [0 0 0 0]; + options.psetOffset double = [0 0 0 0]; + end + + fn = fieldnames(options); + for n = 1:numel(fn) + try + obj.(fn{n}) = options.(fn{n}); + end + end + + end + + function success = set(obj,options) + + arguments + obj + options.active = obj.active; %if not used as arg of .set('active',[0,0,0,0]) method the class values are used... + options.value = obj.value; + end + + try + %connect to device + v = visadev('TCPIP::134.245.243.248::INSTR'); + + writeline(v, '*IDN?;'); + disp(['Active atten: ' readline(v)]); %Keysight Technologies, N7764A, MY49A00696, 1.13.1 + + + slot = 0; %init slot number + %loop to set values for every attenuator slot + for i = 1:length(obj.active) + %convert index number i into slot number + slot = 2*i - 1; + + %check whether attenuator is used or unused + if ~obj.active(i) %atten slot unused + writeline(v, [':OUTP' num2str(slot) ':STAT OFF']); + + else %atten slot used (attenuation/const output power) + + %set attenuation or output power + if obj.active(i) == 1 %attenuation + if obj.value(i) < 0 + obj.value(i) = 0; + hMsgBox = msgbox(['Value for attenuation out of range. Please consider minimum value. Attenuation set to ' num2str(obj.value(i)) 'dB']); + uiwait(hMsgBox); + elseif obj.value(i) > 45 + obj.value(i) = 45; + hMsgBox = msgbox(['Value for attenuation out of range. Please consider maximum value. Attenuation set to ' num2str(obj.value(i)) 'dB']); + uiwait(hMsgBox); + end + writeline(v, [':OUTP' num2str(slot) ':POW:CONTR OFF']); + writeline(v, [':INP' num2str(slot) ':ATT ' num2str(obj.value(i)) 'dB']); + elseif obj.active(i) == 2 %constant output power + if obj.value(i) < -50 + obj.value(i) = -50; + hMsgBox = msgbox(['Value for output power out of range. Please consider minimum value. Output power set to ' num2str(obj.value(i)) 'dBm']); + uiwait(hMsgBox); + elseif obj.value(i) > 20 + obj.value(i) = 20; + hMsgBox = msgbox(['Value for output power out of range. Please consider maximum value. Output set to ' num2str(obj.value(i)) 'dBm']); + uiwait(hMsgBox); + end + writeline(v, [':OUTP' num2str(slot) ':POW:CONTR ON']); + writeline(v, [':OUTP' num2str(slot) ':POW ' num2str(obj.value(i))]); + end + + %set value for speed + if obj.speed(i) < 0.1 + obj.speed(i) = 0.1; + hMsgBox = msgbox(['Value for speed out of range. Please consider minimum value. Speed set to ' num2str(obj.speed(i)) 'dB/s']); + uiwait(hMsgBox); + elseif obj.speed(i) > 1000 + obj.speed(i) = 1000; + hMsgBox = msgbox(['Value for speed out of range. Please consider maximum value. Speed set to ' num2str(obj.speed(i)) 'dB/s']); + uiwait(hMsgBox); + end + writeline(v, [':INP' num2str(slot) ':ATT:SPE ' num2str(obj.speed(i))]); + + %set value for wavelength + if obj.wavelength(i) < 1260 + obj.wavelength(i) = 1260; + hMsgBox = msgbox(['Value for wavelength out of range. Please consider minimum value. Wavelength set to ' num2str(obj.wavelength(i)) 'nm']); + uiwait(hMsgBox); + elseif obj.wavelength(i) > 1640 + obj.wavelength(i) = 1640; + hMsgBox = msgbox(['Value for wavelength out of range. Please consider maximum value. Wavelength set to ' num2str(obj.wavelength(i)) 'nm']); + uiwait(hMsgBox); + end + wavelength_nm = obj.wavelength(i) * 1e-9; + writeline(v, [':INP' num2str(slot) ':WAV ' num2str(wavelength_nm)]); + + %set value for averaging times + if obj.avgTime(i) < 2 + obj.avgTime(i) = 2; + hMsgBox = msgbox(['Value for average time out of range. Please consider minimum value. Average time set to ' num2str(obj.avgTime(i)) 'ms']); + uiwait(hMsgBox); + elseif obj.avgTime(i) > 1000 + obj.avgTime(i) = 1000; + hMsgBox = msgbox(['Value for average time out of range. Please consider maximum value. Average time set to ' num2str(obj.avgTime(i)) 'ms']); + uiwait(hMsgBox); + end + avgTime_sec = obj.avgTime(i) / 1000; + writeline(v, [':OUTP' num2str(slot) ':ATIM ' num2str(avgTime_sec)]); + + %set value for attenuation offset and output power offset + if obj.dbOffset(i) < -200 + obj.dbOffset(i) = -200; + hMsgBox = msgbox(['Value for attenuation offset out of range. Please consider minimum value. Attenuation offset set to ' num2str(obj.dbOffset(i)) 'dB']); + uiwait(hMsgBox); + elseif obj.dbOffset(i) > 200 + obj.dbOffset(i) = 200; + hMsgBox = msgbox(['Value for attenuation out of range. Please consider maximum value. Attenuation offset set to ' num2str(obj.dbOffset(i)) 'dB']); + uiwait(hMsgBox); + end + writeline(v, [':INP' num2str(slot) ':OFFS ' num2str(obj.dbOffset(i)) 'dB']); + + if obj.psetOffset < -60 + obj.psetOffset = -60; + hMsgBox = msgbox(['Value for output power offset out of range. Please consider minimum value. Output power offset set to ' num2str(obj.psetOffset(i)) 'dBm']); + uiwait(hMsgBox); + elseif obj.psetOffset > 60 + obj.psetOffset = 60; + hMsgBox = msgbox(['Value for output power offset out of range. Please consider maximum value. Output power offset set to ' num2str(obj.psetOffset(i)) 'dBm']); + uiwait(hMsgBox); + end + writeline(v, [':OUTP' num2str(slot) ':POW:OFFS ' num2str(obj.psetOffset(i))]); + + %turn on attenuator slot + writeline(v, [':OUTP' num2str(slot) ':STAT ON']); + + end + end + + %pause 2 seconds before measuring output power + pause(2); + + %read output power from attenuator and save to state.ouputpower + for i = 1:length(obj.active) + if obj.active(i) + + %convert index number i into slot number + slot = 2*i - 1; + + %read and save + writeline(v, [':READ' num2str(slot) ':POW?']); + state.outputpower(i) = sscanf(readline(v),'%f'); + end + end + + %create warning message when measured output power differs from + %set output power by 2dB. + for i = 1:length(obj.active) + if obj.active(i) == 2 + + %determine absolute difference + differ = abs(state.outputpower(i)-obj.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); + end + end + end + + %create output + for i = 1:length(obj.active) + if obj.active(i) == 2 + % Output power mode: + %differ = abs(state.outputpower(i)-obj.value(i)); + ch_info_txt = ['OptAtten: ',num2str(i),' -> Desired: ',num2str(obj.value(i)),' dBm; Cur Output: ',num2str(state.outputpower(i)),'dBm']; + disp(ch_info_txt); + else + ch_info_txt = ['OptAtten: ',num2str(i),' -> Attenuated by: ',num2str(obj.value(i)),' dB; Cur Output: ',num2str(state.outputpower(i)),'dBm']; + disp(ch_info_txt); + end + + end + + %disconnect + delete(v); + + catch error + + %disconnect if error occurrs + delete(v); + hMsgBox = msgbox('An error occurred while connecting or writing to the attenuator. Please try again. Otherwise, please restart MATLAB. For more information see error message in command window.'); + uiwait(hMsgBox); + rethrow(error); + + end + + end + + function readvals(obj) + + try + %connect to device + v = visadev('TCPIP::134.245.243.248::INSTR'); + + writeline(v, '*IDN?;'); + disp(['Active atten: ' readline(v)]); %Keysight Technologies, N7764A, MY49A00696, 1.13.1 + + power=[]; + + for s = 1:2:7 + + writeline(v, [':OUTP' num2str(s) ':STAT?']); + line = readline(v); + answer = sscanf(line,'%f'); + disp(['Slot ' num2str(s) ' enabled state: ' num2str(answer)]); + + writeline(v, [':OUTP' num2str(s) ':POW?']); + line = readline(v); + answer = sscanf(line,'%f'); + disp(['P_set (output power) for Slot ' num2str(s) ': ' num2str(answer) ' dBm']); + + writeline(v, [':OUTP' num2str(s) ':POW:CONTR?']); + line = readline(v); + answer = sscanf(line,'%f'); + disp(['Power control for Slot ' num2str(s) ': ' num2str(answer)]); + + if answer == 1 + writeline(v, [':READ' num2str(s) ':POW?']); + line = readline(v); + answer = sscanf(line,'%f'); + disp([' -> Slot ' num2str(s) ' desired outut power: ' num2str(answer) ' dBm']); + else + writeline(v, [':INP' num2str(s) ':ATT?']); + line = readline(v); + answer = sscanf(line,'%f'); + disp([' -> Alpha (fix attenuation) for slot ' num2str(s) ': ' num2str(answer) ' dB']); + end + + writeline(v, [':INP' num2str(s) ':WAV?']); + line = readline(v); + answer = sscanf(line,'%f'); + disp(['Wavelength for Slot ' num2str(s) ': ' num2str(answer) ' nm']); + + writeline(v, [':OUTP' num2str(s) ':ATIM?']); + line = readline(v); + answer = sscanf(line,'%f'); + disp(['Averaging time for Slot ' num2str(s) ': ' num2str(answer) ' s']); + + fprintf('\n') + end + + delete(v); + + catch error + + %disconnect if error occurrs + delete(v); + hMsgBox = msgbox('An error occurred while connecting or writing to the attenuator. Please try again. Otherwise, please restart MATLAB. For more information see error message in command window.'); + uiwait(hMsgBox); + rethrow(error); + + end + end + + end + +end diff --git a/Classes/05_Lab/ScopeKeysight.m b/Classes/05_Lab/ScopeKeysight.m new file mode 100644 index 0000000..31053b9 --- /dev/null +++ b/Classes/05_Lab/ScopeKeysight.m @@ -0,0 +1,279 @@ +classdef ScopeKeysight + %NAME Summary of this class goes here + % Detailed explanation goes here + + properties(Access=public) + model scope_model + fadc scope_fadc + channel + autoScale + extRef + interpolate + recordLen + IPaddress + end + + methods (Access=public) + + function obj = ScopeKeysight(options) + %NAME Construct an instance of this class + % Detailed explanation goes here + + arguments + options.model scope_model = scope_model.DSAZ634A; + options.fadc scope_fadc = scope_fadc.GSa_160; + options.channel logical = [0,0,0,0]; + options.autoScale logical = 0; + options.extRef logical = 1; + options.interpolate logical = 0; + options.recordLen double = 1000000; + options.IPaddress + end + + % + fn = fieldnames(options); + for n = 1:numel(fn) + try + obj.(fn{n}) = options.(fn{n}); + end + end + + switch obj.model + case scope_model.DSAZ634A + + case scope_model.UXR1102A + + end + + + end + + function recordedSignals = read(obj,options) + + arguments + obj + options.channel logical = obj.channel + end + + %%%%%%%%%%%%%%%%%%%%%%%% + %%% CONNECT TO SCOPE %%% + %%%%%%%%%%%%%%%%%%%%%%%% + + if isempty(obj.IPaddress) + switch obj.model + case scope_model.DSAZ634A + 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'); + end + else + try + ressourceName=['TCPIP0::',char(obj.IPaddress),'::inst0::INSTR']; + v=visadev(ressourceName); + catch + error(['Could not connect to instrument using VISA. Check VISA ressourcename -> ', ressourceName]); + end + end + + disp(['Connected to Instrument: ',char(v.Vendor),' ',char(v.Model),' SerNo:',char(v.SerialNumber)]); + + % Check if Scope is ready to go + opdone = 0; + acqdone = 0; + procdone = 0; + 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'); + pause(0.005); + end + armed = sscanf(obj.writeReceiveCheck(v,'aer?'), '%f'); + scope_state = strtrim(obj.writeReceiveCheck(v,'RSTate?')); + + % Set sample rate / activate real edge + if obj.fadc == scope_fadc.GSa_160 + obj.writeNcheck(v,':ACQuire:REDGe 1'); + if length(obj.channel)>1 + if obj.channel(2) == 1 + error("CH 2 is invalid for Real Edge configuration. ONLY use channel 1 and 3 as shown on scope! Simply use [1,0,1,0] as channel config to record 1 and 3.") + end + elseif length(obj.channel)>3 + if obj.channel(2) == 1 || obj.channel(4) == 1 + error("CH 2 is invalid for Real Edge configuration. ONLY use channel 1 and 3 as shown on scope! Simply use [1,0,1,0] as channel config to record 1 and 3.") + end + end + else + obj.writeNcheck(v,':ACQuire:REDGe 0'); + obj.writeNcheck(v,':ACQuire:SRATe %u',double(obj.fadc).*1e9); + end + + % + for n = find(obj.channel == 1) + obj.writeNcheck(v,sprintf(':CHANnel%u:DISPlay ON',n));% display captured data trace + end + + 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)); + end + else + obj.writeNcheck(v,':SINGLE'); + end + + if obj.extRef + switch obj.model + case scope_model.DSAZ634A + obj.writeNcheck(v,':TIMebase:REFClock HFR'); + case scope_model.UXR1102A || scope_model.UXR1104A + obj.writeNcheck(v,':TIMebase:REFClock 1'); + end + + REFclock = strtrim(obj.writeReceiveCheck(v,':TIMebase:REFClock?')); + + if ~strcmp(REFclock,'HFR') + obj.writeNcheck(v,':TIMebase:REFClock ON'); + REFclock = strtrim(obj.writeReceiveCheck(v,':TIMebase:REFClock?')); + assert(~(strcmp(REFclock,'OFF') || strcmp(REFclock,'0')),'No Refclock found: Using internal') + end + + else + obj.writeNcheck(v,':TIMebase:REFClock 0'); + end + + if obj.interpolate + obj.writeNcheck(v,'ACQ:INTerpolate INT16'); + obj.writeNcheck(v,sprintf(':ACQuire:POINts %u',obj.recordLen+100)); + else + obj.writeNcheck(v,sprintf(':ACQuire:POINts %u',obj.recordLen+100)); + obj.writeNcheck(v,'ACQ:INTerpolate OFF'); + end + + obj.writeNcheck(v,':SYSTem:HEADer OFF'); % turn off system headers + obj.writeNcheck(v,':WAVEFORM:FORMAT WORD');% set the data format type to WORD + obj.writeNcheck(v,':WAVEFORM:BYTeorder LSBFirst');% Setup transfer of LSB first + obj.writeNcheck(v,':WAVEFORM:STReaming off');% Turn off waveform streaming + + for n = find(obj.channel == 1) + obj.writeNcheck(v,sprintf(':CHANnel%u:DISPlay ON',n));% display captured data trace + end + + %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% + %%% Acquire Data from Scope %%% + %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% + + obj.writeNcheck(v,':DIGitize'); + + recordedSignals = cell(size(obj.channel)); + for n = find(obj.channel == 1) + recordedSignals{n} = obj.readChannel(v,n); + recordedSignals{n} = recordedSignals{n}(1:obj.recordLen); + obj.writeNcheck(v,sprintf(':CHANnel%u:DISPlay ON',n)); + end + + for n = find(obj.channel == 0) + recordedSignals{n} = []; + end + + obj.writeNcheck(v,sprintf(':%s',scope_state)); + obj.writeNcheck(v,'*OPC?'); + + end + + end + + methods (Access=private) + + function captured_signal = readChannel(obj, v, aCH) + + Format = 'word'; + if strcmpi(Format,'word') + nBytes = 2; + blockReadType = 'int16'; + elseif strcmpi(Format,'byte') + nBytes = 1; + blockReadType = 'int8'; + else + error('Param.Format not recognized'); + end + + maxBlockSize = pow2(20); + + obj.writeNcheck(v,sprintf(':WAVeform:SOURce CHANnel%u',aCH)); % read x origin:WAVeform:SOURce + + xOrg = str2double(obj.writeReceiveCheck(v,':wav:xor?'));% read x origin + + xInc = str2double(obj.writeReceiveCheck(v,':wav:xinc?'));% read x increment => 1/fadc + + xRef = str2double(obj.writeReceiveCheck(v,':wav:xref?')); + + yOrg = str2double(obj.writeReceiveCheck(v,':wav:yor?'));% read y origin + + yInc = str2double(obj.writeReceiveCheck(v,':wav:yInc?'));% read y increment + + yRef = str2double(obj.writeReceiveCheck(v,':wav:yref?')); + + points = str2double(obj.writeReceiveCheck(v,':wav:points?')); + + % timerperiod = max(10,1e-6*points*nBytes); + % set(v,'TimerPeriod',timerperiod); + + timeout = min([max(10,1e-6*points*nBytes) 1000]); + set(v,'Timeout',timeout); + + blockSize = min(pow2(ceil(log2(points))),maxBlockSize); + set(v,'InputBufferSize',nBytes*blockSize); + + blockCount = int8(ceil(max(1,points/blockSize))); + + for block = 0:blockCount - 1 + + startPoint = double(block) * blockSize + 1; + endPoint = min((double(block + 1)) * blockSize,points); + pointsRead = endPoint - startPoint + 1; + + obj.writeNcheck(v, sprintf(':wav:data? %d,%d', startPoint, pointsRead)); + + %Wfm.data(startPoint:endPoint,1) = binblockread(v,blockReadType); %old Matlab + data(startPoint:endPoint,1) = readbinblock(v,blockReadType); + + Junk = fread(v,1,'schar'); %#ok Read out end of file character. + + end + + %scale the signal back to volt + captured_signal = data*yInc+yOrg; + + end + + function writeNcheck(~,visaobj,command) + % 1) send command to intstument + % 2) check if everything is okay + + writeline(visaobj,char(command)); + + % writeline(visaobj,':syst:err?'); + % err_aw=readline(visaobj); + % assert(strfind(err_aw,"No error"),['AWG System Error? after command: ', char(command)]) + + end + + function [aw]=writeReceiveCheck(~,visaobj,command) + % 1) send command to intstument + % 2) receive answer from instrument + % 3) check i f everything is okay + + writeline(visaobj,char(command)); + aw = readline(visaobj); + + % writeline(visaobj,':syst:err?'); + % err_aw=readline(visaobj); + % assert(strfind(err_aw,"No error"),['AWG System Error? after command: ', char(command)]) + + end + + end +end diff --git a/Datatypes/awg_model.m b/Datatypes/awg_model.m new file mode 100644 index 0000000..7487a7b --- /dev/null +++ b/Datatypes/awg_model.m @@ -0,0 +1,10 @@ +classdef awg_model < int32 + + enumeration + M8196A (1) + M8199A (2) + M8199A_ILV (3) + M8199B (4) + end + +end \ No newline at end of file diff --git a/Datatypes/scope_fadc.m b/Datatypes/scope_fadc.m new file mode 100644 index 0000000..cc9c3c6 --- /dev/null +++ b/Datatypes/scope_fadc.m @@ -0,0 +1,22 @@ +classdef scope_fadc < int32 + + + enumeration + GSa_2_5 (2.5) %Gigasamples + GSa_4 (4) %Gigasamples + GSa_5 (5) %Gigasamples + GSa_10 (10) %Gigasamples + GSa_20 (20) + GSa_40 (40) + GSa_80 (80) + GSa_160 (160) + GSa_256 (256) + end + + methods + function val = getValue(obj) + val = double(obj); % Convert int32 to double to get the numeric value + end + end + +end \ No newline at end of file diff --git a/Datatypes/scope_model.m b/Datatypes/scope_model.m new file mode 100644 index 0000000..d5b7587 --- /dev/null +++ b/Datatypes/scope_model.m @@ -0,0 +1,9 @@ +classdef scope_model < int32 + + enumeration + DSAZ634A (1) + UXR1102A (2) + UXR1104A (3) + end + +end \ No newline at end of file diff --git a/projects/Lab_2024/imdd_standard.m b/projects/Lab_2024/imdd_standard.m new file mode 100644 index 0000000..8b13789 --- /dev/null +++ b/projects/Lab_2024/imdd_standard.m @@ -0,0 +1 @@ + diff --git a/projects/standard_system/400G_simulative_setup.mat b/projects/standard_system/400G_simulative_setup.mat index d60da59..c8259fb 100644 Binary files a/projects/standard_system/400G_simulative_setup.mat and b/projects/standard_system/400G_simulative_setup.mat differ diff --git a/projects/standard_system/imdd_minimal_2.m b/projects/standard_system/imdd_minimal_2.m index 50dc18b..ef56fb2 100644 --- a/projects/standard_system/imdd_minimal_2.m +++ b/projects/standard_system/imdd_minimal_2.m @@ -3,15 +3,15 @@ params = struct; params.M = [4]; -params.datarate = [300]; -params.rop = [0]; +params.datarate = [224]; +params.rop = [-5]; -precomp_mode = 0; %0=do nothing ; 1= measure; 2=precomp active +precomp_mode = 2; %0=do nothing ; 1= measure; 2=precomp active 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"; diff --git a/test/awg_test.m b/test/awg_test.m new file mode 100644 index 0000000..98bd341 --- /dev/null +++ b/test/awg_test.m @@ -0,0 +1,15 @@ + +if 0 + + AWG = AwgKeysight("model","M8196A","fdac",92e9,"scaletodac",[1,1,1,1],"skews",[0,0,0,0],"voltages",[0,0,0,0.6]); + + tx_sig = Electricalsignal(clip_out,"fs",92e9); + + tic + AWG.upload("signal4",clip_out); + toc + +end + +SCP = ScopeKeysight("model","DSAZ634A",'autoscale',1,"fadc",'GSa_160',"channel",[1],"recordLen",2000000); +signals = SCP.read();