Уральский государственный экономический университет
Опубликован: 27.05.2014 | Доступ: свободный | Студентов: 467 / 49 | Длительность: 11:44:00
Лекция 1:

Математические методы в моделировании экономики

Лекция 1: 1234 || Лекция 2 >

Вычисление пределов

Для нахождения предела функции Mathcad используется панель Calculus. Ввести функцию под знак \lim_{a \to b} . Обязательно выделить его (клавиша "пробел"), щелкнуть оператор символьного вывода (панель Symbolic или Evaluation).

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

Пример 1.6.

В теории потребления исследуется распределение доходов населения между разными видами расходов. Функции спроса описывают зависимость спроса на товары определенной категории от размера поучаемого дохода. Известны функции Торнквиста, которые описывают зависимость величины спроса на различные группы товаров в зависимости от их роли в процессе потребления. Пределы функций при неограниченном росте доходов показывают насыщение спроса в каждом случае. Пусть Y –спрос, x – доход, a,b,c – фиксированные параметры. Функции Торнквиста:

  • спрос на малоценные товары:
    Y_0=\frac{a\cdot x \cdot (x+c)}{x^2+b} ( 1.7)
  • спрос на товары первой необходимости:
    Y_1=\frac{a\cdot x}{x+с} ( 1.8)
  • спрос на товары второй необходимости (более дорогие):
    Y_2=\frac{a\cdot (x-b)}{x+c} ( 1.9)
  • спрос на товары роскоши:
    Y_3=\frac{a\cdot x \cdot (x-b)}{x+c} ( 1.10)

Примеры вычислений пределов при исследовании спроса на различные виды товаров показаны на рисунке 1.4.

Определение пределов  функций Торнквиста.

Рис. 1.4. Определение пределов функций Торнквиста.

Функции Торнквиста

Спрос на различные виды товаров

a:=20, b:=60, c:=20

Малоценные товары:

Y_0=\frac{a\cdot x \cdot (x+c)}{x^2+b}

\lim_{x \to \infty}Y(x)\to 20

Товары первой необходимости

Y1(x)=\frac{a\cdot x}{x+с}

\lim_{x \to \infty}Y1(x)\to 20 - При увеличении дохода – насыщение

Товары длительного пользования, спрос возникает лишь с некоторого высокого уровня дохода

Y2(x)=\frac{a\cdot (x-b)}{x+c}

\lim_{x \to \infty}Y2(x)\to 20

Y2(x)\; solve,x\to 60

Товары роскоши. Спрос растет неограниченно

Y3(x)=\frac{a\cdot x \cdot (x-b)}{x+c}

\lim_{x \to \infty}Y3(x)\to \infty

Y3(x)\; solve, x \to \begin{pmatrix} 60 \\ 0 \end{pmatrix}

Разложение функции в ряд Тейлора

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

В Mathcad задача разложения функции f(x) в ряд Тейлора в окрестности точки x= x_0 решается символьной операцией с использованием ключевого слова series панели Simbolic. После ключевого слова вводится значение точки разложения с логическим равенством и порядок остаточного члена минус единица.

Если функция f(x) имеет в точке х = x_0 и некоторой ее окрестности производные порядка до (n+1) включительно, то функцию можно разложить в ряд Тейлора для любого значения х этой окрестности

f(x)=f(x_0)+\frac{f'(x_0)}{1!}\cdot (x-x_0)+\frac{f''(x_0)}{2!}\cdot (x-x_0)^2+...+\frac{f^{(n)}(x_0)}{n!}\cdot (x-x_0)^n+\frac{f^{(n)}(x_0)}{n!}\cdot (x-x_0)^n ( 1.11)

(n+1) член разложения:

R_{n+1}(x)=\frac{f^{(n)}(x_0)}{n!}\cdot (x-x_0)^n ( 1.12)
определяет ошибку и называется остаточным членом в форме Лагранжа.

Пример. 1.7.

На рисунке 1.5 показан пример разложения функции f(x)=e^{-x}\cdot \tg (x) в ряд Тейлора в окрестности точки (x=\pi /4) с остаточным членом 3, 4 и 5 порядка. Для сравнения рассчитаны значения функции в точке x=1 по степенным формулам с разным порядком разложения. Точность до второго знака позволяет ограничиться 4 членами разложения.

Следует отметить, запись символьной операции с ключевыми словами serias,x=\pi /4 (разложить в точке (x=\pi /4) и substitute x=1 (подставить x=1) идет в одну строку и не переносится. Рисунок листинга является демонстрационным.

Разложение  функции с остаточным членом 3 и 4 и 5 порядка

Рис. 1.5. Разложение функции с остаточным членом 3 и 4 и 5 порядка

f(x)=e^{-x}\cdot \tg (x)

f(1)=0.573, \frac{\pi}{4}=0.785

f(x)\; series,x=\frac{\pi}{4},3\to e^{-\frac{\pi}{4}}+e^{-\frac{\pi}{4}}\cdot (x-\frac{\pi}{4})+\frac{e^{-\frac{\pi}{4}}\cdot (x-\frac{\pi}{4})^2}{2}substitute,x=1\to \frac{e^{-\frac{\pi}{4}}\cdot (\pi^2-16\cdot \pi+80)}{32}=0.564

f(x)\; series,x=\frac{\pi}{4},4\to e^{-\frac{\pi}{4}}+e^{-\frac{\pi}{4}}\cdot (x-\frac{\pi}{4})+\frac{e^{-\frac{\pi}{4}}\cdot (x-\frac{\pi}{4})^2}{2}+\frac{3\cdot e^{-\frac{\pi}{4}}\cdot (x-\frac{\pi}{4})^3}{2}substitute,x=1\to \frac{e^{-\frac{\pi}{4}}\cdot (200\cdot \pi-40\cdot \pi^2+3\cdot \pi^3-512)}{128}=0.571

f(x)\; series,x=\frac{\pi}{4},5\to e^{-\frac{\pi}{4}}+e^{-\frac{\pi}{4}}\cdot (x-\frac{\pi}{4})+\frac{e^{-\frac{\pi}{4}}\cdot (x-\frac{\pi}{4})^2}{2}+\frac{3\cdot e^{-\frac{\pi}{4}}\cdot (x-\frac{\pi}{4})^3}{2}+\frac{11\cdot e^{-\frac{\pi}{4}}\cdot (x-\frac{\pi}{4})^4}{8}substitute,x=1\to \frac{e^{-\frac{\pi}{4}}\cdot (1696\cdot \pi^2-6144\cdot \pi-224\cdot \pi^3+11\pi^4+11008)}{2048}=0.572

Лекция 1: 1234 || Лекция 2 >