Основы математического аппарата и математические модели линейных ДСАУ. Решетчатые функции, страница 5

Чем замечательны такие уравнения? Тем, что в отличие от дифференциальных уравнений, мы можем найти решение непосредственно через значения исходной управляющей функции g(t), что покажем позднее. А сейчас решим "Пример 2" перехода от формы (15.1.) к (18.1.) и убедимся, что это трудоемкая задача.

________________________________________________________________

5)Примечание5 Доказательства достоверности формул (12.1.), (16.1) и (17.1.) просты, но связаны с огромными рутинными выкладками и расчетами. Например, (12.1) получим, если решим пример 4 дляf[n,0]. Справедливость (14.1) подтвердим, вычислив каждый из коэффициентовi} и {bj} в уравнении (15.1) и просуммировав слагаемые с равными значениями решётчатых функций y[n–i] и e[n–j]. Эти формулы приведены в нескольких источниках, например,[9]. Тем более они не столь - уж важны в нашем курсе.

Пример 2.  Задано неоднородное разностное уравнение в такой форме:

.                  (*)

Чтобы записать данное уравнение в форме (18.1.) рассчитаем его коэффициенты при m=3, k=2 по формулам (17.1.)

При i=0 и ν=0:   c3=a3(-1)0C30=a3=3.

При i=1, ν=0  и ν=1:   c2=a3(-1)1C31+a2(-1)0C20=-3a3+a2=-9+4=-5.

При i=2, ν=0 , ν=1 и ν=2: c1= a3(-1)2C32+a2(-1)1C21+ a1(-1)0C10=3a3 -2a2 +a1=9-8+5=6.

При i=3, ν=0 , ν=1, ν=2 и ν=3: c0= a3(-1)3C33+a2(-1)2C22+ a1(-1)1C11+a0(-1)0C00=-a3 +a2 -a1+a0=

-3+4 -5 +1=-3.

Соответственно, коэффициенты правой части уравнения будут следующими.

При j=0 и ν=0:   d2=b2(-1)0C20=b2=3.

При j=1, ν=0  и ν=1:   d1=b2(-1)1C21+b1(-1)0C10=-2b2+b1=-6+1=-5.

При j=2, ν=0 , ν=1 и ν=2: d0= b2(-1)2C22+b1(-1)1C11+ b0(-1)0C00=b2 -b1 +b0=3-1+1=3.

Итак, вместо исходного выражения (*), получим:

  (**)

Чтобы раскрыть все значения коэффициентов{сi} и {dj}в выражении (**) пришлось вычислить шестнадцать(!!!) сочетаний типа  Сm-ν i-ν и выполнить иные арифметические процедуры. А ведь это было уравнение всего - то третьего порядка! Возможно, для решения задачи использовать пакет "Matcad",но вряд - ли это упростит и ускорит процесс.

Теперь перейдем к описанию процедуры решения уравнения (18.1.). Пусть ДСАУ, модель, которой изображена на рис. 2.1., находилась в покое, и в ней не было запасено энергии. Сигнал на выходе замкнутой ДСАУ . Тогда дискретную выборку (решётчатую функцию) y[n,0] можно, бесспорно, представить в виде  y[0, 0]=y[-1, 0]=y[-2, 0]=   =y[-(m-1), 0]=0,а сигналуправления e [-1,0]=e [-2, 0]=   =e [-(k-1), 0]=0. Это будут исходные нулевые начальные условия (ННУ) решения задачи.

На вход системы можно подать любой аналитически заданный сигнал: . Для упрощения дальнейших рассуждений ограничимся регулярным единичным сигналом: . Тогда в момент времени t=0, соответственно n=0, и g[0, 0]=1.

Запишем разностное уравнение замыкания системы  в обобщенном виде для любого "n"  интервала квантования6).

e[n, 0]=g[n, 0]–y[n, 0],                                                                               (19.1.)

Следовательно, на первом интервале квантования (от момента времени t 0 =0Т до t 1=1Т) выражение (19.1.) будет иметь такое решение:

e*[0, 0]= g[0, 0]– g[0, 0]–y[0, 0]=1–0=1 = const.                                    (19.1.)0

________________________________________________________________

6)Примечание 6. К этим выражениям необходимо дать разъяснения. Уравнение (19.1.) записано для РФ. Но, в силу разных физических причин, которые будут поясняться постепенно уравнение(19.1)0   имеет другой смысл. Прежде всего,  физически  - непрерывный сигнал, а y[n, 0] его фиктивная дискретная выборка. Поэтому y[0, 0] начальное значение непрерывного сигнала  на первом интервале квантования. Сигнал задания g[0, 0]= 1[0, 0] - решётчатая функция. Сигнал ошибки e*[0, 0] непрерывная функция на первом интервале равная РФ g[0, 0]= 1[0, 0]. Особенности формирования функции e*[n, 0] ,будут пояснены в последующих лекциях.

Итак, полученное значение e*[0, 0]= g[0, 0]= 1[0, 0] подставим в уравнение (18. 1.), и рассчитаем y[0, 0], в качестве конечного значения выходного сигнала на первом интервале квантования. В силу ННУ  уравнение (18.1.) будет таким: