write a matlab code that will extract the peak-to-peak amplitude of the m-wave and h-reflex at each stimulation step an create a graph of raw units in matlab

main.m
% Load data
load('data.mat'); % replace with your data filename

% Define constants
fs = 1000; % sampling frequency in Hz
prestim_time = 50; % prestimulus time in ms
poststim_time = 200; % poststimulus time in ms
stim_duration = 1; % stimulus duration in ms
num_stim = length(data); % number of stimulation steps

% Preallocate output variables
mwave_pp = zeros(1,num_stim);
h_reflex_pp = zeros(1,num_stim);

% Loop through each stimulation step
for i = 1:num_stim
    % Get the stimulus onset time
    stim_onset = (i-1)*stim_duration + prestim_time;
    
    % Get the window for prestimulus baseline
    prestim_window = round((stim_onset - prestim_time/2):(stim_onset - prestim_time/4));
    
    % Get the window for m-wave and h-reflex peak-to-peak amplitude
    resp_window = round((stim_onset + stim_duration):(stim_onset + poststim_time));
    
    % Get the data for the windows
    prestim_data = data(prestim_window);
    resp_data = data(resp_window);
    
    % Get the peak-to-peak amplitudes
    mwave_pp(i) = peak2peak(resp_data);
    h_reflex_pp(i) = max(resp_data) - min(prestim_data);
end

% Create a plot of peak-to-peak amplitudes
figure;
plot(1:num_stim, mwave_pp, 'b', 1:num_stim, h_reflex_pp, 'r');
xlabel('Stimulation step');
ylabel('Peak-to-peak amplitude (raw units)');
legend('M-wave', 'H-reflex');
title('M-wave and H-reflex peak-to-peak amplitudes at each stimulation step');
1407 chars
42 lines