Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
метод маткад05(1).pdf
Скачиваний:
37
Добавлен:
06.03.2016
Размер:
889.21 Кб
Скачать

Ермоленко М.В. Введение в систему MathCAD

©Кафедра технической физики, 2008

На Рисунке 8 показано решение системы трех линейных уравнений относительно трех неизвестных

Приближенные решения

Функция Minerr очень похожа на функцию Find (использует тот же алгоритм). Если в результате поиска не может быть получено дальнейшее уточнение текущего приближения к решению, Minerr возвращает это приближение. Функция Find в этом случае возвращает сообщение об ошибке. Правила использования функции Minerr такие же, как и функции Find.

Minerr(z1, z2, . . .)

Возвращает приближенное решение системы уравнений. Число аргументов должно быть равно числу неизвестных.

Если Minerr используется в блоке решения уравнений, необходимо всегда включать дополнительную проверку достоверности результатов.

2.5 Символьное решение уравнений

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

Решение уравнений в символьном виде позволяет найти точные или приближенные корни уравнения:

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

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

Команда Символы®Переменные®Вычислить позволяет решить уравнение относительно некоторой переменной и выразить его корни через остальные параметры уравнения.

Чтобы решить уравнение символьно необходимо:

-Напечатать выражение (для ввода знака равенства используйте комбинацию клавиш [Ctrl]=).

-Выделить переменную, относительно которой нужно решить уравнение, щелкнув на ней мышью.

-Выбрать пункт меню Символы ®Переменные®Вычислить.

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

Чтобы решить систему уравнений в символьном виде, необходимо выполнить следующее:

- Напечатать ключевое слово Given.

24

Ермоленко М.В. Введение в систему MathCAD

©Кафедра технической физики, 2008

-Напечатать уравнения в любом порядке ниже словаGiven. Удостоверьтесь, что для ввода знака = используется [Ctrl]=.

-Напечатать функцию Find, соответствующую системе уравнений.

-Нажать [Ctrl]. (клавиша CTRL, сопровождаемая точкой). MathCAD отобразит символьный знак равенства ®.

-Щелкнуть мышью на функции Find.

Пример 2 Рисунка 7 иллюстрирует символьное решение системы уравне-

ний в MathCAD.

Задания для усвоения темы 2

Упражнение 1

Построить график функции f(x) (Таблица 2) и приблизительно определить один из корней уравнения. Решить уравнение f(x)= 0 с точностью e = 10 – 4 с помощью встроенной функции MathCAD root;

Таблица 2

 

 

 

 

 

f(x)

 

 

 

 

 

 

 

 

 

f(x)

 

 

варианта

 

 

 

 

 

 

 

 

 

 

 

варианта

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

3

 

 

4

 

 

1

ex-1 - x3 - x

 

 

 

 

 

 

9

0.25x3 + x - 2

x Î [ 0, 1]

 

 

 

 

 

 

x Î [ 0, 2 ]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x -

 

 

 

1

 

 

 

 

 

 

 

 

 

arccos

1 - x

2

-x

2

3 + sin( 3.6 x )

 

 

 

 

10

 

 

 

 

 

 

 

 

1 + x 2

 

x Î [ 0, 1]

 

 

 

 

 

 

 

х Î [2, 3]

 

 

 

arccos x -

 

 

 

 

 

 

 

 

 

3x - 4 ln x - 5

3

1 - 0.3x3

11

x Î [ 0, 1]

 

 

 

 

 

 

x Î [ 2, 4 ]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

- arcsin x

 

e x

- e-x - 2

4

 

1 - 0.4 x2

12

x Î [ 0, 1]

 

 

 

 

 

 

x Î [ 0, 1]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3x

- 14

+ ex - e-x

 

 

 

 

 

 

 

 

 

- tg x

5

 

 

 

 

13

 

1 - x

x Î [ 1, 3 ]

 

 

 

 

 

 

x Î [ 0, 1]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 - x + sin x - ln( 1 + x )

6

 

2x2 +

1.2 - cos x - 1

14

x Î [ 0, 1]

 

 

 

 

 

 

x Î [ 0, 2 ]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

æ 2

ö

æ 1

ö

 

1

 

 

 

 

 

 

 

 

 

 

cosç

 

÷

- 2 sinç

 

÷

+

 

 

 

 

х5

х - 0,2

7

 

 

x

15

 

 

è x

ø

 

 

è x

ø

 

х

Î [1, 2]

 

x Î [ 1, 2 ]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8

0.1x 2

- x ln x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x Î[ 1, 2 ]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

25

Ермоленко М.В. Введение в систему MathCAD

©Кафедра технической физики, 2008

Упражнение 2

Для полинома g(x) (Таблица 3) выполнить следующие действия:

1.с помощью команды Символы Þ Коэффициенты полинома создать вектор V, содержащий коэффициенты полинома;

2.решить уравнение g(x) = 0 с помощью функции polyroots;

3.решить уравнение символьно, используя команду Символы Þ Пере-

менные Þ Вычислить.

Таблица 3

g(x)

g(x)

варианта

варианта

 

 

 

 

1.

x4 - 2x3 + x2 - 12x + 20

9.

x4 + x3 - 17x2 - 45x - 100

2.

x4 + 6x3 + x2 - 4x - 60

10.

x4 - 5x3 + x2 - 15x + 50

3.

x4 - 14x2 - 40x - 75

11.

x4 - 4x3 - 2x2 - 20x + 25

4.

x4 - x3 + x2 - 11x + 10

12.

x4 + 5x3 + 7x2 + 7x - 20

5.

x4 - x3 - 29x2 - 71x -140

13.

x4 - 7x3 + 7x2 - 5x + 100

6.

x4 + 7x3 + 9x2 + 13x - 30

14.

x4 + 10x3 +36x2 +70x+ 75

7.

x4 + 3x3 - 23x2 - 55x - 150

15.

x4 + 9x3 + 31x2 + 59x+ 60

8.

x4 - 6x3 + 4x2 + 10x + 75

 

 

Упражнение 3

Решить систему линейных уравнений (Таблица 4):

1.используя функцию Find;

2.матричным способом и используя функцию lsolve.

Таблица 4

Система линейных

Система линейных

варианта

 

уравнений

варианта

 

 

уравнений

1

 

 

 

2

 

 

3

 

 

 

 

 

 

4

 

ì2x + x + 2x + 3x = 8

 

ì2x

+ x

- 5x

+ x

= -4

 

ï

1

2

3

 

4

 

ï

 

1

 

2

3

 

4

1.

ï3x + 3x = 6

 

 

9.

ïx

- 3x

2

- 6x = -7

í

1

3

 

 

 

í

1

 

 

4

 

 

ï2x1 - x2 + 3x4 = 4

 

ï2x2 - x3 + 2x4 = 2

 

ïx + 2x - x + 2x = 4

 

ïx

+ 4x

- 7x

+ 6x = -2

 

î 1

2

3

4

 

 

î

1

 

 

2

3

 

4

 

ìx + 2x + 3x + 4x = 22

 

ìx + 2x + 3x + 4x = 26

 

ï 1

 

2

3

4

 

ï

1

 

2

3

 

4

2.

ï2x + 3x + x + 2x = 17

10.

ï2x + 3x + 4x + x = 34

í

1

2

3

4

í

 

1

 

 

2

3

4

ïx1

+ x2 + x3 - x4 = 8

ï3x1 + 4x2 + x3 + 2x4 = 26

 

ïx - 2x - 3x = -7

 

 

ï4x + x + 2x + 3x = 26

 

î 1

 

3

4

 

 

 

î

 

1

 

2

3

 

4

 

ì9x +10x - 7x - x = 23

 

ì2x - 8x - 3x - 2x = -18

 

ï

1

2

3

4

 

ï

1

 

 

2

3

 

4

3.

ï7x - x - 5x = 37

 

11.

ïx

- 2x + 3x - 2x = 28

í

1

3

4

 

 

í 1

 

 

2

 

3

 

4

ï5x1 - 2x3 + x4 = 22

= 26

ïx2 + x3 + x4 = 10

 

 

ï4x + x + 2x + 3x

4

 

ï11x

2

+ x + 2x = 21

 

î

1

2

3

 

 

î

 

 

3

4

 

 

26

Ермоленко М.В. Введение в систему MathCAD ©Кафедра технической физики, 2008

Продолжение таблицы 4

1

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

3

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

ì6x - x

 

 

+10x - x

 

= 158

 

 

ì2x

- x

+ 4x

 

+ x

4

= 66

 

 

ï

1

2

 

 

 

 

 

 

3

 

 

4

 

= 128

 

ï

1

2

 

 

3

 

 

 

 

 

 

 

 

 

 

4.

ï2x

+ x

 

 

+10x + 7x

 

12.

ï2x

- 6x

 

+ x

 

= -63

 

 

 

 

í

1

 

2

 

 

 

 

 

3

 

 

 

4

7

 

í

2

3

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

ï3x1 - 2x2 - 2x3 - x4 =

 

 

ï8x1 - 3x2 + 6x3 - 5x4 = 146

 

ïx -12x

2

+ 2x - x

4

= 17

 

 

ï2x - 7x + 6x - x

4

= 80

 

î

1

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

î

1

2

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

ìx - 2x

 

+ 6x

+ x

=

88

 

 

ì2x1 - 3x3 - 2x4 = -16

 

 

 

ï

1

 

 

 

 

2

 

 

 

 

3

 

 

 

4

 

 

 

 

 

 

ï5x

+ 2x

 

-

3x

 

 

=

88

 

 

 

ï

 

 

+13x + 4x

 

 

= 213

 

 

 

4

 

 

 

 

ï2x - x

4

5.

 

 

1

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

13.

 

 

1

2

 

 

 

 

3

 

 

 

 

 

 

 

 

í7x

- 3x

 

+

7x

+

2x

= 181

 

í3x1 + x2 + 2x3 + x4 = 72

 

 

ï

1

 

 

 

 

 

2

 

 

 

 

3

 

 

 

4

 

 

 

ï

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ï

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ïx1 -12x3 - 5x4 = -159

 

 

 

- 7x2 + 5x3 + 2x4 = 99

 

 

î

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

î3x1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ìx - 2x

 

- 8x

4

= -7

 

 

 

 

ì7x + 7x - 7x - 2x

4

 

= 5

 

 

ï

1

 

 

 

2

 

 

 

+ 6x

 

= -8

 

 

ï

1

 

2

 

 

 

3

 

 

 

 

 

 

 

6.

 

ïx + 4x

 

- 7x

 

 

14.

 

ï3x + 4x

 

+ 5x + 8x = 60

 

í

1

 

 

 

 

2

 

 

 

3

 

 

 

4

 

 

 

í

1

 

2

 

 

 

3

 

 

 

 

4

 

 

 

 

ïx1 + x2 - 5x3 + x4 = -10

 

 

ï2x1 + 2x2 + 2x3 + x4 = 27

 

 

ï2x

- x

 

+ 2x

4

= 7

 

 

 

 

 

ï2x - 2x - x = -1

 

 

 

 

 

 

î

1

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

î

1

 

3

 

 

4

 

 

 

 

 

 

 

 

 

 

 

ì2x

+ 2x

+ 6x

+ x

 

= 15

 

 

ì6x - 9x

2

+ 5x

 

+ x

 

= 124

 

 

ï

1

 

 

 

 

 

2

 

 

 

 

3

 

 

4

 

 

 

 

ï

1

 

 

 

 

3

 

 

 

 

4

 

 

 

 

7.

 

ï- x + 2x + x

 

= 18

 

 

15.

 

ï7x - 5x - x

 

= -54

 

 

 

í

2

 

 

 

 

 

3

 

 

 

4

 

 

 

 

 

 

 

í

2

 

3

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

ï4x1 - 3x2 + x3 - 5x4 = 37

 

 

ï5x1 - 5x2 + 2x3 + 4x4 = 83

 

 

ï3x - 5x

2

+ x - x = 30

 

 

ï3x - 9x + x + 6x

4

 

= 45

 

 

î

1

 

 

 

 

 

 

 

3

 

 

4

 

 

 

 

 

î

1

 

2

 

 

3

 

 

 

 

 

 

 

 

 

 

ì4x

- 5x

2

+ 7x

+ 5x

= 165

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ï

1

 

 

 

 

 

 

 

 

 

3

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8.

 

ï2x + x

 

- 3x - x

 

= -15

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

í

1

 

 

 

 

2

 

 

 

3

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ï9x1 + 4x3 - x4 = 194

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ïx - x

2

- 2x - 3x = -19

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

î

1

 

 

 

 

 

 

3

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Упражнение 4

 

 

 

 

 

 

 

 

 

Преобразовать

нелинейные

уравнения

системы

из

Таблицы5 к

виду

f 1(x) = y и

f 2 (y)= x. Построить их графики и определить начальное приближе-

ние решения. Решить

систему нелинейных уравнений

с

помощью

функции

Minerr.

 

 

 

 

 

 

 

 

 

 

Таблица 5

 

 

 

 

 

 

 

 

 

 

Система нелинейных

 

 

Система нелинейных

 

варианта

 

уравнений

 

варианта

 

 

уравнений

 

 

1

 

ìsin x + 2 y = 2,

 

9

 

ìsin y + x = -0,4,

 

 

 

ícos(y -1)+ x = 0,7.

 

 

îí2 y - cos(x +1)= 0.

 

 

 

 

î

 

 

 

 

 

 

 

 

 

2

 

ìsin(x + 0,5) - y =1,

 

10

 

ìsin(x + 2) - y =1,5,

 

 

 

îícos(y - 2)+ x = 0.

 

 

îícos(y - 2)+ x = 0,5.

 

 

3

 

ìcos x + y =1,5,

 

11

 

ìcos( x + 0,5) - y = 2,

 

 

 

îí2x - sin(y - 0,5)=1.

 

 

îísin y - 2x =1.

 

 

4

 

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

 

12

 

ìcos( x - 2) + y = 0,

 

 

 

îísin y - 2x =1,6.

 

 

îísin(y + 0,5)- x =1.

 

 

5

 

ìsin(x -1) =1.3 - y,

 

13

 

ìcos( x + 0,5) + y =1,

 

 

 

îx - sin(y +1)= 0.8.

 

 

îsin(y + 0,5)- x =1.

 

 

 

 

í

 

 

 

 

í

 

 

 

 

6

 

ìcos( x + 0,5) + y =1,

 

14

 

ìsin(x) - 2 y =1,

 

 

 

îísin y - 2x = 2.

 

 

îícos( y + 0,5) - x = 2.

 

 

7

 

ì-sin(x +1) + y = 0,8,

 

15

 

ì2 y - sin(x - 0,5) =1,

 

 

 

îísin( y -1) + x =1,3.

 

 

îícos( y) + x =1,5.

 

 

8

 

ìsin(x) - 2 y =1,

 

 

 

 

 

 

 

 

 

í

+ x =1,3.

 

 

 

 

 

 

 

 

 

 

îsin( y -1)

 

 

 

 

 

 

 

 

27

Ермоленко М.В. Введение в систему MathCAD

 

©Кафедра технической физики, 2008

Упражнение 5

 

 

 

Символьно решить системы уравнений:

ì2 y -p z = a,

ì3x + 4 p y = a,

ï

 

í

2x + y = b.

í p z -z = b,

î

ï

3 y + x = c.

 

 

î

Контрольные вопросы

1.Назовите способы нахождения начального приближения.

2.Какие функции для решения одного уравнения вMathCAD вы знаете? В чем их отличие?

3.Какие аргументы функции root не обязательны?

4.В каких случаях MathCAD не может найти корень уравнения?

5.Какая системная переменная отвечает за точность вычислений?

6.Как изменить точность, с которой функция root ищет корень?

7.Как системная переменная TOL влияет на решение уравнения с помощью функции root?

8.Назовите функции для решения систем уравнений вMathCAD и особенности их применения.

9.Опишите структуру блока решения уравнений.

10.Какой знак равенства используется в блоке решения? Какой комбинацией клавиш вставляется в документ?

11.Какие выражения не допустимы внутри блока решения уравнения?

12.Опишите способы использования функции Find.

13.В каких случаях MathCAD не может найти решение системы уравне-

ний?

14.Дайте сравнительную характеристику функциям Find и Minerr.

15.Какие уравнения называются матричными?

16.Как решать матричные уравнения? Назовите способы решения матричных уравнений.

17.Как символьно решить уравнение или систему уравнений Mathв - CAD? Какой знак равенства используется? Какой комбинацией клавиш вставляется в документ?

18.Назовите особенности использования символьного решения уравне-

ний.

28