Если у вас возникли сложности с курсовой, контрольной, дипломной, рефератом, отчетом по практике, научно-исследовательской и любой другой работой - мы готовы помочь.
Предоплата всего
Подписываем
Министерство образования республики Беларусь
Учреждение образования
Белорусский государственный университет информатики и радиоэлектроники
Кафедра РТС
ОТЧЁТ
по лабораторной работе №1
Тема: «Моделирование линейных дискретных систем в среде Matlab»
Вариант 2
Выполнил: Проверил:
ст. гр. 000102 Ходыко Д.Л.
Кажуро А.А.
Рогалевич В.В.
Минск 2014
Цель работы:
Исходные данные: КИХ-фильтр. Частота дискретизации , частота среза ,. Полоса пропускания Вектор воздействия: b=[0.0038 0.016 0.027 0.0164 -0.028 -0.09 -0.122 -0.083 0.022 0.136 0.185 0.136 0.022 -0.083 -0.122 -0.09 -0.09 -0.029 0.016 0.027 0.016 0.0038]. .
Ход работы:
b=[0.0038 0.016 0.027 0.0164 -0.028 -0.09 -0.122 -0.083 0.022 0.136 0.185 0.136 0.022 -0.083 -0.122 -0.09 -0.029 0.016 0.027 0.016 0.0038]
a=[1];
N=length(b);
delta=[1;zeros(N-1,1)];
h=filter(b,a,delta);
stem(0:N-1,h)
grid on
Рис.1. Импульсная характеристика КИХ-фильтра.
b=[0.0038 0.016 0.027 0.0164 -0.028 -0.09 -0.122 -0.083 0.022 0.136 0.185 0.136 0.022 -0.083 -0.122 -0.09 -0.029 0.016 0.027 0.016 0.0038]
a=[1];
N=length(b);
Fs=48000;
[h,nT]=impz(b,a,N,Fs);
stem(nT,h), grid
Рис.2. Импульсная характеристика КИХ-фильтра по коэффициентам РУ.
b=[0.0038 0.016 0.027 0.0164 -0.028 -0.09 -0.122 -0.083 0.022 0.136 0.185 0.136 0.022 -0.083 -0.122 -0.09 -0.029 0.016 0.027 0.016 0.0038]
a=[1]
n=0:32;
x=sin(0.72*n)+sin(0.524*n);
y=filter(b,a,x);
plot(n,x,'-or',n,y,'-ob')
grid
hold on
stem(n,x,'-or')
stem(n,y)
gtext('Output signal')
gtext('input signal')
Рис. 3. Реакция ЛДС на входное воздействие КИХ-фильтра.
b=[0.0038 0.016 0.027 0.0164 -0.028 -0.09 -0.122 -0.083 0.022 0.136 0.185 0.136 0.022 -0.083 -0.122 -0.09 -0.029 0.016 0.027 0.016 0.0038]
h=b;
n=0:32;
x=sin(0.72*n)+sin(0.524*n);
y=conv(h,x);
k=length(y);
hold on
plot(n,x,'-or'), grid
stem(n,x,'-or')
nc=0:(k-1);
plot(nc,y,'-ob')
stem(nc,y)
gtext('Output signal')
gtext('Input signal')
Рис.4.Реакцию ЛДС на основе уравнения свертки.
b=[0.0038 0.016 0.027 0.0164 -0.028 -0.09 -0.122 -0.083 0.022 0.136 0.185 0.136 0.022 -0.083 -0.122 -0.09 -0.029 0.016 0.027 0.016 0.0038]
a=[1];
N=length(b)*10;
h=impz(b,a,N);
Ah=fft(h,N);
H=abs(Ah);
ph=angle(Ah);
dl=round(N/2);
H=H(1:dl);
ph=ph(1:dl);
n=0:dl-1;
plot(n,H,'-or',n,ph,'-ob');
grid
Рис.5.АЧХ и ФЧХ фильтра
Вывод:
В ходе выполнения лабораторной работы были рассмотрены методы математического моделирования обработки сигналов линейной дискретной системы.