plot an half-wave rectified 50hz sine in time and plot its first seven fourier coefficients in frequency in matlab

main.m
% Define the parameters
f = 50; % Frequency of the sine wave in Hz
fs = 5000; % Sampling frequency in Hz
t = 0:1/fs:0.1; % Time vector for one cycle

% Generate the half-wave rectified sine wave
x = abs(sin(2*pi*f*t));
x_rect_freq = fft(x); % Calculate the Fourier coefficients
n_harmonics = 7; % Number of harmonics to consider

% Plot the half-wave rectified sine wave in time
figure;
subplot(2,1,1);
plot(t, x);
xlabel('Time (s)');
ylabel('Amplitude');
title('Half-Wave Rectified 50Hz Sine Wave');

% Plot the Fourier coefficients
frequencies = 0:fs/length(x):fs*(length(x)-1)/length(x);
subplot(2,1,2);
stem(frequencies(1:ceil(length(x)/2)), abs(x_rect_freq(1:ceil(length(x)/2))));
xlabel('Frequency (Hz)');
ylabel('Magnitude');
title('Fourier Coefficients of Half-Wave Rectified 50Hz Sine Wave');
802 chars
26 lines