Выполнение расчетов с использованием пакета программ mathcad

ВЫПОЛНЕНИЕ РАСЧЕТОВ С ИСПОЛЬЗОВАНИЕМ ПАКЕТА ПРОГРАММ MATHCAD

Цель работы:

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

Индивидуальное задание

p1 – последняя цифра в зачетной книжке студента;

р2– предпоследняя цифра в зачетной книжке студента;

р3 – число букв в фамилии студента.

Постановка задачи

На лесной опытной станции были произведены 20 измерений диаметра дерева D на высоте груди (1.3м) и высоты H. Результаты измерений сведены в таблицуТребуется решить задачу интерполяции – для любого значения из интервала от Dmin до Dmax найти значение высоты H.

Интерполяция кубическими сплайнами

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

1. Исходные данные требуется представить в виде матрицы, где каждый из двух столбцов - это вектор значений Di и Hi. Затем, используя функцию V=csort(M,N),отсортировать значения матрицы по столбцу Di в порядке возрастания. В данном формате функции csort(M,N) M– обозначение матрицы, N– номер столбца, по которому производится сортировка.

  1.  
      выполнение расчетов с использованием пакета программ mathcad - student2.ru

    Использование функции s=cspline(x,y) на векторах x и y возвращает вектор s, содержащий значения вторых производных сглаживаемой кривой в заданных точках. В качестве векторов x и yследует задатьзначения Di и Hi, что соответствует столбцам V<0>и V<1>:
  1. Использование функции interpпозволяет найти значения функции q(z) для промежуточных значений zиз интервала от Dmin до Dmax:

выполнение расчетов с использованием пакета программ mathcad - student2.ru
q(z)=interp(s,x,y,z)

Рис. 3.1

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

Приближение функции с использованием метода наименьших квадратов

Время выполнения работы 4 часа.

4.1. Индивидуальное задание

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

Y¢(t) = 2tY(t)

Требуется найти численное решение данного дифференциального уравнения с заданными начальными условиями (задача Коши) Yt=0 = 1 млн. шт.

Задачу решить:

· методом Эйлера Y1(t);

· методом Эйлера-Коши Y2(t);

· с использованием формулы Тейлора второго порядка точности Y3(t).

Сравнить полученные решения с точным (аналитическим) Y(t)=exp(t2) графическим способом и с помощью полученных числовых значений для значений tот0 до 2 (за единицу измерения времени t принято 10 лет).

4.2. Расчетные формулы

Метод Эйлера

Расчетная формула

Yi+1 = Yi + h × f(ti,Yi),

где h – шаг численного решения (h=0.01),

а f(ti,Yi) = Y¢ = 2tY.

Метод Эйлера-Коши

Расчетная формула

Yi+1 = Yi + h/2 × (f(ti,Yi) + f(ti+1, Yi + h × f(ti,Yi))).

где h – шаг численного решения (h=0.01),

Использование формулы Тейлора второго порядка точности

Расчетная формула

Yi+1 = Yi + h × Y¢(t) + h2/2× Y¢¢(t),

где h – шаг численного решения (h=0.01).

Дифференцируя исходное уравнение по t, получим следующее выражение для второй производной:

Y¢¢ = 2Y + 2t Y¢ = 2(1 + 2t2) Y

Принцип золотого сечения

       
  выполнение расчетов с использованием пакета программ mathcad - student2.ru   выполнение расчетов с использованием пакета программ mathcad - student2.ru
 

Основной принцип золотого сечения отражен в следующемсоотношении:

Рис.5.2

выполнение расчетов с использованием пакета программ mathcad - student2.ru
Это правило положено в основу уменьшения отрезка локализации.

Рис.5.3

Исходный отрезок [a0,b0], на котором ищется решение, разбивается двумя точками i0 и j0 по правилу золотого сечения:

       
  выполнение расчетов с использованием пакета программ mathcad - student2.ru   выполнение расчетов с использованием пакета программ mathcad - student2.ru

Правило локализации (уменьшения отрезка) следующее:

если  
 
  выполнение расчетов с использованием пакета программ mathcad - student2.ru

то   выполнение расчетов с использованием пакета программ mathcad - student2.ru
если   выполнение расчетов с использованием пакета программ mathcad - student2.ru то   выполнение расчетов с использованием пакета программ mathcad - student2.ru

В MathCADе данное правило записывается следующим образом:

выполнение расчетов с использованием пакета программ mathcad - student2.ru

Это иллюстрирует следующий рисунок.

выполнение расчетов с использованием пакета программ mathcad - student2.ru

Рис.5.4

выполнение расчетов с использованием пакета программ mathcad - student2.ru Для нахождения точки локального минимума x с заданной точностью необходимо проделать n итераций. Каждая итерация сокращает длину отрезка локализации в раз. Поэтому

 
  выполнение расчетов с использованием пакета программ mathcad - student2.ru

Оценка погрешности определяется:

 
  выполнение расчетов с использованием пакета программ mathcad - student2.ru

При достижении заданной точности (e<0.1) итерации следует прекратить и положить

выполнение расчетов с использованием пакета программ mathcad - student2.ru

ВЫПОЛНЕНИЕ РАСЧЕТОВ С ИСПОЛЬЗОВАНИЕМ ПАКЕТА ПРОГРАММ MATHCAD

Цель работы:

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

Индивидуальное задание

p1 – последняя цифра в зачетной книжке студента;

р2– предпоследняя цифра в зачетной книжке студента;

р3 – число букв в фамилии студента.

Наши рекомендации