create a script for fourier analysis with a sampling frequency of 100, omega of 2pi using elseif functions and plot the data in matlab

main.m
% Create a signal, e.g. a sum of two sine waves
t = 0:0.01:10; % time vector
x = 5*sin(2*pi*5*t) + 2*sin(2*pi*20*t); % signal vector

% Compute the Fourier transform
n = length(x);
f = (0:n-1)*(100/n); % frequency vector
y = fft(x);

% Plot the original signal and its Fourier transform
subplot(2, 1, 1);
plot(t, x);
xlabel('Time (s)');
ylabel('Amplitude');
title('Original Signal');
subplot(2, 1, 2)
plot(f,abs(y));
xlabel('Frequency (Hz)');
ylabel('Amplitude');
title('Fourier Transform');

% Find the dominant frequency component and its amplitude
amp = abs(y(1:n/2+1))*2/n; % normalized amplitude
[~, idx] = max(amp);
f_dom = f(idx);
amp_dom = amp(idx);
disp(['Dominant frequency component: ', num2str(f_dom), ' Hz']);
disp(['Amplitude of dominant frequency component: ', num2str(amp_dom)]);
796 chars
29 lines