Спектр сигнала

clear all;clc;
time=10; % Длительность
N=1024; % Количество линий Фурье спектра
Fd=1024; % Частота дискретизации
 
Ak=0.5; % Постоянная составляющая
A1=1; % Амплитуда первой синусоиды
f=10; % Гц
w=2*pi*f; % Рад/с
 
t=0:1/Fd:time; % Дискретное время
dt=time/(Fd-1); % Шаг по времени
 
Signal=Ak+A1*sin(w.*t);
 
Z=fft(Signal,N); % Преобразование Фурье
S=fftshift(Z)*2/N; % Сдвиг 
F=abs(S(N/2+1:N-1)); % Модуль преобразования
F(1)=F(1)/2; % Амплитуда постоянной составляющей
df=Fd/(N*time); % Шаг по частоте в спектре, Гц
Nyq=Fd/(2*time); % Частота Найквиста, Гц
f=-Nyq+df*(0:N-1); % Задание значений частот
f=f(N/2+1:N-1);
figure;
plot(f,F);

clear all;clc;
time=10^-7; % Длительность
N=2048; % Количество линий Фурье спектра
Fd=2048; % Число отсчетов по времени
 
dt=time/(Fd-1); % Шаг по времени
t=0:dt:time; % Дискретное время
 
 
AM=10; % Амплитуда
f=10^9; % Гц
w=2*pi*f; % Рад/с
 
nt=AM.*sin(w.*t); 
%figure;
%plot(t,nt);
Z=fft(nt,N); % Преобразование Фурье
S=fftshift(Z)/N; % Сдвиг 
S=S(N/2+1:N-1)*2;
S(1)=S(1)/2;  % Амплитуда постоянной составляющей
F=abs(S); % Модуль преобразования
df=Fd/(N*time); % Шаг по частоте в спектре, Гц
Nyq=Fd/(2*time); % Частота Найквиста, Гц
f=-Nyq+df*(0:N-1); % Задание значений частот
f=f(N/2+1:N-1);
figure;
plot(f,F);