1. Численный поиск корня уравнения.
Любой численный метод решения уравнения основан на уточнении какого-либо начального значения корня до заданной точности. Точность вычислений в MathCAD задается встроенной переменной TOL и по умолчанию равна 0,001. Для задания начального значения корня удобно сначала построить график функции, задающей уравнение. Например, найдем корни уравнения . Для задания начального значения x построим график функции .
Для нахождения правого корня уравнения зададим точность TOL=0,000001, начальное значение x=4 и используем функцию root(f(x),x).
Можно локализовать корень не с помощью начального значения, а с помощью интервала, которому принадлежит корень, при этом границы интервала нужно задать как аргументы функции root. Например, найдем левый корень исходного уравнения как корень, принадлежащий интервалу [0.01, 1].
Можно задать функцию r(x)=root(f(x),x), которая возвращает корень уравнения, полученный из начального приближения x. Например:
Пользуясь этой функцией можно получить вектор корней уравнения:
Для получения комплексного корня начальное приближение следует задавать комплексным.
2. Нахождение корней полиномов.
Для нахождения корней полинома в MathCAD имеется встроенная функция polyroots(a), аргументом которой является вектор коэффициентов полинома . Например, для уравнения вектор а имеет вид:
Коэффициенты полинома и корни могут быть и комплексными. Например для уравнения
Для того, чтобы выделить из многочлена вектор коэффициентов можно воспользоваться символьными преобразованиями. Например решим уравнение :
3. Нахождение корней уравнений аналитически (путем символьных преобразований).
|
Во многих случаях MathCAD позволяет найти аналитическое решение. Для того, чтобы найти решение уравнения необходимо записать выражение и выделить в нем переменную (поставить указатель курсора возле переменной), и воспользоваться пунктом Solve for Variable (Переменная®Решение) из пункта меню Symbolic (Символика).
Можно также воспользоваться функцией root, поставив вместо знака «=» знак «®».
Например:
При этом для первого уравнения найден только 1 корень, хотя их бесчисленно много, а для полинома найдены все 3 корня, в том числе и комплексные.
4. Численный поиск решения системы уравнений и неравенств.
Системы линейных и нелинейных уравнений и неравенств позволяет решать в MathCAD блок given в сочетании с функцией Find.
После слова given записывается система уравнений и неравенств, подлежащих решению, при этом вместо знака «=» следует набирать Ctrl+=. Перед блоком given необходимо указывать начальные приближения для всех переменных, если нужно найти комплексный корень, то следует задавать комплексное начальное приближение. Признаком окончания системы служит функция Find, если надо найти точное решение или функция Minerr, если система не может быть решена точно, и требуется найти наилучшее приближение, обеспечивающее минимальную погрешность.
Функции Minerr и Find должны иметь столько же или меньше аргументов, сколько уравнений и неравенств содержит блок given. Если окажется, что блок содержит слишком мало уравнений или неравенств, то его можно дополнить тождествами или повторяющимися выражениями.
Например,
|
Если необходимо найти решение при различных начальных приближениях, имеет смысл определить новую функцию. Например:
Если система не имеет решения, но возникает необходимость найти значения переменных, при которых уравнения, входящие в систему удовлетворяются хотя бы приближенно, то используется функция miner. Например:
5. Решение систем линейных уравнений.
Систему линейных уравнений можно решать численным методом, описанном в предыдущем пункте, но если определитель матрицы из коэффициентов при неизвестных в уравнениях системы будет равен 0 (система не имеет решения либо имеет множество решений), то численный метод не даст результата. Также численный метод может дать приближенный результат вместо точного. Поэтому для решения систем линейных уравнений можно воспользоваться методом Гаусса, матричным методом или формулами Крамера. Однородную же систему линейных уравнений можно решить только методом Гаусса. Рассмотрим решение систем линейных уравнений всеми этими способами:
а) численный метод
б) метод Гаусса
в) матричный метод
г) метод Крамера
Рассмотрим решение однородной системы методом Гаусса:
6. Решение систем уравнений в символьном виде.
Во многих случаях решение системы уравнений может быть найдено не только численно, но и аналитически. Для этого также используется блок given и функция Find, но вместо знака равенства после функции следует поставить знак символического преобразования ®.
Например:
Решение записано в вид матрицы, каждый столбец которой является решением.
|
Решить систему уравнений в символьном виде можно также с помощью оператора solve. Для этого систему уравнений записывают в виде матрицы, состоящей из 1 столбца и нужного количества строк. Например:
В данном случае решением является каждая строка матрицы.
Контрольные задания.
1. Дано уравнение. Построить график для определения начальных значений корня. Определить один из корней с помощью задания начального значения, другой с помощью задания интервала, которому принадлежит корень. Найти вектор корней уравнения. Все вычисления производить с точностью 0.00001.
1), 10)
2), 11)
3), 12)
4), 13)
5), 14)
6), 15)
7), 16)
8), 17)
9), 18)
2. Найти все корни полиномиального уравнения:
1)
2)
3)
4)
5)
6)
7)
8)
9)
10)
11)
12)
13)
14)
3. Решить аналитически уравнение, заданное в первом задании.
4. Найти одно из решений системы уравнений численно:
5. Найти приближенное решение системы уравнений, не имеющей точного решения, построить график, показав на нем найденное решение.
1), 7), 13)
2), 8), 14)
3), 9), 15)
4), 10), 16)
5), 11), 17)
6), 12), 18)
6.Решить систему линейных уравнений численно, методом Гаусса, матричным методом, по формулам Крамера и с помощью символьных преобразований:
7. Решить однородную систему линейных уравнений:
1), 7), 13)
2), 8), 14)
3), 9), 15)
4), 10), 16)
5), 11), 17)
6), 12), 18)