Как найти минимальный элемент в маткаде

Mathcad: программой или одной формулой?

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

процедурный способ поиска минимального по модулю значения в массиве

(код написан так, чтобы он работал независимо от значения системной переменной ORIGIN).

Но куда быстрее сделать это с помощью оператора векторизации:

быстрый способ с оператором векторизации

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

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

а можно просто узнать длину вектора, составленного из вхождений минимального элемента в исходный вектор:

количество минимальных элементов в массиве одним вызовом

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

Современная тенденция языков и пакетов — готовые алгоритмы вместо их программирования, например, библиотека <algorithm> в C++, функции обработки и фильтрации массивов в PHP или встроенные операторы суммирования, перемножения, векторизации и т.п. в Mathcad.

Общие правила расчётов по типовым алгоритмам можно сформулировать так:

• для каждой искомой величины создаётся по переменной того же типа данных, что исследуемые значения;
• до цикла обработки данных всем искомым в нём переменным присваиваются начальные значения. В простейших случаях достаточно присвоить:
  • ноль – количеству, сумму или арифметическому среднему;
  • единицу – произведению;
  • заведомо малое для анализируемых данных значение или первый элемент последовательности данных – максимуму;
  • заведомо большое для анализируемых данных значение или первый элемент последовательности данных – минимуму;
• в цикле обработки, если очередной элемент данных t (в качестве t может выступать элемент массива, выражение и т.д.) должен быть обработан алгоритмом, применяются операторы вида:
  • k := k + 1; при поиске количества элементов k;
  • s := s + t; при поиске суммы (или арифметического среднего) s;
  • p := p * t; при поиске произведения элементов p;
  • if t<min then min := t; при поиске минимума min из элементов данных t;
  • if t>max then max := t; при поиске максимума max из элементов данных t;

где «:=» — оператор присваивания. Эти правила не зависят от языка или пакета, в которых вы работаете.

Однако если в языке, как и в Mathcad, есть встроенные циклические операторы и встроенная условная функция, все 3 шага типового алгоритма можно выполнить одним выражением.

Приведём несколько примеров, иллюстрирующих это утверждение, найдя для массива A одним оператором:

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

Вот скриншоты расчётов, сделанных в том же порядке:

некоторые расчёты для типовых алгоритмов операторами суммы и произведения

Последний пример показывает, что во всём нужно знать меру. В самом деле, нам пришлось трижды фильтровать массив прежде, чем его стало можно обработать стандартной функцией min.
Трудно придумать аргументы, чем это лучше и нагляднее классического кода

минимальный из положительных элементов массива — классический способ лучше

Также при поиске произведений следует иметь в виду вот это замечание (2).

 Скачать документ .xmcd Mathcad 14/15 с этими расчётами в архиве .zip (20 Кб)

29.03.2017, 15:21 [9382 просмотра]

---

К этой статье пока нет комментариев, Ваш будет первым

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

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

• Minimize
  (f, x1, … ,хn)
  – вектор значений аргументов, при
  которых функция f
  достигает минимума;

• Maximize
  (f, x1, … ,хn)
  – вектор значений аргументов, при
  которых функция f
  достигает максимума;

• f(x1,
  … ,хn)–
  заданная целевая функция;

• x1,
  … ,хn
  – аргументы, по которым производится
  минимизация(максимизация).

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

Пример
7. Поиск локального экстремума в
окрестности заданной точки.

Найти
максимум функции
 в
окрестности точки (4;5).

Ответ:
функция имеет максимум, равный 4, в
точке(1;1).

Пример
8. Поиск условного экстремума функции.

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

1. Задаем
   целевую функцию, матрицу системы
   ограничений и вектор правой части этой
   системы

2. Задаем
   начальное приближение решения.

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

Ответ:
минимум функции равен 32.155 и достигается
в точке (1,0.623,0.343,1,0.048,1).
Контрольные
вопросы.

1. Что
   значит отделить корень
    уравнения
   ?

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

3. Опишите
   конструкцию вычислительного блока.

4. В
   чем состоит градиентный метод?

5. В
   чем различие между функциями Find
   и Minner
   для решения систем нелинейных уравнений?

6. Где
   необходимо расположить ограничительные
   условия при решении задачи оптимизации?

7. Как
   ограничено число ограничительных
   условий для решения задачи оптимизации?

Варианты заданий

Вариант
1

1. Решить
   уравнение
   ,
   используя встроенные функции root
   и
   Find.
   Сравнить полученные решения.

2. Найти
   все корни полинома
   .
   Проиллюстрировать решение графически.

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

4. Найти
   максимум функции .

Вариант
2

1. Решить
   уравнение 
   ,
   ,
   используя встроенные функции root
   и
   Find.
   Сравнить полученные решения.

2. Найти
   все корни полинома
   .
   Проиллюстрировать решение графически.

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

4. Найти
   максимум функции
    при
   ограничении
   .

Вариант
3

1. Решить
   уравнение 
   ,
   ,
   используя встроенные функции root
   и
   Find.
   Сравнить полученные решения.

2. Найти
   все корни полинома
   .
   Проиллюстрировать решение графически.

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

4. Найти
   максимум функции при
   ограничении.

Вариант
4

1. Решить
   уравнение
   ,
   ,
   используя встроенные функции root
   и
   Find.
   Сравнить полученные решения.

2. Найти
   все корни полинома
   .
   Проиллюстрировать решение графически.

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

4. Найти
   максимум функции
   .

Вариант
5

1. Решить
   уравнение
   ,
   ,
   используя встроенные функции root
   и
   Find.
   Сравнить полученные решения.

2. Найти
   все корни полинома
   .
   Проиллюстрировать решение графически.

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

4. Найти
   минимальное и максимальное значения
   функции
   .

Вариант
6

1. Решить
   уравнение
   ,
   ,
   используя встроенные функции root
   и
   Find.
   Сравнить полученные решения.

2. Найти
   все корни полинома 
   .
   Проиллюстрировать решение графически.

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

4. Найти
   максимум функции
    при
   условиях
   ,
   , .

Вариант
7

1. Решить
   уравнение
   ,
   ,
   используя встроенные функции root
   и
   Find.
   Сравнить полученные решения.

2. Найти
   все корни полинома
   .
   Проиллюстрировать решение графически.

3.  Решить
   систему нелинейных уравнений:
   .
   Выполнить
   проверку.

4. Найти
   минимум функции
    при
   условиях
   ,
   ,
   .

Вариант
8

1. Решить
   уравнение
   ,
   ,
   используя встроенные функции root
   и
   Find.
   Сравнить полученные решения.

2. Найти
   все корни полинома
   .
   Проиллюстрировать решение графически.

3. Решить
   систему нелинейных уравнений:
   .
   Выполнить
   проверку.

4. Найти
   минимум функции
    при
   условиях
   ,
   , .

Вариант
9

1. Решить
   уравнение
   ,
   ,
   используя встроенные функции root
   и
   Find.
   Сравнить полученные решения.

2. Найти
   все корни полинома
   .
   Проиллюстрировать решение графически.

3. Решить
   систему нелинейных уравнений:
   .
   Выполнить
   проверку

4. Найти
   минимум функции
    при
   условиях 
   ,
   ,
   .

Вариант
10

1. Решить
   уравнение, предварительно оделив корни
   7,
   ,
   используя встроенные функции root
   и
   Find.
   Сравнить полученные решения.

2. Найти
   все корни полинома
   .
   Проиллюстрировать решение графически.

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

.
Выполнить проверку.

1. Найти
   минимум функции
    при
   условиях
   ,
   ,
   .

Тема
4:
Элементы программирования в пакете
инженерных расчетов MathCAD.
Цель
работы:
Изучение  возможностей символьного
пакета MathCAD для программирования условных
и циклических выражений. Приобретение
навыков написания простейших
программ.
Используемые
программные средства:
Пакет MathCAD.

Теоретические
сведения.

Для
повышения гибкости Mathcad в системе
предусмотрена возможность написания
небольших программ для решения тех
проблем, которые не могут быть реализованы
стандартными средствами. Обычно прибегать
к программированию приходится в тех
случаях, когда стандартные средства
либо не могут решить задачу, либо
неэффективны. Для написания программ
используется программная палитра,
которая вызывается кнопкой панели
управления. Как видно, всего имеется 10
операторов, из которых и строится
программа. Причём операторы должны
вводиться только из палитры, писать их
«вручную» не рекомендуется.

Практическая
часть
Для
примера приведём простую программу
возвращающую 1, если число чётное, и 0 в
противном случае.
 –
начинаем создание программы с кнопки
Add
Line.
Вертикальная линия играет роль операторных
скобок.
 –
оператор локального присваивания. В
программе нельзя использовать оператор
присваивания «:=», вместо него
используется оператор локального
присваивания, отличие которого заключается
в том, что локальная переменная определена
только внутри своего блока и при выходе
из программы теряет своё значение.
Пример:

–
условный
оператор создаёт конструкцию видa:
где
первый операнд выполняется, если
справедливо условие являющееся вторым
операндом. Из всех программных операторов
оператор условия является, пожалуй,
наиболее важным. Его приходится
использовать практически во всех
создаваемых алгоритмах. Как уже показано,
условный оператор if
имеет два маркера. В правый маркер
вводится условие, в левый – операция,
которая должна быть проделана в случае,
если условие выполнится (если же оно не
выполнится, система просчитывает
программу, пропуская данный фрагмент).
В маркер оператора может быть внесено
несколько выражений условий или операций.
В случае задания комплекса условий
будьте предельно внимательны и всегда
помните, чем отличаются формы его
определения через программный блок и
с использованием логического умножения.
Неверное задание формы комплекса условий
– самая распространенная ошибка при
работе с данным оператором.
Пример:
Функция возвращает остаток от деления
на 2:

После
того, как функция определена, она может
использоваться наравне со встроенными
функциями.
С оператором if
возможны более сложные конструкции при
использовании ещё одного оператора
–
реализующего альтернативу. Аналог
традиционной конструкции Если
… То … Иначе … .
Т.е. оператор предназначен для определения
того
действия, которое должно быть выполнено,
если условие оператора if
окажется истинным. Одновременно может
быть использовано несколько условных
операторов if.
Оператор otherwise
в таком случае будет задействован, если
не выполняется условия всех операторов
if.

 –
оператор
цикла обеспечивает повторяющиеся
вычисления, если известно количество
шагов.
Первый операнд – переменная
цикла, это интервальная переменная и
её
значения определены во втором операнде.
Третий операнд – тело цикла может
состоять из блока операторов и выполняется
пока не исчерпаются все значения
переменной цикла.
Пример:

Рассмотрим
программу упорядочения чисел по убыванию
в одномерном массиве. Пусть дан массив
чисел:  Мы воспользуемся вложенными
циклами и в качестве тела цикла по i
используем ещё один цикл по j.
Здесь
реализован простейший алгоритм
сортировки, так называемый «метод
пузырька». На самом деле, большие
числа как бы всплывают наверх при каждом
шаге цикла по i, в то время как в цикле
по j на каждом шаге происходит сравнение
пары чисел и замена, если большее число
находится ниже, причём эта замена
осуществляется снизу.
Отметим,
что в системе имеется стандартная
функция сортировки sort().
Примечание:
второй цикл мы организовали с отрицательным
шагом от конечного значения к начальному.
Можно использовать программные
возможности Mathcad просто для задания
функций более сложного вида.

 –
оператор
break
введён для повышения гибкости
программирования и позволяет закончить
цикл досрочно, не исчерпав всего списка.
Как правило, ввиду того, что цикл бывает
нужным остановить при выполнении
некоторого условия, данный оператор
почти всегда используется с условным
оператором if.

Рассмотрим
следующую задачу: нам необходимо найти
первое вхождение 0 в числовом массиве
и вернуть его индекс. Мы приводим
работающую программу, где введена
функция last(M)
которая возвращает последний индекс
массива. Возвращаемым значением программы
является последний выполняемый оператор
— k.
 –
ещё один полезный оператор позволяющий
прервать выполнение текущей итерации
и перейти к следующей.
Рассмотрим
например задачу нахождения максимального
и минимального элемента
массива. 

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

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

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

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

 –
оператор
on
error
служит для обработки ошибочных ситуаций
типа «деления на нуль». Данный
оператор по своему синтаксису полностью
соответствует условному оператору if.
В правый его маркер следует внести
величину или выражение, ошибка в
вычислении которого должнв быть
зарегистрирована. В левом маркере
следует прописать условие, которое
должно быть выполнено в случае ошибочной
ситуации.
Например, нам необходимо
описать функцию
,
чтобы не было особенности в 0.

 тогда:

Примечание.
Оператор
on
error
может
использоваться в арифметических
выражениях.

Варианты
заданий

Вариант
1

1. Составить
   программу которая будет менять местами
   2 строки матрицы.

2. Используя
   оператор on error для предотвращения
   появления ошибки «деление на нуль»,
   вычислить функцию
   .

3. Написать
   программу, где функции, возвращающая
   –1, 0 или 1 в зависимости от знака аргумента
   (соответственно «-», 0, «+»).

Вариант
2

1. Составить
   программу, которая будет выводить сумму
   геометрической прогрессии, при заданных
   значениях
   .

2. Написать
   программу по выводу матрицы размером
   3х3, где на второстепенной диагонали
   стояли бы 1.

3. Вычислить
   значение функции
   .

Вариант
3

1. Вычислить
   сумму чисел 1+2+3+….+n.

2. Для
   x изменяющего от -2 до 2 вычислить значение
   .

3. Составить
   программу по вычислению длины вектора.

Вариант
4

1.Создайте
программу pr(n)
для вычисления произведения
чисел
1*2*3*….*n.
2.Создайте
программу для вычисления и вывода двух
корней квадратного уравнения
.
3.Используя
оператор on error для предотвращения
появления ошибки «деление на нуль»,
вычислить функцию
.

Вариант
5

1. Создайте
   программу для вычисления и вывода двух
   корней квадратного уравнения
   f(x)=ax2+bx+c.

2. Дано
   натуральное число n, действительное x.
   Вычислить
   .

3.  Написать
   программу по выводу на экран знака «+»
   («-»), если значение
   .

Вариант
6

1. Написать
   программу, которая будет выводить
   единичную матрицу.

2. Вычислить
   значение
   ,
   используя итерационную процедуру
    В
   качестве приближенного значения корня
   квадратного берется такое значение
   ,
   которое удовлетворяет условию
   ,
   где
    –
   заданная точность вычисления.

3. Вычислить
   ,
   где
    меняется
   от 0 до 10.

Вариант
7

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

2. Составим
   программу для вычисления переменной
   z по формуле
   .

3. Написать
   программу, которая выводит абсолютное
   значение функции
   .

Вариант
8

1. Определить
   функцию, которая равна 1, если аргумент
   размещён между чётным и нечётным числом,
   и 0 в противном случае.

2. Вычислить
   сумму бесконечной геометрической
   прогрессии со знаменателем 0,4 и начальным
   элементом -8.

3. Написать
   программу по выводу на экран матрицы,
   каждый элемент которой будет вычисляться
   по правилу:
   ,
   где
    и
    –
   элементы матриц соответственно
   

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

• #
• #
• #
• #
• #
• #
• #
• #
• #
• #
• #

Часто используемые функции

augment(A,B) — Возвращает матрицу, сформированную путем размещения массива В справа от массива А. А и В должны иметь одинаковое число строк:

ceil(y) – возвращает наименьшее целое, большее или равное y. Значение y должно быть вещественным числом:

cols(A) – число столбцов матрицы А.

csort(B,n) – сортирует строки матрицы В таким образом, чтобы расположить элементы столбца n в порядке возрастания. Нумерация столбцов по умолчанию начинается с нуля:

eigenvals(A) – определяет вектор собственных значений для квадратной матрицы А:

find(x,y,…) – возвращает значения x,y,…, удовлетворяющие ограничениям: равенствам или неравенствам, заданным в блоке given решения уравнений. Число уравнений должно равняться числу неизвестных.

floor(y) – возвращает наибольшее целое, меньшее или равное y. Значение y должно быть вещественным числом:

identity(n) – создает единичную матрицу размером

, в которой диагональные элементы равны 1, а остальные элементы равны 0:

length(v) – длина вектора v.

max(v) – максимальный по значению элемент вектора v.

min(v) – минимальный по значению элемент вектора v.

maximize(f,v) – возвращает вектор размерности n, который обеспечивает функции f максимальное значение. Функция f – функция n переменных; вектор v – вектор начальных приближений ответа;

minimize(f,v) — возвращает вектор размерности n, который обеспечивает функции f минимальное значение. Обращение аналогично функции Maximize.

mean(M) – среднее арифметическое элементов вектора или матрицы.

median(v) – медиана вектора. Элементы вектора должны быть заданы в порядке возрастания:

rbeta (m, s ₁, s₂)Возвращает вектор m случайных чисел, имеющих бэта-распределение. есть параметры формы.

rbinom (m, n, p) Возвращает вектор m случайных чисел, имеющих биномиальное распределение. . n есть натуральное число.

rgamma (m, s) Возвращает вектор m случайных чисел, имеющих гамма- распределение, есть параметр формы.

rgeom (m, p) Возвращает вектор m случайных чисел, имеющих геометрическое распределение. .

rlogis (m, l, s) Возвращает вектор m случайных чисел, имеющих логистическое распределение, в котором l является п а раметром расположения, а есть параметр масштаба.

rnorm (m,l, s) Возвращает вектор m случайных чисел, имеющих нормальное распределение. .

rpois (m,d) Возвращает вектор m случайных чисел, имеющих распределение Пуассона. .

rt (m, d) Возвращает вектор m случайных чисел, имеющих t-распределение Стьюдента..

runif (m, a, b) Возвращает вектор m случайных чисел, имеющих равномерное распределение, в котором b и a являются граничными точками интервала. .

rnd (x) Возвращает равномерно распределенное случайное число между 0 и x. Эквивалент runif (1, 0, x).

round(y,n) – округляет вещественное число y до n знаков справа от десятичной точки. Если n отсутствует, то y округляется до ближайшего целого числа. Если n<0, то y округляется до n знаков слева от десятичной точки.

rows(A) – число строк матрицы А.

rsort(B,n) – сортирует столбцы матрицы В таким образом, чтобы расположить элементы строки n в порядке возрастания:

Нумерация строк по умолчанию начинается с нуля.

submatrix(M,ir,jr,ic,jc) – подматрица, состоящая из элементов матрицы М, содержащихся в строках от ir до jr и столбцах от ic до jc: