Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Машинный анализ и моделирование электрических цепей

..pdf
Скачиваний:
4
Добавлен:
12.11.2023
Размер:
12.84 Mб
Скачать

увеличении быстродействия АЦВК (без использования суперпо­ зиции) по сравнении с ЭЦВМ, включенной в состав аналогоцифрового комплекса, поэтому кратность применения принципа суперпозиции должна быть 3 -4 и уж во всяком случае'не пре­ вышать 10.

Проведены исследования, целью которых было выяснение зависимости точности решения задачи от кратности примене­ ния принципа суперпозиции, а также влияния порядна подклю­ чения кодоуправляемых элементов к узлам сетки (ломодульно, в шахматном порядке, через ряд й т.п.). Исследования показа­ ли, что подключение кодоуправляемых источников тока или на­ пряжения (через КС) к узлам АП можно производить в. любой последовательности (это не влияет на точность решения). Что касается кратности применения принципа суперпозиции, то с ее увеличением растет ошибка, определяемая приборной по­ грешностью измерительного устройства;

ànp *■ яДлр '

где ц - * кратность применения принципа суперпозиции; Ар* - приборная ошибка.

Поскольку применение принципа суперпозиции позволяет более точно произвести конечно-разностную аппроксимацию исходного уравнения (благодаря использованию более мелкой пространственной сетки), проведено сравнение величины âpp с ошибкой конечно-разностной аппроксимации. Сравнивались результаты решения двумерной нестационарной задачи тепло­ проводности на АП с двукратным использованием принципа су­ перпозиции с решением той же задачи без применения супер­

позиции на АП, имеющем вдвое меньшее число узлов. Установ­ лено, что ошибки, возникающие в результате использования принципа суперпозиции, оказываются на порядок .меньше по­ грешности пространственной аппроксимации. Таким образом, основным фактором, ограничивающим кратность применения принципа суперпозиции, является снижение быстродействия, и ебли ограничить кратность указанным выше пределом (3 )- (4 ), то использование, суперпозиции оказывается целесообразной.

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

икодоуправляемым источником тока, включенного между узлом

и"землей", показана на рис. 1. Между тем представляет инте­

рес использование принципа суперпозиции в АП с кодоуправляе-

Рис, 2,

мыми источниками напряжения (К Н ), подключаемыми к уэлак

.сетки через кодоуправляемые проводимости (К П ), во-первых, потому что в этом АП исключен кодоуправляемый источник тока, который более сложен, чем КН, а во-вторых, потому что принцип суперпозиции реализуется на АЦВК [ ЗУ, включа» шем указанный выше АП. Упрощая схему кодоуправляемых элементов, мы усложняем схему коммутатора (вдвое увеличя вается количество двухпозиционных ключей). Это объясняете» тем, что при использовании принципа суперпозиции источники тока (или напряжения) заменяются эквивалентами их внутрев

них сопротивлений.

Если

в случае КТ, внутренняя проводимое

которого близка

нулю,

исключение из цепи КТ равносильно

разрыву в цепи, то отключение КН, внутреннее сопротивление которого равно нулю (в случае идеального источника), озна­ чает короткое замыкание в цепи, т.е, КП долмсна быть подки чека второй клеммой к 'зе м л е ' (первая клемма соединена с узлом сетки). Переход от обычной схемы замещения (схема Либмана) к схеме с использованием принципа суперпозиции для обоих видов АП показан на рис. 2 (на примере одномер* ной задачи нестационарной теплопроводности). Надо заметит^ что поскольку указанные значения внутреннего 'сопротивлении имеют идеальные источники, для уменьшения ошибок, возни­ кающих в результате неучета их действительных значений, необходимо номиналы 'временных' сопротивлений в АП с ко* неуправляемыми источниками напряжений делать достаточно большими (при Ду500 Ом и внутренним сопротивлением КН

~ 1 Ом погрешность составит всего 0,2%), а при испольэ(

вашш КТ сопротивления R должны быть намного меньше внутреннего сопротивления кодоуправляемого ис­ точника тока.

Таким образом, при использо­ вании принципа суперпозиции комму­ татор, обслуживающий АЦВК с КН и КП в аналоговом процессоре, дол­ жен иметь структуру, обеспечиваю-; щую подключение КН к одной из

двух кодоуправляемых проводимостей (в создаваемом АЦВК предусмотрено двукратное применение принципа су­ перпозиции) с одновременным соеди­ нением второй КП с ‘'землёй1 (рис.З). В качестве ключей коммутатора ис­ пользуются двухпозиционные реле с. тремя контактными группами типа РЭС-32 (одна контактная группа ра­

ботает на переключение КН, другая - на подключение КП к шине "земля", третья служит для блокировки обмотки реле).

В заключение отметим, что применение метода суперпо­ зиции, сокращающее количество кодоуправляемых источников тока (напряжения), особенно оказывается ^эффективным в ана­ логовых процессорах с большим количеством узлов, поскольку в этом случае уменьшение габаритов АЦВК и его стоимости особенно существенно.

Литература

1, Николаев И.С., Козлов Э.С., Максимов М.М, Назначение и принципы построения аналого-цифрового вычислительного комплекса "Сатурн". - В кн.: Средства аналоговой и анало­ го-цифровой вычислительной техники. - М,: Машинострое­ ние, 1068, с. 71-79.

2, А.с, 497603 (С С С Р), Устройство для моделирования задач теории поля, ' Ю,М,Мацевитый, В.А.Маляренко, С,Ф,Лущпен-

ко и др.-Опубл. в Б,И., 1975, № 48.

3, À.c. 491963 (С С С Р), Устройство для моделирования тепло­ проводности, Ю.М.Мацевитый, В,А,Маляренко, В,А;Пал.еЙ и др.-Опубл, в Б?И,,1975| № 2,

УДК 881.323:518.5

О.Н.Пьявченко

РАЗНОСТНО-КВАНТОВАННЫЕ СХЕМЫ ИНТЕГРИРОВАНН НА Ш А С БАЗОВОЙ СИСТЕМЫ ДИФФЕРЕНЦИАЛЬНЫХ

УРАВНЕНИЙ

Для решения на цифровых интегро-арифметических стр] турах (ЦИАС) [ \ ] системы обыкновенных дифференциальны уравнений представляются в базовой форме

 

Если

( f - v a f o ,

то у ^ М 0,уп£ ,- -,уггУ(м;г'',у11)1

«

•А —2,31 /

р —0, f t ••• t в » &

 

;

ÿp«■ const

(р*а+1,Си-2,,..,Ь);

I

/ # - '7

W

*

J>/

Базовая ситема ЦИАС отличается от решаемой на цифровых интегрирующих структурах системы К.Шеннона [ 2 ] наличиек, в её составе суперпозиций алгебраических и трансцендентш функций. При этом включению подлежат функции, вычислений значений которых методом интегрирования порождающих ура нений либо невозможно из-за их гипертрансцендентного харе тера, либо нецелесообразно по каким-то причинам, наприме; из-за высокой алгоритмической сложности, узкого допустим го диапазона изменения переменных и т.п. Для интегрировав ния дифференциальных уравнений базовой системы (1 ) испол эуются формулы на основе значений переменных f S j , При этом вычисления в ЦИАС строятся по иным раэностно-квай ванным схемам, чем соответствующие схемы интегрировани) уравнений Шеннона. Отсюда вытекает необходимость оёосно ванной разработки интерполяционных и экстраполяционных р ностно-квантованных схем моделирования решения базовой системы (1 ), пригодных для построения наборов макроопера ций и алгоритмов реализации в Ш А С .

Для решения этой задачи применительно к ЦИАС, числ в которых представляются с фиксированной запятой, сделав)

,ряд допущений, Будем считать, ч.то для облегчения подбора

' supIjjpWU /

svp 1ц (x)l* J;

 

-

™ Р ! Ир (*>Ы I

sup I Ц ( х ) Ы 1 ;

(2 )

_

P ~ 0t f f •••* П j ft 1"? г ’—t to f & — 1,21 ... ftt ,

 

Запятая фиксируется перед первым'цифровым разрядом непо­ средственно после разряда, в котором записан знак числа. Ко­ личество разрядов, отводимых для представления хранимых в

памяти интегро-арифметических модулей

(М ИА)

ЦИАС цифро­

вых значений переменных у^ (X) ( $=1,2,3,...,и)ъ

констант' у^

(It ар=(Г1-1,а.+2г ..} Ъ ) , ограничим числом N . Обозначим

квантованные N -разрядные значения величин

 

(t

= 0,

m-f).

у^ (х )

 

Ур

Квантование значений переменных

и констант

производится оператором

 

 

 

*Р-КЕУPdH-çrf ~ ?41ЕHfid+l-ç)+Of?Я

(3)

 

(it = 1,2,.,.,6 ;

<

 

 

 

 

 

 

позволяющим с

помощью функции

расчленения

P-v [

2] огра­

ничить количество разрядов

числа

y p ( i t j - ç )

до разряда -

включительно. Если

квантуется приращение

^ Ук (i+tŸ^ftkd+th

~ H i > то

? H k c iH P ^ V - N ^ h d t n ^

где оператор

4>.м в

общем случае отличается от оператора

<р-р

 

 

 

Исходя из изложенного, введем квантование чисел и при­

ращений в сцстему интерполяционных разностных уравнений,

моделирующую решение базовой системы (1 )

 

 

 

Ук(Ш ) “ Jfki

ф ^Уk(i t?) »

 

 

 

 

 

 

V.SkiiH) =

лЛ} Vsn u i) ,

/ e

 

 

 

 

 

V 9j ( i + i ) ~ ÿ p r f t t )

 

 

 

 

*

 

 

’‘(Ураtt-fl)~ypd-fi))(9a(£+t-у) ~ÿ afi-v) ^ ’

$к}е f °>Ц‘

I Если

Е ц = 0

( j = fa r) ,

то

 

 

 

 

 

Уш и) = ^к (9o’9j(éH)tÿ2(iti)>-"fÿadd),y(at!),’",^to

/( —293,

р ff /

P •’'О f î / • • • t Я ft fk tt

^ Л

7о 9

У р -û on st

. S r ’ -

/ г * ;

UP (ftto^ m Урв *

i , . r f i t

,

yp(i+f-jg>

В уравнениях (4 ) приращение интеграла V i j d t ! ) выщ

сляется по формуле численного интегрирования по Стилгьесу построенной на основе значений переменных

H tH iH -'O (ft - > т; , а значения функций (|

 

 

 

/ а ~.Рк ^9p>9ll92>-:>9a> !/(а+{)>

'• 9 Р

 

рассчитываются с

помощью операторов

(

4 ),

реализую­

щих методы интерполирования, разложения

в ряды

и т.п.

. При ограничении разрядной сетки величин в разностной

системе (4 ) неквангованные значения

9k (i+ i) > 9p fi+ t -ft)'•

9й (6 * 1 - У ? »

^fiY

заменяются квантованными

fl

-разрядны­

ми 9 k ( i H ) f ‘9P (t+ / - fi)

• f f ( P H ~ v ) •

Âfiv

 

В

этом

случае

приращения интеграла

V s j r f + f j рассчитываются

не

по раз­

ностной

(4 ), а по разностно-квантованной формуле

 

 

дг-/ т-#

^sj ( i t î )

9p(i+1) ( 9 f ( t + l > ~ '9 f i ^

 

 

 

 

 

 

 

 

 

 

 

Л

 

 

; (f y a 't i - v ) " 9f a - y p -

В результате

такой замены значение интеграла

^ j ( i t i ) вы

числяется в М ИА

с погрешностью квантования второго рода

 

 

L

 

 

 

i

 

 

 

 

 

 

 

 

9 p (e+ i)(°9 ti(e+ v ~°9$е>*е^0 ü9p(e+f??9f(e+i> *

*Л Л Л

 

^Ур(е +*~<Р} F

 

 

 

 

 

 

*ЛЛЛо ***{0i P ^

f l - j t r 09 p (C 'fiP ?9 f (e t i-

 

 

 

J

t

T

АР у Fy p fe * t -/ »(0 9 tK e H -i» ~ 09 f (e - y P '

(61

e*o/ > *fiX

 

В этом выражении погрешностями квантования первого рода

являются остатки

0ÿp ( e H -fi>,

^

отбрас#

 

 

ол ру

 

ваемще при квантовании чисел ^

, 9 t ( e + l - 0 > * f i *

(fi e 0,1, ...,т ; V

, яг) :

 

 

^9p(e+l-'fi) 9p(e+t-fi) ~9p(e-fi) •

" 09$ie+1-v> ~ 9q.(e*1-y)~ 9 f(e ~ V •.

 

L OAfly

**Ajgy ~^py’

 

 

Для orip ^деления .характера изменения

и оценки

величиныпо­

грешности квантования

второго рода

/3 i ( i t j )

( е ) найдем

выражение ее

среднеквадратичного значения б ; • Если счи­

тать ординаты

Цр(е . п

(г*~/,2,...,а, е ‘= -ю ;-n t,1 ,.„ rW rP

независимыми величинами, то погрешности

O f& (e+f? ( е=х

- -

т

, - т

i ? при округлении по 1 /2 могут принимать

любые одинаково возможные значения от

до +1/2

(f t

-

основание

системы счисления). При этих условиях сред­

неквадратичная погрешность квантования первого рода состав­ ляет

G

Г

*

(7 )

2]f3

 

 

 

 

Среднеквадратичная погрешность квантования второго рода на основании выражения (в ) записывается в виде

2

Уа V 9 i(* t »69(S+r>

«wW ■■■—

~

Т'

с~0А=0'*жи

г

______г

 

ш—шшш

■— ШЁШШШттшяш——

 

^

я-i m -fi-î

 

 

 

 

у ^

 

 

6 p ci+ t-jn

 

 

ê

 

 

 

 

 

^

h < e + ï'V f i

Gp (e - p ) + ( У У р а + I - f l f i G9 (i+ 1 "*> *

* ( ^Ура-ру) Gçc-У) * 2^ ( P Sp(e+1- p fi G if,(9-7?J

Отбрасывая в подкоренном выражении величины четвертого, порядка малости и полагая G p e ^ t - G (7 ) ( е=-1Я+1,1. - , О,/..., i +1 ) , перепишем последнее выражение в виде

 

i

-t

 

ai * °>2 Ш -~f iyi

> lt + h J ,

t

*

{Очевидно, что

Qj <0,2Ш~ЫJ / 2sup Цр(Х) t h 2(ù-J) sgPftfp W)JZ+

+ У (ь + 1 ) sup [tfç (x>]Z+/t zf 2 sup [ y pl (X )]* t

m-J n-p-1

t ! stip[ ÿ '( X > ] !} 2 2 A .,

* r Jj=û ?*>0 *

В соответствии с условиями (2 ) это неравенство npeoôpaayei ся к виду

/т-Ут-/>-1 г

Giг,- *0,407R

у 1 *Ь!( i.tZ 2 2

А0

)

_ .

 

 

А -0 * - •

 

 

 

т-1 Я'Р-?

Ал» < /

 

 

 

Так как

2

2

Г 3 7, то

окончательно

получаем

р*0

У“0

J

 

 

 

 

 

Gj < 0,407 R ~Nÿ l + h zi

(8)

Из этого соотношения следует, что погрешность квантования второго рода (6 ) имеет тенденцию к накоплению. Однако ско рость накопления довольно низкая. Например, условие 6i<lOR соблюдается! при ингегрировании с шагом f}= 0,2 и погрешно­ стью d -1СГ4 на протяжении /< 15000 шагов.

Недостатком разностно-квантованной формулы (5 ) я в л я й ся ее значительная трудоемкость из-за наличия операций ум­

ножения . 20 -разрядных чисел

на О -разрядные и суммирова­

ния J0 -разрядных чисел. Для оценки сложности, благодаря

свойствам коэффициентов AfiV : Аоа =-1/2 ;

AfiV = -Ayfi (fit t>D

и f i * у ) f Ару~0 (fi= v > О)

t формула (5 )

переписывается

в более удобной форме

 

 

(ÿpd+1) + îp i

^ *

*4* ( (ip (*+ t rfi) 'f p (i - f i) n fy d *h fi)

Подсчет доказывает, что сложность этой формулы составляет

Q - ( 9т - 7 2С t ( * т - 3 ) У ( f 4 / n ^ 3 )

операций суммирования С и умножения У Д/ -разрядных чисел,

88

Сложность вычислений существенно уменьшается при ин­ тегрировании раэностно-квантованной формулы

V3j(iff> ~ М2 ( ÿprffV +Уp i > ( / " * » f y i) +

-V-Jt f ( §p(itf-fi) - Уp ci -p) Mÿa (itf -p ) "fa.( i ~р)Я] '( 11)

вкоторой используются произведения только Jlf -разрядных

сомножителей и суммируются слагаемые, имеющие не более 2N разрядов. Сложность формулы (11) составляет всего

Q = (7т - 5)С 1 - (Зт~ 2)У (Г*яг*3).

(12)

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

 

d я-р-1

d X-A-I

(13)

G j< 0,407R

(22 Z л‘+л*)+г+22 2

4

 

 

р -о р*р+1 хр

ржеp *p ti

W

 

Это обстоятельство служит серьезным аргументом против практического применения формулы (11).

Если в формуле (9 ) постоянные коэффициенты Apj± ввести под квадратные скобки и проквантовать произведения этих ко­ эффициентов на разности ординат, то получим следующую раз­ ностно-квантованную формулу:

^ sj ( i + f ? ~ ^

(ÿ p (ù tt) +У'pi X ÿ ÿ rfti? ~УÇ* ) +

 

*pmû

t

( F p (iH -p )~ ÿ p fi-p p jc / ^ tl-p ÿ ^ 'y i^

~ Ÿ -ii £ Afifl (ifP(£+7-/i) ~У p c i - ( HtfrCit!~ p ) ~if f c i -p ) V *

( 1 4 )

Среднеквадратичная погрешность квантования второго рода

этой формулы удовлетворяет условию

 

ffj<ü,407R~M]/f -с h2mi

(IS)

а сложность .составляет

 

 

Q *

(pff, - f i) e t (4m - 3)7

(16)

Сравнение фррмул (9 ), (14) показывает, что формула (9 ) не, много точнее формулы (14), но содержит на (т - 1) больше операций суммирования N -разрядных чисел.

Однако вычислительную сложность формулы (9 ) можно снизить до уровня (16), С этой целью оценим возможность сокращения разрядности 3N -разрядных слагаемых до М*2Н разрядов. Количество разрядов М должно быть таким, чтобы вводимая погрешность квантования

G j{ М ) = 0,208R ~Mf in T

( 171

не сказывалась на точности вычисления значения интеграла,

Т 6

G j { M ) ~ K G j ( К

« 1 )

 

(18)

Подставляя в это уравнение выражения

(1

7 ),.(8 ),

находим,

что погрешность вычисления интеграла

S jd + ig практически

не изменяется, если после округления слагаемых оставлять

М

> j ÿ

j [~1д 8 -0 ,1 5 +Ш д8 +0,51дт - О,SIg (f*/r Zi ) ] (19)

раэряаов.

Для

практических оценок воспользуемся упрошенной

формулой

 

М - —

(0,05+ R lg R + 0 ,5 lg m j ) .

 

( 20)

 

 

 

Расчет показывает,

что при R -

10, /V=6, f f l- 3, /=1Сг

количест­

во разрядов М можно взять равным 10, т.е. меньше

2 N . Ест

же выбрать

М = 2N

$ то погрешность G j ( M )

удовлетворяет

условию (18)

при

# 4 2* 1 о ! Т а к и м образом, на практике

можно выбрать величину

М в пределах

 

 

 

 

 

N < М * 2 t i .

 

 

 

При этом сложность разностно-квантованной формулы

 

^ s] ( i +1? ~ ¥-м

[ 2

U fp d + V fÿpi K f y t i + I ) ~Fgi

 

 

 

 

 

 

 

 

( 21)

~ ( И

р

 

<

*

+J-J3) ~ ÿ g (c -jt p ]}

 

составит не более

0, (16) операций суммирования и умноже­

ния N -разрядных чисел. Отсюда следует, что алгоритмы опе­

раций интегрирования МИА ЦИАС можно строить

на основе

9.0

Соседние файлы в папке книги