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

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

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

вновь корректировать значение правой части соотношения ( 8 )

ит.д. до установления решения.

Вкачестве примера рассмотрим интегрирование уравнения

ÊJL

-

i l l

( 12)

â t

~

â x2

 

с начальными и граничными условиями следующего вида £8]х

U(x,О) - 4х ( f - X) , 0*х$ 1,

 

 

 

( 1 9 )

U ( 0 , t ) = 0, V (t , t )= 0 ,

Схему интегрирования

(9 )

с точностью до коэффициентов за­

пишем в виде

 

 

 

Л

ÛU

д ги

д ги

âx2 r

f t 2

â t

( 1 4 )

* â t 2

с граничными и начальными условиями (13). Полагаем, что на оставшейся части границы выполняется условие

М - . л

(|3)

i t

 

которое верно при достаточно медленном изменении величины

U либо при достаточно

больших t и асимптотическом поведе­

нии решения.

 

Для уравнения (12)

с граничными условиями (13) можио

определить максимальное и минимальное значения нормы ре­ шения. Действительно, значение нормы максимальное в том случае, если начальное распределение функции lf(X, Ç ) с те­ чением времени не изменяется. Среднеарифметическое значе­ ние функции для фиксированной дискретной системы точек обо­ значим через 100%,

Значение нормы минимально, если во вс$й области инте­ грирования значение функции равно нулю. Естественно! что максимальная и минимальная оценки, нормы решения дают до­ вольно широкий интервал изменения решения. После решения •уравнения

£ J L

+

£ I . 0 .

(>в)

Фл*

*

â t *

 

при граничных условней (19), ( 15) значение норми решения будет равно 60%< Таким образом, кроме максимальной и ми­ нимальной оценок нормы решения, получена приближенная оцен­

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

лений максимальной оценки. Можно заключить, что

искомое

решение

находится между решением уравнения* ( 16)

и макси*

мальной

оценкой исходного уравнения ( 12).

 

Действительно, при проведении итерационного

процесса

(14) решение стремится в сторону максимальной оценки реше­ ния, Если в какой-либо итерации решение окажется близким максимальной оценке, то правая часть полученного решения будет близка нулю и. в последующей: итерации оно приблизит» ся к решению уравнения (1 6 ). Следовательно, искомое реше­ ние должно лежать в интервале между первым приближением! максимальной оценкой,

При практической реализации итерационной процедуры (14) в правую часть уравнения вводится такой множитель У, чтобы от итерации к итерации значение нормы решения моно* тонно возрастало на одну и ту же величину при стремлении к некоторому постоянному, значению, соответствующему иско­ мому решению. При этом значения у изменяются от 0,1 до!.

и,%

Распределение .функции U при / =0,6 с изложенным, вы­ ше управлением сходимостью решения представлено не рисун­ ке, где 1 - точное, 2 - приближенное решение; в квадратах указаны номера итерации. Вычисления, выполнены .В.Н.Курпеней на ЦВМ " £ £ - 1020"; уравнение Лапласа решено методом релаксации [ 1]\ значение нормы равно 93,5%.

Литература

1. Козлов Э.С., Сергеев Н.П*, Николаев Н,С. Автоматизация процессов решения краевых задач с помощью сеточных АЦВМ. - М.: Энергия, 1974. - 111 с.

2. Пухов Г.Е. Методы анализа и синтеза квазианалоговых электрических цепей, - Киев: Наук.думка, 1967, - 608 с.

3.Толмачев В.П. Решение краевых задач методом дифферен­ циальных приближений с использованием АВМ. - В кк: Методы и средства решения краевых задач. ' !Рига: Изд-во РПИ, 1975, с. 60 - 62.

4.Красносельский М .А., Вайникко Г.М., Забрейко П.П. и пр. Приближенное решение операторных уравнений. - М.: Наука» 1Ô69. - 456 с.

5.Коллатц Л. Функциональный анализ и вычислительная Мате­

матика. - М .: Мир, 1967. - 185 с.

6. Демидович Б.П., Марон И.А.. Шувалова Э.З. Численные ме­ тоды анализа. - М.: Наука, 1983. - 368 с.

7.Вазов В., Форсайт Дж. Разностные методы решения •диффе­ ренциальных уравнений в частных производных. - М.: Изд-во иностр.пит., 1963. - 484 с.

УДК 621.372:519.85

Е.Д.Михайлова

МАШИННАЯ РЕАЛИЗАЦИЯ ПАРАМЕТРИЧЕСКОГО МЕТОДА СИНТЕЗА НА БАЗЕ КОМПОНЕНТНЫХ УРАВНЕНИЙ

RL C -ЦЕПЕЙ

Всистемах машинного проектирования параметрический синтез может быть осуществлен путем решения компонентны'; уравнений

■fj ( X j t

X 2

* С] ~ O.t

( 1 )

связывающих коэффициенты

C j

- {

М, схем­

ных функций

 

 

 

и значения параметров X; элементов цепи с заданной тополо­ гией.

Здесь

 

 

 

 

 

 

( 3 )

где

- степень, ÇX: » Q, 1; N -

число

элементов в

цени;

It -

число сомножителей в слагаемом,

имеющих степень сг^-Ч;'

в х

-, число слагаемых, равное числу деревьев

/ -й

степени

в соответствующих графах

цепи»

j

; Лг, I

- степень one-

ратора р .

 

 

 

 

 

 

Решение системы (1 )

можно найти численными

(итерато

ными) методами оптимизации /"27 и прямым решением на ба­ зе теории исилючения [ \ ] , Обычно применяют различные ва­ рианты первого способа, которые эффективны, если известно начальное приближение '“'достаточно близкое" к решению или жестко обусловленное [ 2 ], В противном случае гарантия решениг отсутствует, что при несовместности системы ( 1) мо­ жет привести к ненужным затратам времени.

Второй путь решения не требует знания начального пррб№ женяя, -относится к точным методам и дает однозначный ответ на вопрос о существовании решения (1 ). Недостатком, этого метода является сложность формирования систем результантов на каждом щаге исключения и большой объем аналитических вычислений.

Возможность применения метода исключения к решению системы (1 ) исследована в работе f â j . В результате обосно­ вана методика формирования отдельных уравнений в системах результантов непосредственно из исходной для каждого щага исключения системы уравнений; решение систем результантов расчленено на части, в которых используется ограниченное число исходных уравнений ( j ) , соответственно ограничено •число уравнен,1й в системах результантов;

Ограничение -числа уравнений в системах результантов ^о некоторого минимального числа R * 3 позволяет довести процедуру исключения неизвестных в выбранных уравнениях К до конца. Решения совместных систем уравнений (\ ) в этом

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

(1) уравнений, значит не были учтены ограничения, вносимые неиспользованными уравнениями из системы (1 ) . Истинность или ложность решения определяется подстановкой его во все уравнения исходной системы (1 ). Такой путь более эффекти­ вен, чем решение полной системы, так как обратный ход под­ становки решений основывается на численных методах, в то время как формирование полных систем результантов - это комбинаторные .задачи с преобразованием сложных аналитиче­ ских выражений. К тому же число уравнений в полных систе­ мах результантов значительно превышает число исходных, уравнений ( 1 ), поэтому обратный ход решения полных систем даже для истинных решений (корней) связан с большим .объе­ мом вычислений.

Преимущество методики [ Ъ ] заключается в возможности

продолжать "усеченные" системы результантов до поянЬех на любом этапе, если в этом есть необходимость или позволяют возможности ЭВМ .

Основной круг задач, требующих решения при применении методики Г 37, состоит из нескольких известных задач: рас-* крытие определителя с аналитическими и числовыми элемента­ ми, определение корней многочленов и вычисление значений многочленов, а также оригинальных комбинатопных (логиче­ ских) задач по формированию результантов. Задача раскрытия определителя с произвольным алгебраическим многочленом от многих неизвестных в произвольных степенях может быть отнесена к задачам, которые ждут оптимального решения.

Принципиальное решение указанных задач не представля­ ет трудностей. При практическом применении метода возника­ ют вычислительные трудности, обусловленные большими поряд­ ками определителей и высокими степенями неизвестной в мно­ гочлене. Ограничения связаны, в основном с объемом памяти ЭВМ и временем решения, которые катастрофически возраста­ ют с увеличением числа неизвестных и увеличением максималь­ ных степеней неизвестных в исходных уравнениях.

Оценим возможные степени многочленов, к которым При-

135

В (Х$ )

водятся уравнения в системах результангов от исключаемой на данном шаге неизвестной, и порядки результантов при еле-,

дующих предположениях.

Предположи, что система (1 ) инвариантна к последова­

тельности . исключения неизвестных, хотя реально ход решения существенно зависит от порядка исключения неизвестных, и

что

на каждом

f -м этапе исключения получается многочлен

от

Jfj+j со свободным членом, т.е. слагаемым, в

котором

неизвестная

j имеет нулевую степень. Оценка

в этом

случае максимальна, Действительно, при выполнении этого предположении степень Sf многочлена от и 6 (

Ипорядок результанта ДЛИ некоторой пары уравнений

иВ2(Х() ~ Û при St=<?i2 всегда будут связаны. следующей

зависимостью:, наибольший порядок

( j j

t j

результанта

равен

удвоенному значении) максимальной степени

âf

неизвестной

в исходной системе, а

максимальная

степень

S j + j

неизве­

стной

X i + i

исключаемой на следующем этапе, равна

максимальной .степени

д{

неизвестной

 

Xj

увеличешюй в.

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

 

6 f +i

т.е.

 

 

= 2 Si ,

â i+ i « S i Cf

= 2 âf.

 

(4)

Отсюда, на

I

шаге

исключения

неизвестных степень

/-А. неизвестной,

 

 

 

определится

числом

 

 

 

 

 

А

- г 2 "

' ' ,

 

 

 

 

 

что при

/=? 1,2,3,4,5,..,

составляет

для

 

— 1,2,8,128,22768....

а для порядков определителей соответственно составит

2,4,16,266,45536,...

Числа

эти настолько

велики,

что

при

4 метод неприемлем. Однако, если

минимальная оценка осу*

ществляетсн в предположении» что свободный член в многочле­ нах относительно xf- при i > 2 отсутствует и что производит­ ся. сокращение на общие множители в виде произведений неиз­

вестных x f é x f i j 1- * ПРИ Условии Xj Ф 0 , , то на каждом шаге исключения получаются уравнения первой степени относительно X j и соответственно результанты второ­ го порядка, f i i = î, 2 , ...

Наиболее вероятной является оценка более близкая к ми­ нимальной границе, особенно для i $ 5. Как показывает опыт, степени неизвестных не превышают 4, а порядок результанта & Это объясняется наличием в исходной системе ( 1) уравнений,

не имеющих полного набора неизвестных, а также тем, что хотя каждое слагаемое в ( 1) содержит ограниченное число сомножителей п < N , но уже к третьему шагу исключения кортеж сомножителей в слагаемом представляет полный или почти полный набор неисключенных неизвестных, что на по­ следующих шагах исключения приводит к появлению общих множителей. Полагая Xj т4 0 и сокращая на общий множитель

на каждом шаге исключения,

получим уравнения значительно

более низкой степени, чем определяемое по

(5 ). Например,

если считать, что общие множители появились только на тре­

тьем шаге исключения в виде

х £ Х $ Х.£ »

то <£*=1,2,6,70,

9800, 6 j+f ** 2,4,12,140,19600.

Это объясняет

механизм умень­

шения максимальной оценки. Возможность получения общих

множителей дли системы ( 1)

становится реальной уже при ис­

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

На основании изложенного о возможных реальных значе­ ниях степеней неизвестных в уравнениях и порядке результан­ тов, учитывая, что эти значения при А/$ 10 приемлемы в смыс­ ле затрат машинного времени при использовании ЭВМ,' разра­

ботан алгоритм и программы на языке ФОРТРАН

.и АНАЛИ­

ТИК для ЭВМ "М И Р-2" [ 4 J. Алгоритм содержит

несколько

легкоформализуемых операций:

 

 

 

1. Преобразование алгебраического уравнения от несколь­

ких неизвестных в уравнение от неизвестной

X j

в виде за­

писи массивов коэффициентов при степенях

Xj

 

 

2. Определение наличия общего множителя и сокращение

на

него.

âj неизвестной

 

3.

Определение максимальной степени

Xj

в уравнениях.

 

 

 

4.

Выбор из системы ураьиений пар (сочетаний по два)

уравнений для построения результантов

порядка, 6j+j «=

=> â if +

â (2 ^

 

 

 

5.

Раскрытие результанта для каждой пары уравнений:

а) выбор элементов для члена результанта, б) определение знака члена, в) определение значения члена результанта после перемножения всех его элементов и приведения подобных чле­ нов, г ) суммирование членов результанта с приведением по­ добных.

6,

 

Переход к операции 1 для исключения следующей не.

известной*

 

 

 

 

Дйя э к о н о м и и

памяти ЭВМ и времени решения алгоритм

реализовал в виде программ, учитывающих возможные огранн.

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

тов до некоторого значения R

вместо всех возможных попар.,

ных сочетаний уравнений [Z J .

При этом

отпадает необходи­

мость в занесении в память ЭВМ всех

M>R исходных урав­

нений из

(1 ), так как они не используются. Учитывая, что а

системе

( 1)

возможны уравнения, независящие от некоторых

Xу , т.е. имеющие

неполный набор: неизвестных, из всех

уравнений ( 1) обязательно должны быть введены в новую

систему те

Q < М

уравнений, которые содержат полное мно­

жество переменных

(не менее

Q& 4 при

N > А). При слабой

связи уравнений { М &Q ) необходимо и число R принимать

не менее

0.

( R Ь й

), что гарантирует достаточность R } 3

для решения задачи. Чтобы не нарушать необходимую полноту множеств, неизвестных' в процессе исключений из-за малых К, в следующую систему результантов сначала записываются (I уравнений из предыдущей системы, которые не зависели от исключаемой на этом этапе неизвестной, и только после них R или Q ~в при R>Q~tf результантов, построенных для пар уравнений предыдущей системы.

Ход решения по такому алгоритму тесно связан с чис­ лом R результантов, которые берутся для построения урав» ний в промежуточных системах, и числом тех уравнений в системе, которые зависят от исключаемой неизвестной. Поэто му алгоритмом предусмотрен переход к исключению следую­ щей неизвестной +j при условии, что неизвестная X^ ос* тается неисключенной в последующих системах, если оказыва­ ется, что имеющихся уравнений относительно неизвестной Ц недостаточно для построения R результантов, а этих же урав

нений относительно неизвестной

Xj , /

> ï f

достаточно.

Если при любых Ху взятое

число

R

результантов недо­

статочно для проведения исключения всех неизвестных, необ­ ходимо увеличить R и повторить решение.

Если процесс исключения идет до конца и в последней системе результантов все результанты равны нулю, значит исходная система уравнений совместна, в противном случаенесовместна. При "недоопределенных" системах ( 1) получение нулевых результантов в случае совместных уравнений проис­

ходит уже в системах результантов, номер которых превыша­ ет, число заданных в (1 ) уравнений. Тогда оставшиеся нерсключенными неизвестные Xj могут принимать любые значения. Номера этих неизвестных выдаются на печать.

При совместности системы ( 1 ) осуществляется переход

ко второму этапу решения - определению корчей,

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

щих уравнения (1 ). Для этого находят все

или оддн на корней

Xjf из предпоследней системы результантов

/?/

Д л я

RlC -цепей интерес представляют в основном положительные

действительные корни. Поэтому, получив первый такой корень для КN , подставляют его в уравнения Bj (*#-/> Хц ) а 9

( ^ - 2 )-й системы, результантов, откуда'находят корень. Хи-)* Затем в уравнения В£ (Х#-2>*л/-1>х# )=й (/V- 3)-й системы ре­ зультантов подставляют значении корней Хц и Хц- t Я опре­ деляют корень Хц~2 и т.д., пока не будут найдены рее кррни. Этот этап также легко алгоритмизуем и реализован В виде программы на языке ФОРТРАН и АНАЛИТИК для ЭВМ "МИР-2*. Для поиска корней полинома использован метод Хичкока.

Программы, реализующие алгоритм на языке АНАЛИТИК и ФОРТРАН', различаются способом записи исходной информа­ ции и обработки промежуточных результантов, имеют при реа­ лизации различные временные характеристики и разные требо­ вания к объему памяти, определяемые в основном числом Неиз­ вестных и числом слагаемых в уравнениях (1 ). Г?рн раскрытии числовых результантов, нулевые значения которых определяют совместность'системы ( 1), важным является значение нормы вектора коэффициентов уравнений / / И так как при ее малых значениях, несовместность системы ( I ) может выра­ жаться ненулевыми значениями результантов очень малых зна­ чений, равных приближенно Ц 4 [ 1 1 Во избежание этого при раскрытии результантов на последнем шаге исключения нужно вводить поправочный множитель ffÿfj//*, S 9 i » - норма векто­ ра коэффициентов / -й системы результантов.

Основное время решения задачи О ) по предложенным программам затрачивается на раскрытие результантов. Это ог­ раничивает использование метода. Усовершенствование этого блока программы позволит решать уравнения с большим чис­ лом неизвестных. Поскольку наиболее высокие порядки резуль­ тантов бывают на последних шагах йсключе|шя( для раскрытия предпоследнего шага (результант имеет вид полиномиального определителя) и последнего (результант является числовым

определителем) можно в виде подпрограмм применять извест- | ные алгоритмы, которые значительно эффективнее ислользуе- , мого в программе универсального алгоритма. Однако это свя­ зано со специальным формированием информации об определи­ теле и хранением указанных подпрограмм. В предложенной же программе вопрос о формировании информации об определителе решен простым заданием степеней неизвестных во взятой паре, уравнений и записью коэффициентов при степенях хf в каждом из уравнений, что вполне достаточно для раскрытия результан­ тов, структура которых имеет конкретный вид f \ J .

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

Литература

1, Ван-дер-Варден Б.А. Современная алгебра. - М .: ОГИЗ, Щ Т . - Ч. 2, с. 8-15.

2, /1аннэ А,А,, Саркисян Б.С. Нелинейные задачи в теории реализации электронных схем. - В кн,: Пятая Всесоюз, ‘межвуз. конференция по теории и методам расчета нелиней­ ных электрических цепей и систем (2 2 -2 4 октября 1976 г.): Теэ.докл. Ташкент, 1975, с. 21-23.

3.Михайлова Е.Д. Решение систем нелинейных уравнений при определении параметров элементов электрической цепи, - Теоретическая электротехника, 1977, вып. 2 2 , с. 34-41.

4.Михайлова Е.Д., Дерябина А.Г. Решение систем нелинейны! алгебраических уравнений по методу Ван-дер-Вардена. - Республиканский фонд алгоритмов и программ. Киев, ИК АН УССР, Регистрац. № 188 от 22.11.76 г., с. 15.

УДК 621.37/39

В. В.Первое

ОПТИМАЛЬНОЕ РАСПИСАНИЕ ТЕСТОВОГО ЦИКЛИЧЕСКОГО КОНТРОЛЯ РАБОТОСПОСОБНОСТИ ЭЛЕМЕНТОВ СИСТЕМЫ

Вопрос рационального использования ресурса, в част ноет! временной избыточности, отводимого на тестовый контроль те*'; нического состояния аппаратуры системы, представляет несоМ' ненный интерес. Подобные задачи при различных условиях ис-

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