рефераты конспекты курсовые дипломные лекции шпоры

Реферат Курсовая Конспект

Численные методы.

Численные методы. - раздел Математика, Обзор возможностей математических пакетов MathCAD 2000, MathLAB 5.0, Mathematica Значения Элементарных Функций И Многочисленных Специальных Математических Фун...

Значения элементарных функций и многочисленных специальных математических функций в вещественных и комплексных точках с вещественными координатами, можно найти, просто вычислив соответствующие выражения. Найдем, например, значение функции синус в точке 0.7 и значение функции Бесселя (x) первого рода нулевого порядка в комплексной точке 1.+0.5I:

Sin[0.7]

BesselJ[0, 1. + 0.5 I]

Команды символьных вычислений Sum,Integrate, DSolve и т.п. имеют своих двойников NSum, NIntegrate, NDSolve и т.д., выполняющих приближенные численные вычисления.

Integrate[Sin[x^4 + x], {x, 1, 2}]

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

NIntegrate[Sin[x^4 + x], {x, 1, 2}, WorkingPrecision -> 25]

Приближение дискретных данных с наименьшим среднеквадратичным отклонением с помощью задаваемой пользователем системы функций осуществляет команда Fit. В рассмотренном ниже примере порождается 10 случайных точек на координатной плоскости и находятся наилучшие приближения с помощью полинома третьей степени или показательных функций 1, E^x, E^(2x), E^(3x)

data = Table[{Random[], Random[]}, {i, 10}] poly = Fit[data, {1, x, x^2, x^3}, x] expo = Fit[data, {1, E^x, E^(2x), E^(3x)}, x]

С помощью команды FindMinimum можно численно найти минимум функции на заданном множестве. Найдем, например, первый минимум функции Бесселя BesselJ[0,x] на отрезке [0,5]. Прежде всего, найдем грубо положение минимума, нарисовав график функции

Plot[BesselJ[0, x], {x, 0, 5}];

Мы видим, что минимум расположен вблизи точки x = 4. Уточним его положение

FindMinimum[BesselJ[0, x], {x, 4}]

ФункцияConstrainedMin находит минимум линейной функции от нескольких аргументов в области, ограниченной системой линейных неравенств

ConstrainedMin[2x - 3y, {x > y, y > 1, 2x + y < 5}, {x, y}]

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

slv = NDSolve[{x'[t] == -3(x[t] - y[t]), y'[t] == -x[t]z[t] + 30x[t] - y[t], z'[t] == x[t]y[t] - z[t], x[0] == 0, y[0] == 1, z[0] == 0}, {x[t], y[t], z[t]}, {t, 0, 30}, MaxSteps -> 10000, Method -> RungeKutta, MaxStepSize -> 0.004]; ParametricPlot3D[Evaluate[{x[t], y[t], z[t]} /. slv], {t, 0, 30}];

Фазовый портрет иллюстрирует наличие так называемого "странного аттрактора" у траекторий системы Лоренца [13].

2. Разработка заданий для выполнения лабораторных и расчётно-графических работ по курсу «Вычислительная математика»

 

2.1. Лабораторная работа №1.

Знакомство с пакетом MathCAD 2000

2.1.1. Простейшие приёмы работы.

– Конец работы –

Эта тема принадлежит разделу:

Обзор возможностей математических пакетов MathCAD 2000, MathLAB 5.0, Mathematica

Обзор пакета MathCAD... Введение в MathCAD...

Если Вам нужно дополнительный материал на эту тему, или Вы не нашли то, что искали, рекомендуем воспользоваться поиском по нашей базе работ: Численные методы.

Что будем делать с полученным материалом:

Если этот материал оказался полезным ля Вас, Вы можете сохранить его на свою страничку в социальных сетях:

Все темы данного раздела:

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

Формульный редактор.
Фактически система MathCAD интегрирует три редактора: формульный, текстовый и графический. Для запуска формульного редактора достаточно установить курсор мыши в любом свободном месте окна редактиро

Наборные панели и шаблоны.
Подготовка вычислительных блоков облегчается благодаря выводу шаблона при задании того или иного оператора. Для этого в MathCAD служат наборные панели с шаблонами различных математических символов

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

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

Инструментальная панель
Инструментальная панель командного окна системы MATLAB позволяет обеспечить простой доступ к операциям над М-файлами (рис. 1.4)

Программирование в среде Matlab 5.
Файлы, которые содержат коды языка MATLAB, называются M-файлами. Для создания M-файла используется текстовый редактор; вызову М-файла предшествует присваивание значений входным аргументам; результа

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

Анализ.
Наряду с алгебраическими преобразованиями "Математика'' позволяет выполнять операции математического анализа. Базовыми являются операции интегрирования Integrate и дифференцирования D. Проинте

Формульный редактор.
Фактически система MathCAD интегрирует три редактора: формульный, текстовый и графический. Для запуска формульного редактора достаточно установить курсор мыши в любом свободном месте окна редактиро

Наборные панели и шаблоны.
Подготовка вычислительных блоков облегчается благодаря выводу шаблона при задании того или иного оператора. Для этого в MathCAD служат наборные панели с шаблонами различных математических символов.

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

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

Графы вычислительных процессов.
Рассмотрим более удобный способ подсчёта распространения ошибки в каком-либо арифметическом вычислении [6]. С этой целью мы будем, изображать последовательность операций в вычислении с пом

Деление
Если выполняется деление a1/a2, то стрелка от a1 к косой черте в кружке получает коэффициент +1, а стрелка от a2 к косой черте в кру

Методы решения алгебраических и трансцендентных уравнений.
  2.3.1. Приближённое решение уравнения f(x) = 0 методом деления пополам (методом бисекций). Пусть задана непрерывная функция f(x) и требуется найти корень у

Метод простых итераций.
Метод простых итераций (метод последовательных приближений) решения уравнения f(x) = 0 состоит в замене исходного уравнения эквивалентным ему уравнением x = j(x) и построении последовательности x

Порядок выполнения лабораторной работы с помощью метода простых итераций.
Графически или аналитически отделить корень уравнения f(x) = 0. 1. Преобразовать уравнение f(x) = 0 к виду x = j(x) так, чтобы в некоторой окрестности [a,b] корня x производная j¢(x)

Порядок выполнения лабораторной работы с помощью метода Ньютона
1. Графически или аналитически отделить корень уравнения f(x) = 0. Убедиться, что на найденном отрезке [a,b] функция f(x) удовлетворяет условиям сходимости метода Ньютона. 2. Выбрать начал

Метод Вегстейна.
При выводе метода Вегстейна решения задачи о неподвижной точке x=φ(x) будем использовать как аналитические, так и геометрические соображения [1]. Пусть уже найдены:

Алгоритм Вегстейна.
Шаг 0. Ввод x0 (начального приближения), φ(x) (исходной функции), q (оценки модуля производной), ε (допустимой абсолютной погрешности). Шаг 1. Вычислить

Метод Чебышева.
Требуется найти вещественный корень уравнения f(x) = 0, изолированный в интервале (a, b). Функция f(x) предполагается непрерывной вместе с производными до n-го порядка включительно, причём в интерв

Метод Данко.
Для отыскания действительного корня уравнения f(x)=0, изолированного в интервале (a, b), рассматривается кривая [5]

Метод простых итераций.
Рассмотрим произвольную нелинейную систему уравнений в Rn [10]. или в более к

Метод Ньютона.
Рассмотрим систему n нелинейных уравнений с n неизвестными или в векторной форме f

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

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

Тригонометрическая интерполяция.
Пусть функция f(х) задана на отрезке [0,2p] таблицей значении f(xi) в равноотстоящих узлах (i=1, 2, ..., 2N

Метод градиентного спуска.
Общая задача нелинейного программирования без ограничений состоит в минимизации функции f(x)=f(x1, x2, ..., xп}, заданной во всем n-мерном евклидовом

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

Порядок выполнения лабораторной работы.
1. Составить подпрограмму-функцию для вычисления значений целевой функции. 2. Составить программу-функцию для нахождения экстремума целевой функции методом Нелдера-Мида. 3. Провес

Метод Эйлера.
Пусть требуется найти приближенное решение дифференци­ального уравнения y'=f(x,у), удовлетворяющее начальному усло­вию у(х0)=у0. Численное решение задачи состоит

Метод Рунге-Кутта.
Пусть требуется найти решение дифференциального уравнения y¢=f(x, y), удовлетворяющее начальному условию y(x0) = y0. Численное решение задачи состоит в построени

Решение одной из задач.
В качестве примера возьмём задачу разложения функции в ряд Тейлора (задача №1). Поскольку ЭВМ способна непосредственно вычислять только функции, содержащие арифметические операции (полином

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

К лабораторной работе №3. Методы решения алгебраических и трансцендентных уравнений.
Найти корень уравнения с точностью e = 0.0001. · Отделяем корень графически.

Моделирование сплайн-интерполяции.
Для исследования сплайн-интерполяции составим программу, вычисляющую сплайн-коэффициенты по граничному условию А. Функция Spline реализует алгоритм вычисления сплайн-коэффициентов с учётом

Хотите получать на электронную почту самые свежие новости?
Education Insider Sample
Подпишитесь на Нашу рассылку
Наша политика приватности обеспечивает 100% безопасность и анонимность Ваших E-Mail
Реклама
Соответствующий теме материал
  • Похожее
  • Популярное
  • Облако тегов
  • Здесь
  • Временно
  • Пусто
Теги