fouriertransform a function in matlab

main.m
% Define function
Fs = 200; % Sampling frequency
t = 0:1/Fs:1-1/Fs; % Time vector
x = sin(2*pi*10*t) + cos(2*pi*20*t); % Function to be transformed

% Perform Fourier Transform
N = length(x); % Length of signal
X = fft(x); % Compute Fourier Transform
f = Fs*(0:(N/2))/N; % Create frequency vector
P = abs(X/N); % Compute Power of Signal

% Plot results
subplot(2,1,1);
plot(t,x);
xlabel('Time (s)');
ylabel('Amplitude');
title('Original Signal');
grid on;
subplot(2,1,2);
plot(f,P(1:N/2+1));
xlabel('Frequency (Hz)');
ylabel('Power');
title('Power Spectrum');
grid on;
569 chars
25 lines