книги / Микропроцессорное управление технологическими процессами в радиоэлектронике
..pdfуправления технологическими процессами намотки и навивки. Наибольшими возможностями для повыше ния быстродействия и точности обладают частотные методы измерения. Несмотря на большое количество литературы, посвященной анализу этих методов, до сих пор не решены вопросы оценки в динамических режимах точности воспроизведения непрерывных диф ференцируемых функций натяжения упругого мате риала.
Рассмотрим |
простейшую |
классическую |
упругую |
систему — струну. Поперечные колебания |
такой си |
||
стемы описываются следующими уравнениями: |
|||
неподвижная |
струна |
|
|
|
д2у |
д2у |
(2. 1) |
|
dt2 |
дх2 |
|
|
|
||
неподвижная струна с учетом диссипативных |
|||
потерь |
|
|
|
|
= |
|
(2.2) |
движущаяся струна
движущаяся в упругой среде струна
где у — поперечное смещение струны; х — координа та точки в продольном движении струны; т — масса
струны; а = |
К 7 7 т — скорость распространения бе |
||||
гущих. волн; |
Т — натяжение |
струны; |
2р = R/m — |
||
коэффициент |
сопротивления |
среды; |
v — линейная |
||
скорость движения упругого |
материала. |
|
|||
Параметры |
поперечных |
колебаний |
упругого |
||
материала |
определяются коэффициентами |
v, а, р. |
21
Натяжение Т может изменяться как под воздей ствием случайных величин, так и параметров движе ния упругого материала (ускорения, частоты, ампли туды продольных колебаний, неточности изготовле ния механических узлов технологической системы намотки), и в общем случае описывается некоторой непрерывной функцией Т (0. которую можно разло жить в ряд Фурье,
1 |
00 |
(ak cos /гоу + Ьк sin kco0t) == |
Т (0 = -о- aQ+ |
Е |
|
z |
fe=i |
|
|
= |
£ с * 1* * '. |
|
|
k = —00 |
Для управления |
движением упругого материала |
и его параметрами необходимо воспроизвести эту функцию в виде частоты поперечных колебаний упру гого материала. Частота колебаний и параметр Т в простейшем случае, описываемом уравнением (2 1)» связаны зависимостью, которую и используют в час тотном методе измерения параметра Т в вибропреобра зователях
|
|
|
h ^ - h - V T /m , |
|
(2.5) |
где |
I — длина струны. |
|
|
||
Однако это уравнение применимо для статических |
|||||
(Т = |
const) |
или |
динамических (Т = |
var) режимов |
|
при |
условии |
f0 |
fn, где fn = л сор/2л — я-я |
гармо |
|
ника функции Т (0, при которой ап |
amln, Ъп |
bmwt |
0min и frmin — любые, наперед заданные минимальные амплитуды i-й гармоники функции Т (t), величины ко торых выбираются из заданной точности воспроизве дения функции Т ({).
Для динамических режимов, при которых частоты /о и fn соизмеримы (но /0 > fn)> частотный метод, бази-
22
руюшийся на уравнении (2.5), неприменим вследствие низкой точности. В этом случае необходимо иметь информацию о поперечных колебаниях упругого мате риала, т. е. знать меру достоверности информации о предельной величине отношения f0/fnt при которой от сутствует потеря информации, а также длительность минимального информативного временного интервала. Поскольку частота колебаний является переменной при Т — = Т (/), то можно измерить лишь ее усредненную величи ну (Лф) за определенный ин тервал времени (4, 4+0» на пример за период (рис. 1).
Этот интервал времени опре деляет усредненную частоту поперечных колебаний
/ср = l/(4+i — 40*
а следовательно, усредненное значение функции Т (t) на этом интервале. Очевидно, что чем меньше этот интер вал, тем точнее определяется
(воспроизводится, восстанавливается) функция Т (t). Временные интервалы (4, 4-и) можно формировать за полный периохц_тогда У Т = fcp (функция иг), или за
полпериода К Г = 2/ср (функция и2). Эти интервалы соответствуют нулям функции у (t) (поперечных коле баний упругого материала). Для увеличения точности определения функции Г (0 интервал (4, 4+i) можно сократить за счет формирования его в экстремальных точках функции у (0 и ее производных, т.. е. в нулях производных функции у (0 (функция «,). Из рис. 1 видно, как уменьшается временной интервал (4, 4+0 и, как следствие, увеличивается точность определения
23
функции Т (i). Например, при использовании эк стремальных точек самой функции y (f), т. е. при значе
нии у' (0 = |
0 ]Л7, = |
4/ср. Отношение частоты следова |
||
ния импульсов /ср = |
l/(/*+i — tk) |
к частоте |
изменения |
|
/„ функции |
Т (/) увеличится в |
4 раза |
и, следова |
тельно, во столько же раз увеличится точность вос произведения функции Т (0- Однако вопрос о выборе минимального интервала, соответствующего неиска женной информации о частоте поперечных колебаний упругого материала, так же, как и вопрос о достовер ности получаемой информации в теории колебаний, не решен. Исследования частотного метода, основан ного на применении формулы (2.5), не содержат оцен ки точности воспроизведения исходной непрерывной функции Т {t) в динамических режимах (при соизме римых частотах fn и /ср). Это объясняется отсутствием
критерия и методов |
воспроизведения функции Т (О |
в динамических режимах. |
|
При разработке |
обобщенного критерия точности |
методов воспроизведения функции Т (f) возникают следующие задачи: общая задача об информативности поперечных колебаний упругой системы, которая сводится к определению информативности функции
колебаний |
у (0 и ее производных yin) (/), где п = |
= 1, 2, 3, |
...» и исследованию точности воспроизведе |
ния Т {f) в зависимости от отношения /ср//„; частные задачи 1) определения минимального информативно го временного интервала функции колебаний у (/), величина которого однозначно связана с параметром колебаний Т; 2) повышения точности воспроизведения функции Т (/) при соизмеримых fn и /ср за счет допол нительной информации колебаний упругой системы, возникающей в динамических режимах.
Для решения этих задач необходимо: доказать информативность интервалов времени (tk, 4+i), обра зованных нулями функции поперечных колебаний
24
материала у (О, и информативность (зависимость от параметра Т) интервалов времени (4, 4 + 0 , образо ванных нулями производных функции у (f) для всех упругих колебательных систем, описываемыхуравне ниями (2.1) — (2.4); определить погрешность восстанов ления функции Т (t) и минимальный информативный интервал (4, 4-j-i)- Эти задачи являются актуальными для исследования движения не только упругого ма териала, но и обширного класса виброчастотных преобразователей, работающих в динамических ре жимах.
2.2. МАТЕМАТИЧЕСКОЕ ОБОСНОВАНИЕ МЕТОДА ПРЕДЕЛЬНЫХ ИНТЕРВАЛОВ
Простейшая колебательная система. Свободные по перечные колебания ограниченной струны с закреп ленными концами описываются уравнением
___ г |
* > 0 |
(2.6) |
“ 0 дх2 |
|
|
и краевыми условиями. |
|
|
У(О» /) = */(/, 0 = 0 . |
|
(2.7) |
Считая, что при t <z 0 натяжение постоянно и стру на колеблется с постоянной частотой, можно задать начальные условия в следующем виде:
|
у(х, |
0) = A sin ах; |
|
(2 .8) |
|
|
|
|
|
|
|
|
• | г |
<*> |
°) = °- |
|
|
Учитывая |
краевые |
условия |
(2.7), |
находим |
|
A sin аI = 0, |
откуда |
а |
= nkll (k = |
± 1 , |
± 2 , ...). |
Уравнение (2.6) с краевыми условиями (2.7) и началь ными условиями (2.8) решим методом разделения пе ременных. Найдем решение уравнения (2.6) в виде
25
у (x, t) — X (*) Ф (0- |
После |
подстановки значения |
|||
у (.х, f) в уравнение (2.6) определяем |
|
||||
X” {х) |
|
Ф" (0 |
1 |
(2.9) |
|
Х{Х) |
~ |
Ф(<) <2а (О |
|||
|
|||||
где X — некоторая |
константа. |
|
|
Из выражения (2.9) и краевых условий (2.7) полу чаем следующую краевую задачу для определения
Х(х):
X" (JC) + WX {х) = 0, |
X (0) = X (/) = |
0. |
(2.10) |
|||||
Общее решение уравнения (2.10) имеет вид |
|
|||||||
X (х) = |
Сг sin Хх 4- С2 cos Хх. |
|
|
|||||
Решая уравнение в точках х = 0 и х = |
/, записываем |
|||||||
X (0) = |
С2 = 0, |
X (/) = |
Cj sin X/ = |
0. |
|
|||
Поэтому Я, = |
шг//, |
п — ± 1 , |
± 2 , ± |
«•• • |
Из вы |
|||
ражения (2.9) для |
каждого |
п |
получаем |
следующее |
||||
уравнение: |
|
|
|
|
|
|
|
|
Ф"п (0 + |
(n n /lf а* (/) Фп (0 = 0. |
|
(2.11) |
Обозначим через Ф„ (/) (i = 1, 2) решение уравне ния (2.11), удовлетворяющее следующим краевым условиям:
dt
= 1 ; Фп( 0 ) |
II о |
|
с |
¥ |
___1 |
FII1с |
а |
1 О |
о |
|
|
(2. 12)
Тогда |
решение |
начально-краевой задачи (2.6) — |
(2.8) имеет вид |
|
|
у (х, |
о = Е |
(АпФ'п(<) + В Л (/)) sin (япхП), |
где Ап, Вп — неизвестные постоянные.
26
Учитывая начальные условия (2.8) и условия (2.12), получаем следующие равенства:
-рСО |
|
|
|
|
|
|
|
2 |
|
Апsin (nnx/l) = A sin {лЫ 1)\ |
|||||
П——оо |
|
|
|
|
|
||
|
|
+ ° о |
|
|
|
|
|
|
|
2 |
Bnsin(ntix/l) = 0. |
|
|
||
|
|
—оо |
|
|
|
|
|
Отсюда и |
из полноты |
функций |
(sin (:nnxlt)} на |
||||
интервале 0 ^ |
х ^ |
/ делаем вывод, |
что Ап — 0 при |
||||
п Ф k, Вп = |
0 для всех л и |
= А . Поэтому решение |
|||||
начально-краевой |
задачи |
(2.6) — (2.8) |
запишем в |
||||
следующем |
виде: |
|
|
|
|
|
|
|
у {х%t) = |
АФ\ (t) sin (nkxll), |
|
||||
где Ф1 (0 — решение |
уравнения (2.11) |
при п — k и |
|||||
краевых условиях |
(2.12). |
|
|
|
|||
Определим натяжение струны с помощью частоты |
|||||||
ее колебаний при k = |
1 |
|
|
|
|||
|
у (*, /).= |
ЛФ[ (/) sin [пхИ). |
|
Частота колебаний струны совпадает с частотой функции и (t) = Ф| (t). Функция и (t) удовлетворяет следующему уравнению:
u"(0 + (JL)3 Q2(< )u (i)= o. |
(2.13) |
Докажем несколько теорем, необходимых для ана |
|
лиза функции и (t). |
|
Теорема 1. Об информативности нулей |
функции |
и (t) (или у (0 — поперечных колебаний |
упругого |
материала). Пусть t0 и tx два последовательных нуля функции u(t). Положим / = 1/[2 (^ — /0)]. ‘ Тогда на интервале t0 < t с tx существует точка ? такая, что / = а
27
Доказательство. Рассмотрим .функцию |
z (t) — |
||||
= sin 2л/ (t |
— /0). Эта функция |
обращается |
в нуль |
||
в точках t0 и /i и удовлетворяет уравнению |
|
||||
|
|
z" (/) + 4л2/2г ( 0 |
= 0 . |
|
(2.14) |
Теорему докажем методом от противного. Пред |
|||||
положим сначала, что при всех t £ [t0, /х] / > |
a (0/(2/). |
||||
Тогда при |
всех t £ U0, / J выполняется |
неравенство |
|||
(л//)2 а2 (t) |
< |
4л2/2. Применим |
теорему |
Штурма к |
уравнениям (2.13) и (2.14). Между любыми последо
вательными нулями |
функции |
и (/) на |
интервале |
|||
U0, |
/ J |
должен находиться нуль |
функции |
г (/). |
Но |
|
г (/) |
на |
интервале [/0, |
/ J нулей |
не имеет, |
тогда |
как |
/0 и tx — последовательные нули |
функции |
и (/). |
По |
лученное противоречие показывает, что наше пред
положение / > |
a (t)l(2J) при всех |
t £ [/0, |
tx] |
неверно. |
|||
Покажем, что при всех / £ [(0, |
не может выпол |
||||||
няться |
и |
неравенство / < |
а (0/(2/). Действительно, |
||||
тогда 4л2/2 < |
а2 (0 (л/02 и |
по теореме Штурма пуля |
|||||
ми /0 и 0 |
функции г (0 должен |
быть нуль функции |
|||||
и (0, |
что |
противоречит |
предположению |
теоремы. |
|||
Поскольку |
функция а (0 |
неперерывна, |
и |
функция |
£а (0
/----- щ - принимает как положительные, так и от
рицательные значения, то по теореме Вейерштрасса
существует точка |
в которой |
/ = |
a {t')/2l, что и |
||||||
требовалось |
доказать. |
|
|
|
|
|
|
||
Теорема 2. Об информативности нулей первой про |
|||||||||
изводной |
функции и (/). Пусть |
/0 — нуль |
функции |
||||||
и (/), а |
наименьшее |
такое |
число, что |
и' (tj) |
= 0. |
||||
Положим / = |
1/[4 (/х — /„)], |
тогда на |
интервале |
/0 < |
|||||
< t < . tx существует точка Г, в которой |
справедливо |
||||||||
равенство / |
= а (/")/(2/). |
|
|
|
|
|
|
||
Доказательство. |
Поскольку |
/х > |
/0 — наимень |
||||||
ший корень функции и' (/) на полуоси t > |
t0, то при |
||||||||
t £ U0> tx] функция и (/) не изменяет знак. |
|
|
28
Рассмотрим функцию г (?) = sin 2nf (t — ?0). Тог да г" (?) -f 4л2/2z (?) = 0.
Умножив это равенство на и (?) и вычтя из него равенство (2.13), умноженное на г (?), получим
z" (?) и (?) -|- 4л2/2г (?) и (?) — z (?) и" (?) —
— л2//2а2 (?) z (?) и (?) = 0.
Откуда
-А . (г' (t) u(t) — z (t) и’ (t)) = г" (0 и (t) — z (t) и" (<) =
= 4я2 (/2 — а2 (<)/(4/2)) «(*) г (<).
Интегрируя это равенство в пределах от ?0 до ?lf записываем
г' (У и (У — г (У и' (У — г' (У и (У +
+ г (У и' (У = |
4я“ J (Р _ |
д 2 (t)/(iP)) и (t) z (0 dt. |
||||
|
|
h |
|
|
|
|
Так как |
и (?0) |
= и' (?х) |
= 2 (?0) |
= |
г' |
(?х) = 0, то |
} |
( Р - о ^ ) /( 4 Л » «< *)*№ |
я |
= |
0. |
||
to |
|
|
|
|
|
Из неотрицательности функции z (?) на интервале [?0, ?J, знакопостоянства функции « (?) на этом интер вале и полученного выше равенства находим, что функция f — а (?)/(21) не может быть на интервале [?0, tj] знакопостоянной. Поэтому по теореме Вейер-
штрасса |
существует |
такое ?" £ [?0, |
?J, |
что / = |
= а (п |
/ т . |
|
|
|
Теорема 3. О неинформативности нулей второй |
||||
производной функции |
и (?). Каждый |
нуль |
функции |
и (?) является нулем функции и" (?), и каждый нуль функции и" (?) является нулем функции и (?).
29
Доказательство. Пусть и (У = 0. Подставив в уравнение (2.13) вместо Означение /0, получим а" (У = = 0. Пусть и" ( у = 0. Подставляя в уравнение (2.13) значение t = У находим
(я//)2 a2(*i) ы(У = 0.
Но а (У ^ 0, поэтому w (У = 0. Следовательно, теорема доказана.
Докажем теперь, что нули функции и(п) при п > 2 неинформативны. Покажем вначале, что нули функ
ций |
и(п) |
неинформативны при /г > |
2 и |
7' = const. |
|||||||
Теорема 4. О неинформативности нулей функций |
|||||||||||
и{п)• |
Пусть |
Т (0 = |
const. |
Тогда |
нули |
функции |
|||||
M(2W(& = |
0, |
1, |
...) |
совпадают |
с |
нулями |
функции |
||||
и (f)у а нули функции u{2k+l) (k = |
0, |
1, 2..,) |
с нулями |
||||||||
функции |
« '(/). |
|
|
|
|
|
|
|
|
||
Доказательство. Дифференцируя п раз равенство |
|||||||||||
U" (0 |
+ с&и (0 |
= 0, |
получаем |
и(п+2) (f) + |
alu{n) (t) = |
||||||
= 0, |
где |
|
с& — — |
Таким |
образом, |
имеем, что |
|||||
и{п+2) (t) = |
|
—do и(п> {t) для |
любого |
/г |
0. |
Отсюда |
|||||
определяем |
|
|
|
|
|
|
|
|
|
||
|
м<2*> (t) = |
- а2ои^~2) (t) = |
|
aW 2k~ A) (t) = |
|||||||
|
|
|
|
= . . . |
= ( _ 1 ) V |
0*« (0 ; |
|
|
|||
|
a <“ +,) (if) = — |
|
(0 |
= a40a (2ft_3) (0 |
= |
= . . . = ( - l ) kabku'(t).
Поскольку a0 Ф 0, то из этих равенств видно, что
нули и(п) (0 при л > 2 совпадают с нулями и (t) или и' (f) в зависимости от четности я, и поэтому дополни тельной информации о Т (t) по сравнению с нулями и (t) или и' (t) не несут, если Т (t) является посто янным.
30