Решение обыкновенных линейных разностных уравнений с постоянными коэффициентами. Разностные уравнения и их применение в экономике Линейные разностные уравнения первого порядка примеры

Решение обыкновенных линейных разностных уравнений

с постоянными коэффициентами

Связь выхода и входа линейной дискретной системы может быть описана обыкновенным линейным разностным уравнением с постоянными коэффициентами

где y[ n] - выходной сигнал в момент n ,

x[ n] - входной сигнал в момент n ,

a i , b k – постоянные коэффициенты.

Для решения таких уравнений могут использоваться два метода

Прямой метод,
Метод Z – преобразования.

Вначале рассмотрим решение линейного разностного уравнения с помощью прямого метода.

Общее решение неоднородного (с отличной от нуля правой частью) линейного разностного уравнения равно сумме общего решения линейного однородного разностного уравнения и частного решения неоднородного уравнения

Общее решение однородного разностного уравнения (zero- input response ) y h [ n]

определяется в виде

Подставляя это решение в однородное уравнение, получаем

Такой полином называют характеристическим полиномом системы. Он имеет N корней . Корни могут быть действительными или комплексными и некоторые корни - совпадающими (кратными).

Если корни являются действительными и разными, то решение однородного уравнения имеет вид

где коэффициенты

Если некоторый корень, например, λ 1 имеет кратность m , то соответствующий ему член решения приобретает форму

Если все коэффициенты однородного уравнения и соответственно характеристического многочлена действительны, то два члена решения, соответствующие простым комплексно сопряженным корням можно представить (записать) в виде , при этом коэффициенты A, B определяются по начальным условиям.

Вид частного решения y p [ n] уравнения зависит от правой части (входного сигнала) и определяется согласно нижеприведенной таблице

Таблица 1. Вид частного решения для различного характера правой части

Входной сигнал x[n]

Частное решение y p [n]

A (константа)

Решение линейного разностного уравнения методом Z – преобразования заключается в применении Z – преобразования к уравнению с использованием свойств линейности и временного сдвига. В результате получается линейное алгебраическое уравнение относительно Z - изображения искомой функции. Обратное Z – преобразование дает искомое решение во временной области. Для получения обратного Z – преобразования чаще всего используется разложение рационального выражения на простые (элементарные) дроби, так как обратное преобразование от отдельной элементарной дроби имеет простой вид.

Заметим, что для перехода во временную область могут использоваться и другие методы вычисления обратного Z – преобразования.

Пример . Определим отклик (выходной сигнал) системы, описываемой линейным разностным уравнением , на входной сигнал

Решение .

1. Прямой метод решения уравнения.

Однородное уравнение . Его характеристический полином .

Корни полинома .

Решение однородного уравнения .

Поскольку,то частное решение определяем в виде .

Подставляем его в уравнение

Для нахождения константы К примем n = 2 . Тогда

Или , К=2,33

Отсюда частное решение и общее решение разностного уравнения (1)

Найдем константы С 1 и С 2 . Для этого положим n = 0 , тогда из исходного разностного уравнения получаем . Для данного уравнения

Поэтому . Из выражения (1)

Следовательно,

Из выражения (1) для n = 1 имеем .
Получаем следующие два уравнения для С 1 и С 2

Решение этой системы дает следующие значения: С 1 =0,486 и С 2 = -0,816.

Следовательно, общее решение данного уравнения

2. Решение методом Z – преобразования.

Возьмем Z – преобразование от исходного разностного уравнения , учитывая свойство (теорему) временного сдвига . Получаем

На практике простейшие разностные уравнения возникают при исследовании например величины банковского вклада. Эта величина является переменной Y x , представляющей сумму, которая накапливается по установленному закону при целочисленных значениях аргумента x . Пусть сумма Y o положена в банк при условии начисления 100 r сложных процентов в год. Пусть начисление процентов производится один раз в год и x обозначает число лет с момента помещения вклада (x = 0, 1, 2,...). Обозначим величину вклада по истечении x лет через Y x . Мы получаем

Y x = (1+r)Y x-1.

Если начальная сумма составляет Y o , мы приходим к задаче отыскания решения полученного разностного уравнения, подчиненного начальному условию Y x = Y o при x = 0. Полученное разностное уравнение содержит Y x и значение этой переменной на один год раньше, т.е. Y x-1; в данном случае аргумент x явно не входит в разностное уравнение.

Вообще говоря, обыкновенное разностное уравнение устанавливает связь между значениями функции Y = Y(x ), рассматриваемой для ряда равноотстоящих значений аргумента x , но можно без ограничения общности считать, что искомая функция определена для равноотстоящих значений аргумента с шагом, равным единице. Таким образом, если начальное значение аргумента есть x , то ряд его равноотстоящих значений будет x , x+1, x+2,... и в обратном направлении: x , x-1, x-2,.... Соответствующие значения функции будем обозначать Y x, Y x+1, Y x+2, ... или Y x , Y x-1, Y x-2, .... Определим так называемые разности различных порядков функции Y x с помощью следующих формул:

Разности первого порядка

D Y x = Y x+1 - Y x ,

D Y x+1 =Y x+2 - Y x+1,

D Y x+2 = Y x+3 - Y x+2,

... ... ... ... ...

Разности второго порядка

D 2 Y x = D Y x+1 - D Y x ,

D 2 Y x+1 = D Y x+2 - D Y x+1 ,

D 2 Y x+2 = D Y x+3 - D Y x+2 ,

... ... ... ... ...

Разности третьего порядка

D 3 Y x = D 2 Y x+1 - D 2 Y x ,

D 3 Y x+1 = D 2 Y x+2 - D 2 Y x+1 ,

... ... ... ... ...

Обыкновенным разностным уравнением называется уравнение, связывающее значения одного независимого аргумента x , его функцииY x и разностей различных порядков этой функции D Y x , D 2 Y x, D 3 Y x, .... Такое уравнение можно записать в общем виде следующим образом:

j (x , Y x , D Y x , D 2 Y x D 3 Y x , D n Y x ) = 0, (10.1)

которое по форме аналогично дифференциальному уравнению.

Порядком разностного уравнения называется порядок наивысшей разности, входящей в это уравнение. Разностное уравнение (10.1) часто удобнее записать, пользуясь не разностями неизвестной функции, а ее значениями при последовательных значениях аргумента, то есть выразить D Y x , D 2 Y x, D 3 Y x ,... через Y x , Y x+1 , Y x+2, .... Уравнение (10.1) можно привести к одной из двух форм:

y (x , Y x , Y x+1, ...,Y x+n ) = 0, (10.2)

x (x , Y x , Y x-1, ...,Y x -n) = 0. (10.3)

Общее дискретное решение Y x обыкновенного разностного уравнения n -го порядка представляет функцию x (x = 0, 1. 2,...), содержащую ровно n произвольных постоянных:

Y x = Y(x, C 1 , C 2 ,..., C n ).

Паутинообразная модель

Пусть рынок какого-либо отдельного товара характеризуется следующими функциями спроса и предложения:

D = D(P), S = S(P).

Для существования равновесия цена должна быть такой, чтобы товар на рынке был распродан, или

D( P) = S(P).

Цена равновесия задается этим уравнением (которое может иметь множество решений), а соответствующий объем покупок-продаж, обозначаемый через , - следующим уравнением:

D () = S( ).

Динамическая модель получается при наличии запаздывания спроса или предложения. Простейшая модель в дискретном анализе включает неизменное запаздывание или отставание предложения на один интервал:

D t = D (P t) и S t = S (P t-1).

Это может случиться, если для производства рассматриваемого товара требуется определенный период времени, выбранный за интервал. Действие модели таково: при заданном P t-1 предшествующего периода объем предложения на рынке в текущем периоде будет S (P t-1), и величина P t должна установиться так, чтобы был куплен весь объем предложенного товара. Иными словами, P t и объем покупок-продаж X t характеризуются уравнением:

X t = D (P t) = S (P t-1).

Итак, зная исходную цену P o , с помощью этих уравнений мы можем получить значения P 1 и X 1. Затем, используя имеющуюся цену P 1, из соответствующих уравнений получим значения P 2 и X 2 и т.д. В общем изменение P t характеризуется разностным уравнением первого порядка (одноинтервальное отставание):

D (P t) = S (P t-1).

Решение можно проиллюстрировать диаграммой, представленной на рис.5, где D и S - соответственно кривые спроса и предложения, а положение равновесия (со значениями и ) соответствует точке их пересечения Q. Цена в начальный момент времени равна P o . Соответствующая точка Q o на кривой S дает объем предложения в период 1. Весь этот предложенный объем товара раскупается при цене P 1 , заданной точкой Q 1 на кривой D с той же ординатой (X 1), что и Q o . Во второй период времени движение происходит сначала по вертикали от точки Q 1 к точке на кривой S, дающей X 2, а затем по горизонтали - к точке Q 2 на кривой D. Последняя точка характеризует P 2 . Продолжение этого процесса и дает график паутины , показанный на рис. 5. Цены и объемы (покупок - продаж) в последовательные периоды времени являются соответственно координатами точек Q 1 , Q 2 , Q 3 ,... на кривой спроса D. В рассматриваемом случае последовательность точек стремится к Q. При этом точки поочередно располагаются на левой и правой стороне от Q. Следовательно, и значения цены P t стремятся к , располагаясь поочередно по обе стороны от . Точно так же обстоит дело и с объемами покупок - продаж (X t ).

Решение можно получить алгебраически для случая линейных функций спроса и предложения: D = a +aP , S = b +bP . Значения равновесия и будут заданы уравнениями

A +a = b +b ,

то есть

= (a - b )/(b - a), = (b a - a b )/(b - a). (10.4) . р t-1. (10.7)

Уравнения (10.7) аналогичны (10.5), за исключением того, что они описывают отклонения от уровней равновесия (теперь уже известно, что таковые существуют). Оба эти уравнения являются разностными уравнениями первого порядка. Положим c = b /a и подставим его в уравнение (10.7), так что разностное уравнение относительно р t будет

р t = c р t-1 . (10.8)

При данном значении р o в момент t = 0 из (10.8) получаем решение:

р t = р o c t,

или

P t = + (P o - ) c t .

Уравнение вида

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

Обычно вместо уравнения (1) рассматривается уравнение, которое получается из (1) путем перехода от конечных разностей к значению функции, т. е. уравнение вида

Если в уравнении (2) функция, то такое уравнение называется однородным.

Рассмотрим однородное уравнение

Теория линейных разностных уравнений аналогична теории линейных дифференциальных уравнений.

Теорема 1.

Если функции являются решениями однородного уравнения (3), то функция

также является решением уравнения (3).

Доказательство.

Подставим функции в (3)

т. к. функция является решением уравнения (3).

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

(4)

Если (4) имеет место только при то функции , называются линейно независимыми.

Любое k линейно независимымых решений уравнения (3) образуют фундаментальную систему решений.

Пусть линейно независимымые решения уравнения (3), тогда

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

Будем искать решение уравнения (3) в виде:

Подставим в уравнение (3)

Поделим уравнение (5) на

Характеристическое уравнение. (6)

Положим, что (6) имеет только простые корни Нетрудно убедиться, что являются линейно независимыми. Общее решение однородного уравнения (3) имеет вид

Пример.

Рассмотрим уравнение

Характеристическое уравнение имеет вид

Решение имеет вид

Пусть корень имеет кратность r. Этому корню соответствует решение

Если предположить, что остальные корни не являются кратными, то общее решение уравнения (3) имеет вид

Рассмотрим общее решение неоднородного уравнения (2).

Частное решение неоднородного уравнения (2), тогда общее решение

ЛЕКЦИЯ 16

План лекции

1. Понятие о D и Z - преобразованиях.

2. Область применения D и Z - преобразований.

3. Обратные D и Z - преобразования.

ДИСКРЕТНОЕ ПРЕОБРАЗОВАНИЕ ЛАПЛАСА.

Z – ПРЕОБРАЗОВАНИЕ.

В прикладных исследованиях, связанных с использованием решетчатых функций, широко применяется дискретное преобразование Лапласа (Д – преобразование) и Z – преобразование. По аналогии с обычным преобразованием Лапласа дискретное задается в виде

где (1)

Символически Д – преобразование записывается в виде

Для смещенных решетчатых функций

где - смещение.

Z – преобразование получается из Д – преобразования подстановкой и задается соотношением

(3)

Для смещенной функции

Функция называется оригиналом, если

2) существует показатель роста, т. е. найдутся такие и , что

(4)

Наименьшее из чисел (или предел, к которому стремится наименьшее число), для которого справедливо неравенство (4), называется абсциссой абсолютной сходимости и обозначается

Теорема.

Если функция является оригиналом, то изображение определено в области Re p > и является в этой области аналитической функцией.

Покажем, что при Re p > ряд (1) абсолютно сходится. Имеем

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

Вторую часть теоремы примем без доказательств.

Изображение является периодической функцией с мнимым периодом

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

и является в этой полу полосе аналитической функцией.

Найдем область определения и аналитичности функции F(z), положив . Покажем, что полу полоса плоскости p преобразованием переводится в область на плоскости z: .

Действительно, отрезок , ограничивающий полу полосу на плоскости p, переводится на плоскости z в окрестность: .

Обозначим через линию, в которую преобразование переводит отрезок . Тогда

Окрестность .

Т. о. Z – преобразование F(z) определено в области и является в этой области аналитической функцией.

Обратное Д – преобразование позволяет по изображению восстановить решетчатую функцию

(5)

Докажем справедливость равенства.

Лежат внутри окрестности.

(7)

(8)

В равенствах (7) и (8) вычеты берутся по всем особым точкам функции F(s).

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

\[Г (z+1) = zГ(z)\]

Для разностных уравнений с постоянными коэффициентами существуют детально разработанные методы нахождения решения в замкнутой форме. Неоднородное и однородное разностные уравнения n-го порядка задаются соответственно уравнениями, где \ постоянные коэффициенты.

Однородные разностные уравнения.

Рассмотрим уравнение n-го порядка

\[(a_nE^n +a{n-1}E^n1 + \cdots +a_1Е + a_1)y(k) = 0 \]

Предлагаемое решение следует искать в виде:

где \ - подлежащая определению постоянная величина. Вид предполагаемого решения, задаваемый уравнением, не является наиболее распространенным. Допустимые значения \ служат корнями многочлена от \[ е^r.\] При\[ \beta = е^r \]предполагаемое решение становится таким:

где \[\beta\] - подлежащая определению постоянная величина. Подставляя уравнение и учитывая \, получим следующее характеристическое уравнение:

Неоднородные разностные уравнения. Метод неопределенных коэффициентов. Рассмотрим разностное уравнение n-го порядка

\[ (a_nЕn +а_{n-1}Еn^-1+\cdots+ а_1Е +a_1)y(k) =F(k) \]

Ответ имеет следующий вид:

Где можно решить разностное уравнение онлайн?

Решить уравнение вы можете на нашем сайте https://сайт. Бесплатный онлайн решатель позволит решить уравнение онлайн любой сложности за считанные секунды. Все, что вам необходимо сделать - это просто ввести свои данные в решателе. Так же вы можете посмотреть видео инструкцию и узнать, как решить уравнение на нашем сайте. А если у вас остались вопросы, то вы можете задать их в нашей групе Вконтакте http://vk.com/pocketteacher. Вступайте в нашу группу, мы всегда рады помочь вам.