книги / Машинный анализ и моделирование электрических цепей
..pdfвновь корректировать значение правой части соотношения ( 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 многочлена от X£ и 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
В. В.Первое
ОПТИМАЛЬНОЕ РАСПИСАНИЕ ТЕСТОВОГО ЦИКЛИЧЕСКОГО КОНТРОЛЯ РАБОТОСПОСОБНОСТИ ЭЛЕМЕНТОВ СИСТЕМЫ
Вопрос рационального использования ресурса, в част ноет! временной избыточности, отводимого на тестовый контроль те*'; нического состояния аппаратуры системы, представляет несоМ' ненный интерес. Подобные задачи при различных условиях ис-