clear all, close all, clc
N=8192; % Ухнплйкьт бсйимьт дейгмЬфщн уфзн чспнйкЮ дйЬскейб 3.4 seec
f=10; % ухчньфзфб 10Hz
T=3.4; % чспнйкЮ дйЬскейб уЮмбфпт. РспупчЮ Мзн фп мресдеэефбй ме фп Ф фзт ресйьдпх
dt=T/N; % вЮмб Ьопнб чсьнпх
t=0:dt:T;

y=sin(2*pi*f*t);
figure(1)
subplot(2,1,1)
plot(t,y)
title('Signal')
xlabel('time (seconds)')
subplot(2,1,2)
plot(t(1:1024),y(1:1024))
title('Signal')
xlabel('time (seconds)')


z=fft(y); % Discrete Fourier Transform (with FFT)
s=abs(z);
s=s/(N/2);
s=s.^2;

df=1/T;
fs=1/dt;
f=0:df:fs;
figure(2)
semilogy(f,s)
axis([0 50 10^-1 10^1])
title('Fasma Isxyos')
xlabel('Frequency(Hz)')
ylabel('Isxys')









% N=10000;  деЯгмбфб
%f=10; % ухчньфзфб 100Hz
%T=3.4; % ресЯпдпт 3.4s
%vx=(1/f)/N; % вЮмб Ьопнб чсьнпх
%t=0:0.0002:1;

%y=sin(2*pi*f*t);
%figure(1)
%plot(t(1:200),y(1:200))

% Discrete Fourier transform (via the Fast Fourier Transform)
%Y = fft(y);

% Compute the amplitude of the spectrum and center it
%S = abs(fftshift(Y));
%f=0:1:5000;
%figure(2)
%plot(f,S)

%z=fft(y,2048); % Discrete Fourier Transform (with FFT)
%s=abs(z);
%i=[s(1:N/2)/N/2];


%f2=[0:N/2-1]*1/T;
%figure(2);
%semilogy(f2, i.^2)
%axis([0 20 10^-8 10^3])
%title('Fasma Isxyos')
%xlabel('Frequency(Hz)')
%ylabel('Isxys')
