Расчет цифрового фильтра

Автор работы: Пользователь скрыл имя, 03 Декабря 2012 в 19:45, курсовая работа

Краткое описание

В данной курсовой работе рассчитывается цифровой фильтр с ЛФЧХ по заданной амплитудно-частотной характеристике. Вычисляются аналетические выражения импульсных характеристик фильтра, далее строится график ФЧХ, АЧХ. Также в данной работе исследуется прохождение через фильтр различных сигналов, а именно единичного импульса и сигнала состоящего из суммы трех гармоник.

Содержание

1 Анализ ТЗ и выбор метода синтеза фильтра…………………………………....5
2 Импульсная характеристика..................................................................................6
3 Графики АЧХ и ФЧХ..............................................................................................7
4 Моделирование
4.1 Единичный импульс..........................................................................................10
4.2 Сумма трех гармоник........................................................................................12
5 Структурная схема КИХ-Фильтра…………………………………………..…14
6 Проверка оптимального порядка………………………………………………15
ЗАКЛЮЧЕНИЕ.......................................................................................................16
СПИСОК ЛИТЕРАТУРЫ.......................................................................................17

Скачать в ZIP архиве (272.66 Кб) Сколько стоит заказать работу?
Вложенные файлы: 1 файл

Kursa4 новый.doc

— 436.00 Кб (Скачать файл)

 

 

 

[HY,wY]=freqz(Y,1,2^18);

HY1=20*log10(abs(HY)/max(abs(HY)));

 

grid on;

 

figure(4);                          subplot(2,1,1);stem(xn.*2*N,X0);xlabel('N');ylabel('X');title ('Единичный импульс');

grid on;

subplot(2,1,2);stem(xn.*2*N,Y);xlabel('N');ylabel('Y');title ('Импульсная харакнеристика');

grid on;

 

figure(5);

subplot(3,1,1);plot(n,sin(2*pi*0.004*n));ylabel('X');xlabel('n');title('Гармоника с частотой 50 Гц')

grid on;axis([0 1400 -1.1 1.1]);

subplot(3,1,2);plot(n,sin(2*pi*0.3*n));ylabel('X');xlabel('n');title('Гармоника с частотой 3,6 КГц')

grid on;axis([0 1000 -1.1 1.1]);

subplot(3,1,3);plot(n,sin(2*pi*0.45*n));ylabel('X');xlabel('n');title('Гармоника с частотой 5,4 КГц')

grid on;axis([0 1000 -1.1 1.1]);

 

 

X2=sin(2*pi*0.004*n)+sin(2*pi*0.3*n)+sin(2*pi*0.45*n);

spx2=fft(X2);

 

Y01=filter(h5_1,1,X2);

Y1=filter(h6_2,1,Y01);

spy2=fft(Y1);

figure(6);

subplot(2,2,1);plot(n,X2);ylabel('X');xlabel('n');title('Сумма трёх гармоник')

grid on;axis([0 1000 -3.1 3.1]);

subplot(2,2,3);plot(n,Y1);ylabel('Y');xlabel('n');title('Сигнал на выходе')

grid on;axis([0 1000 -1.1 1.1]);

subplot(2,2,2);plot(n/N.*Fc+6000,abs(spx2)./max(abs(spx2)));xlabel('Частота, Гц');ylabel('X');title('Спектр входного сигнала')

axis([0 6000 0 1])

grid on;

subplot(2,2,4);plot(n/N.*Fc+6000,abs(spy2)./max(abs(spy2)));xlabel('Частота, Гц');ylabel('Y');title('Спектр выходного сигнала')

axis([0 6000 0 1])

grid on;

figure(7);

plot(wY./(2*pi).*Fc,HY1,'r');xlabel('Частота, Гц');ylabel('dB');title ('Проверка АЧХ Фильтра');grid on;

grid on;

 

 

Ne =14873;

fe1=0.003;

fe2=0.0045;

fe3=0.0055;

fe4=0.007;

we01=(fe1*2*pi+fe2*2*pi)/2;

we02=(fe3*2*pi+fe4*2*pi)/2;

ne = -(Ne-1)/2:1:(Ne-1)/2;

he = (sin(we01*ne)-sin(we02*ne))./(ne*pi);

he((Ne+1)/2)=1+((we01-we02)./pi);

 

 Ne_2 =37;

fe1_2=0.25;

fe2_2=0.4167;

we0_2=(fe1_2*2*pi+fe2_2*2*pi)/2;

ne_2 = -(Ne_2-1)/2:1:(Ne_2-1)/2;

he_2 = sin(w0_2*ne_2)./(ne_2*pi);

he_2((Ne_2+1)/2)=we0_2./pi;

 

he5_1=he.*BLACKMANHARRIS(Ne)';

he6_2=he_2.*GAUSSWIN(Ne_2)';

 

 

[He5_1,we5_1]=freqz(he5_1,1,2^14);

[He6_2,we6_2]=freqz(he6_2,1,2^14);

 

He1=20*log10(abs(He6_2.*He5_1)/max(abs(He6_2.*He5_1)));

 

figure(8);

subplot(2,2,1);plot(w6_2./(2*pi).*Fc,H1,'r');xlabel('Частота, Гц');ylabel('dB');title ('АЧХ Фильтра при порядке 39');grid on;axis([3000 5500 -80 2]);

subplot(2,2,3);plot(w6_2./(2*pi).*Fc,H1,'r');xlabel('Частота, Гц');ylabel('dB');title ('АЧХ Фильтра при порядке 14875');grid on;axis([35 65 -80 2]);

subplot(2,2,2);plot(we6_2./(2*pi).*Fc,He1,'r');xlabel('Частота, Гц');ylabel('dB');title ('АЧХ Фильтра при порядке 37');grid on;axis([3000 5500 -80 2]);

subplot(2,2,4);plot(we6_2./(2*pi).*Fc,He1,'r');xlabel('Частота, Гц');ylabel('dB');title ('АЧХ Фильтра при порядке 14873');grid on;axis([35 65 -80 2]);

 

end


Информация о работе Расчет цифрового фильтра