Как найти максимум в массиве маткад

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. Написать
   программу по выводу на экран матрицы,
   каждый элемент которой будет вычисляться
   по правилу:
   ,
   где
    и
    –
   элементы матриц соответственно
   

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

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

Slawa писал(а):А что делает функция match? (не нахожу я ее в help).
Наравне с мин/макс значением меня интересует НОМЕР его элемента (координаты) в матрице. Есть что- то подобное в маткаде?

Кроме функции match есть функции для построчного/”постолбцового” поиска:

lookup(z, A, B) Looks in a vector or matrix, A, for a given value, z, and returns the value(s) in the same position(s) (i.e., with the same row and column numbers) in another matrix, B. When multiple values are returned, they appear in a vector in row-wise order, starting with the top left corner of B and sweeping to the right.

hlookup(z, A, r) Looks in the first row of a matrix, A, for a given value, z, and returns the value(s) in the same column(s) in the row specified, r. When multiple values are returned, they appear in a vector.

vlookup(z, A, c) Looks in the first column of a matrix, A, for a given value, z, and returns the value(s) in the same row(s) in the column specified, c. When multiple values are returned, they appear in a vector.

В М11 есть, по идее должны быть и М2000/2001

MathCAD: cделай это по-разному…

Просто пара-тройка “весёлых картинок” с занятий по MathCAD, чтоб не выкидывать, версия – Prime II или III, или 15-я, какая попалась:

4 способа извлечь строку из матрицы

2 способа найти скалярное произведение векторов

2 способа узнать длину вектора

Перестановка местами максимального и минимального элемента матрицы. На самом деле, условие не совсем корректно. Если максимальный и минимальный элементы единственны или нужны только первые их вхождения (или последние, если поменять знаки “<“, “>” на “≤”, “≥”), то достаточно следующего “школьного” и однопроходного кода:

Перестановка местами максимального и минимального элемента матрицы (только одна)

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

Перестановка местами максимального и минимального элемента матрицы (все вхождения)

Конечно, значения максимума и минимума можно было найти стандартными функциями max, min, но они всё равно бы делали такой же проход.

Тем не менее, обмен местами минимальных и максимальных значений легко было бы сделать и без программирования (приведён пример для вектора):

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

То же самое – с определением номера максимального элемента вектора. Корректно было бы вернуть вектор, состоящий из всех номеров максимальных элементов, на случай, если их несколько:

список номеров элементов вектора, которые совпадают с максимальным

В матрице найти номера строк, суммы элементов которых максимальны, учесть, что строк может быть больше одной:

Mathcad, номера строк, суммы элементов которых максимальны

Решаем символьно квадратное уравнение и сразу же упрощаем вид полученных mathcad формул:

Решаем символьно квадратное уравнение

09.06.2015, 16:19 [12829 просмотров]

---