Часть I

Текст m-файла к лабораторной работе №4

clear; close all;

b = [0.1 -0.11 0]; %коэффициенты числителя

a = [1 -1.55 0.75]; %коэффициенты знаменателя

N = 5; %число двоичных разрядов

am = size(a);

bm = size(b);

nmax = 2^N;

%квантование коэффициентов

bq = round(b*nmax)/nmax;

aq = round(a*nmax)/nmax;

%последовательное программирование

sos = tf2sos(b,a);

sosq = round(sos*nmax)/nmax;

[bsosq, asosq] = sos2tf(sos);

%построение АЧХ

[h, f] = freqz(b,a);

hq = freqz(bq,aq);

hsosq = freqz(bsosq,asosq);

figure;

plot(f/pi,abs(h),'r',f/pi,abs(hq),f/pi,abs(hsosq),'g-.');

xlabel('Frequency');

ylabel('Magnitude');

legend('original','quantized, direct form','quantized, second-order sections form');

%Расчет импульсных характеристик

plottime=100;

x2(plottime)=0;

x21(plottime)=0;

%импульсная характеристика фильтра без округления

%результатов вычислений

for k=1:plottime

if k<=bm(2)

x21(k) = b(k);

end;

for i=2:am(2)

if k>i-1

x2k(i) = x21(k-i+1);

else x2k(i) = 0;

end;

x21(k) = x21(k) - a(i)*x2k(i);

end;

end;

%импульсная характеристика фильтра с округлением

%результатов вычислений

for k=1:plottime

if k<=bm(2)

x2(k) = round(b(k)*nmax)/nmax;

end;

for i=2:am(2)

if k>i-1

x2k(i) = x2(k-i+1);

else x2k(i) = 0;

end;

x2(k) = x2(k) - round(a(i)*x2k(i)*nmax)/nmax;

end;

end;

figure;

plot(x21,'r--')

hold on;

plot(x2);

xlabel('Time');

ylabel('Impulse response');

legend('original','quantized');

Е.Л. Собакин

ЦИФРОВАЯ СХЕМОТЕХНИКА

Часть I

Учебное пособие

Томск 2002

УДК 681.325.6

Собакин Е.Л. Цифровая схемотехника. Учеб. пособие. Ч.I. - Томск: Изд. ТПУ, 2002. - 160с.

В пособии изложены основные вопросы курса лекций для студентов специальности 210100 - Управление и информатика в технических системах. Пособие подготовлено на кафедре автоматики и компьютерных систем ТПУ, соответствует учебной программе дисциплины и предназначено для студентов института дистанционного образования.

Печатается по постановлению Редакционно-издательского Совета Томского политехнического университета

Рецензенты:

В.М. Дмитриев - профессор, доктор технических наук, заведующий кафедрой теоретических основ электротехники Томского университета систем управления и радиоэлектроники;

С.И. Королёв - директор ТОО НПО «Спецтехаудитсервис»,

кандидат технических наук, старший научный сотрудник.

Темплан 2002

ã Томский политехнический университет, 2002