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

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

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

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

Рис.44 АЧХ

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

Рис.45 ФЧХ

2 Для полученных частотных характеристик реальной и мнимой осуществить обратное преобразование Фурье

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

x=rectpuls(t-0.035,0.07);

plot(t,x)

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

x1=ifft(x,201);

x1=x1/0.0007;

plot(t,real(x1))

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

plot(t,imag(x1))

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

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

x=tripuls(t-0.035,0.07);

plot(t,x)

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

x1=ifft(x,201);

x1=x1/0.0007;

 plot(t,real(x1))

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

 plot(t,imag(x1))

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

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

x=exp(-t/0.028);

plot(t,x)

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

x1=ifft(x,201);

x1=x1/0.0007;

plot(t,real(x1))

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

plot(t,imag(x1))

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

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

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

plot(t,x)

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

x1=ifft(x,201);

x1=x1/0.0007;

plot(t,real(x1))

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

plot(t,imag(x1))

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

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

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

plot(t,x)

Рис.58 периодическая последовательность

x1=ifft(x,201);

x1=x1/0.0007;

plot(t,real(x1))

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

plot(t,imag(x1))

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

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

x=sawtooth(2*pi*t/0.028);

plot(t,x)

Рис.61 пилообразный импульс

x1=ifft(x,201);

x1=x1/0.0007;

plot(t,real(x1))

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

plot(t,imag(x1))

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

Реакция на синусоидальные сигналы в зоне пропускания и задерживания