Спектр фаз для прямоугольного импульса

@bronikovskii · Регистрация: 02.10.2020

Author24 — интернет-сервис помощи студентам

Здесь (https://www.cyberforum.ru/post14902653.html) уже задавал вопрос по поводу спектра фаз периодической последовательности прямоугольных импульсов. Ответа не получил. Покопался в теме, научился строить графики спектра амплитуд. Использую следующий код для построения спектра фаз меандра. Но результат мне не нравится. Помогите исправить код.

Matlab M

tau =0.5; % Длительность импульсов
q = 2; % Скважность
T=tau*q; % Период следования импульса
N=11; % Число импульсов
d=(0:N-1)+tau/2; % Вектор задержки для pulstran
A=1; % Амплитуда импульсов
Fs=1000; % Частота дискретизации
t=0:1/Fs:10; % Длительность сигнала
S=A.*pulstran(t,d,'rectpuls', tau); % Синтез меандра
n = length(S); % Длина самого большого измерения массива
FFTsig = fft(S, n); % Быстрое преобразование Фурье
theta = angle(FFTsig);
theta = theta(1:n/2); % Нормировка
F = 0:Fs/n:Fs/2-Fs/n; % Частоты гармоник
plot(F,theta); % Построение фазового спектра

@grisha_lu · 19.10.2020, 03:56

Честно, не очень разбираюсь в этом и возможно не правильно

Кликните здесь для просмотра всего текста

Matlab M

tau =0.5; % Длительность импульсов
q = 2; % Скважность
T=tau*q; % Период следования импульса
N=11; % Число импульсов
d=(0:N-1)+tau/2; % Вектор задержки для pulstran
A=1; % Амплитуда импульсов
Fs=100; % Частота дискретизации
t=0:1/(2*Fs):10-1/(2*Fs); % Длительность сигнала
S=A.*pulstran(t,d,'rectpuls', tau); % Синтез меандра
FFTsig = fft(S); % Быстрое преобразование Фурье
theta = angle(FFTsig);
%theta = theta(1:n/2); % Нормировка
f=-Fs:1/t(end):Fs;
figure
hold on
title('Фазовый спектр')
plot(f,theta); % Построение фазового спектра
hold off
 
re=real(fft(S));
im=imag(fft(S));
phase=asin(im);
 
% Непрерывно фазовый спектр
uphase=unwrap(abs(phase));
figure 
hold on
title('Непрерывно фазовый спектр')
plot(f, fftshift(uphase))
hold off
 
figure 
hold on
title('Real and Imag')
plot(f, fftshift(im), 'color', 'black')
plot(f, fftshift(re),'color', 'b')
hold off
 
 
% Фаза 1
figure 
hold on
title('Phase, способ 1')
plot(f, fftshift(phase))
hold off
 
% Фаза 2
spi=imag(fft(S));
sp=real(fft(S));
phs=fftshift(spi./abs(sp));
   for k=1:1:length(spi)
      if (abs(spi(k))<1)
          phs(k)=0;
      else
          phs(k)=spi(k)./abs(spi(k));
      end
end
phs=fftshift(phs*pi/2.);
 
figure
hold on
title('Phase, способ 2')
stem(f,phs);
xlabel('Частота'); 
zoom on;
hold off

	19.10.2020, 03:56