IviInfo = instrhwinfo('ivi');
installedDrivers = IviInfo.Modules

% Create the MATLAB driver
makemid('AgMD1', 'AgMD1.mdd', 'ivi-c');

initOptions = 'Simulate=true, DriverSetup= Cal=0, Trace=false, model=M9210A';
visaAddress = 'PXI17::10::0::INSTR';
myDigitizer = icdevice('AgMD1.mdd', visaAddress, 'optionstring', initOptions);

% Connect to the digitizer using the device object created above
connect(myDigitizer);

% Abort present acquisition if any
invoke(myDigitizer.Waveformacquisitionlowlevelacquisition, 'abort');

% Set the input impedance values of the individual channels
myDigitizer.RepCapIdentifier = 'Channel1';
myDigitizer.Channel.Input_Impedance = 50; % (Ohms)

myDigitizer.RepCapIdentifier = 'Channel2';
myDigitizer.Channel.Input_Impedance = 50; % (Ohms)

% Set the acquisition parameters
numberOfRecords = 1;
ptsPerRecord = 1e4;
samplingRate = 2e9;

invoke(myDigitizer.Configurationacquisition, 'configureacquisition',...
       numberOfRecords, ptsPerRecord, samplingRate);

% Set the individual channel parameters
Range = 0.2;
Offset = 0.0;
Coupling = 1;
Enabled = true;
invoke(myDigitizer.Configurationchannel, 'configurechannel', 'Channel1',...
       Range, Offset, Coupling, Enabled);
invoke(myDigitizer.Configurationchannel, 'configurechannel', 'Channel2',...
       Range, Offset, Coupling, Enabled);

% Set the trigger source, and trigger type
myDigitizer.RepCapIdentifier = 'Channel1';
myDigitizer.Trigger.Active_Trigger_Source = 'External1';

% The hex value can be found in the MD1 driver documentation
AGMD1_VAL_IMMEDIATE_TRIGGER = hex2dec('000003E9');
myDigitizer.Trigger.Trigger_Type = AGMD1_VAL_IMMEDIATE_TRIGGER;

% Determine the minimum amount of memory needed to fetch or read data from
% the digitizer for maximum performance
dataWidth = 64;
numRecords = 1;
offsetWithinRecord = 0;
numPointsPerRecord = ptsPerRecord;

arrayElements = invoke(myDigitizer.Waveformacquisitionlowlevelacquisition,...
                       'queryminwaveformmemory', dataWidth, numRecords,...
                       offsetWithinRecord, numPointsPerRecord);
WaveformArray = zeros(arrayElements, 1);

maxTimeMilliseconds = 50;

% Initiate an acquisition on all enabled channels, wait (up to
% |maxTimeMilliseconds|) for the acquisition to complete, and return the
% waveform for this channel

[sig1, ActualPoints, FirstValidPoint, ~, ~, ~, XIncrement] = ...
    invoke(myDigitizer.Waveformacquisition, 'readwaveformreal64', ...
    'Channel1', maxTimeMilliseconds, arrayElements, WaveformArray);
sig1 = sig1(FirstValidPoint+1:FirstValidPoint+ActualPoints);

% |FETCHWAVEFORMREAL64| returns a previously acquired waveform for the selected
% channel. An acquisition must be made prior to calling this method. For
% this case the previous call to |READWAVEFORMREAL64| has performed the waveform
% acquisition already. Call this method separately for each channel

[sig2, ActualPoints, FirstValidPoint, ~, ~, ~, ~] = ...
    invoke(myDigitizer.Waveformacquisitionlowlevelacquisition, ...
    'fetchwaveformreal64', 'Channel2', arrayElements, WaveformArray);
sig2 = sig2(FirstValidPoint+1:FirstValidPoint+ActualPoints);

% Create the time vector using the XIncrement value returned by the
% |READWAVEFORMREAL64| method
t = 0:XIncrement:(length(sig1)-1)*XIncrement;

figure;
% Plot the channel 1 waveform
subplot(2, 1, 1);
plot(t, sig1);
xlabel('Time (s)');
ylabel('Voltage (V)');
ylim([-.2 0.2]);
title('Digitizer Channel 1');

% Plot the channel 2 waveform
subplot(2, 1, 2);
plot(t, sig2);
xlabel('Time (s)');
ylabel('Voltage (V)');
ylim([-.2 0.2]);
title('Digitizer Channel 2');

disconnect(myDigitizer);
delete(myDigitizer);
clear myDigitizer;