Файл: Лабораторная работа № 8.doc

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

ЛАБОРАТОРНАЯ РАБОТА №8.

ТЕМА: РАСПОЗНАВАНИЕ ОБРАЗОВ.


Создание модели бензина.

В предыдущих лабораторных работах проводилось распознавание класса неизвестного бензина. Но после такого распознавания желательно определить свойства этого бензина. Каждая марка бензина должна по целому ряду показателей удовлетворять требованиям соответствующего ГОСТа или ОСТа.

В данной лабораторной работе мы будем определять три таких показателя: октановое число, полученное по исследовательскому методу, октановое число, полученное по моторному методу и фракционный состав бензина в начале кипения. В таблице 1 приведены требуемые ГОСТом значения этих показателей для трех рассматриваемых нами марок бензина. Таблица1

№№

Показатель

Марка бензина

А-76

АИ-95

АИ-92

1

Октановое число по исследоват. методу

>76

>95

>92

2

Октановое число по моторн.методу

>76

>85

>83

3

Фракционный состав при начале кипения

>35

>30

>35

Постановка задачи.

Имеется 12 спектров трех известных нам марок бензина ( по четыре для каждой марки).

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

Метод решения задачи.

  1. На основании спектров и паспортных данных (см. таблицу 2) известных бензинов необходимо построить математические модели определения каждого показателя.

  2. В эти математические модели вводятся спектры «неизвестных» бензинов и определяются показатели.

  3. Полученные показатели сравниваются с известными нам на самом деле показателями исследуемых бензинов (см. таблицу 3). Таким образом работоспособность моделей.

Таблица 2. Паспортные данные известных бензинов



Показатели

А-76

АИ-95

АИ-92

734

736

811

902

285

286

289

778

441

462

504

533

1

Октановое число по исследов. методу, ед.

-

81.0

80.7

79.8

96.2

96.1

95.3

95.0

93.4

92.6

93

92.6

2

Октановое число по моторн. методу, ед.

76

76.9

77.2

76.2

85.2

85.8

85.5

85.1

85.4

83.5

85

84.3

3

Фракционный состав: начало кипения, град C

37

38

39

35

32

32

32

40

37

36

35

45

Таблица 3 Паспортные данные «неизвестных» бензинов

№№

Показатели

А-76

АИ-95

АИ-92

1

Октановое число по исследов. методу, ед.

81.3

96.1

93.2

2

Октановое число по моторн. методу, ед.

76.8

85.8

85.1

3

Фракционный состав: начало кипения, град C

39

31

41


Решение задачи.

Первые шаги полностью повторяют решение лабораторных работ 5-7:

Перед началом решения командой ФАЙЛ – НАСТРОЙКА СТРАНИЦЫ зададим альбомный формат ( LANDSCAPE).

Все спектры расположены в папке МАТКАД – ДАННЫЕ ДЛЯ ЛАБ. РАБОТЫ. В ней созданы три папки для известных марок бензина и папка « неизвестных» спектров.

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

2.Сначала необходимо ввести все известные ( опорные) спектры в маткад. Это производится с помощью команд меню ВСТАВИТЬ- ДАННЫЕ – ВВОД ФАЙЛА. После ввода этих команд открывается окно FILE OPTIONS ( опции файла), в котором имеется кнопка « BROWSE» ( искать).Нажав на эту кнопку, откроем окно READ FROM FILE ( читай из файла) и укажем путь: МАТКАД –ДАННЫЕ ДЛЯ ЛАБ.РАБОТЫ – БЕНЗИН А- 76- А-76спектр1.Потом нажмем кнопку ОТКРЫТЬ.

Затем нажмем два раза кнопку ГОТОВО в окне FILE OPTIONS. В маткаде появится рамка с надписью А-76 …….txt. Присвоим ему имя F11.

Аналогично введем остальные спектры бензина А76, присвоив им имена F12, F13, F14, и все спектры бензинов АИ-95, присвоив им имена F2i, и все спектры бензина АИ-98, присвоив им имена F3i.Здесь индексы i меняются от 1 до 4.

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

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

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

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


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





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


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



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


  1. Сформируем векторы паспортных данных известных бензинов (см. рис.3)


Паспортные данные известных бензинов





























Рис.3. Векторы паспортных данных известных бензинов.

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

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





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

6.Мы будем строить математические модели, используя одну из возможных мер сходства – коэффициенты корреляции. Вычислим матрицы этих коэффициентов для спектров известных бензинов ( см. рис.5).



Корреляционные матрицы спектров








Рис.5. Матрицы коэффициентов корреляции.

7 . Математические модели для всех показателей и для всех марок бензина будем искать в виде




Здесь сorr(Y,Sj) – коэффициент корреляции между спектрами известных бензинов и спектром «неизвестного» бензина, а коэффициенты b нужно определить заранее в процессе обучения математической модели. Для каждого показателя и для каждого сорта бензина эти коэффициенты будут различными. Так как мы считаем, что принадлежность к той или иной марке для «неизвестных» бензинов уже определена в предыдущих лабораторных работах, то, чтобы не вводить каждый раз множество различных исходных данных в программу, составим несколько вспомогательных программ( см.рис.6). В них, в зависимости от марки исследуемого бензина мы автоматически формируем для подстановки в обучающие программы матрицы коэффициентов корреляции, паспортных данных, обучающих ( т.е. известных нам) спектров и «неизвестных» спектров.






















Рис.6 Вспомогательные программы формирования исходных данных

8 . Обучение будем проводить методом наименьших квадратов с помощью вычислительного блока GIVEN и встроенной функции MINERR( минимум ошибки). Для любого численного метода, как известно, должны быть заданы начальные приближения. У нас в процессе обучения ищутся коэффициенты b, которым мы задаем начальные значения, равные 1 ( для линейных задач, а наша задача линейна, величина начальных приближений не имеет значения). На рис. 7 показана программа задания начальных приближений, а на рис.8 приведен один из вычислительных блоков.






Рис.7. Программа и матрица начальных приближений

















Рис.8 Один из вычислительных блоков.


Здесь Z – вектор всех паспортных данных данного показателя для всех известных бензинов,

 - Матрица коэффициентов корреляции между всеми известными бензинами данной марки,

b – искомые коэффициенты.

Выражение



Является теоретической функцией, которой мы аппроксимируем заданные значения вектора Z. Аппроксимация производится методом наименьших квадратов, о чем говорит выражение







Компьютер подбирает такие значения коэффициентов b, при которых сумма квадратов разностей между «теоретическими» и практическими значениями показателя минимальна ( minerr- минимум ошибки). На рис.8 показаны вычисленные из этих условий значения коэффициентов b.

Так как мы вычисляем три показателя, то таких вычислительных блоков в нашей задаче должно быть три.

9. После обучения математической модели, т.е. после определения значений коэффициентов b, подставив в формулу




з начения найденных ранее коэффициентов b и матрицы коэффициентов корреляции между известными и неизвестными бензинами, определим значение искомого показателя. В данном примере мы искали для «неизвестного» бензина марки АИ-95 октановое число, определяемое по исследовательскому методу. Ниже приведен ответ 95,928. Легко убедиться, что он соответствует ГОСТу.



Ниже, на рис. 9 приведен полный листинг программы решения задачи.

Задание студентам.

  1. Разобраться в приведенном выше материале.

  2. Набрать программу по листингу.

  3. Определить значения трех показателей: октанового числа по исследовательскому методу;

Октанового числа по моторному методу;

Фракционного состава при начале кипения

Для трех марок «неизвестных» бензинов: А-76,АИ-95,АИ-92.

4. Сравнить полученные значения с паспортными данными для этих трех бензинов.