15.3.3. Устранение тренда

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

• Вычислить регрессию f(x), например линейную, исходя из априорной информации о тренде (предпоследняя строка листинга).

• Вычесть из данных у (х) тренд f (x) (последняя строка листинга).

Листинг 15.18. Устранение тренда

На рис. 15.22 показаны исходные данные (кружками), выделенный с помощью регрессии линейный тренд (сплошной прямой линией) и результат устранения тренда (пунктир, соединяющий крестики).

Рис. 15.22. Устранение тренда (листинг 15.18)

15.3.4. Полосовая фильтрация

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

Рис. 15.23. Результат полосовой фильтрации (листинг 15.19)

Алгоритм полосовой фильтрации приведен в листинге 15.19, а результат его применения показан на рис. 15.23 сплошной кривой. Алгоритм реализует такую последовательность операций:

• Приведение массива данных y к нулевому среднему значению путем его вычитания из каждого элемента у (третья и четвертая строки листинга).

• Устранение из сигнала у высокочастотной составляющей, имеющее целью получить сглаженный сигнал middle, например с помощью скользящего усреднения с малым окном w (в листинге 15.19 w=3).

• Выделение из сигнала middle низкочастотной составляющей slow, например, путем скользящего усреднения с большим окном w (в листинге 15.19 w=7) либо с помощью снятия тренда (см. разд. 15.3.3).

• Вычитание из сигнала middle тренда slow (последняя строка листинга), тем самым выделяя среднемасштабную составляющую исходного сигнала у.

Листинг 15.19. Полосовая фильтрация

Лекция 5. Числовой ввод-вывод данных. Создание двух- и трехмерных графиков.

16.1. Числовой ввод-вывод

Наиболее простой и распространенный ввод-вывод данных в Mathcad реализован присваиванием и (либо численным, либо символьным) выводом непосредственно в документе. Фактически документ Mathcad является одновременно и кодом программы и результатом ее выполнения. Поэтому самый простой и распространенный способ ввода-вывода — это непосредственное присвоение и вывод вычисленных значений в документах.

Листинг 16.1. Числовой ввод данных

Листинг 16.2. Числовой вывод данных (продолжение листинга 16.1)

16.2. Создание графиков

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

• Двумерные графики:

• XY (декартовый) график (XY Plot);

• полярный график (Polar Plot).

• Трехмерные графики:

• график трехмерной поверхности (Surface Plot);

• график линий уровня (Contour Plot);

• трехмерная гистограмма (3D Bar Plot);

• трехмерное множество точек (3D Scatter Plot);

• векторное поле (Vector Field Plot).

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

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

Некорректное определение данных приводит, вместо построения графика, к выдаче сообщения об ошибке.

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

• Поместите курсор ввода в то место документа, куда требуется вставить график.

• Если на экране нет панели Graph (График), вызовите ее нажатием кнопки с изображением графиков на панели Math (Математика).

• Нажмите на панели Graph (График) кнопку X-Y Plot для создания Декартового графика (рис. 16.1) или другую кнопку для иного желаемого типа графика.

• В результате в обозначенном месте документа появится пустая область графика с одним или несколькими местозаполкителями (рис. 16.1, слева). Введите в местозаполнители имена переменных или функций, которые должны быть изображены на графике. В случае Декартова графика это два местозаполнителя данных, откладываемых по осям х и Y.

Рис. 16.1. Создание Декартового графика при помощи панели Graph

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

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

Самый наглядный способ создания графика — с помощью панели инструментов Graph (График). Однако точно так же создаются графики путем выбора соответствующего элемента подменю Insert / Graph (Вставка / График), показанного на рис. 16.2, либо нажатием соответствующей типу графика горячей клавиши.

Чтобы удалить график, щелкните в его пределах и выберите в верхнем меню Edit (Правка) пункт Cut (Вырезать) или Delete (Удалить).

Рис. 16.2. Создание графика посредством меню