На практике простейшие разностные уравнения возникают при исследовании например величины банковского вклада. Эта величина является переменной 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. Основные понятия и примеры разностных уравнений

Разностные уравнения играют большую роль в экономической теории. Многие экономические законы доказывают с помощью именно этих уравнений. Разберем основные понятия разностных уравнений.

Пусть время t выступает как независимая переменная, а зависимая переменная определяется для времени t, t-1, t-2 и т.д.

Обозначим через значение в момент времени t; через - значение функции в момент, сдвинутый назад на единицу (например, в предыдущем часу, на предыдущей неделе и т.д.); через - значение функции y в момент, сдвинутый на две единицы назад, и т.д.

Уравнение

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

Уравнение

В котором =0, называется разностным однородным уравнением n-го порядка с постоянными коэффициентами. Решить разностное уравнение n-го порядка - значит найти функцию, которая обращает это уравнение в верное тождество.

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

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

где и - произвольные постоянные.

Теорема 2. Если - частное решение неоднородного разностного уравнения (1) и - общее решение однородного уравнения (2), то общим решением неоднородного уравнения (1) будет функция

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

где - вектор из неизвестных функций, - вектор из известных функций.

Есть матрица размера nn.

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

§ 2. Решение разностных уравнений

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

Соответствующее однородное уравнение есть

Проверим, будет ли функция

решением уравнения (3).

Подставляя в уравнение (4), получаем

Следовательно, есть решение уравнения (4).

Общее решение уравнения (4) есть функция

где C - произвольная постоянная.

Пусть - частное решение неоднородного уравнения (3). Тогда общее решение разностного уравнения (3) есть функция

Найдем частное решение разностного уравнения (3), если f(t)=c, где c - некоторая переменная.

Будем искать решение в виде постоянной m. Имеем

Подставив эти постоянные в уравнение

получаем

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

Пример1 . Найти с помощью разностного уравнения формулу прироста денежного вклада А в сбербанке, положенного под p % годовых.

Решение . Если некоторая сумма положена в банк под сложный процент p, то к концу года t её размер составит

Это однородное разностное уравнение первого порядка. Его решение

где C - некоторая постоянная, которую можно рассчитать по начальным условиям.

Если принять, то C=A, откуда

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

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

и соответствующее однородное уравнение

Если k является корнем уравнения

есть решение однородного уравнения (6).

Действительно, подставляя в левую часть уравнения (6) и учитывая (7), получаем

Таким образом, если k - корень уравнения (7), то - решение уравнения (6). Уравнение (7) называется характеристическим уравнением для уравнения (6). Если дискриминант характеристическое уравнение (7) больше нуля, то уравнение (7) имеет два разных действительных корня и, а общее решение однородного уравнения (6) имеет следующий вид.

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

Однако мы постараемся вам показать, что диффуры – это не так сложно, как кажется.

Основные понятия теории дифференциальных уравнений

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

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

Дифференциальное уравнение (ДУ ) – это уравнение, содержащее производные функции y(х), саму функцию, независимые переменные и иные параметры в различных комбинациях.

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

Решением дифференциального уравнения является функция, которая обращает его в тождество. Существуют общие и частные решения ДУ.

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

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

Обыкновенные дифференциальные уравнения

Обыкновенные дифференциальные уравнения – это уравнения, содержащие одну независимую переменную.

Рассмотрим простейшее обыкновенное дифференциальное уравнение первого порядка. Оно имеет вид:

Решить такое уравнение можно, просто проинтегрировав его правую часть.

Примеры таких уравнений:

Уравнения с разделяющимися переменными

В общем виде этот тип уравнений выглядит так:

Приведем пример:

Решая такое уравнение, нужно разделить переменные, приведя его к виду:

После этого останется проинтегрировать обе части и получить решение.

Линейные дифференциальные уравнения первого порядка

Такие уравнения имеют вид:

Здесь p(x) и q(x) – некоторые функции независимой переменной, а y=y(x) – искомая функция. Приведем пример такого уравнения:

Решая такое уравнение, чаще всего используют метод вариации произвольной постоянной либо представляют искомую функцию в виде произведения двух других функций y(x)=u(x)v(x).

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

Пример решения ДУ с разделяющимися переменными

Вот мы и рассмотрели простейшие типы ДУ. Теперь разберем решение одного из них. Пусть это будет уравнение с разделяющимися переменными.

Сначала перепишем производную в более привычном виде:

Затем разделим переменные, то есть в одной части уравнения соберем все "игреки", а в другой – "иксы":

Теперь осталось проинтегрировать обе части:

Интегрируем и получаем общее решение данного уравнения:

Конечно, решение дифференциальных уравнений – своего рода искусство. Нужно уметь понимать, к какому типу относится уравнение, а также научиться видеть, какие преобразования нужно с ним совершить, чтобы привести к тому или иному виду, не говоря уже просто об умении дифференцировать и интегрировать. И чтобы преуспеть в решении ДУ, нужна практика (как и во всем). А если у Вас в данный момент нет времени разбираться с тем, как решаются дифференциальные уравнения или задача Коши встала как кость в горле или вы не знаете, как правильно оформить презентацию , обратитесь к нашим авторам. В сжатые сроки мы предоставим Вам готовое и подробное решение, разобраться в подробностях которого Вы сможете в любое удобное для Вас время. А пока предлагаем посмотреть видео на тему "Как решать дифференциальные уравнения":

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

y(k) = F(k) (92)

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

Вспомогательное решение.

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

a 1 (k)y(k+1) + a 0 (k)y(k) = 0, (93)

где a 0 (k)≠0 и a 1 (k)≠0. Его можно переписать в виде

y(k+1) = a(k)y(k). (94)

при k=0,1,2...

у(1)=а(0)у(0),

у(2)=а(1)а(0)у(0)

у(3)=а(2)а(1)а(0)у(0)

или, в общем случае,

так что общее решение уравнения (94) равно

Нижний предел произведения произволен, так как любое фик­сированное число множителей а(0), а(1), и а(2), ... можно объединить с произвольной постоянной С.

Решение однородного уравнения выше первого порядка в общем случае не выражается в виде элементарных функций, так как процедура, основанная на уравнениях (81) и (82), при зависящих от k коэффициентах перестает быть справедливой. Если известны все, кроме одного, независимые решения уравнения, то можно определить и ос­тавшееся решение. Как и для дифференциальных уравнений, в ряде отдель­ных случаев удается получить, решение в явном виде. Уравнение вида

a n f(k + n)y(k + n) + ... + a 1 f (k + 1)y(k + 1) + a n f(k)y(k) = 0,

где коэффициенты а i - постоянные величины, путем подстановки z(k)=f(k)y(k) сводится к разностному уравнению с постоянными коэффициентами. Проце­дура отчасти сходна с используемой для дифференциального уравнения Эй­лера, но замене в данном случае подлежит зависимая (а не независимая) пе­ременная. Этот метод широко используется при решении уравнений с пе­ременными коэффициентами.

Дифференциальные уравнения систем автоматического регулирования. Методика составления дифференциальных уравнений систем автоматического регулирования.

Общие замечания.

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

Любая САР представляет совокупность отдельных взаимодействующих друг с другом элементов, соединенных между собой связями. Первым этапом при составлении дифференциальных уравнений САР является разделение системы на отдельные элементы и составление дифференциальных уравнений этих элементов. Уравнения элементов и уравнения связей между отдельными элементами описывают процесс в системе регулирования, т.е. изменение по времени всех координат системы. Зная уравнения элементов и уравнения связей, можно составить структурную схему САР.

Структурная схема САР характеризует геометрию системы, т.е. показывает, из каких элементов состоит САР и как эти элементы связаны между собой. Состояние САР, а также каждого входящего в нее элемента характеризуется некоторым числом независимых переменных. Этими переменными могут быть как электрические величины (ток, напряжение, и т.д.), так и механические (скорость, угол поворота, перемещение и т.д.). Обычно, чтобы характеризовать состояние системы или ее элемента, выбирают одну обобщенную координату на входе системы или элемента (g(t)) и одну – на выходе (x(t)). В ряде случаев такое представление невозможно, так как система или ее элемент могут иметь несколько входных и выходных величин. В многомерных системах можно рассматривать векторные входную и выходную величины с размерностями, совпадающими соответственно с числом входных и выходных величин САР.

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

При составлении дифференциальных уравнений САР основной задачей является составление дифференциальных уравнений отдельных элементов системы. Уравнение отдельных элементов составляются на основе тех физических законов, которые характеризуют поведение элемента.

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

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

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

Если дифференциальное уравнение элемента нелинейно из-за нелинейности его статической характеристики, то линеаризация уравнения сводится к замене нелинейной характеристики элемента x =φ(g ) некоторой линейной функцией x = ag + b . Аналитически эта замена производится с помощью разложения в ряд Тейлора функции x =φ(g ) в окрестности точки, соответствующей установившемуся состоянию и отбрасывания всех членов содержащих отклонение ∆g входной величины элемента в степени выше первой. Геометрически это означает замену кривой x =φ(g ) касательной, проведенной к кривой в точке (x 0 , g 0), соответствующей установившемуся состоянию работы элемента (рис. 29). В других случаях линеаризация производится путем проведения секущей, мало отклоняющейся от функции x =φ(g ) в требуемом диапазоне изменения входной величины элемента.

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

Рассмотрим процесс линеаризации нелинейного уравнения элемента с помощью ряда Тейлора. Пусть поведение элемента описывается нелинейным дифференциальным уравнением

F(x n , x ’ , x, g) = 0 (1). Тогда установившееся состояние элемента характеризуется уравнением F(0, 0, x, g) = 0 (2). пусть g 0 и х 0 – значения установившегося состояния. Тогда координаты g и х можно записать в виде x = x 0 + ∆x, g = g 0 + ∆g, где ∆g и ∆x – отклонение координат g и х от установившегося состояния. Уравнение (1) в отклонениях имеет вид:

F(∆x ’’ , ∆x ’ , x 0 + ∆x, g 0 + ∆g) = 0 (3).

Разложим левую часть уравнения (3) в ряд Тейлора относительно точки установившегося состояния (0, 0, x 0 , g 0):

Частные производные в левой части уравнения (4) представляют собой некоторые числа, величины которых зависят от вида функции F(x ’’ , x ’ , x, g) и значений координат x 0 и g 0 .

Считая отклонения ∆g, ∆x от установившегося состояния, а также их производные по времени малыми и пологая, что функция F(x ’’ , x ’ , x, g) достаточно гладкая по всем аргументам в окрестности точки, соответствующей установившемуся состоянию, отбросим в уравнении (4) все члены, которые содержат отклонения ∆g и ∆x, а также их производные в степени выше первой. Полученное уравнение (5) является линейным дифференциальным уравнением с постоянными коэффициентами,,,и представляет собой результат линеаризации уравнения (1).

Очевидно, что необходимым условием линеаризации является возможность разложения в ряд Тейлора функции F(x ’’ , x ’ , x, g) в окрестности точки, соответствующей установившемуся состоянию.

Процесс линеаризации уравнения (1) может быть геометрически интерпретирован следующим образом. В пространстве переменных x ’’ , x ’ , x, g уравнение (1) задает некоторую поверхность. Переход от уравнения (1) к линейному уравнению (5) означает замену поверхности некоторой касательной плоскостью, проведенной к поверхности в точке, соответствующей установившемуся состоянию. Естественно, что ошибка при такой замене тем меньше, чем меньше отличаются друг от друга точки поверхности и точки плоскости. Это справедливо лишь в некоторой малой окрестности установившегося состояния.

Понятие управляемости и наблюдаемости.

Процесс или объект принято называть полностью управляемым, если он может быть переведен из некоторого состояния х(t 0) в желаемое состояние равновесия х(t 1) за конечный интервал времени t 1 – t 0 . Другими словами, процесс является полностью управляемым, если существует управляющее воздействие m(t), определенное на конечном интервале времени t 0 ≤ t ≤ t 1 , которое переводит процесс из начального состояния х(t 0) в желаемое состояние равновесия х(t 1) за время t 1 – t 0 .

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

Линейный дискретный процесс n-го порядка является полностью управляемым тогда и только тогда, когда векторы

s 1 = φ(-T)h(T),

s 2 = φ(-T)h(T),

s n = φ(-T)h(T)

линейно независимы.

Эти векторы возникают в связи со следующими преобразованиями.

(t) = Ax(t) + d m(t),

в котором m(t) – единственное управляющее воздействие. Случай единственного управляющего воздействия рассматривается ради упрощения интерпретации получаемых выражений. Уравнение переходных состояний процесса имеет вид

где φ(Т) – матрица перехода процесса и
.

Понятию управляемости можно дать еще и другое толкование, способствующее лучшему его пониманию. Пусть линейный многомерный процесс описывается векторным дифференциальным уравнением (t) = Ax(t) + D m(t), где х – n-мерный вектор состояния;

m – r-мерный вектор, представляющий управляющие воздействия;

А – квадратичная матрица коэффициентов n-го порядка;

D – матрица управления размера n×r.

Матрица А может быть приведена к диагональной форме

,

где λ i – собственные значения матрицы А линейного процесса, которые предполагаются все различными.

Применяя подстановку x=Tz, уравнение запишем в канонической форме

(t) = Λz(t) + ∆ m(t),

где
. Векторz будем называть каноническим вектором состояния.

Процесс, описываемый уравнением (t) = Ax(t) + D m(t), является управляемым, если матрица ∆ не содержит строк, все элементы которых равны нулю; координаты, соответствующие ненулевым строкам ∆, считаются управляемыми.

Пример:

Вывести дифференциальное уравнение центробежного маятника, который применяется в качестве чувствительного элемента в некоторых САР. Схема маятника изображена на рисунке. Входной величиной является угловая скорость ω, а выходной величиной – перемещение х платформы. При увеличении скорости вращения шары под действием центробежной силы расходятся и перемещают платформу. На платформу воздействует также сила упругости пружины, силы демпфирования и силы инерции.

Введем обозначения: с – коэффициент жесткости пружины; k – коэффициент вязкого трения; m – масса шара; М – масса частей, участвующих в поступательном движении вдоль оси ОХ; ω – угловая скорость вала; f 0 – сила предварительного поджатия пружины.

Для составления дифференциального уравнения центробежного маятника используем уравнение Лагранжа второго рода:
(I = 1, 2,…, n) (*). В качестве обобщенной координаты x i выберем выходную координату – перемещение платформы х. Найдем выражение для кинетической энергии Т, потенциальной энергии П и диссипативной функции R центробежного маятника. Из рисунка видно, что

ρ = r + l sin α, x = 2a(1 – cos α).

Кинетическая энергия системы Т = Т 1 + Т 2 + Т 3 , где Т 1 – кинетическая энергия во вращательном движении вокруг оси ОХ; Т 2 – кинетическая энергия шаров во вращательном движении вокруг точек А и А’; Т 3 – кинетическая энергия масс в поступательном движении вдоль оси ОХ. Имеем:

,

,
. (*1)

Потенциальная энергия маятника П = П 1 + П 2 + П 3 , где П 1 – потенциальная энергия масс, движущихся параллельно оси ОХ; П 2 – потенциальная энергия; П 3 – потенциальная энергия пружины. Для рассматриваемого случая имеем:

,
,
. (*2)

Найдем обобщенную диссипативную силу Q R . Благодаря наличию демпфера сила сухого трения мала по сравнению с силой вязкого трения и ею можно пренебречь. Согласно формуле
будем иметь

. (*3)

Вычислим значение отдельных слагаемых, входящих в уравнение Лагранжа (*):

,

,

.

Подставим полученные выражения в уравнение Лагранжа второго рода (*), тогда

Введем следующее обозначения:

,
,

; (*5)

. (*6)

С учетом принятых обозначений уравнение центробежного маятника запишется в виде

Уравнение (*7) представляет собой нелинейное дифференциальное уравнение. Состояние равновесия (х 0 , ω 0) является решением уравнения

Рассмотрим малые колебания маятника относительно состояния равновесия

х = х 0 + ∆х, ω = ω 0 + ∆ω. (*9)

Разложим функции f 1 (x), f 2 (x), f 3 (x, ω) в ряд Тейлора в окрестности состояния равновесия (х 0 , ω 0).

где функции F 1 (∆x), F 2 (∆x), F 3 (∆x, ∆ω) имеют более высокий порядок малости по сравнению с ∆x и ∆ω. Учитывая, что x’ = ∆x’ и x” = ∆x”, и принимая во внимание выражения (*8), (*9), (*10), уравнение (*7) можно переписать в виде

где функция

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

, (*11)

,

(*12)

.

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

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

Обычно вместо уравнения (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).