Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Курсач черновик.docx
Скачиваний:
0
Добавлен:
17.09.2019
Размер:
204.81 Кб
Скачать

Министерство образования и науки Российской Федерации

Филиал ФГБОУ ВПО «УГАТУ» в г.Кумертау

Кафедра _____Промышленной автоматики______________________

Курсовая работа

По курсу «Программирование и основы алгоритмизации»

Руководитель

канд. техн. наук, доцент

______________Л. Н. Грачева

«______»_________20___г.

Исполнитель

Студент группы АТП-124

___________________Н. В. Васильев

«______»_________20___г.

Кумертау 2012

Аннотация

ВВЕДЕНИЕ

2. Решение задания №1

2.1 Исходные данные и условия задачи.

Условия задачи:

1)Решить систему линейных алгебраических уравнений с помощью матричных операторов MATLAB.

2)Решить систему линейных алгебраических уравнений методом итераций.

Исходные данные:

Уравнение типа АХ=В, которое необходимо решить.

Матрица А:

2 16 -3 0

-8 5 0 40

25 0 -2 3

0 -3 20 0

Матрица В:

9

98

5

-7

Метод исключения Гаусса является одним из наиболее простых и эффективных методов. Алгоритм метода Гаусса основан на приведении матрицы А к треугольному виду (прямой ход) и последовательном вычислении неизвестных (обратный ход). Эти процедуры можно выполнять над невыраженными матрицами, в противном случае метод Гаусса неприменим. Недостатком метода является накапливание погрешностей в процессе округления, поэтому метод Гаусса без выбора главных элементов используется обычно для решения сравнительно небольших (n < 100 - 200) систем уравнений с плотно заполненной матрицей и не близким к нулю определителем.

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

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

В данной задаче входной информацией являются:

A и В - матрица системы уравнений

norma - норма матрицы А

Название величины

Тип величины

Объем памяти, байт

Назначение величины

Примечание

1

А

массив

64

Коэф. неизвестных сис-м ур-ий

Матрица 4*4

2

В

массив

16

Коэф. неизвестных сис-м ур-ий

Матрица 4*1

Таблица 1

Промежуточные параметры:

Параметры выводимые в процессе решения задачи.

В данной задаче промежуточными параметрами являются:

i – размерность матрицы

Таблица 2

Название величины

Тип величины

Объем памяти, байт

Назначение величины

Примечание

1

i

Короткое целое число

4

Размерность матрицы

Количество строк и столбцов в данной матрице

Выходные параметры:

Выодными параметрами является результат полученный в результате выполнения программы.

В данной задаче выходными параметрами являются:

а и b - вектор столбец свободных членов

mn и me – норма вектора

Таблица 3

Название величины

Тип величины

Объем памяти, байт

Назначение величины

Примечание

1

a

массив

64

Нижний придел интегрирования

2

b

массив

16

Верхний придел интегрирования

3

mn

переменная

4

Норма вектора

4

me

переменная

4

Норма вектора

2.3 Формализация задачи.

Систему уравнений в MATLAB можно решить различными методами:

1) Использованием обратной матрицы.

X = A-1B. Данное выражение непосредственно вычисляется в MATLAB. Получить обратную матрицу можно с помощью функции inv(A), либо с помощью оператора A^-1 , что означает A в степени -1.

2) Решение системы уравнений методом Гаусса.

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

Недостатком метода является накапливание погрешностей в процессе округления, поэтому метод Гаусса без выбора главных элементов используется обычно для решения сравнительно небольших (n < 100 - 200) систем уравнений с плотно заполненной матрицей и не близким к нулю определителем.

В MATLAB для решения систем уравнений (1) методом исключения Гаусса применяются следующие операторы:

X = A\B

- правое деление, для систем AX = B;

X = B/A

- левое деление, для систем XA = B. Очевидно, что в этом случае X и B есть не векторы-столбцы, а векторы-строки. Для решения данным способом необходимо преобразовать векторы-столбцы в векторы-строки.

2.4 Разработка программы решения задачи в МАТЛАБ.

1)М-файл:

2) М-файл:

2.5 Контрольный пример.

1)

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]