Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Лабораторные работы 1-7

.pdf
Скачиваний:
16
Добавлен:
27.02.2016
Размер:
610.74 Кб
Скачать

Совместная система, имеющая единственное решение, называется определённой; в противном случае неопределённой.

Система называется крамеровской, если в ней число уравнений равно числу неизвестных и определитель матрицы коэффициентов А отличен от нуля.

Крамеровская система имеет единственное решение, которое можно

найти по правилу Крамера: xi = i , где = det A , i - определитель,

получаемый из заменой i - го столбца на столбец свободных членов.

MathCAD решает уравнения при помощи итерационных методов, поэтому требуется задавать начальные приближения, на основе которых строится последовательность, сходящаяся к искомому решению. Различные начальные приближения приводят к различным решениям.

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

3.2. Решение систем уравнений

Для решения системы уравнений необходимо:

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

2.Напечатать ключевое слово Given. Оно указывает MathCAD, что далее следует система уравнений.

3.Ввести уравнения системы ниже ключевого слова Given. Между левыми и правыми частями уравнений должен стоять логический знак

24

равенства. Для этого используйте сочетание клавиш [Ctrl] и [=] или его символ с палитры отношений.

4. Ввести выражение, включающее функцию Find.

Функция Find (х1, х2, х3, … ) возвращает решение системы уравнений. Число аргументов должно быть равно числу неизвестных.

Функция Find возвращает ответ в виде вектора, содержащего значения переменных х1, х2, х3, …, являющиеся решением системы уравнений .

Пример 1. Решить систему уравнений

sin( x + y ) 1,1 x = 0,1

x 2 + y 2 =1

Решение.

x := 1

y := 0.5

Given

sin(x + y) 1.1 x 0.1

x2 + y2 1

Find(x,y) = 0.8050.594

3.3. Решение систем линейных уравнений

Рассмотрим способы решения крамеровской системы линейных уравнений.

Крамеровскую систему линейных уравнений можно решить четырьмя способами:

25

1. Используя блок решения уравнений Given – Find .

2.С помощью функции lsolve .

Функция lsolve (A, B) возвращает вектор решения X такой, что

AX = B .

3.Используя обращение матрицы коэффициентов системы.

4.Методом Крамера.

Пример 2. Решить систему уравнений

3x1 + 6x2 = 9,2x1 + 0,5x2 = 4.

Решение. Найдем решение системы различными способами. 1. Используя блок Given – Find:

x1:= 1

x2:= 0

{задаем начальные приближения}

Given

3

x1

+ 6 x2

 

9

 

 

 

 

 

 

 

 

 

 

2

x1

+ 0.54 x2

 

4

 

 

 

 

 

 

X1 := Find(x1,x2)

 

 

 

 

 

X1 =

 

1.844

 

 

 

 

 

 

0.578

 

 

 

 

 

 

 

 

 

 

 

2. С помощью функции lsolve.

Создаем матрицу коэффициентов системы и вектор свободных членов:

ORIGIN:= 1

A :=

3

6

B :=

9

 

 

 

2

0.54

4

 

 

 

 

 

 

 

Находим решение системы:

 

 

 

 

 

 

X2 := lsolve(A ,B)

 

X2 =

 

1.844

 

 

 

 

 

 

0.578

 

 

 

 

 

 

26

3. Для нахождения решения используем формулу Х=А-1В:

X3 := A

1

B

 

 

X3 =

 

1.844

 

 

 

 

 

 

 

 

 

0.578 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4. Методом Крамера.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а) Составляем главный определитель системы

и убеждаемся, что он не

равен нулю:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

:=

 

A

 

 

= −10.38

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

б) Составляем вспомогательные определители

i:

 

 

 

 

 

d1 := augment(B,A

 

 

)

d1 =

9

 

6

 

 

1

:=

 

d1

 

 

 

 

 

 

 

2

4

0.54

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(

A

1

 

 

)

d2 =

3

9

 

 

 

2 :=

 

d2

 

 

 

 

 

 

 

d2 := augment

 

 

 

,B

2

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

в) Находим хi= i/ . Формируем матрицу решения Х и выводим ответ:

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1.844

 

 

 

 

 

X4 :=

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

X4 =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0.578

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3.4.Задания для самостоятельного решения

1.Найти решение системы уравнений (1), задавая различные приближения (не менее трех).

2.Получить решение линейной системы (2) четырьмя указанными способами. Если найти решение тем или иным способом не удается, указать причину.

27

 

Система (1)

 

 

 

 

 

Система (2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sin( x +1) y =1,2

х

1

+ 7х

2

 

2х

3

= 3

1.

 

 

 

 

 

 

 

 

 

 

 

 

 

3х1 + 5х2 + х3 = 5

 

2x + cos y = 2

 

 

 

 

 

 

 

 

 

 

 

 

2х3 = −4

 

 

2х1 + 5х2

 

 

 

 

 

 

 

 

 

 

 

cos( x 1) + y = 0,5

 

2х

 

+

х

3х

= −4

2.

 

 

 

 

1

 

 

 

2

 

 

 

 

3

 

 

4х1 + 7х2 2х3 = −6

 

x cos y = 3

 

х

 

8х

2

+5х

 

 

=1

 

 

 

 

1

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sin x + 2 y = 2

 

 

2х

 

+

3х

 

+

 

х

= 4

3.

 

 

 

 

1

 

 

 

 

2

 

 

 

3

 

 

 

4х1 х2 +5х3 = 6

 

cos( y 1) + x = 0,7

 

 

х 2х + 4х

= 9

 

 

 

 

 

 

1

 

 

 

 

2

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

cos x + y =1,5

 

х

 

+3х

2х

 

 

= −5

4.

 

 

1

 

 

 

 

2

 

 

3

 

 

 

 

х1 +9х2 4х3 = −1

 

2x sin( y 0,5) =1

 

2х

 

 

+ 6х

2

 

3х = 6

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

sin( x +0,5)y =1

 

 

2

х

х +

х

 

 

= −4

5.

 

 

 

 

1

 

 

 

 

2

 

3

 

 

 

 

 

3х1 9х2 +8х3 = 5

 

cos( y 2) +x =0

 

 

 

 

 

5х2

+5х3 = 4

 

 

 

 

2х1

 

cos( x +0,5 ) + y = 0,8

 

3х

+

5х

6х = −5

6.

 

 

 

 

1

 

 

 

 

2

 

 

 

 

 

3

 

 

2х1 3х2 +5х3 = 8

 

sin y 2x =1,6

 

х

 

+ 4х

 

х

=1

 

 

 

1

 

 

 

 

2

 

3

 

 

 

 

 

 

 

 

 

 

 

 

sin x( x 1) =1,3 y

 

2х

+

5х

6х = −8

7.

 

 

 

 

1

 

 

 

 

2

 

 

 

 

 

3

 

 

х1

+7х2 5х3 = −9

 

x sin( y +1) = 0,8

 

4

х

+

 

2х

х

 

 

= −12

 

 

 

 

 

1

 

 

 

 

 

2

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

28

 

2 y cos( x +1) = 0

2х

+ 4х

 

 

 

3х

 

= −10

8.

 

1

 

 

 

2

 

 

 

 

3

 

 

х1 +5х2 2х3 = 5

 

x + sin y = −0,4

 

 

2х2 + 4х3 = 3

 

 

3х1

 

 

 

 

 

 

 

cos( x +0,5) y = 2

х

5х +

3х

= −1

9.

1

 

 

2

 

 

 

 

 

3

 

 

 

 

2х1 + 4х2 + х3 = 6

 

sin y 2x =1

3х

+3х

2

 

7х

= −13

 

 

 

 

1

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

sin (x + 2) y =1,5

4х

+ 7

х

 

 

 

3х

 

= −10

10.

 

1

 

 

 

2

 

 

 

 

3

 

 

2х1 +9х2 х3

= 8

 

x + cos (y 2) = 0,5

х

 

+6х

 

 

 

3х

= 3

 

 

 

1

 

 

2

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

sin (y +1) x =1,2

2х

+3х

 

 

+ 2

х

= 9

11.

 

 

1

 

 

 

2

 

 

3

 

х1

+ 2х2 3х3

=14

 

2y + cos x = 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3х1 + 4х2 + х3 =16

 

 

 

 

 

 

 

 

 

 

 

cos (y 1) + x = 0,5

х

+ 2

х

 

 

+3х

 

=14

12.

 

1

 

 

 

2

 

 

 

3

 

 

х1 + х2 + х3 = 6

 

y cos x = 3

х

+ х

 

= 3

 

 

 

 

 

 

1

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sin y + 2x = 2

х

+5х

 

+ 4

х

=1

13

 

1

 

 

2

 

 

 

 

3

 

 

 

2х1 10х2 +8х3 = 0

 

cos (x 1) + y = 0,7

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3х1 +15х2 2х3 = 5

 

 

 

 

 

 

 

 

 

 

cos y + x =1,5

х

3х

 

 

+ 2х

 

= −1

14.

 

1

 

 

2

 

 

 

 

3

= 3

2y sin (x 0,5) =1

х1 +9х2 + 6х3

 

 

 

+3х2 + 4х3 =1

 

 

х1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

29

 

sin ( y +0,5) x =1

3х

+ 2

х

+

х = 5

15.

 

 

1

 

 

 

2

 

 

 

3

cos (x 2) + y

= 0

х1

+ х2 х3

= 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4х1 х2 +5х3 = 3

 

 

 

 

 

 

 

 

 

 

 

 

cos( y +0,5 ) + x = 0,8

2

х

 

+ х

 

х

 

= 5

16.

 

1

2

 

 

 

3

 

 

= −5

sin x 2 y =1,6

 

х1 2х2 + 2х3

 

 

 

 

 

 

 

х3 =10

 

 

 

7х1 + х2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sin( y 1) + x =1,3

х

 

+ х

 

 

х

 

= 0

17.

 

1

2

 

 

3

 

 

y sin( x +1)

=0,8

3х1 4х2 +3х3 = −1

 

 

 

х2 3х3 = −8

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

2x cos( y +1) = 0

2х

х

 

 

+ х

 

 

= −1

18.

 

 

1

 

2

 

3

 

 

y + sin x = −0,4

х1 +3х3

= 7

 

 

 

+ х2

 

+3х3 = 6

 

 

 

х1

 

 

 

 

 

 

 

 

 

cos( y +0,5 ) x = 2

3х

2

х

= −5

19.

 

 

1

 

 

 

2

 

 

 

= −1

sin x 2 y =1

 

х1

2х2 + х3

 

 

 

 

 

3х2 х3 = 0

 

 

 

х1

 

 

 

 

 

 

 

 

sin ( y +2) x =1,5

х

3х

+ х

= −2

20.

1

 

 

2

 

 

 

3

 

 

= −11

y +cos (x 2)

= 0,5

х1 2х2 4х3

 

 

 

 

х1 х2 =1

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

30

ЛАБОРАТОРНАЯ РАБОТА № 4 Тема: Построение графиков в декартовой системе координат

Цели работы:

1.Научиться строить графики явно заданных функций в декартовой системе координат.

2.Научиться строить графики параметрически заданных функций в декартовой системе координат.

4.1. Построение графиков явно заданных функций

Если зависимость между переменными х и у выражена уравнением, разрешенным относительно у, то у называется явно заданной функцией и записывается в виде у=f(х) .

Для того чтобы построить график функции у=f(х), необходимо:

1.Определить х как дискретную переменную (в пределах области определения).

2.Задать функцию f(х) .

3.Щелкнуть мышью в свободном месте. Выбрать из меню «Графика» XYPlot (Декартов график).

4.В появившемся шаблоне напечатать х в среднем поле по оси абсцисс, напечатать f(х) в среднем поле по оси ординат.

Остающиеся поля предназначены для ввода границ на осях. Если их оставить незаполненными, MathCAD заполнит их автоматически при создании графика.

5.Щелкнуть мышью вне графика.

Пример 1. Построим график функции у=-2,3х2+5,1.

Решение.

31

x:= −3,−2.8..3

 

 

 

 

f(x) := −2.3 x2 + 5.1

 

 

 

 

10

 

 

 

4

2

0

2

4

f(x)

 

 

 

 

 

10

 

 

 

 

20

 

 

 

 

 

x

 

 

 

рис.1

 

 

 

График можно форматировать. Для этого необходимо:

1.Щелкнуть мышью по чертежу, чтобы выделить его (появляется синяя выделяющая рамка). Дважды щелкнуть мышью по выделенному чертежу. Появится окно форматирования.

2.В верхней части окна форматирования расположены закладки: «ХУ-

оси», «Графики», «Надписи».

«ХУ-оси» позволяет форматировать оси. В частности, можно задать стиль оформления осей: рамка, пересечение, ничего. Рамка окружает график наружной рамкой. Пересечение покажет оси, пересекающиеся в начале координат. Ничего – оси не будут отражены вообще.

«Графики» позволяет форматировать отдельные кривые: указать цвет, тип, толщину и т.п. для каждой кривой.

«Надписи» позволяет оформить график: указать заголовок, подписать название осей.

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

32

зависит от своего аргумента, то на одном рисунке можно построить только 10 графиков.

Чтобы построить несколько графиков функций одной и той же переменной, необходимо:

1.Определить независимую переменную и все функции, графики которых будем строить.

2.Создать шаблон декартового графика.

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

4.Щелкнуть мышью вне графика.

Пример 2. Построим на одном чертеже графики функций у=1-х и у=0,3х2+0,66.

Решение.

x := −5,−4.8.. 5

 

 

y1(x) := 1 x

y2(x) := 0.3 x2 + 0.66

 

 

10

 

y1(x)

5

 

 

 

y2(x)

 

 

5

0

5

 

5

 

 

x

 

 

рис.2

 

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

1. Определить все независимые переменные и все функции, графики которых будем строить.

33