% fft abs lenght n/2 22
clear all
clc
Fs=1000       % فرکانس نمونه برداری
L=1            % ثانیه
t=0:1/Fs:L-1/Fs
f1=10       % فرکانس موج سینوسی
f2=30       % فرکانس موج سینوسی
f3=50       % فرکانس موج سینوسی

y1=3*sin(2*pi*f1*t)
y2=1*sin(2*pi*f2*t)
y3=2*sin(2*pi*f3*t)
 
figure          % حفظ نمودار زیر
subplot(311)        % از سه نمودار بعد در یک ستون، اولین نمودار بشرح زیر
plot(t,y1)          % رسم نمودار با محور افقی تی و محور عمودی وای یک
subplot(312)        % از سه نمودار بعد در یک ستون،  دومین نمودار بشرح زیر
plot(t,y2)
subplot(313)
plot(t,y3)
pause(2)                % دو ثانیه توقف


% رسم جمع سه تابع در حوزه زمان

yt=y1+y2+y3
figure
plot(t,yt)
pause(2)

% رسم بر حسب حوزه فرکانس - فوریه
figure
ff1=fft(yt)         % تبدیل حوزه زمان وای تی به حوزه فرکانس به نام اف اف یک
subplot(411)        % ازچهار نمودار بعد در یک ستون، اولین نمودار بشرح زیر
plot(ff1)
pause(2)

% رسم قدر مطلق همراه  رسم  تقارن آن 
ff2=abs(ff1)        % جذر جمع مجذور حقیقی و موهومی اف اف یک
subplot(412)        % ازچهار نمودار بعد در یک ستون، دومین  نمودار بشرح زیر
plot(ff2)
pause(2)

%رسم نصفه اول
N=length(ff2)       % طول ماتریس اف اف 2 یعنی اف اف دو چند تا عدد دادرد
ff3=ff2(1:N/2)      % نصف اعداد اف اف دو را در اف اف 3 بریز
subplot(413)
plot(ff3)
pause(2)
% عدد فرکانسها صحیح نیست یکی بیشتر است

% تصحیح محور ایکس
x=(0:N/2-1).*Fs/N   % به تعداد ان دوم در اف اس تقسیم بر ان و ضرب یک به یک کن
subplot(414)
plot(x,ff3)
pause(2)

% رسم آخرین نمودار بشکل مجزا
figure
plot(x,ff3)     % یک بار دیگر نمودار را به تنهایی رسم کن


% it is better Fs mazrab 2 bashad Fs=1024

% mishod tabdil ff ba tool delkhah anjam dad  ff1=fft(yt,200) deghat kamtar
% it is better baraye signalhaye bozorg ein addad 200 nazdiktarin tavan 2 baad az Fs bashad ke
% ba in dastoor dar zeer khat 30
% n=nextpow2(length(yt))
% ff1= fft(yt,2^n)