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

маткад учебник

.pdf
Скачиваний:
1817
Добавлен:
10.02.2015
Размер:
565.34 Кб
Скачать

3. ОПЕРАТОРЫ

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

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

Векторные и матричные операторы были рассмотрены раньше. Ниже приводится таблица основных операторов для переменных любого типа. Здесь Х и Y обозначают любую переменную; f – функцию; m и n – целые числа; i – дискретный аргумент; z и w обозначают вещественные или комплексные числа.

 

 

 

Таблица 3

 

Основные операторы

 

 

 

 

Операция

Обозначение

Клавиши

Описание

 

 

 

 

Круглые

(Х)

'

Группирование операторов

скобки

 

 

 

 

 

 

 

 

Факториал

n!

!

Возвращает значение, равное

 

 

 

1·2·…·(n – 1)·n, n – неотрицатель-

 

 

 

ное целое число

31

 

 

 

 

 

 

 

 

 

Продолжение табл. 3

 

 

 

 

 

 

 

 

 

Операция

Обозначение

Клавиши

 

 

 

Описание

 

 

 

 

 

 

 

Комплексное

 

 

 

 

"

Меняет знак мнимой части

X

 

сопряжение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Степень

 

zw

^

Возводит z в степень w

Изменение

–X

Умножает Х на –1

 

знака

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Квадратный

 

 

 

z

\

Возвращает положительный квад-

корень

 

 

 

 

ратный корень для положитель-

 

 

 

 

 

 

 

 

 

 

 

ного z; главное значение для

 

 

 

 

 

 

невещественных z

 

 

 

 

 

 

Корень n-й

n

z

Ctrl+\

Возвращает корень n-й степени

степени

 

 

 

 

 

z; возвращает вещественный ко-

 

 

 

 

 

 

рень, когда возможно

 

 

 

 

 

 

Абсолютное

 

|z|

|

Возвращает Re z

2 + Im z 2

значение

 

 

 

 

 

 

 

 

 

 

Деление

 

X

/

Делит выражение Х на скаляр z.

 

 

z

 

Если X – массив, делит каждый

 

 

 

 

 

 

элемент на z

 

Умножение

X·Y

*

Возвращает произведение Х и Y,

 

 

 

 

 

 

если Х и Y – скаляры. Умножает

 

 

 

 

 

 

каждый элемент Y на Х, если Y –

 

 

 

 

 

 

массив, а Х – скаляр. Возвращает

 

 

 

 

 

 

скалярное произведение, если Х

 

 

 

 

 

 

и Y – векторы одного размера.

 

 

 

 

 

 

Выполняет умножение матриц,

 

 

 

 

 

 

если Х и Y – матрицы соответст-

 

 

 

 

 

 

вующих размеров

 

 

 

 

 

 

Суммирование

n

X

Ctrl+Shift+4

Выполняет суммирование Х по

 

 

i

=

m m

n X

может быть

 

i = m

 

 

 

,

+ 1, ..., .

 

 

 

любым выражением, m и n долж-

 

 

 

 

 

 

 

 

 

 

 

 

ны быть целыми числами

 

 

 

 

 

 

 

Произведение

n

 

Ctrl+Shift+3

Выполняет перемножение Х по

 

X

 

i = m, m + 1, ..., n. X может быть

 

i = m

 

любым выражением, m и n долж-

 

 

 

 

 

 

ны быть целыми числами

32

 

 

 

 

 

 

 

 

 

 

 

 

Окончание табл. 3

 

 

 

 

 

 

 

 

 

Операция

 

Обозначение

 

Клавиши

 

 

 

Описание

 

 

 

 

 

 

 

 

Суммирование

 

X

 

$

 

Возвращает суммирование Х по

по дискретному

i

 

 

 

 

 

дискретному аргументу i. X мо-

аргументу

 

 

 

 

 

 

 

жет быть любым выражением

 

 

 

 

 

 

 

 

Перемножение

 

X

 

#

 

Возвращает произведение Х по

по дискретному

i

 

 

 

 

 

дискретному аргументу i. X может

аргументу

 

 

 

 

 

 

 

быть любым выражением

20

n = 210

20

(n + 1) = 5.109× 1019

20

10

n m = 2.554× 1013

n =0

 

 

n = 0

 

 

 

 

 

n =0 m =0

 

i := 0.. 20

xi := sin(0.1 i π )

20

xn = 0

 

20

xn n = −63.138

 

 

 

 

 

 

 

n =0

 

 

 

n =0

j := 1.. 10

y j :=

 

ij y1 = 210

 

 

ij = 2.554× 1013

 

 

 

i

 

y10 = 2.416×

1013

i

j

 

 

 

 

 

 

 

 

 

Переменный верхний предел суммирования

 

 

 

i := 0.. 10

f(x) := i2 (i x)

g(n) := n

j

m

 

 

 

 

i

 

 

 

 

j

=1 m =1

 

 

f(0) = 0

g(7) = 84

f(2) = 5

 

 

 

 

 

f( 30) = 55

 

g(20) = 1.54 ×

103

 

 

 

 

 

Рис. 11. Суммы и произведения

На рис. 11 приведены примеры вычисления различных сумм и произведений.

33

x := 2

 

 

 

d

x3 = 12

 

 

 

y := 10 t := 0

g(t) := 5t4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

Численное значение производной

Точный результат

 

d

x5 = 79.9999999999999

g(x) = 80

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

d

x5y = 799.999999999999

g(x)y = 800

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

d

x5 y = 32

 

 

 

 

 

x5 = 32

 

 

 

 

 

 

 

 

dy

 

 

 

 

 

 

 

 

 

 

 

 

Поскольку выражение не зависит от t,

d x5y = 0

его производная по t равна нулю:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

Вычисление производной в последовательности точек:

 

 

 

f(t) := d g(t)

i := −2.. 3

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

i

 

 

 

 

g

 

 

=

f

 

 

=

f(2) = −160

 

 

 

 

 

 

 

 

 

2

 

 

2

 

 

 

 

 

5

 

 

 

–20

 

f(2) = 160

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0.313

 

 

 

–2.5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f(3) = 540

 

 

 

 

 

0

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0.313

 

 

 

2.5

 

f(2 + i) = 40 + 220i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

 

 

 

20

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

25.313

 

 

67.5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 12. Дифференцирование с помощью Mathcad

34

Для вычисления производных в Mathcad необходимо выполнить следующее:

сначала определите точку, в которой необходимо найти производную. Например, чтобы найти производную x3 по x в точке x = 2, наберите [x:2];

щелкните ниже или справа от определения x. Затем наберите [?]. Появляется оператор производной с двумя пустыми полями;

щелкните на поле в знаменателе и наберите [x]. Это имя переменной, по которой проводится дифференцирование;

щелкните на поле справа от d/dx и наберите [x^3]. Это выражение, которое нужно продифференцировать;

нажмите знак [=], чтобы увидеть результат.

На рис. 12 приведены примеры использования операций дифференцирования.

Выражение, которое нужно продифференцировать, может быть вещественным или комплексным. Переменная дифференцирования должна быть простой неиндексированной переменной. От алгоритма вычисления производной, который используется в Mathcad, можно ожидать, что первая производная будет вычислена с точностью 7 или 8 значащих цифр, если точка, в которой ищется производная, удалена от особенностей функции. Точность этого алгоритма уменьшается на одну значащую цифру при каждом увеличении порядка производной. Необходимо помнить, что результат дифференцирования есть не функция, а число – значение производной в указанной точке переменной дифференцирования. В предыдущем примере производная от x3 не есть выражение 3x2, а значение 3x2, вычисленное в точке x = 2. Информацию о символьных вычислениях производных см. в разделе "Символьные вычисления".

Использование оператора интегрирования в Mathcad очень похоже на использование оператора дифференцирования. Для вызова оператора интегрирования можно набрать знак [&] или использовать соответствующую палитру. Далее заполняются все свободные поля. Для приближенного вычисления определенного интеграла Mathcad использует численный алгоритм интегрирования Ромберга. Пределы интегрирования должны быть вещественными числами. Выражение, которое нужно интегрировать, может быть, однако, вещественным или комплексным. Кроме переменной интегрирования, все переменные в подынтегральном выражении должны быть определены ранее в другом месте рабоче-

35

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sin(x)2 dx = 0.143

i := 0.. 5

 

f(x) := x2 + 3 x + 2

 

 

 

 

0

 

i

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

gi =

fi =

gi :=

f(x) dx

fi :=

f(x) dx

 

 

 

 

 

0

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3.833

 

0.833

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

12.667

 

0.667

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

28.5

 

1.5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

53.333

 

5.333

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

89.167

 

14.167

 

 

 

 

 

 

 

 

 

 

 

 

Влияние значения встроенной переменной TOL

 

 

 

 

 

 

 

 

на вычисление интеграла:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

e5

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

Точность по умолчанию

TOL = 1 × 10

 

 

 

 

dt = 5.00000000006814

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

6

 

e5

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

Высокая точность

 

TOL := 10

 

 

 

 

 

dt = 5.00000000000001

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

e5

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

Низкая точность

 

TOL := 10

 

 

 

 

 

dt = 5.00000005437786

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

Криволинейный интеграл в комплексной плоскости

 

 

 

 

x(t) := 2 cos (t)

y(t) := 2 sin(t)

z(t) := x(t) + i y(t)

f(z) :=

1

 

z

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f(z(t)) d z(t) dt = 3.142i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 13. Вычисление интегралов

36

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

Подобно всем численным методам, точность алгоритма интегрирования в Mathcad зависит от особенностей подынтегрального выражения. Если выражение, которое нужно интегрировать, имеет разрывы или быстро осциллирует, численное решение, найденное Mathcad, может быть неточно. Поскольку метод интегрирования Mathcad делит интервал на четыре подынтервала, а затем удваивает число точек разбиения, это может привести к неправильным результатам для периодических функций с периодом 1/2n от длины интервала. Чтобы обойти эту проблему, делите интервал на два подынтервала, не кратных периоду функции и интегрируйте по каждому подынтервалу отдельно. Для вычисления точного значения интеграла могут быть использованы возможности символьных вычислений Mathcad.

На рис. 13 приведены примеры использования оператора интегрирования.

37

4.ВСТРОЕННЫЕ ФУНКЦИИ

Вэтом разделе описаны основные функции Mathcad и даны примеры по их использованию.

Функции, используемые для работы с векторами и матрицами, описаны выше.

Чтобы просмотреть список всех встроенных функций, доступных в вашей версии Mathcad, и их описание, выберите команду Функция из меню Вставка.

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

Следует отметить, что существуют Пакеты Функций Mathcad, приспособленные к специфической области применения. Описания этих функций содержатся в самом Пакете Функций. Другим источником дополнительных функций являются Электронные Книги. Функции и описания к ним можно непосредственно копировать в рабочий документ. И, наконец, можно создавать собственные встроенные функции [1]. Для этого нужно иметь 32-разрядный компилятор языка С.

4.1. Тригонометрические функции

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

sin z = eiz eiz ; 2i

cos z = eiz + eiz ; 2

38

eiz = cos z + i sin z .

Все эти тригонометрические функции используют аргумент, выраженный в радианах. Чтобы перейти к градусам, используется встроенная единица deg. Например, чтобы вычислить sin45°, введите [sin(45*deg)=].

Обратные тригонометрические функции asin(z), acos(z), atan(z) возвращают угол в радианах между 0 и 2π . Можно отобразить результат в градусах или в единицах π . Это делается следующим образом. Набрать с помощью клавиатуры [asin(0.5)=], на экране получим

asin(0.5)=0.524

Щелкнем левой кнопкой мыши в конце этого выражения и введем в

появившееся окно ввода "deg", получим asin(0.5)=30 deg

Теперь вместо deg введем из палитры греческих букв π , получим asin(0.5)=0.167π

4.2. Логарифмические и показательные функции

Логарифмические и показательные функции в Mathcad могут использовать комплексный аргумент и возвращать комплексные значения. Значения экспоненциальной функции для комплексного аргумента вычисляются с применением формулы

ex+iy = ex (cos y + i sin y) .

В общем случае значения натурального логарифма определяются формулой

ln(x + iy) = ln x + iy + i arg( y x) + 2nπ i .

В Mathcad функция возвращает основное значение логарифма, соответствующее n = 0, а именно:

ln(x + iy) = ln x + iy + i arg( y

x) = ln x + iy +

+ i arctg( y x) + iπ (1− Φ

(x))sign y,

где Ф(х) – единичная функция Хевисайда (см. подразд. 4.7).

ВMathcad используются следующие встроенные логарифмические

ипоказательные функции:

39

exp(z) - возвращает e в степени z;

ln(z) - возвращает натуральный логарифм z (z ≠ 0); log(z) - возвращает логарифм z по основанию 10 (z ≠ 0).

Для вычисления логарифма по основанию n можно ввести функцию

logn (z , n ) := ln (z ) ln (n )

В этом случае получим logn (1 6 , 2 ) = 4

4.3. Специальные функции и функции усечения

Функция erf(x) возвращает значение интеграла ошибок

 

2

x

 

erf(x) =

et2

dt ,

π

 

0

 

 

 

 

где x должен быть вещественным.

Эйлерова гамма-функция Г(z) для вещественного z совпадает со следующим интегралом:

Γ (z)= t z1etdt .

0

Гамма-функция Эйлера удовлетворяет рекуррентному соотношению

Γ (z+ 1)= zΓ (z) .

Откуда следует, что для положительных целых z

Γ (z+ 1)= z!

Для усечения и округления функций могут использоваться следующие встроенные функции:

Re(z) – вещественная часть; Im(z) – мнимая часть;

arg(z) – аргумент, заключенный между –π и π ;

floor(x) – наибольшее целое число, не превосходящее x (x вещественное);

ceil(x) – наименьшее целое число, превышающее x (x вещественное);

40

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