Использование дискретного преобразования Фурье для вычисления частотных характеристик и выходных сигналов линейных систем

Использование дискретного преобразования Фурье для вычисления частотных характеристик и выходных сигналов линейных систем.

1  Для заданных сигналов получить частотные характеристики с помощью дискретного преобразования Фурье:

Прямое преобразование Фурье:

t=0:0.0007:0.14;

Реакция на прямоугольный импульс

x=rectpuls(t-0.035,0.07);

plot(t,x)

Рис.1 прямоугольный импульс

xf=fft(x,201);

xf=xf*0.0007;

f=0:1/(0.0007*201):1/(0.0007*2);

plot(f,real(xf(1:101)))

Рис.2 реальная часть

plot(f,imag(xf(1:101)))

Рис.3 мнимая часть

plot(f,abs(xf(1:101)))

Рис.4 АЧХ

 plot(f,angle(xf(1:101)))

Рис.5 ФЧХ

Реакция на треугольный импульс

x=tripuls(t-0.035,0.07);

plot(t,x)

Рис.6 треугольный импульс

xf=fft(x,201);

xf=xf*0.0007;

f=0:1/(0.0007*201):1/(0.0007*2);

plot(f,real(xf(1:101)))

Рис.7 реальная часть

plot(f,imag(xf(1:101)))

Рис.8 мнимая часть

plot(f,abs(xf(1:101)))

Рис.9 АЧХ

plot(f,angle(xf(1:101)))

Рис.10 ФЧХ

Реакция на экспоненциальный импульс

x=exp(-t/0.028);

plot(t,x)

Рис.11 экспоненциальный импульс

xf=fft(x,201);

xf=xf*0.0007;

f=0:1/(0.0007*201):1/(0.0007*2);

plot(f,real(xf(1:101)))

Рис12. реальная часть

plot(f,imag(xf(1:101)))

Рис13. мнимая часть

plot(f,abs(xf(1:101)))

Рис.14 АЧХ

plot(f,angle(xf(1:101)))

Рис.15 ФЧХ

Реакция на синусоидальный сигнал

x=sin(pi*t/0.07).*(t<=0.07);

plot(t,x)

Рис.16 синусоидальный сигнал

xf=fft(x,201);

xf=xf*0.0007;

f=0:1/(0.0007*201):1/(0.0007*2);

plot(f,real(xf(1:101)))

Рис.17 реальная часть

plot(f,imag(xf(1:101)))

Рис.18 мнимая часть

plot(f,abs(xf(1:101)))

Рис.19 АЧХ

plot(f,angle(xf(1:101)))

Рис.20 ФЧХ

Реакция на периодическую последовательность

x=0.5*(1+square(2*pi*t/0.028,0.014/0.028*100));