автореферат диссертации по процессам и машинам агроинженерных систем, 05.20.01, диссертация на тему:Повышение эффективности функционирования МТА путем оптимизации скорости движения при выполнении сельскохозяйственных работ

ВВЕДЕНИЕ.

1 СОСТОЯНИЕ ВОПРОСА И ЗАДАЧИ ИССЛЕДОВАНИЙ

1.1 Управляемость и её влияние на режимы работы сельскохозяйственных МТА

1.2 Условия функционирования сельскохозяйственных МТА.

1.3 Влияние управляемости на агротехнические и эксплуатационно-технологические показатели МТА.

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

При работе трактора в агрегате с сельскохозяйственными машинами на первый план выдвигаются ограничения, связанные с агротехникой возделывания сельскохозяйственных культур, условиями труда водителя, воздействиями на окружающую среду. Эти ограничения напрямую зависят от свойства трактора - управляемости.

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

Поэтому в соответствии с нормативно-технической документацией на испытания тракторов предусматривается определение их управляемости [1, 2, 3]. Однако управляемость трактора при проведении испытаний оценивается субъективно. Количественная оценка этого важного показателя не определяется, что вызвано отсутствием метода и аппаратуры, позволяющих оценивать управляемость трактора при его предварительных или приёмочных испытаниях [4] инженерами-испытателями.

Разработанные в настоящее время методы оценки управляемости тракторов при проведении испытаний новых моделей неприемлемы в силу их сложности и наличия большого числа показателей, оценивающих отдельные стороны управляемости. Единого показателя управляемости, достаточного для оценки этого свойства трактора при испытаниях, не предложено. Оценка управляемости, изложенная в [5, 6], может использоваться только для частного случая, при представлении трактора как участника дорожного движения и не учитывает его основное назначение - выполнение полевых сельскохозяйственных работ в составе машинно-тракторного агрегата.

Таким образом, целью данной работы является повышение эффективности функционирования МТА путём оптимизации скорости движения с учётом управляемости при выполнении сельскохозяйственных работ.

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

Для реализации указанной цели необходимо было:

- предложить метод оценки управляемости трактора величиной удельной информации связи, рассчитанной по данным реализаций двух процессов: воздействие водителя на рулевое управление трактора (вход), угловое перемещение рамы трактора (выход);

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

Актуальным вопросом является также получение количественной оценки управляемости трактора МТЗ-80. Этого можно достигнуть с помощью расчёта допустимых значений колебаний рамы трактора с учётом агротехнических требований, например, при работе трактора МТЗ-80 в агрегате с культиватором КРН-5,6 и обосновать ограничение скорости его движения при междурядной культивации кукурузы.

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

Данная работа выполнялась в соответствии с планом научно-исследовательских и опытно-конструкторских работ КубНИИТиМ в 1998 г, отчёт № 26-98 по научно-исследовательской работе № 960132.7, УДК 628.114.2.075.001.4, гос. регистрация № 01.9.80007774.

На защиту выносятся следующие основные положения:

- показатель, оценивающий управляемость трактора и метод его определения;

- метод определения допустимых колебаний рамы трактора;

- оптимизация режима работы МТА;

- приборы и аппаратура.

Научную новизну работы представляют:

- новый информационный параметр (бит/с), оценивающий управляемость трактора - удельная информация связи системы: орган управления -рама трактора;

- метод определения этого информационного параметра;

- функциональная зависимость допустимой величины среднеквадрати-ческого отклонения угла поворота рамы трактора МТЗ-80 от величины защитной зоны при междурядной культивации с культиватором КРН-5,6;

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

- метод определения допустимых колебаний рамы трактора.

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

Разработана методика для определения оптимальной скорости движения МТА с учётом ограничений связанных с управляемостью трактора.

Результаты исследований автора внедрены в КубНИИТиМе (приложение 1) с СКОС ВЙМ (приложение 3).

ОБЩИЕ ВЫВОДЫ

1. Управляемость, как свойство трактора, должна определяться при проведении испытаний новой модели трактора с учётом его основного назначения - выполнение полевых сельскохозяйственных работ в составе машинно-тракторного агрегата. Наиболее существенной особенностью трактора как звена динамической системы МТА является его работа под воздействием управляющих и возмущающих сил, которые можно отнести к категории случайных в вероятностно-статистическом смысле. Исследования показали, что в качестве показателя, оценивающего управляемость трактора, может быть принят информационный параметр - удельная информация связи системы: орган управления - рама трактора.

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

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

3. Для регистрации процессов на входе и выходе для колёсного трактора разработана аппаратура, описание которой приводится в работе. Эта аппаратура позволяет за время контрольного опыта при эксплуатационно-технологических испытаниях МТА осуществить регистрацию процессов движения рулевого колеса и рамы трактора.

4. В результате проведения экспериментальных исследований с применением разработанной аппаратуры получены значения удельной информации связи. Эти значения составляют для трактора МТЗ - 80 в агрегате с культиватором КРН - 5,6 при междурядной обработке кукурузы 1,149 бит/с при скорости движения 7,63 км/ч и 1,190 при скорости движения 9,18 км/ч, и в агрегате с культиватором КПС-4,0 для сплошной обработке почвы 1,421 бит/с при скорости движения 6,12 км/ч и 1,452 бит/с при скорости движения 8,78 км/ч.

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

5. Для выбора максимальной скорости движения МТА с учётом ограничений, связанных с условиями труда водителя, трактор, как объект управления, следует представлять в виде звена замкнутой динамической системы.

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

Наиболее высокие требования к управляемости трактора предъявляются при работе его на междурядной обработке сельскохозяйственных культур, в частности, при работе трактора МТЗ-80 с культиватором КРН-5,6 на междурядной обработке кукурузы.

6. В результате теоретических и экспериментальных исследований получены значения допустимой величины среднеквадратического отклонения угла поворота рамы трактора МТЗ - 80 в зависимости от защитной зоны, принимаемой при междурядной культивации кукурузы МТА, состоящего из трактора МТЗ - 80 с культиватором КРН - 5,6.

При защитной зоне 10 см допустимое значение среднего квадратического отклонения рамы трактора, подсчитанное по формулам (2.35 - 2.37), исходя из агротехнических требований должно быть не более 0,905 град, при допустимой величине психофизической нагрузки на водителя. Допустимую психофизическую нагрузку на водителя при оптимальной частоте воздействий на рулевое управление трактора можно выразить значением среднеквад-ратического отклонения движения рулевого колеса, которое должно быть не более 28,9 град.

7. Установлено, что на междурядной культивации кукурузы оптимальной скоростью движения МТА является 7,63 км/ч. На этой скорости движения среднеквадратическое значение отклонения угла поворота рамы составляет 0,954 град, что определяет количество срезанных культурных растений не более 3% при допустимой величине психофизической нагрузки на водителя.

Повышение скорости движения МТА до 9,18 км/ч приводит к возрастанию колебаний рамы трактора. При максимальном напряжении труда водителя (среднеквадратическая величина колебаний руля составила 33,8 град) среднеквадратическая величина колебаний рамы трактора составила 1,592 град, что вызывает подрезание культурных растений, составляющее 13%.

8. Расчёт экономической эффективности работы МТА при выполнении междурядной культивации кукурузы показал, что эффект от повышения производительности МТА за счёт увеличения скорости движения получен значительно меньшим, чем ущерб, наносимый увеличением количества срезанных растений.

Экономический эффект от выбора оптимальной скорости движения МТА с учётом ограничений, связанных с условиями труда водителя, составляет 622,85 руб. на 1 га в год.

1. ГОСТ 7057-81. Тракторы сельскохозяйственные. Методы испытаний. Взамен ГОСТ 7057-73; Введ. 29.04.81. - М.: Изд-во стандартов, 1981. - 24 с. УДК 629.114.2.001.4. Группа Д29. СССР.

2. ГОСТ 25836-83. Тракторы. Виды и программы испытаний.

3. Коробейников А.Т. Испытания сельскохозяйственных тракторов / А.Т. Коробейников, B.C. Лихачёв, В.Ф. Шолохов. М.: Машиностроение, 1985. - 243 с.

4. Левитану с А. Д. Ускоренные доводочные испытания тракторов / А. Д. Левитанус. -М.: Машиностроение, 1983. 151 с.

5. ОСТ 37.001.051-86. Управляемость и устойчивость автомобилей. Термины и определения.

6. ОСТ 37.001.054-86. Управляемость и устойчивость автомобилей. Методы испытаний.

7. ГОСТ 24055-88. Методы эксплуатационно-технологической оценки. Введ. 01.01.89. М.: Изд-во стандартов, 1988. - 15 с. УДК 631.001.4:006.354. Группа Т51. СССР.

8. ГОСТ 24056-88. Методы эксплуатационно-технологической оценки машин на этапе проектирования. Введ. 01.01.89. М.: Изд-во стандартов, 1988. - 9 с. УДК 631.3.001.4:006.354. Группа Т51-СССР.

9. ГОСТ 24057-88. Методы эксплуатационно-технологической оценки машинных комплексов, специализированных и универсальных машин на этапе испытаний. Введ. 01.01.89. М.: Изд-во стандартов, 1988. - 8 с. УДК 631.3.001.4:006.354. Группа Т51. СССР.

10. ГОСТ 24059-88. Методы эксплуатационно-технологической оценки транспортных средств на этапе испытаний. Введ. 01.01.89. М.: Изд-во стандартов, 1988. - 9 с. УДК 631.3.001.4:006.354. Группа Т51. СССР.

11. Калоев А.В. Основы проектирования систем автоматического вождения самоходных машин / А.В. Калоев. Орджоникидзе: Ир, 1975. - 295 с.

12. А. с. 398197 МКИ В25 J 14/01. Способ ориентации самоходных машин в индукционном поле токонесущего провода / Зайченко В.В., Драновский А.И. (СССР). № 5267383/22 08; Заявл. 24.02.76; Опубл. 21.04.79, Бюл. № 6. - 44с.

13. Бахмутский М.М. Количественная оценка влияния передаточного числа механизма поворота на курсовую устойчивость автопоезда / М.М. Бахмутский, В.К. Шеронов // Труды НАМИ. М.: 1971. - вып. 129.-С. 32-35.

14. А. с. 138098 МКИ В25 J 18/02. Электроконтактная головка копирующего устройства для автоматического вождения трактора / И.Ч Логинов (СССР). № 8465123/11 02; Заявл. 21.10.73; Опубл. 21.08.75, Бюлл. № 1.-61 с.

15. Яровой А.Ч. Исследование бокового увода тракторных шин и его влияние на некоторые показатели криволинейного движения агрегата: Дис. канд. техн. наук. Харьков, 1970

16. Бахмутский М.М. Количественная оценка влияния передаточного числа механизма поворота прицепа на курсовую устойчивость автопоезда / М.М. Бахмутский, В.К. Шарапов // Труды НАМИ. М., 1971.-вып. 129.-С. 36-38.

17. Скотников В.А. Основы теории расчёта трактора и автомобиля / В. А. Скотников, А. А. Мещенский, А. С. Сол омский. М.: Агро-промиздат, 1986. -451 с.

18. Гуськов В.В. Тракторы. Теория / В.В. Гуськов, Н.Н. Велев, Ю.Е. Атаманов, Н.Ф. Бочаров, И.П, Ксиневич, А.С. Солонский. М.: Машиностроение, 1988.-611 с.

19. Фарабин Я.Е. Теория поворота тракторных машин / Я.Е. Фарабин. -М.: Машиностроение, 1970. 251 с.

20. Брянский Ю.А. Управляемость большегрузных автомобилей / Ю.А. Брянский. М.: Машиностроение, 1982. - 349 с.

21. Коновалов В.Ф. Динамическая устойчивость тракторов / В.Ф. Коновалов. М.: Машиностроение, 1981. - 674 с.

22. Адамович Н.В. Управляемость машин. Эргономические основы оптимизации рабочего места человека-оператора / Н.В. Адамович. -М.: Машиностроение, 1977. 277 с.

23. Калинов А.Я., Подсистема САПР для расчётной оценки и управляемости МТА / А.Я. Калинов, Ф.В. Кальянов, А.В. Черницер // Тракторы и сельскохозяйственные машины, №2, 1988. С 67.-71

24. Котик М.А. Саморегуляция и надёжность человека-оператора / М.А. Котик. Таллин: Валгус. 1974. - 470 с.

25. Старожук И.А. Влияние типа механизма поворота трактора на напряжённость труда тракториста / И.А. Старожук, Г.С. Цейтлина, B.C. Сафронов // Тракторы и сельхозмашины, № 2, 1998. С. 11-22.

26. Гетьман Н.И. Оценка условий труда механизаторов / Н.И. Гетьман, А.В. Калюжный, И.Э. Липковйч // Механизация и электрификация сельского хозяйства, № 7, 2001. С. 45 - 49.

27. Хаматуров А.А. Динамика системы: дорога, шина, автомобиль, водитель / А.А. Хаматуров, В.А. Афанасьев, B.C. Васильев. М.: Машиностроение, 1976. - 239 с.

28. Лавров Т.М. Некоторые вопросы управляемости тракторных агрегатов / Т.М. Лавров // Труды НАТИ. М., 1975. - вып. 243. - С. 40 -46.

29. Лурье А.Б. Статистическая динамика сельскохозяйственных агрегатов / А.Б. Лурье. -М.: Колос, 1981. 382 с.

30. Лурье А.Б. Расчёт и конструирование сельскохозяйственных машин /А.Б. Лурье, А.А. Громбчевский. -Л.: Машиностроение, 1977.-528 с.

31. Гельденбейн С.П, Основы автоматизации сельскохозяйственных агрегатов / С.П. Гельденбейн. М.: Колос, 1975. - 490 с.

32. Кербер В.Н. Повышение уровня функционирования сельскохозяйственных агрегатов на основе их моделирования: Дисс. д-ра техн. Наук. Санкт-Петербург, 1993

33. Борский И.Б., Динамика трактора / И.Б. Борский, В.Я. Анилович, Г.М. Кутьков. М.: Машиностроение, 1975. - 390 с.

34. Жалнин Э.В. Математическое моделирование процессов земледельческой механики / Э.В. Жалнин // Тракторы и сельхозмашины, № 1,2000.-С. 123-130.

35. Бурков Л.Н. Концепция общей теории с/х машин / Л.Н. Бурков // Тракторы и сельхозмашины, № 8, 1998. С. 34-38.

36. Ксиневич И.П. Системы автоматического управления ступенчатыми трансмиссиями тракторов / И.П. Ксиневич, В.П. Тарасик. М.: Машиностроение, 1979. - 581 с.

37. Свирщевский А.Б. Технологические основы автоматизации сельскохозяйственного производства. / А.Б. Свирщевский, С.П. Гель-денбейн. М.: Колос, 1966. - 446 с.

38. Хачатрян Х.А. Стабильность работы почвообрабатывающих агрегатов / Х.А. Хачатрян. М.: Машиностроение, 1974. - 601 с.

39. Фортуна В.И. Основы внедрения прогрессивных технологий возделывания и уборки сельскохозяйственных культур и их обеспечение / В.И. Фортуна, Г.Г. Маслов, Е.И. Трубилин. Краснодар, 1995. - 140 с.

40. Лисютин В.П. Дифференциальная обработка междурядий / В.П. Лисютин // Кукуруза и сорго. № 3, 1988. С. 23-29.

41. ГОСТ 1114-84. Культиваторы пропашные. Типы и основные параметры.

42. ОСТ 10.4.3-99. Машины и орудия для обработки пропашных культур. Методы оценки функциональных показателей.

43. Саакян Д.Н. Система показателей комплексной оценки мобильных машин / Д.Н. Саакян. М.: Агропромиздат, 1988. - 415 с.

44. Гячев Л.В. Устойчивость движения сельскохозяйственных машин и агрегатов / Л.В. Гячев. М.: Машиностроение, 1981. - 358 с.

45. Царапкин О.И. Оптимизация режимов преобразования случайных процессов при экспериментальных исследованиях тракторов / О.И. Царапкин // Труды НАТИ. Мм 1976. - Вып. 244, С. - 23-25.

46. Каспаров Б.М. О методики исследования устойчивости прямолинейного движения сельскохозяйственного агрегата / Б.М. Каспаров // Труды ВИМ. М., 1970. - № 50. - С. 34-40.

47. Каспаров Б.М. Об устойчивости прямолинейного движения управляемых колёсных навесных сельскохозяйственных агрегатов / Б.М. Каспаров // Труды ВИМ. М., 1972. - № 60. - С.15-21.

48. Каспаров Б.М. О повышении устойчивости прямолинейного движения машинотракторного агрегата / Б.М. Каспаров // Труды ВИМ.- М., 1984. № 101. - С. 51-60 .

49. Стратонович P.J1. Теория информации / P.JI. Стратонович. М.: Советское радио, 1975. - 422 с.

50. Драновский А .И. Исследование систем управления тракторами методами информации / А.И. Драновский, JI.P. Беспрозванный // Труды НАТИ.-М., 1972.-вып. 29.-С. 11-19

51. Драновский А.И. Автоматическое вождение колёсных и гусеничных машин по постоянным трассам / А.И. Драновский, И.И. Кри-нецкий. М.: Машиностроение, 1971. - 341 с.

52. Драновский А.И. Новые методы исследования систем вождения тракторных агрегатов / А.И. Драновский // Труды НАТИ. М., 1975.- вып. 243. С. 24-26

53. Драновский А.И. Анализ управления в системе «водитель-трактор» / А.И. Драновский, B.C. Сафронов // Труды НАТИ. М., 1972. -вып. 219.-С. 31-38

54. Коган И.М. Прикладная теория информации / И.М. Коган. М.: Радио и связь, 1981. - 283 с.

55. Востроухов Н.Ч. Информационно-измерительная техника / Н.Ч. Востроухов, Н.Н. Евтихеев. -М.: Высшая школа, 1977. 561 с.

56. Nicolaos S. Trannes et al. Entropy and data compression IEEE Trans. Of Jnf. Theory, IT-16, № 1, January, 1970

57. Кендалл M. Теория распределений / M. Кендалл, А. Стьюарт. М.: Наука, 1966.-860 с.

58. Фрейнбруд Г.Н. Синтез оптимальной системы автоматического вождения тракторного агрегата со следовой программой движения / Г.Н. Фрейнбруд // Труды НАТИ. М., 1982. - вып. 5. - С. 15-18

59. Кодоев Р.В. Динамика пропашеого культиваторного агрегата как объекта автоматизации: Дие. канд. техн. наук. Орджоникидзе, 1973

60. Калаев С.С. Исследование и разработка системы автоматического вождения пропашного культиваторного агрегата с коррекцией движения по уклону: Дис. канд. техн. наук. Орджоникидзе, 1973

61. Надыкто В.Т. Управляемость и устойчивость движения агрегата на основе МЭС / В.Т. Надыкто // Тракторы и сельхозмашины, № 7, 1998

62. Стаюжук И.А. Влияние типа механизма поворота трактора на напряжённость труда тракториста / И.А. Стаюжук, Г.С. Цейтлина, B.C. Сафронов // Тракторы и сельхозмашины, № 2,1999

63. Галимжанов В.Н. Анализ формирования люфта рулевого управления / В.Н. Галимжанов, А,П. Егоров И Тракторы и сельхозмашины, № 12, 2000

64. Бурдин В.В. Координирование на местности МТА для междурядной обработки пропашных культур / В.В. Бурдин, А.Н. Площадинов // Тракторы и сельхозмашины, № 5, 2000

65. Рославцев А.В. Средства исследования движения МТА / А.В. Ро-славцев // Тракторы и сельхозмашины, № 3, 1999

66. Площадинов А.Н. Улучшение управляемости МТА на базе универсально пропашного трактора класса 2 путём рационального выбора параметров направляющих колёс: Автореферат дис. канд. техн. наук. М: 1987

67. Протокол № 07-52-00 приёмочных испытаний прибор для определения свободного хода рулевого колеса тракторов и самоходных с/х машин. Новокубанск: КубНИИТиМ, 2000

68. Инструкция по эксплуатации аппаратуры ИП 238. - Новокубанск: КубНИИТиМ, 2001

69. Вентцель Е.С. Теория вероятности / Е.С. Вентцель. М.: Высшая школа, 1969.-590 с.

70. Смирнов Н.В. Курс теории вероятности и математической статистики / Н.В. Смирнов, И.В. Дунин-Барковский. М.: Наука, 1965. -461 с.

71. Райбман Н.С. Построение моделей процессов производства / Н.С. Райбман, В.М. Чадеев. М.: Советское радио, 1975. - 280 с.

72. Разработка методики оценки управляемости движения мобильных энергетических средств: Отчёт о НИР / МИИСП. М., 1989

73. Агеев JI.E. Основы расчёта оптимальных и допустимых режимов работы машинотракторных агрегатов / J1.E. Агеев. JL: Колос, 1978. - 295 с.

74. Бендат Д. Измерение и анализ случайных процессов / Д. Бендат, А. Пирсон. М.: Мир, 1974. - 465 с.

75. Лившиц Н.А. Вероятностный анализ систем автоматического управления / Н.А. Лившиц, В.Н. Пугачёв. М.: Советское радио, 1969. - 896 с.

76. Мирский Г.Я. Аппаратурное определение характеристик случайных процессов / Г.Я. Мирский. М.: Советское радио, 1972. - 485 с.

77. Finkelstein J. Fundamental concepts of measurement //Acta IMEKO. -1973. Vol. 1 - Budapest, 1974.

78. Чеголин П.М. Автоматизация спектрального и корреляционного анализа / П.М. Чеголин. М.: 1969. - 311 с.

79. Свешников А.А. Прикладные методы теории случайных функций / А.А. Свешников. М.: Наука, 1968. - 463 с.

80. Рославцев А.В. Методы исследования движения машинотракторных агрегатов / А.В. Рославцев, В.А. Хоустов, В.М. Авдеев // Тракторы и сельхозмашины № 6, 1998. С. 13-16

81. Хаймис Г.А. Исследование сельскохозяйственной техники и обработка опытных данных / Г.А. Хаймис, М.М. Ковалёв. М.: Колос, 1994.-341 с.

82. Доспехов Б.А. Методика полевого опыта / Б.А. Доспехов. М.: Аг-ропромиздат, 1985. - 289с.

83. Сборник агротребований на тракторы и сельскохозяйственные машины, том 31, М: ЦНИИТЭИ, 1982

84. Нкимов В.В. Вероятностные методы в строительной механике корабля / В.В. Нкимов. М.: Мир, 1974. - 290 с.

85. Рубинов A.M. Элементы математического анализа / A.M. Рубинов, К.Ш. Шапиев. М.: Просвещение, 1972. - 277 с.

86. Маркушевич А.И. Исследования по современным проблемам теории функции комплексного переменного / А.И. Маркушевич. М.: Наука, 1961.-514 с.

87. Джонсон Н. Статистика и планирование эксперимента в технике и науке. Методы обработки данных / Н. Джонсон, Ф. Лион. М.: Мир, 1980.-456 с.

88. Костюк А.Ф. Методы и аппаратура для анализа характеристик случайных процессов /А.Ф. Костюк, В.В. Ольшевский, Э.И. Цветков. -М.: Мир, 1967.-751 с.

89. Мирский Т.Я. Аппаратура для определения характеристик случайных процессов / Т.Я. Мирский. М.: 1972. - 401 с.

90. Колмогоров А.Н. Основные понятия теории вероятностей / А.Н. Колмогоров. М.: Наука, 1974. - 119 с.

91. Колмогоров А.Н. Элементы теории функций и функционального анализа / А.Н. Колмогоров, С.В. Фомин. М.: Наука, 1989. - 621 с.

92. Тьюки Д. Анализ результатов наблюдений. Пер. с англ. / Д. Тьюки. М.: Мир, 1981.-684 с.

93. Марчук Г.И. Методы вычислительной математики / Г.И. Марчук. -М.: Наука, 1980.-535 с.

94. Турчак Л.И. Основы численных методов: Учеб. Пособие / Л.И. Тур-чак. М.: Наука, 1987. - 320 с.

95. Цветков Э.И. Процессорные измерительные средства / Э.И. Цветков. Л.: Энергоатомиздат, 1989. - 224 с.

96. Феллер В. Введение в теорию вероятностей и её приложения. В 2-х томах. Т. 1: Пер. с англ. / В. Феллер. М.: Мир, 1984. - 528 с.

97. Перевёрткин С.М. МикроЭВМ в информационно-измерительных системах / С.М. Перевёрткин, НИ. Гаранин, Ю.Н. Костин, И.И. Миронов. М.: Машиностроение, 1987. - 245 с.

98. Чемоданов Б.К. Математические основы теории автоматического регулирования / Б.К. Чемоданов. М.: Высшая школа, 1977. - 452 с.

99. Гнеденко Б.В. Курс теории вероятности / Б.В. Гнеденко. М.: Наука, 1965.-350 с.

100. Кендалл М. Статистические выводы и связи. Пер. с англ. / М. Кендалл, А. Стьюарт М.: Наука, 1973. - 587 с.

101. Курош А.Г. Курс высшей алгебры / А.Г. Курош. М.: Наука, 1962. -431 с.

102. Сергеев В.A. Visual Basic 6.0 / В.А. Сергеев. СПб.: БХВ - Санкт-Петербург, 1999. - 720 с.

103. Дьяконов В.В. MathCAD 8/2000 / В.В. Дьяконов. СПб.: Питер, 2000. - 300 с.

104. Джеффри М. М. Обработка баз данных на Visual Basic 6.0. Пер. с англ. / М.М. Джеффри. М.: Вильяме, 1999. - 911 с.

105. Михеева В.Д. Microsoft Access 2000 / В.Д. Михеева, И.А. Харитонова. СПб.: БХВ - Санкт-Петербург, 1999. - 789 с.

106. Фролов А.В., Фролов Г.В. Операционная система Windows 95 для программистов / А.В. Фролов, Г.В. Фролов. М.: Диалог - МИФИ, 1995.-267 с.

107. Макашарипов С.В. Программирование баз данных на Visual Basic 5.0 в примерах / С.В. Макашарипов. СПб.: Питер, 1997. - 305 с.

108. DAQCard Е Series User Manual. Austin. National Instruments, 1996

109. The Measurement and Automation 2000. Austin. National Instruments, 2000

110. Measurement Studio. Austin. National Instruments, 2000

111. Direct Data Acquisition and Instrument Control with Microsoft Excel. Austin. National Instruments, 1998

112. Агропромышленный комплекс Кубани. Статистический сборник 2000 г. Краснодар, 2001

113. ВМ 13.025 84. Временная методика определения сравнительной экономической эффективности мероприятий, направленных на улучшение санитарно-гигиенических показателей условий труда на новой (модернизированной) сельскохозяйственной технике.

114. ГОСТ 12.2.019 86. Тракторы и машины самоходные сельскохозяйственные. Общие требования безопасности. - М.: Изд-во стандартов, 1986.- 11 с.

115. Антошкевич B.C. Совершенствование методик и нормативов определения экономической эффективности новой техники / B.C. Антошкевич. М.: ВИСХОМ, 1983

116. ГОСТ 23728 88 - ГОСТ 23730- 88. Техника сельскохозяйственная. Методы экономической оценки. - М.: Изд-во стандартов, 1988. - 25 с. УДК 631.3.004.15:006.354. Группа Т51. СССР.

117. Самойленко Е.М. Метод оценки новой сельскохозяйственной техники по социально-экономическому критерию / Е.М. Самойленко //

118. Испытание и исследование новой сельскохозяйственной техники и технологий: Сб. трудов КубНИИТиМ. Новокубанск, 1997. - С 67-74.

119. ОСТ 10.18 2001. Стандарт отрасли. Испытания сельскохозяйственной техники. Методы экономической оценки.

120. Методы определения экономической эффективности технологий и сельскохозяйственной техники. М., 1998.1. Утверждаю

121. Дире££рр КубНИИТиМ . , прцф: ^>А.Т. Табашмиков2001 г.1. Акт

122. Предложенная методика позволяет устанавливать ранее не определяемые функциональные показатели тракторов и сельскохозяйственных машин, определяющие условия выбора функциональных режимов эксплуатации сельскохозяйственной техники.

123. Председатель Члены комиссии:1. Ю.Г. Коваленко1. А.Б. Иванов1. Утверждаю

124. Директор КубНИИТиМ t Д/пн., Проф. -A А.Т. Табашникоиc&^ribSy* 2001 г." ; . У " г1. Акт

125. В 2002 году начат выпуск опытной партии ИИС ИП-238 по заявкам машиноиспы тательных и научно-исследовательских станций стран СНГ.

126. Председатель Члены комиссии:

127. Я. Любашин А.Н. Красиков Н.И. Ладоша С.М. Деняк1. АКТо производственной проверке законченных научно-исследовательских и опытно-конструкторских работ1. НИР и ОКР)

128. Наименование научно-исследовательского, проектио-конструкторской организации-разработчика: Кубанский ордена Ленина Научно-исследовательский институт по испытанию тракторов и сельскохозяйственных машин

129. Наименование законченных НИР и ОКР, поставленную па производственную проверку: Аппаратура для оценки управляемости МТА, метод выбора оптимальной скорости движения с учётом ограничений, связанных с управляемостью трактора

130. Авторы законченных НИР и ОКР: зав. лабораторией Таркивский В.Е. зав. отделом, канд. техн. наук. Поляков О.А., зав. отделом Ладоша Н.И.

131. Законченные НИР и ОКР рекомендованы к производственной проверке решением Учёного совета (научно-методической комиссии):

132. Кубанского научно-исследовательского института по испытанию тракторов и сельскохозяйственных машин. Протокол № 12 от 5 декабря 2000 г.

133. Производственная проверка проводилась в ГНУ Северо-Кавказская опытная станция ВИМ

134. Ответственные за проведение производственной проверки: канд. техн. наук. Хваталин Ю. А.

135. Условия проведения проверки: Работа трактора МТЗ-ЙО с культиватором КРН-5.6 при междурядной обработке кукурузы на Северо-Кавказской опытной станции ВИМ

136. Объём производственной проверки25 га• ■

137. Сроки проведения производственной проверки: май-август 2001 г.

138. Ме тодика производственной проверкиопределена ГОСТ 2455 2457 - 81, ГОСТ 7057 - 81. ГОСТ 25836 — 83

139. Что рекомендуется для внедрения в производство Методика для определения оптимальной скорости движения МТА и аппаратура для 9 определения управляемости трактора „.

140. Ответет ве н и ы й п редставител ь гваетров Е. В.

141. Ответственный представиселг> научной организации1. Таркивский В.Г 2001 г.

142. Результаты, полученные при проведении полевого опыта на междурядной культивации.

143. Рабочая скорость 7,63 км/ч Рабочая скорость 9,18 км/ч

144. Время, Исходные дан- Данные после остацио- Исходные дан- Данные после оста-с ные наривания и центриро- ные ционаривания и цен-вания трирования

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

146. Рабочая скорость 6,12 км/ч Рабочая скорость 8,78 км/ч

147. Время, с Исходные дан- Остационаренные и Исходные дан- Остационаренные иные центрированные дан- ные центрированные дан-ные ные

148. Программа для ЭВМ для расчёта удельной информации связи. Программа написана на языке Visual Basic for Application для среды Microsoft Excel.

149. Dim Ег(ЮОО) As Boolean If Where = 1 Then Sheets.Add Else

150. Workbooks.Add End If Sum = 0 E = 0 XX = 41. KT = к / w1. Com Then

151. Mes = " между (" + MX + " и " + MY + "): " Else1. Mes = " End If

152. For i = 0 To KT x = 1 ((DataV(i) л 2) / (DataX(i) * DataY(i))) If x <= 0 Then x = l

153. Er(i) = True E = E + 1 End If

154. Range("A" & XX + 1 + i).Value = -((Log(x) / 0.693147180559945) * W) / 6.283185307

155. Sum = Sum + (Log(x) / 0.693147180559945) * W Next

156. UIS =-Sum/6.283185307 Range("B2:K" & E + 4).Select With Selection1.terior. Co lorlndex = 36 End With

157. With Selection.Borders(xlEdgeBottom) .LineStyle = xlContinuous .Weight = xlMedium .Colorlndex = xlAutomatic End With

158. With Selection.Borders(xlEdgeTop) .LineStyle = xlContinuous .Weight = xlMedium .Colorlndex = xlAutomatic End With

159. With Selection.Borders(xffidgeRight) .LineStyle = xlContinuous .Weight = xlMedium .Colorlndex = xlAutomatic End With

160. With Selection.Borders(xffidgeLeft) .LineStyle = xlContinuous .Weight = xlMedium .Colorlndex = xlAutomatic End With

161. Range("B2:K2").Borders(xlEdgeBottom). Weight = xlThin Range("B2"). Select With Selection .Font. Size = 12 .Font.Italic = True

162. FormulaRlCl = "Удельная информация связи" & Mes & UIS End With

163. Range("B4").FormulaR 1С1 = "Пропущено точек: " & E Range("B3 "). FormulaR 1С1 = "Окно спектра от 0 до " & К & " рад"

164. For i = 0 То КТ If Er(i) Then XX = XX + 1 Range("B" & XX).Select With Selection .Font.Colorlndex 3

165. FormulaRlCl = "Пропущена точка на частоте " & i * W & "." End With End If Next1. End Sub

166. Программа для ЭВМ для расчёта автокорреляционной функции процесса. Программа написана на языке Visual Basic for Application для среды Microsoft Excel.

167. Sub Корфункция(Graf, Where, Com, DataSource, K, KtK, T) ' Graf рисовать график (boolean)

168. Where расположение выходных данных (1-новый лист,2-новая книга) ' Com - в первой ячейке - описание (boolean) ' DataSource - первая ячейка данных ' К - количество ординат ' KtK - количество точек корфункции ' Т - смещение ординат процесса

169. Dim А(1000), KF(1000), KFN(IOOO) As Variant Sum = 0считывание значений ординат If Com Then Mes = Range(DataSource)

170. M = Va.(Right$(DataSource, Len(DataSource) 1)) Else

171. M = Val(Right$(DataSource, Len(DataSource) -1)) -1 End If

172. For i = 1 To К A(i) = Range(Left$(DataSource, 1) & i + M) Sum = Sum + A(i) Next

173. Mat = Sum / К 'Математическое ожидание процесса If Where = 1 Then Sheets.Add Else1. Workbooks. Add End If

174. SName = ActiveSheet.Name ' Построение таблицы значений корфункции Columns("C").ColumnWidth= 11.57 Columns("D:E").Column Width = 26.71 Range("C2"). Select With Selection .Font.Size = 12 .Font.Bold = True .Font.Italic = True

175. FormulaRlCl = "Таблица значений автокорреляционной функции" End With

176. Range("C3").FormulaRlCl = "Смещение, с" Range("D3").FormulaRlCl = "Ненормированные значения" Range("E3").FormulaRlCl = "Нормированные значения"

177. Range("C3:E" & KtK + 4).Select With Selection.Borders(xlEdgeLeft) .LineStyle = xlContinuous .Weight = xlMedium .Colorlndex = xlAutomatic End With

178. With Selection.Borders(xlEdgeTop) .LineStyle = xlContinuous .Weight = xlMedium .Colorlndex = xlAutomatic End With

179. With Selection.Borders(xlEdgeBottom) .LineStyle = xlContinuous .Weight = xlMedium .Colorlndex = xlAutomatic End With

180. With Selection.Borders(xlEdgeRight) .LineStyle = xlContinuous .Weight = xlMedium .Colorlndex = xlAutomatic End With

181. Range("D3:D" & KtK + 3).Select With Selection.Borders(xlEdgeRight) .LineStyle = xlContinuous .Weight = xlThin .Colorlndex = xlAutomatic End With

182. Range("C3:C" & KtK + 3).Select With Selection.Borders(xlEdgeRight) .LineStyle = xlContinuous .Weight = xlThin .Colorlndex = xlAutomatic End With

183. Range("C3 :E3 ") .Select With Selection.Borders(xlEdgeBottom) .LineStyle = xlContinuous .Weight = xlThin .Colorlndex = xlAutomatic End With

184. Range("C" & KtK + 4 & & "E" & KtK + 4).Select With Selection.Borders(xlEdgeTop) .LineStyle = xlContinuous .Weight = xlThin .Colorlndex = xlAutomatic End With

185. Range("C" & KtK + 4).FormulaRlCl = "Математическое ожидание процесса: " &1. Mat

186. Расчёт значений автокорфункции For i = 1 То KtKх = О

187. For j = 1 То (К i + 1) х — х + ((АО) - Mat) * (AG + i -1) - Mat»1. Next1. KF(i) = x / (К i + 1)

188. Range("C" & i + 3).FormulaRlCl = T * (j -1) Range("D" & i + 3).FormulaRlCl = KF(i) Range("E" & i + 3).FomiulaRl CI = KF(i) / KF(1) § Next

189. Построение графика корфункции If Graf Then Charts.Add

190. ActiveChart.ChartType = xlLineMarkers

191. ActiveChart.SetSourceData Source:=Sheets(SName).Range("E4:E" & KtK + 3),1. PlotBy:=xlColumns

192. ActiveChart. Location Where :=xlLocationAsObject, Name:=SName ► ActiveChart.SeriesCollection( 1).XValues = & SName & "!R4C3:R" & KtK + 3"C3"

193. With ActiveChart .HasTitle = True

194. ChartTitle.Characters.Text = "График нормированной корфункции (" & Mes &liyi

195. Axes(xlCategory, xlPrimary).HasTitle = True

196. Axes(xlCategory, xIPrimary).AxisTitle.Characters.Text = "Смещение, с" .Axes{xlValue, xlPrimary) .HasTitle = False End With• With ActiveChart.Axes(xlCategory)

197. HasMajorGridlines = True .HasMinorGridlines = False End With

198. With ActiveChart.Axes(xlValue) .HasMajorGridlines = True .HasMinorGridlines = False End With

199. ActiveChart.HasLegend = False ActiveChart.HasDataTable = False ' ActiveSheet.Shapes('^Harp. 1 ").ScaleWidth 1.3, msoFalse,msoScaleFromBottomRiglit

200. ActiveSheet.Shapes("Диагр. l").ScaleWidth 1.5, msoFalse, msoScaleFromTopLeft ActiveSheet.Shapes("flnarp. 1 ").IncrertientLeft -2.25 ActiveSheet.Shapes("Диагр. l").IncrementTop 108 End If End Sub

201. Программа для ЭВМ для расчёта спектральной плотности процесса. » Программа написана на языке Visual Basic for Application для среды Microsoft1. Excel.

202. К максимальное значение спектральной плотности

203. W приращение частоты спектральной плотности

204. F Количество точек корфункции, используемых для расчётов

205. Т частота квантования корфункции

206. Dim KF(1000), S(1000), А(1000), KFN(IOOO) As Variant

207. Считывание массива значений ординат If Com Then

208. Mes = Range(DataSource) » M = Val(Right$(DataSource, Len(DataSource) -1))1. Ms = Mes

209. Mess = "(" + Mes + ")" Else

210. M = Val(Right$(DataSource, Len(DataSource) -1» -1 Ms = "" Mess = "" End If

211. For i = 1 To KK A(i) = Range(Left$(DataSource, 1) & i + M) Sum = Sum + A(i) Next

212. Mat = Sum / КК 'Математическое ожидание процесса KT = К / W 'Количество точек спектральной плотности If Mode Then If Where = 1 Then1. Sheets Add Else1. Workbooks.Add End If

213. SName = ActiveSheet.Name Columns("D:E").Column Width = 19 Range("D2").Select With Selection .Font. Size = 14 .Font.Bold = True .Font.Italic True

214. FormulaRlCl = "Таблица значений нормированной" End With

215. Range("D3"). Select With Selection .Font.Size = 14 .Font.Bold = True .Font. Italic = True

216. FormulaRlCl = " спектральной плотности" End With

217. Range("D4").FormulaRlCl = "Частота" Range("E4").Formu.aRlCl "Значения" Range("D4:E" & KT + 6).Select With Selection.Borders(xlEdgeLeft) .LineStyle = xlContinuous .Weight = xlMedium .Colorlndex = xIAutomatic End With

218. With Selection.Borders(xlEdgeTop) .LineStyle = xlContinuous .Weight = xlMedium .Colorlndex = xIAutomatic End With

219. With Selection.Borders(xlEdgeBottom) .LineStyle = xlContinuous .Weight = xlMedium .Colorlndex = xIAutomatic End With

220. With Selection.Borders(xlEdgeRight) .LineStyle = xlContinuous .Weight = xlMedium .Colorlndex = xIAutomatic End With

221. Range("D4:D" & KT + 5).Select With Selection.Borders(xlEdgeRight) .LineStyle = xlContinuous Weight = xlThin .Colorlndex = xIAutomatic End With

222. Range("D4"). Select With Selection.Borders(xlEdgeBottom) .LineStyle = xlContinuous .Weight xlThin .Colorlndex = xIAutomatic End With1. Range("E4") .Select

223. With Selection.Borders(xlEdgeBottom) .LineStyle = xlContinuous .Weight = xlThin .Colorlndex = xIAutomatic End With

224. Range(MD" & KT + 6 & ":EM & KT + 6).Select With Selection.Borders(xlEdgeTop) .LineStyle = xlContinuous .Weight = xlThin

225. Colorlndex = xlAutomatic End With

226. Range("D" & KT + 6).FormulaRlCl = "Использовано " & F & " точек корфункции"1. End If

227. Расчёт значений автокорфункции For i = 1 To F x = 01. For j = 1 To (KK i + 1)x = x + ((AG) Mat) * (AG -H - 1) - Mat)) Next

228. KF(i) = x / (KK i + 1) KFN(i) = KF(i) / KF(1) Next i

229. Data(i) = Sum * Т/ 3.14159254 If Mode Then Range("D" & i + 5).FormulaRlCl = W * i Range("E" & i + 5).FormulaRlCl = Data(i) End If Next

230. Построение графика нормированной спектральной плотности If Graf And Mode Then Charts. Add

231. ActiveChart.ChartType = xlLineMarkers

232. ActiveChart.SetSourceData Source :=Sheets(SName).Range("E5:E" & KT + 5),1. PlotBy:=xlColumns

233. ActiveChart.Location Where:=xlLocationAsObject, Name:=SName ActiveChart.SeriesCollection(l).XValues = "=" & SName & "!R5C4:R" & KT + 5 &1. C4"

234. With ActiveChart .HasTitle = True

235. ChartTitle.Characters.Text = "График нормированной спектральной плотности& Mess

236. Axes(xlCategory, xlPrimary).HasTitle = True

237. Axes(xlCategory, xlPrimary).AxisTitle.Characters.Text = "Частота, рад" .Axes(xlValue, xlPrimary).HasTitle = False End With

238. With ActiveChart.Axes(xlCategory) .HasMajorGridlines = True .HasMinorGridlines = False1. End With

239. With ActiveChart.Axes(xlValue) .HasMajorGridlines = True .HasMinorGridlines = False End With

240. ActiveChart.HasLegend = False ActiveChart.HasDataTable = False

241. Программа для ЭВМ для расчёта взаимной корреляционной функции процесса. Программа написана на языке Visual Basic for Application для среды Microsoft Excel.

242. Sub Взаимнаякорфункция(ОгаП, Graf2, Where, Com, DataSourcel, DataSource2, K, KtK, T)1. К количество ординат

243. KtK количество точек взаимной корфункции

244. Ml = Val(Right$(DataSourcel, Len(DataSourcel) -1)) 1 M2 = Val(Right$(DataSource2, Len(DataSource2) -1)) -1 End If

245. For i = 1 To К A(i) Range(Left$(DataSourcel, 1) & i + Ml) B(i) = Range(Left$(DataSource2, 1) & i + M2) Suml = Suml + A(i) Sum2 = Sum2 + B(i) Next

246. Matl = Suml / К 'Математическое ожидание процессов Mat2 = Sum2 / К If Where = 1 Then

247. Sheets. Add Else Workbooks.Add End If

248. SName = ActiveSheet.Name ' Построение таблицы значений корфункции Columns("C").ColumnWidth = 11.57 Columns("D:E").Column Width = 26.71 Range("A2"). Select With Selection .Font.Si2e = 14 .Font.Bold = True .Font.Italic = True

249. FormulaRlCl = "Таблица значений взаимной нормированной корреляционнойфункции"1. End With

250. With Selection.Borders(xffidgeTop) .LineStyle = xlContinuous ► .Weight = xlMedium

251. Colorlndex = xlAutomatic End With

252. With Selection.Borders(xlEdgeBottom) .LineStyle = xlContinuous .Weight = xlMedium .Colorlndex = xlAutomatic End With

253. With Selection.Borders(xlEdgeRight) .LineStyle = xlContinuous * .Weight = xlMedium

254. Colorlndex = xlAutomatic End With

255. Range("C3:C" & KtK + 3).Select With Selection.Borders(xlEdgeRight) .LineStyle = xlContinuous .Weight = xlThin .Colorlndex = xlAutomatic End With

256. Range("D3:D" & KtK + 3).Select With Selection.Borders(xlEdgeRight) .LineStyle = xlContinuous .Weight = xlThin .Colorlndex = xlAutomatic End With

257. Range("C3:E3").Select With Selection.Borders(xlEdgeBottom) .LineStyle ~ xlContinuous .Weight = xlThin .Colorlndex = xlAutomatic End With

258. Range("C" & KtK + 4 & ":" & "E" & KtK + 4).Select With Selection.Borders(xlEdgeTop) .LineStyle = xlContinuous .Weight = xlThin .Colorlndex = xlAutomatic End With1. Расчёт взаимной дисперсии1. XI =0 Х2 = 0

259. Fori = 1 То К XI = XI + ((A(i) Matl) л 2) Х2 = Х2 + ((B(i) - Mat2) л 2) Next

260. Disl =Х1/К Dis2 = Х2 / К DisV = Sqr(Disl * Dis2)

261. Range("D" & KtK + 4).FormulaRlCl = "Взаимная дисперсия процессов: " & DisV ' Расчёт значений корфункции For i = 1 То KtK XI = 0 Х2 = 0

262. Forj = 1 То (K-i+1) XI = XI + ((A(j) Matl) * (B(j + i -1) - Mat2)) X2 = X2 + ((B(j) - Mat2) * (A(j + i - 1) - Matl)) Next

263. KFP(i) = 2 * (XI / (K i + 1)) / DisV KFO(i) = 2 * (X2 / (K - i + 1)) / DisV Range("C" & i + 3).FormulaRlCl = T * (i - 1) Range("D" & i + 3).FormulaRlCl = KFP(i) Range("E" & i + 3).FonnulaRlCl = KFO(i) Next1. Graft Then

264. Построение графика взаимной прямой корфункции Charts. Add

265. ActiveChart.ChartType = xlLineMarkers

266. ActiveChart. SetSourceData Souxce:=Sheets(SName).Range("D4:D" & KtK + 3),1. PIotBy:=xlColumns

267. ActiveChart. Location Where:=xlLocationAsObject, Name:=SName

268. ActiveChart.SeriesCoIIection( 1).XValues = "=" & SName & "!R4C3:R" & KtK + 3"C3"

269. With ActiveChart .HasTitle = True

270. ChartTitle.Characters.Text = 'Трафик взимной прямой нормированной корфункции" & Mes

271. Axes(xlCategory, xlPrimary).HasTitle = True

272. Axes(xlCategory, xlPrimary). AxisTitle.Characters.Text = "Смещение, с" .Axes(xlValue, xlPrimary) .HasTitle = False End With

273. With ActiveChart.Axes(xlCategory) .HasMajorGridlines = True .HasMinorGridlines = False End With

274. With ActiveChart.Axes(xlValue) .HasMajorGridlines = True .HasMinorGridlines = False End With

275. ActiveChart.HasLegend = False ActiveChart.HasDataTable = False

276. ActiveSheet.Shapes(",nnarp. l").ScaleWidth 1.2, msoFalse,^ msoScaleFromBottomRight ActiveSheet.Shapes("flHarp. l").ScaleWidth 1.5, msoFalse, msoScaleFromTopLeft End If1. Graf2 Then

277. Построение трафика взаимной обратной корфункции Charts.Add

278. ActiveChart.ChartType = xlLineMarkers

279. ActiveChart.SetSourceData Source: =Sheets(SName).Range("E4:E" & KtK + 3),1. PlotBy:=xlColumns

280. ActiveChart.Location Where :=xlLocationAsObject, Name:=SName

281. ActiveChart. SeriesCollection( 1 ).XValues = & SName & "!R4C3:R" & KtK + 3"C3"

282. With ActiveChart .HasTitle = True

283. ChartTitle.Characters.Text = "График взимной обратной нормированной корфункции" & Mes

284. Axes(xlCategory, xlPrimary).HasTitle = True

285. Axes(xlCategory, xlPrimary). AxisTitle.Characters.Text = "Смещение, с" .Axes(xlValue, xlPrimary) .HasTitle = False End With

286. With ActiveChart. Axes(xICategory) .HasMajorGridlines = True .HasMinorGridlines = False End With

287. With ActiveChart.Axes(xlValue) .HasMajorGridlines = True .HasMinorGridlines = False End With

288. ActiveChart.HasLegend = False ActiveChart.HasDataTable = False

289. ActiveSheet.Shapes('^Harp. 2").ScaleWidth 1.2, msoFalse, msoScaleFromBottomRight

290. ActiveSheet.Shapes("fliiarp. 2"). Scale Width 1.5, msoFalse, msoScaleFromTopLeft End If1. Grafl And Graf2 Then

291. ActiveSheet.Shapes("flnarp. 2").IncrementTop 300 ActiveSheet.Shapes("flHarp. l").IncrementTop 100 End If

292. Range("A" & KtK + 20).Select End Sub

293. Программа для ЭВМ для расчёта взаимной спектральной плотности ► процесса. Программа написана на языке Visual Basic for Application для среды Microsoft Excel.

294. К Предельная частота взаимной спектральной плотности 'W - Приращение частоты

295. F Количество точек корфункции, используемых для расчётов

296. Т Частота квантования корфункции

297. Dim А(1000), В(1000), KFP(IOOO), KFO(IOOO) As Variant

298. Считывание массива значений ординат1. Com Then

299. Mesl = Range(DataSourcel) Mes2 = Range(DataSource2) Mes = " (" + Mesl + " и " + Mes2 + ")" Ml = Val(Right$(DataSource 1, Len(DataSourcel) -1)) M2 Val(Right$(DataSource2, Len(DataSource2) - 1)) Else

300. Ml = Val(Right$(DataSourcel, Len(DataSourcel) 1)) - 1 M2 = Val(Right$(DataSource2, Len(DataSource2) - 1)) - 1 End If

301. For i = 1 To KK A(i) = Range(Left$(DataSourcel, 1) & i + Ml) B(i) = Range(Left$(DataSource2, 1) & i + M2) Suml = Suml + A(i) Sum2 = Sum2 + B(i) Next

302. Matl = Suml / KK 'Математическое ожидание процессов Mat2 = Sum2 / KK

303. KT = К / W 'Количество точек спектральной плотности ' i номер точки спектральной плотности ' m - счётчик цикла суммы для корфункции If Model Then If Where = 1 Then Sheets.Add Else Workbooks.Add End If

304. SName = ActiveSheet.Name Columns("C:D").Column Width = 22

305. Columns("E").ColumnWidth = 19 Range("C 1 ").Select With Selection .Font.Size = 14 .Font.Bold = True .Font.Italic = True

306. FormulaRlCl =" Таблица значений нормированной" End With

307. Range("C2"). Select With Selection .Font.Size = 14 .Font.Bold = True .Font.Italic = True

308. FormulaRlCl = "взаимной спектральной плотности" End With

309. Range("C3").Select With Selection .Font.Size = 14 .Font.Bold = False .Font.Italic = True FormulaRlCl = Mes End With

310. With Selection.Borders(xlEdgeTop) .LineStyle = xlContinuous .Weight = xlMedium .Colorlndex = xlAutomatic End With

311. With Selection.Borders(xlEdgeBottom) .LineStyle = xlContinuous .Weight = xlMedium .Colorlndex = xlAutomatic End With

312. With Selection.Borders(xlEdgeRight) .LineStyle = xlContinuous .Weight = xlMedium .Colorlndex = xlAutomatic End With

313. Range("B4:D" & KT + 5).Select

314. Selection.Borders(xlInsideVertical).LineStyle = xlContinuous With Selection.Borders(xlEdgeRight) .LineStyle = xlContinuous .Weight = xlThin

315. Colorlndex = xIAutomatic End With

316. Range("B4:E4").Select With Selection. Borders(xlEdgeBottom) .LineStyle = xlContinuous , .Weight = xlThin

317. Colorlndex = xIAutomatic End With

318. Range("B" & KT + 6 & ":E" & KT + 6).Select Selection.Borders(xlInsideVertical).LineStyle = xlNone With Selection.Borders(xlEdgeTop) .LineStyle = xlContinuous .Weight = xlThin .Colorlndex = xIAutomatic ► End With

319. Range("C" & KT + 6).FormulaRlCl = "Использовано " & F & " точек взаимных кор-функций" End If

320. Расчёт взаимной дисперсии для расчёта взаимной корфункции XI =0 Х2 = 0

321. For i = 1 То КК XI = XI + ((A(i) Matl) Л 2) * Х2 = Х2 + ((B(i) - Mat2) л 2)1. Next

322. Disl = XI / КК Dis2 = Х2 / КК DisV = Sqr(Disl * Dis2) ' Расчёт значений корфункции For i = 1 То F XI =0 Х2 = 0

323. For j = 1 То (КК i + 1) XI = XI + ((АО) - Matl) * (B(j + i -1) - Mat2)) X2 = X2 + ((B(j) - Mat2) * (AG + i - 1) - Matl)) Next

324. KFP(i) = 2 * (XI / (KK i + 1)) / DisV KFO(i) = 2 * (X2 / (KK - i + 1)) / DisV Next

325. Расчёт значений нормированной взаимной спектральной плотности For i = 0 То КТ SumD = (KFO(l) + KFP(l)) / 2 SumM = 0 SP = 01. For M = 2 To F

326. SumD SumD + (((KFO(M) + KFP(M)) / 2) * Cos((M - 1) * W * i * T)) SumM - SumM + (((KFO(M) - KFP(M)) / 2) * Sin((M - 1) * W * i * T)) Next

327. SD = SumD * T / 3.141592654 SM SumM * T / 3.141592654 Data(i) = Sqr(SD Л 2 + SM A 2)1.Model Then

328. Range("B" & i + 5).FormulaRlCl = W * i Range("C" & i + 5).FormulaRlCl = SD Range("D" & i + 5).FormulaRlCl = SM Range("E" & i + 5).FormulaRlCl = Data(i) End If Next1. Graf And Model Then

329. Построение графика нормированной взаимной спектральной плотности Charts Add

330. ActiveChart.ChartType = xlLineMarkers

331. ActiveChart.SetSourceData Source:=Sheets(SName).Range("E5:E" & KT + 5), PlotBy:= xlColumns

332. ActiveChart.Location Where :=xlLocationAsObject, Name:=SName ActiveChart.SeriesCollection( 1).XValues = "=" & SName & "!R5C2:R" & KT + 5 & "C2" With ActiveChart .HasTitle = True

333. ChartTitle.Characters.Text = "График нормированной взаимной спектральной плотности" & Mes

334. Axes(xlCategory, xlPrimary).HasTitle = True

335. Axes(xlCategory, xlPrimary).AxisTitle.Characters.Text = "Частота, рад" .Axes(xlValue, xlPrimary).HasTitle = False End With

336. With ActiveChart.Axes(xlCategory) .HasMajorGridlines = True .HasMinorGridlines = False End With

337. With ActiveChart.Axes(xlValue) .HasMajorGridlines = True .HasMinorGridlines = False End With

338. ActiveChart.HasLegend = False ActiveChart.HasDataTable = False

339. Программа для ЭВМ для остационаривания, центрирования и приведе-» ния процесса к нормальному. Программа написана на языке Visual Basic for

340. Application для среды Microsoft Excel.

341. Sub Station(Grafl, Graf2, Graf3, Mode, Where, Com, DataSource, DataNumber, CFilter, WhereC)

342. Mode режим работы (1-по окружности, 2- потеугольнику,3-по прямоугольнику ' CFilter - ёмкость фильтра в количестве ординат(не чётное) ' WhereC - адрес первой ячейки столбца выходного диапазона Dim А(1000), АС(ЮОО), MatF(lOOO) As Variant If Com Then

343. Mes = Range(DataSource) M = Val(Right$(DataSource, Len(DataSource) -1)) Ms " (" + Mes + ")" Else

344. M = Val(Right$(DataSource, Len(DataSource) -1)) -1 Ms ="" End If

345. For i = 1 To DataNumber A(i) = Range(Left$(DataSource, 1) & i + M) Next

346. XX = DataNumber CFilter + 1 'Количество обрабатываемых ординат Select Case Where Case 1 Sheets. Add HereSold = "C" Herelold = 3 HereStrend = "D" Hereltrend = 2 HereSc = "E" Herelc = 3 Case 2

347. Workbooks.Add HereSold = "C" Herelold = 3 HereStrend = "D" Hereltrend = 2 HereSc = "E" Herelc = 3 Case 3

348. HereSold = LeftS(DataSource, 1)

349. Herelold = Val(Right$(DataSource, Len(DataSource) -1))

350. HereStrend = Left$(WhereC, 1)1. HereSc = Left$(WhereC, 1)

351. Herelc = Val(Right$(WhereC, Len(WhereC) -1))

352. Hereltrend = Herelc + XX 11. End Select

353. SName = ActiveSheet.Name If Mode = 1 Then For i = 1 To XX Mat = 0 L = -CFilter / 2 For j = i To (i + CFilter -1) Mat = Mat + (0.75 / CFilter) * (1 L / CFilter) л 2 * A(j) L = L + 1 Next

354. AC(i) = A(i) Mat MatF(i) = Mat Next End If Suml =0 Sum2 = 0 Sum3 = 0 For i = 1 To XX Suml = Suml + A(i) Sum2 = Sum2 + MatF(i) Sum3 = Sum3 + AC(i) Next

355. Matl = Suml / XX 'Матожидание исходного процесса

356. Mat2 = Sum2 / XX 'Матожидание тренда

357. Mat3 = Sum3 / XX 'Матожидание стационарного процесса

358. Расчёт дисперсий исходного и стационарного процессов1. Ml =01. М2 = 0

359. For i = 1 То XX Ml = Ml + (A(i) Matl) л 2 M2 = M2 + (AC(i) - Mat3) л 2 Next

360. Disl = Ml / XX 'Дисперсия исходного Dis2 = M2 / XX 'Дисперсия стационарного ' Создание таблиц If Where < 3 Then Range("C2").Select With Selection .Font.Size = 14 .Font.Bold = True

361. FormulaRlCl = "Приведение процесса к стационарному виду" End With

362. Range(HereSold & Herelold & ":" & HereSc & Herelc + XX l).Select With Selection .Borders(xlEdgeTop).Weight = xlMedium .Borders(xlEdgeBottom). Weight = xlMedium .Borders(xlEdgeLeft). Weight = xlMedium .Borders(xlEdgeRight). Weight = xlMedium End With

363. Range("G4:K9").Select With Selection .Borders(xlEdgeTop). Weight = xlMedium

364. Borders(xlEdgeBottom).Weight = xlMedium .Borders(xlEdgeLeft). Weight = xlMedium .Borders(xlEdgeRight).Weight = xlMedium End With

365. For i = 1 To XX If Where < 3 Then

366. Range{HereSold & Herelold + i l).FormulaRlCl = A(i) End If

367. Range(HereSc & Herelc + i l).FormulaRlCl = AC(i) - Mat3 Range(HereStrend & Hereltrend + i).FormulaRlCl = MatF(i) - Mat2 + Matl Next1. Grafl Then

368. Построение графика исходного(действительного) процесса Range("Al 5").Select Charts. Add

369. ActiveChart.ChartType = xlLine 'Markers

370. ActiveChart.SetSourceData Source:=Sheets(SName).Range(HereSold & Herelold & ":" & HereSold & (Herelold + XX -1», PlotBy:= xlColumns If Grafi Then

371. ActiveChart. SeriesCollection.NewSeries

372. ActiveChart.SeriesCollection(2).Values = "=" & SName & "!R" & Hereltrend + I & "C4:R" & Hereltrend + XX & "C4"1. End If ---

373. ActiveChart.Location Where:=xlLocationAsObject, Name:=SName With ActiveChart .HasTitle = True

374. ChartTitle.Characters.Text = "Исходный процесс" & Ms .Axes(xlCategory, xlPrimary).HasTitle = False .Axes(xl Value, xlPrimary).HasTitle = False End With

375. With ActiveChart. Axes(xlCategory) .HasMajorGridlines = True .HasMinorGridlines = False End With

376. With ActiveChart.Axes(xlValue) .HasMajorGridlines = True .HasMinorGridlines = False End With

377. ActiveChart.HasLegend = False ActiveChart.HasDataTable = False

378. ActiveSheet.Shapes("flnarp. l").ScaleWidth 2.5, msoFalse, msoScaleFromBottomRight

379. ActiveSheet.Shapes(" Диагр. 1").Scale Width 1.5, msoFalse, msoScaleFromTopLeft End If

380. Graf2 Then ' Построение графика остационаренного процесса Charts.Add

381. ActiveChart.ChartType = xlLineMarkers

382. ActiveChart.SetSourceData Source :=Sheets(SName).Range(HereSc & Herelc & & HereSc & (Herelc + XX 1», PlotBy:= xlColumns

383. ActiveChart. Location Where :=xlLocationAsObject, Name:=SName With ActiveChart .HasTitle = True

384. ChartTitle.Characters.Text = "Стационарный процесс" & Ms .Axes(xlCategory, xlPrimary).HasTitle = False .Axes(xlValue, xlPrimary).HasTitle = False End With

385. With ActiveChart.Axes(xlCategory) .HasMajorGridlines = True .HasMinorGridlines = False End With

386. With ActiveChart.Axes(xlValue) .HasMajorGridlines = True .HasMinorGridlines = False End With

387. ActiveChart.HasLegend = False ActiveChart.HasDataTable = False

388. ActiveSheet.Shapes("flHarp. 2").ScaleWidth 2.5, msoFalse, msoScaleFromBottomRight ActiveSheet.ShapesCTlHarp. 2"). Scale Width 1.5, msoFalse, msoScaleFromTopLeft End If

389. Grafl And Graf2 Then ActiveSheet.Shapes("flnarp. 2").IncrementTop ISO End If End Sub