Maple 9.5/10 в математике, физике и образовании

Дьяконов Владимир Павлович

Теперь приступим к тестированию фильтра. Зададим входной сигнал в виде зашумленного меандра с частотой 500 Гц и размахом напряжения 2 В:

> l := round(fs/2/500):

> for n from 0 by 2*1 to T do

> for n2 from 0 to l-1 do

> if n+n2 <= T then

> x[n+n2] := evalf(-l+rand/10^12-0.5);

> fi;

> if n+n2+1 <= T then

> x[n+n2+1] := evalf(1+rand/10^12-0.5);

> fi;

> od:

> od:

Временная

зависимость синтезированного входного сигнала представлена на рис. 11.35.

Рис. 11.35. Синтезированный входной сигнал

Вычислим реакцию фильтра на входной сигнал:

> for n from 0 to Т do

> y[n] : = sum(h[k]*х[n-k], k=0..N);

> od:

Построим график выходного сигнала.

> р := [seq([j/fs, x[j]], j=0..T)]: q:=[seq([j/fs, y[j]], j=0..T)]:

> plot(p, time=0..T/fs/4, labels=[time,volts], title= `Входной сигнал`, color=black);

> plot(q, time=0..T/fs/4, labels=[time, volts], title=`Выгодной сигнал`, color=black);

Временна́я зависимость выходного сигнала показана на рис. 11.36. Нетрудно заметить, что, в конце концов, выходной сигнал вырождается в пятую гармонику входного сигнала, но этому предшествует довольно заметный переходной процесс. Он связан с узкополосностью данного фильтра.

Рис. 11.36. Временна́я зависимость выходного сигнала цифрового фильтра

Вычислим спектры входного и выходного сигналов, подготовив массивы выборок сигналов и применив прямое преобразование Фурье с помощью функции FFT:

> ri := array(1..Т+1): ii := array(1..Т+1):

> for n from 0 to T do

> ri[n+1] := x[n]*2/T; ii[n+1] := 0;

> ro[n+1] := y[n]*2/T; io[n+1] := 0;

> od:

> FFT(m, ri, ii): FFT(m,ro,io):

Построим график спектра входного сигнала, ограничив масштаб по амплитуде значением 0.5 В:

> р := [seq([j*fs/(T+1), abs(ri[j+1]+ii[j+1]*I)], j=0..T/2)]:

> q := [seq([j*fs/(T+1), abs(ro[j+1]+io[j+1]*I)], j=0..T/2)] :

> plot(p, frequency=0..fs/2, y=0..0.5, labels=[частота, V],

titles=`Частотный спектр входного сигнала`, color=black);

Этот график представлен на рис. 11.37. Из него хорошо видно, что спектральный состав входного сигнала представлен только нечетными гармониками, амплитуда которых убывает по мере роста номера гармоники. Пятая гармоника на частоте 2500 Гц находится посередине полосы пропускания фильтра, ограниченной граничными частотами фильтра 2300 и 2700 Гц. Заметны также беспорядочные спектральные линии шума сигнала

в пределах полосы прозрачности фильтра.

Рис. 11.37. Спектрограмма входного сигнала

Теперь построим график спектра выходного сигнала:

> plot(q, frequency=0..fs/2, y=0..0.5, labels=[частота, V],

title=`Частотный спектр выходного сигнала`, color=black);

Он представлен на рис. 11.38. Хорошо видно эффективное выделение пятой гармоники сигнала и прилегающей к ней узкой полосы шумового спектра.

Рис. 11.38. Спектрограмма выходного сигнала цифрового фильтра

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

11.3.5. Моделирование цепи на туннельном диоде

А теперь займемся моделированием явно нелинейной цепи. Выполним его для цепи, которая состоит из последовательно включенных источника напряжения Es, резистора Rs, индуктивности L и туннельного диода, имеющий N-образную вольт-амперную характеристику (ВАХ) — см. файл tdc. Туннельный диод обладает емкостью С, что имитируется конденсатором С, подключенным параллельно туннельному диоду.

Пусть ВАХ реального туннельного диода задана выражением:

> restart;

> А:=-3: а:=10: В:=1*10^(-8): b:=20:

> Id:=Ud->A*Ud*exp(-a*Ud)+B*(exp(b*Ud-1));

Id := Ud →A Udе^(-aUd) + Вe^(bUd-1)

Построим график ВАХ:

> plot(Id(Ud), Ud=-.02..0.76, color=black);

Этот график представлен на рис. 11.39. Нетрудно заметить, что ВАХ туннельного диода не только резко нелинейна, но и содержит протяженный участок отрицательной дифференциальной проводимости, на котором ток падает с ростом напряжения на диоде. Это является признаком того, что такая цепь способна на переменном токе отдавать энергию во внешнюю цепь и приводить к возникновению колебаний в ней различного типа.

Рис. 11.39. ВАХ туннельного диода

Работа цепи описывается системой из двух дифференциальных уравнений:

di/dt= (Es-i(t) * Rs-u(t))/L

du/dt= (i(t)-Id(u(t))/C

Пусть задано Es=0,35 В, Rs=15 Ом, C=10∙10^{– 12}, L=30∙10^{– 9} и максимальное время моделирования tm=10∙10^{– 9}. Итак, задаем исходные данные:

> Es: = .35: Rs:=15: C:=10*10^(-12): L:=30*10^(-6): tm:=10*10^(-9):