Ввод данных закончен

Решение.

Перед началом решения командой файл- настройка страницы зададим альбомный формат (landscape).

Все спектры расположены в папке Маткад- данные для лаб. работы. В ней созданы три папки для известных марок бензина и папка «неизвестных» спектров.

1. Для начала счета с единицы введем ORIGIN: =1

2.Сначала необходимо ввести все известные (опорные) спектры в Маткад. Это производится с помощью команд меню- вставить- данные-ввод файла. После ввода этих команд открывается окно file options (опции файла), в котором имеется кнопка

«browse» (искать). Нажав на эту кнопку, откроем окно read from file (читай из файла) и укажем путь: Маткад-данные для лаб. Работы- бензин А-76- спектр1. Потом нажмем кнопку открыть.

Затем нажмем два раза кнопку готово в окне file options. В Маткаде появится рамка с надписью А-76 …….txt. Присвоим ему имя F1₁.

Аналогично введем остальные спектры бензина А76, присвоив им имена F1₂, F1₃, F1₄, и все спектры бензинов АИ-95, присвоив им имена F2_i, и все спектры бензина АИ-98, присвоив им имена F3_i.Здесь индексы i меняются от 1 до 4.

3. Введенные спектры представляют собой матрицу из 7462 строк и 3 столбцов: первый столбец – порядковый номер, второй – частота и третий – амплитуда спектра (см. рис.2)

1: 399.95 -0.000131

2: 400.44 0.001253

3: 400.92 0.001301

4: 401.40 0.000805

5: 401.88 0.000681

6: 402.37 0.001579

7: 402.85 0.003724

8: 403.33 0.006809

9: 403.81 0.010327

Рис.2 Вид введенного спектра.

Затем сформируем сами матрицы как показано на рис.2. Цифра или буква в уголках в показателе матрицы обозначает соответствующий столбец данной матрицы.

Нас интересует только амплитуда. Поэтому сформируем из амплитуд всех введенных спектров три матрицы из четырех столбцов каждая: для бензина А-76 матрицу S1, для бензина АИ-95 матрицу S2 и для бензина АИ-98 матрицу S3.

Cначала введем заголовок «Формирование матриц». Для этого нажав одновременно кнопки shift и Э (в английском шрифте), откроем окно надписей. Перейдя на русский и введя в меню шрифт Arial cyr, запишем этот заголовок.

Затем сформируем сами матрицы как показано на рис.3

Формирование матриц

Рис.3 Формирование матриц.

4.Спектры введены в формате текстов. Их следует перевести в цифровой формат с помощью встроенной функции str2num. (string – строка, num- число). Перевод показан на рис.4

Рис.4. Перевод текстовых файлов в цифровые.

5. Подготовим «неизвестные» спектры. На самом деле спектры, которые мы называем «неизвестными» нам известны. Первый из этих спектров принадлежит бензину А-76, второй – бензину АИ-95, третий – бензину АИ-98. Мы должны проверить, правильно ли распознает их формируемая нами программа.

Введем «неизвестные» спектры, как мы это делали для известных спектров, дав им имена X₁, X₂, X₃,выделим третьи столбцы и переведем их в цифровую форму (см.рис.4).

Рис.5. Ввод исследуемых спектров.

6. Составим программу поиска максимального коэффициента корреляции каждым из «неизвестных» и всеми известными спектрами (см рис.6).

7.Так как у нас три «неизвестных» спектра, то мы организуем цикл по j от 1 до 3.

8. Переменными max₁,max₂,max₃ мы обозначим максимальные коэффициенты корреляции между спектрами бензинов А-76, АИ-95 и АИ-98, соответственно. В начале мы присвоим им нулевые значения.

9. В каждой марке бензинов мы имеем четыре известных спектра. Поэтому организуем цикл по «к» от 1 до 4-х. В этом цикле мы ищем максимальные коэффициенты корреляции между каждым из «неизвестных» спектров и спектрами каждой марки.

10. Эти вычисленные максимальные коэффициенты корреляции циклом по «i» мы помещаем в вектор MAX.

Следующим циклом по «i» мы помещаем в вектор «К» условный номер марки бензина. Группа 1 – бензин А-76, группа 2 – бензин АИ-95, группа 3 – бензин АИ-98.

12.Наконец, для большей наглядности присваиваем группе 1 имя «76», а остальным, соответственно, имена «95», «98».

13. ниже приведен вектор ответа. Видим, что распознавание произведено правильно.

Рис.6. Маткад - программа определения марки бензина по максимальному коэффициенту корреляции.

Дата добавления: 2015-06-27