книги из ГПНТБ / Измайлов А.Ф. Аналоговые вычислительные машины и их применение в ВМФ учебное пособие с элементами программированного обучения
.pdfxg(0)^px(0)+alfx(0) i- a2x(0) = ^ 2 0 + |
аг х ао . (3 . 7 ) |
Существует бесчисленное множество вариантов преоб разования одного уравнения в систему. Каждому варианту соответствует определенный набор на АВМ. Может слу читься, что переход от уравнения к системе приведет к оптимальной схеме набора, что даст возможность получить решение с наименьоиыи погрешностями. Так, например, пе реход от уравнения (3 .1 ) к системе (З .б ) значительно умень&ает количество сигналов, подаваемых на входы суммирущих блоков (сравните схемы рис. 3.7 и З .П ) .
Однако при решении на АВМ систем уравнений общего вида в схеме набора могут появиться контуры, склонные к самовозбуждению (рис. 3 .7 ), что вызовет либо мгно венную перегрузку операционных блоков, либо усилит вли яние различных помех. Чтобы этого не случилось, любая система уравнений должна быть в ходе подготовки пре образована, по крайней мере, к виду (3 .8 ).
138
Дальнейвее исключение неизвестных из левой части уравнений (3 .8 ), например преобразование к виду ( З .б ), сохраняет указанное свойство приведенной системы. С методикой подобных преобразований можно познакомиться при выполнении задания 13.
Переход от уравнения к системе связан с затратами ручного труда на вычисление ее коэффициентов и началь ных условий. Поэтому о целесообразности перехода судят, сопоставляя предполагаемое повышение точности с объе мом дополнительных вычислений. Возможности рассмотрен ной оптимизации изучены мало.
Задание 13. Систему уравнений
(p2+2p-h3)x1+(f+3p+2)oc2+ ( f + 2p+l)*:3=fi ;
(p 2f f 4 ) x i +(2p2+2р+3)хг+(3р*+р+2)хэ=/г ; .
(р г+Зр Ч )х+ (зр + р Н )х г +(2 р 2+3р +3)x=f3
-
преобразуйте к виду (9 ) . После изучения методов набора уравнений ( § 5 гл. Ш) составьте блок-схемы, отвечаю щие этим двум системам.
139
Решение,
(р+2р+3)х + (р+3р+2)эс2 +(p2t 2 f i--l')x3=fll 7
(-р + ^ х /(-р г+рЧ)ос+(- 2р2+р~1)х3 ~f, ~ /2 ;
( р +Ь ) х Г З х г + ( ~ 3 p ^ 3 p ) x = f r 2 f 2+ f3 .
Характер проделанных преобразований нетрудно уста
новить по виду правых частей уравнений.
§ 3 . Замена аргумента
Предварительно выясним закономерности преобразова
ния функций |
cc(t) |
и ее производных р х ( { ) f p zx ( t ) , , . . 7 |
|
в результате |
замены |
в них аргумента t |
пропорциональ |
ным аргументом *Z |
: |
|
|
|
t = |
rn*"t - |
|
|
|
(3 .9 ) |
||
|
|
|
|
|
|
|
|
|
|
Коэффициенты пропорциональности mt = co n st |
при |
||||||||
нято называть масштабом аргумента. |
|
|
|
||||||
Произведя указанную систему, получаем функцию |
|||||||||
v(m t Т) |
. |
Если теперь аргументом |
считать |
не произведе |
|||||
ние m t (T ) , |
а только |
сомножитель |
Т |
, приходим к |
|||||
функции X (?) . |
|
|
|
|
|
|
|
||
Пример, |
x(t) =0,(Ht +i |
|
|
|
|
||||
|
|
|
|
|
|
||||
|
|
|
m t =2. |
|
|
|
|
|
|
|
x(™i '?)=0,C H (2Z)4 ; Х (7 ) = 0 ,0 8 Т 3Н . |
||||||||
Связь |
между функциями |
x ( t ) |
и Х (Т ) |
такова: |
|||||
X ( t ) = X ( ? ) |
при условии |
t = nit Z |
, |
(3 .10) |
|||||
|
|
|
|
|
|
|
|
|
|
что иллюстрируется графиком (рис. 3 .3 ). |
|
|
|||||||
Дифференцируя |
выражение (1 0 ), |
определяем соотношение |
|||||||
между производными: |
|
|
|
|
|
|
|||
d x ( t) |
olх (mf Т ) |
_ |
/ |
с(Х(Z ) |
|
t = m t Z j |
|
||
d t |
cl(mt 7) |
|
m t |
d Z |
’ |
|
|||
|
|
|
|
( З . И )
IM
**h
Рис, 303
а в общем случае |
(рис. 3 .4 ; |
|
|
|
|
|
|
|
|
|
(3 .12) |
Продолжая пример, можно проверить соотношения |
|||||
(3 .1 0 )-(3 .1 2 ). Эти формулы, |
разумеется, справедливы |
||||
для известной функции - f ( t ) |
и всех ее производных. |
||||
Переход от заданного |
( I ) к машинному уравнению |
||||
(2) частично заключается |
в |
замене функций x ( t) ,p x (t) ... |
|||
f ( t ) переменными |
Х (Л ) |
, |
р Х ( 7 ) , . . . , |
Г * (7 ) |
соглас |
но формулам ( I I ) |
и (12). |
В результате |
замены |
получаем |
или |
|
|
|
|
(3 .13) |
|
|
|
|
|
|
|
|
|
|
р х ю |
+ |
|
|
|
|
|
(3 .14) |
Начальные условия |
преобразуются |
аналогично |
при |
||
t =T = 0 |
. Соотношения ( I I ) и (12) |
выполняются безу |
|||
словно, |
отсюда |
|
|
|
|
|
X (О) - эс(О) |
рХ(О) = mt р х (О) р |
|
||
р |
X( 0) ~mt f |
x ( 0 ) f ... } р *X(0) = m “ *px(0) |
|||
Уравнение (14) |
является машинным и после введения |
||||
обозначений |
|
|
|
|
143
M l
Рис. 3.4
г |
а г |
- - |
А |
|
— " А , ; |
|
2 Т ’ ' 7 |
||
« о |
|
|
||
|
|
|
|
п ч а п . |
л |
- |
, а а - |
пЧ |
|||
/77, |
|
|
|
а„
II ib
7 |
п |
F (7)-Г(Т) |
|
m r , ОГ„ |
|||
|
полностью соответствует форме записи (2 ):
Р * Х *А , f " ' ' x * • • +АЯ., f X + A „ X = Г .
Замена аргумента влияет на коэффициенты уравнения, на диапазон изменения переменный, на частоту и длитель ность, и величина масштаба rnt определяется из стрем ления к оптимизации перечисленных характеристик машин ного уравнения и процесса, развивающегося в модели*
Было бы желательным добиться равенства
(3 .15)
i а I у 2 ,3 , • • m/7
и одинаковых диапазонов изменения переменных
l * L KC Ip ^ I mokc I f 3 ^ Iм а к с ~ |
|
|
^2 |
||||
Каждому коэффициенту |
уравнения |
Ai |
соответствует |
||||
произведение |
коэффициентов передачи |
блоков |
( |
) |
|||
(см. гл . |
Ш, |
§ 4 ) , и условие(3.15)может |
быть реализова |
||||
но путем |
установки |
= I . В этом |
случае |
отпадает не |
обходимость введения потенциометров в схему, а блоки моделируют математические операции с погреоностями, близкими к минимуму. Если все переменные изменяются в одном и том же диапазоне, то пределы выходных напряже ний всех интегрирующих блоков можно выбрать максималь ными ( - 100 в ) . Этим будет достигнуто наилучиее сог
1 0 |
145 |
ласование пределов шкалы вольтметра и измеряемого на пряжения и уменьшения влияния "дрейфа нуля".
Условия (15) и (16) - лишь ориентир для подготовки уравнения, вероятность строгого их выполнения ничтожна, и практика вполне допускает значительные отступления от них. Безусловно, обязательными являются условия, вы текающие из допустимых значений коэффициентов передачи блоков, частот и длительностей. Практика показывает, что реально достижимому выравниванию коэффициентов уравнения или переменных сопутствуют вполне приемлемые частоты и длительности процесса - модели.
Для правильного определения величины масштаба нуж но располагать некоторой информацией о характеристиках процесса-оригинала: диапазонах изменения параметров, скоростях, ускорениях, частотах, времени затухания пе реходных режимов и т .д . Такую информацию могут дать анализ физической картины исследуемого процесса, сопо ставление с аналогичными процессами и, наконец, проб ное решение на АВМ.
Буравнениях и в системах уравнений, встречающихся
винженерной практике, часто наблюдаются тенденции к монотонному изменению коэффициентов и пределов измене ния производных, т .е . тенденции к одновременному выпол
нению неравенств (3.17) и (3.18) или (3.19) и (3 .2 0 ): (3.17)
макс
(3.19)
146
|
V |
|
|
J |
(3.20) |
Несмотря на весьма приближенный характер этих не |
||
равенств, выбор масштаба аргумента на их основе |
обычно |
|
дает положительный результат. |
|
|
Для того, чтобы коэффициенты |
стали примерно |
одинаковыми, целесообразно добиться равенства коэффи циентов в крайних членах левой части уравнения (3 .1 3 ):
откуда следует формула для расчета величины масштаба:
(3.21)
Имея целью в первую очередь добиться выравнивания диапазонов переменных, производят расчет по формулам (3 .2 2 ), которые вытекают из соотношений (3 .1 0 ),(3 .1 2 ), и в качестве масштаба их принимают наименьший резуль тат :
> |
(3 .22) |
|
1V7