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

plot(t,x)

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

xf=fft(x,201);

xf=xf*0.0007;

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

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

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

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

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

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

Рис.24 АЧХ

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

Рис.25 ФЧХ

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

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

plot(t,x)

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

xf=fft(x,201);

xf=xf*0.0007;

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

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

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

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

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

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

Рис.29 АЧХ

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

Рис.30 ФЧХ

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

x=sin(2*pi*335*t);

plot(t,x)

Рис.31 синусоидальный сигнал в зоне пропускания

xf=fft(x,201);

xf=xf*0.0007;

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

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

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

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

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

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

Рис.34 АЧХ

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

Рис.35 ФЧХ

x=sin(2*pi*455*t);

plot(t,x)

Рис.36 синусоидальный сигнал в зоне задерживания

xf=fft(x,201);

xf=xf*0.0007;

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

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

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

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

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

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

Рис.39 АЧХ

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

Рис.40 ФЧХ

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

x=sin(2*pi*335*t)+sin(2*pi*455*t);

plot(t,x)

Рис.41 бигармонический сигнал

xf=fft(x,201);

xf=xf*0.0007;

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

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

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