С целью устранения недостатков и сохранения ценной сущности теории проф. Н.Г.Новикова попытаемся сначала установить закономерность изменения величины отношения высот волн изгиба нитей в ткани по мере перехода от рассматриваемой фазы к последующей на основе данных Табл.1.
Анализ различных вариантов изыскания математического отображения этого закона показал целесообразность использования основного метода аналитической геометрии – метода координат. В качестве значений аргумента по оси абсцисс принимаем X= NF, а значения функции по оси ординат - Y = KhN = f (NF). С помощью координат точек можно будет алгебраически решить поставленную геометрическую задачу: определить вид функции f (NF) и форму кривой как множество точек, объединенных общим и характерным для них геометрическим свойством.
Предположим, что данная интерполирующая функция f(X) может быть представлена в виде многочлена n -ной степени по отношению к X (NF), Яковлев, [7]:
Y = f (X) = ao Xn + a1 Xn-1 + a2 Xn-2 +…+an (1)
Такой многочлен графически изображается кривой n -го порядка параболического типа.
|
|
Чтобы найти искомую функцию, необходимо определить степень многочлена n и величину его коэффициентов ai. Теоретически значения коэффициентов аi можно найти из решения системы линейных уравнений, однако на практике эта процедура является сложной вычислительной задачей при неизвестной степени многочлена n. В нашем случае, ориентируясь на Табл.1, можно заметить, что значения аргумента X образуют арифметическую прогрессию с постоянным интервалом ∆NF = 1. Это позволяет использовать случай равноотстоящих значений аргумента методом оперирования разностями функции f(X) различных порядков [7].
Для определения степени этого многочлена заполним диагоналевую таблицу разностей (Табл.2), где Xo=1; Yo=0; h =1; n = 7; DYo=Y1-Yo = 0,143; DY1=Y2-Y1 = 0,190; D2Yo= DY1- D Yo = 0,047; D3Yo=D2Y1-D2Yo = 0,030; и т.д.
Таблица 2
Таблица разностей DnY
X (NF) | Y (Kh) | Δy (∆Kh) | Δ2 y | Δ 3y | Δ4 y | Δ5 y | Δ6 y | Δ7 y |
0,143 | 0,143 | |||||||
0,333 | 0,190 | 0,047 | ||||||
0,600 | 0,267 | 0,077 | 0,030 | |||||
1,0 | 0,4 | 0,133 | 0,056 | 0,026 | ||||
1,666 | 0,666 | 0,266 | 0,133 | 0,077 | 0,051 | |||
3,0 | 1,334 | 0,668 | 0,402 | 0,269 | 0,192 | 0,141 | ||
7,0 | 4,0 | 2,666 | 1,998 | 1,596 | 1,327 | 1,135 | 0,994 |
Вычисляя последовательно разности всё более высокого порядка, к сожалению, мы не находим разности одинаковой величины, что указывало бы на окончание поиска степени многочлена. Поэтому пришлось заполнить всю таблицу для всех возможных порядков разностей DnY. В результате из Табл. 2 можно определить искомую степень нашего многочлена n = 7.
Для определения значений коэффициентов аi нашей интерполирующей функции воспользуемся формулой Ньютона [7]:
|
|
Y= Yo + (X - Xo) + (X - Xo)(X –X1) + (X –X0)(X – X1)(X – X2) +…
…+ (X -Xo)(X –X1)(X – X2)…(X – Xn-1) (2)
После подстановки численных значений рассматриваемых параметров из Табл.2 в формулу (2) находим искомый вид функции с конкретными значениями коэффициентов аi. В итоге,применительно к нашим обозначениям получаем формулу для определения коэффициента отношения высот волн изгиба нитей в однослойной ткани по мере перехода от данной фазы к последующей в соответствии с теорией проф. Н.Г.Новикова:
КhN = 0,143(NF -1) + 0,0235(NF -1)(NF-2) + 0,005(NF -1)(NF –2)(NF -3) +
+ 0,00108(NF -1)(NF -2)(NF -3)(NF -4) + 0,000425 (NF -1)(NF -2)(NF -3)(NF -4)(NF -5) + + 0,0001958(NF -1)(NF -2)(NF -3)(NF -4)(NF -5)(NF -6) +
+ 0,0001972(NF -1)(NF -2)(NF -3)(NF -4)(NF -5)(NF -6)(NF -7) (3)
Анализ громоздкой формулы (3) показывает на существенную ограниченность её использования.
Геометрическая аппроксимация функции (3) представлена на Фиг.2 в виде пунктирной кривой линии KhN, которая в качестве образца позволяет провести изыскание более простого закона Кh = f (NF).
ФИГ.2
- ТАНГЕНЦИАЛЬНЫЙ ЗАКОН РАСПРЕДЕЛЕНИЯ ПОРЯДКА ФАЗ СТРОЕНИЯ ОДНОСЛОЙНОЙ ТКАНИ 2D ПОЛОТНЯНОГО ПЕРЕПЛЕТЕНИЯ
Анализ различных видов функциональных зависимостей обнаружил очень близкий аналог геометрической интерпретации уравнения (3). Оказалось, что замечательную точку М (Кh = 1, NF = 5), соответствующую средней (пятой) фазе строения ткани на кривой КhN = f(NF), можно обнаружить и на кривой Khtg графика, апроксимирующего функцию тангенса. Известно, что тангенс 450 также равен единице в точке M, соответствующей середине диапазона изменения аргумента NF.