- •Требования к выполнению лабораторних работ
- •Постановка задачи
- •Содержание отчета
- •Теоретические сведения
- •Примеры выполнения заданий
- •Расчет в пакете mathcad
- •Постановка задачи
- •Содержание отчета
- •Постановка задачи.
- •Выводы. Теоретические сведения
- •Примеры выполнения заданий
- •Уточнение корня комбинированным методом.
- •Пример вычислений по методу хорд в пакете mathcad.
- •Пример реализации модифицированного метода
- •Варианты лабораторных работ
- •Лабораторная работа № 3 решение алгебраических и трансцендентных уравнений методом простой итерации
- •Постановка задачи
- •Содержание отчета
- •Постановка задачи.
- •Теоретические сведения
- •Примеры выполнения заданий
- •П ример решения трансцендентного уравнения в пакете mathcad.
- •Пример расчета в пакете matlab.
- •Встроенные функции пакетов mathcad, matlab для приближенного решения уравнений Пакет mathcad
- •Пакет matlab
- •Постановка задачи
- •Содержание отчета
- •Теоретические сведения.
- •Выводы. Теоретические сведения
- •Примеры выполнения заданий
- •Пример вычислений в пакете mathcad.
- •Программная реализация метода Гаусса для решения слау в пакете matlab.
- •Варианты лабораторных работ
- •Лабораторная работа № 5 итерационные методы решения системы линейных алгебраических уравнений (слау)
- •Постановка задачи
- •Содержание отчета
- •Постановка задачи.
- •Теоретические сведения.
- •Выводы. Теоретические сведения
- •Пример выполнения заданий
- •Решение системы методом простой итерации.
- •Пример расчета в пакете mathcad.
- •Программная реализация метода простой итерации для решения слау в пакете matlab.
- •Встроенные функции пакетов mathcad и matlab
- •Для приближенного решения систем линейных
- •Алгебраических уравнений
- •Пакет mathcad
- •Пакет matlab
- •Тема 3. Интерполяция и приближение полиномами Лабораторная работа № 6 интерполяционный многочлен лагранжа
- •Постановка задачи
- •Теоретические сведения.
- •Выводы. Теоретические сведения
- •Пример выполнения заданий
- •Построение полинома Лагранжа для таблично заданной функции.
- •Построение полинома Лагранжа в пакете mathcad.
- •Пример построения полинома Лагранжа в пакете matlab.
- •Пример использования узлов Чебышева.
- •Пример вычисления погрешности интерполирования.
- •Варианты лабораторных работ
- •Лабораторная работа № 7 интерполирование для таблиц с постоянным шагом. Численное дифференцирование. Обратное интерполирование
- •Постановка задачи
- •Содержание отчета
- •Постановка задачи.
- •Теоретические сведения.
- •Выводы. Теоретические сведения
- •Пример выполнения заданий
- •Построение первой формулы Ньютона в пакете mathcad.
- •Пример программ для вычислений по первой формуле Ньютона в matlab.
- •Обратное интерполирование с помощью полинома Ньютона в пакете mathcad.
- •Решение задачи обратного интерполирования
- •Варианты заданий
- •Лабораторная работа № 8 интерполирование кубическими сплайнами
- •Постановка задачи
- •Содержание отчета
- •Постановка задачи
- •Теоретические сведения.
- •Выводы. Теоретические сведения
- •Пример выполнения заданий
- •Построение кубического сплайна в пакете matlab.
- •Построение кубического сплайна в пакете mathcad.
- •Варианты лабораторных работ
- •Варианты лабораторных работ (окончание)
- •Встроенные функции интерполирования Пакет mathcad
- •Пакет matlab
- •Библиографический список
- •Содержание
- •Тема 1. Приближенное решение алгебраических
- •Тема 2. Приближенное решение систем линейных
- •Тема 3. Интерполяция и приближение полиномами……50
Примеры выполнения заданий
П ример решения трансцендентного уравнения в пакете mathcad.
Исходное уравнение:
– производная
В качестве параметра можно положить
Функция, реализующая метод итераций, может быть такой:
Пример расчета в пакете matlab.
% приведенное алгебраическое уравнение для метода итераций
Function res = p_iter_func(x,lambda);
res = x - lambda*(x^3+5*x+11);
return
% метод итераций
function res = p_iter(a,lambda,eps);
x0=a;
x1=p_iter_func(a,lambda);
i=0;
while (abs(x0-x1)>eps)
x0=x1;
x1=p_iter_func(x0,lambda);
i=i+1;
end;
i;
res=x1;
return
Варианты заданий (трансцендентное уравнение)
Номер варианта |
Уравнение |
Номер варианта |
Уравнение |
1 |
|
16 |
|
2 |
|
17 |
|
3 |
|
18 |
|
4 |
|
19 |
|
5 |
|
20 |
|
6 |
|
21 |
|
7 |
|
22 |
|
8 |
|
23 |
|
9 |
|
24 |
|
10 |
|
25 |
|
11 |
|
26 |
|
12 |
|
27 |
|
13 |
|
28 |
|
14 |
|
29 |
|
15 |
|
30 |
|
Варианты для выполнения лабораторной работы № 3 (алгебраическое уравнение) взять из лабораторной работы № 2.
Встроенные функции пакетов mathcad, matlab для приближенного решения уравнений Пакет mathcad
Для поиска корней обычного многочлена p(x) степени n MATHCAD имеет очень удобную функцию: polyroot(V). Она возвращает вектор корней многочлена степени n, коэффициенты которого находятся в векторе v, имеющем длину, равную n+1.
Заметим, что корни многочлена могут быть как вещественными, так и комплексными числами. Не рекомендуется пользоваться этой функцией, если степень полинома выше пятой-шестой, так как тогда трудно получить малую погрешность вычисления корней.
Для решения уравнений с одним неизвестным MATHCAD имеет встроенную функцию root, которая, в зависимости от типа задачи, может содержать либо два, либо четыре аргумента и, соответственно, работает несколько по-разному: root(f(х), х); root(f(х), х ,а ,b),
где f(х) – скалярная функция, определяющая уравнение, х – скалярная переменная, относительно которой решается уравнение; а,b – границы интервала, внутри которого происходит поиск корня.
Первый тип функции root требует дополнительного задания начального значения (guess value) переменной х. Для этого нужно просто предварительно присвоить х некоторое число. Поиск корня будет производиться вблизи этого числа. Таким образом, присвоение начального значения требует априорной информации о примерной локализации корня.
Если правильное завершение работы функции root отсутствует, то это значит, что :
– нет корней;
– задание начальных условий некорректно (начальное условие слишком далеко от корней, несоответствие типов);
– на пути к корню встретился локальный экстремум;
– встретилась точка разрыва функции.
Для задания начальных условий удобнее всего построить график функции. Если в окрестности корня функция имеет малую кривизну, то экстремум рекомендуется искать в виде
.
Для нахождения второго корня необходимо либо задать другое начальное приближение, либо использовать вызов
,
где а – уже вычисленный корень.
Примеры использования функции root:
1) x:=nu root (F(x),x)=...;
2) если далее корень понадобится в выражениях, то, например, xv:=root (F(x),x).
Если точность мала, то изменением значения TOL-погрешности вычислений следует задать требуемую точность, например, TOL:=0,0001.
С помощью символьных преобразований также можно решать уравнения, неравенства.
Пример1. Найти корни уравнения .
Решение.
Ввести заданное уравнение. В качестве знака равенства следует использовать оператор отношения = с панели инструментов Boolean.
На панели инструментов Symbolic выбрать оператор Solve ■ →.. В качестве требуемой переменной ввести x.
Щелкнуть левой кнопкой мыши за пределами введенной формулы. Справа от оператора → появится результирующее значение:
.
Заметим, что данная последовательность действий остается верной и при наличии в уравнении параметра, например
.