Главная Цепи и сигналы [0] [1] [2] [3] [4] [5] [6] [7] [8] [9] [10] [11] [12] [13] [14] [15] [16] [17] [18] [19] [20] [21] [22] [23] [24] [25] [26] [27] [28] [29] [30] [31] [32] [33] [34] [35] [36] [37] [38] [39] [40] [41] [42] [43] [44] [45] [46] [47] [48] [49] [50] [51] [52] [53] [54] [55] [56] [57] [58] [59] [60] [61] [62] [63] [64] [65] [66] [67] [68] [69] [70] [71] [72] [73] [74] [75] [76] [77] [78] [79] [80] [81] [82] [83] [84] [85] [86] [87] [ 88 ] [89] [90] [91] [92] [93] [94] [95] [96] [97] [98] [99] [100] [101] [102] [103] [104] [105] [106] [107] [108] [109] [110] [111] [112] [113] [114] [115] [116] [117] [118] [119] [120] [121] [122] [123] [124] [125] [126] [127] [128] [129] [130] [131] [132] [133] [134] [135] [136] [137] [138] [139] [140] [141] [142] [143] [144] [145] [146] [147] [148] [149] [150] [151] [152] [153] [154] [155] [156] [157] [158] [159] [160] [161] [162] [163] [164] [165] [166] [167] [168] [169] nrv-v-i e(t) -O-r- Интегратор Интегратор Рис. 8.51. Линейный четырехполюсник Рнс. 8.52. Структурная схема моделирования уравнения (8.109) уравнения системы результаты аналогового моделирования нетрудно затем преобразовать в алгоритм моделирования на цифровой ЭВМ. Моделирование дифференциального уравнения проиллюстрируем сначала на примере линейной цепи, показанной на рис. 8.51. Напряжение (в операторной форме) на выходе цепи \1Ср Lp4-r+\/Cp /.Cp2-f-rC/)-f 1 Перепишем это выражение в форме Е (р) = LCp и (р) + rCpU (р) -j-U(p). (8.107) Учитывая, что умножение U (р) на р эквивалентно дифференцированию и (t), а умножение на р эквивалентно двухкратному дифференцированию, приходим к следующему дифференциальному уравнению: e{t)LC d + rC + uit)LCu + rCu-u. dt dt (8.108) Решим это уравнение относительно наивысшей производной: ЬСи - е - гСи- (8.109) Уравнению (8.109) соответствует схема на рис. 8.52. Основными эле" ментами этой схемы являются интеграторы на операционных усилителях, сумматоры и умножители на число. Выходы интеграторов ими характеризуют состояние системы в рассматриваемый момент времени / с учетом начальных условий, т. е. состояния в момент времени о. а также с учетом действия входного сигнала на интервале от /„ до t. В этом смысле выходы интеграторов и и и называют переменными состояния, а основанный на них способ интегрирования дифференциальных уравнений называют методом переменных состояв и я (МПС). По существу, представленное на рис. 8.52 устройство, основанное на применении интеграторов, решает заданное дифференциальное уравнение (8.108) и определяет реакцию цепи (см. рис. 8.51) на воздействие е (t). Принцип построения схемы, описанный выше, можно применить и для моделирования нелинейной системы. В качестве иллюстрации построим 1 Этот метод получил особенно широкое распространение при анализе линейных систем со многими входами и выходами, когда эффективно применение матричного исчисления и векторного анализа для операций с большим числом переменных, входящих в исследуемые задачи [8]. В данном случае нас будет интересовать применение МПС к нелинейным задачам. (См. Деруссо П., Рой Р., Клоуз Ч. Пространство состояний в теории управления: Пер. с англ. - М.: Наука, 1970. схему, соответствующую нелинейному дифференциальному уравнению (8.81), которое запишем в форме (7 = - cos u>t --2aq~u)b q + bgq-. (8.110) Схема представлена на рис. 8.53. Дополнительно к двум интеграторам и сумматору потребовался функциональный блок, осуществляющий операцию возведения в куб функции q (t). Определение реакции нелинейных цепей на заданное воздействие решением дифференциального уравнения требует, как правило, больших вычислительных затрат. Задача упрощается, если рассматриваемая нелинейная система допускает дискретное представление процессов, происходящих в отдельных ее элементах, как линейных, так и нелинейных. Такое положение имеет место, в частности, если нелинейные элементы являются безынерционными. В подобных случаях открывается путь к эффективному цифровому моделированию. Покажем это на примере моделирования простого амплитудного детектора (см. рис. 8.24). Вольт-амперную характеристику диода представим в виде нелинейной функции д=/("д) = /(-"вых). (8.111) где "д = е - - напряжение, действующее на диоде. Соотношение между токами в линейной части схемы определяется выражениями 1«(0 + с(0=гд(0. hit)- откуда вытекает следующее нелинейное уравнение: ~ + "вых it)=f\e it) - «еых (01 Ul д С"яых=/(е-"ьых) - (8.112) Это уравнение моделируется аналоговой структурной схемой, представленной на рис. 8.54. На выходе сумматора 1 образуется разность е - «вых> которая после нелинейного преобразования в диоде дает ток /д, т. е. первое слагаемое в правой части (8.112). Нелинейное безынерционное преобразование (8.111) легко реализуется в ЭВМ. Функция / (Ujj) может быть задана либо в виде таблицы (например, по экспериментальным данным), либо аналитически. В первом случае таблица должна храниться в памяти ЭВМ, во втором случае вычисление значений 1д производится по программе, соответствующей аналитическому выражению. Переходя к дискретному моделированию, задаем шаг Т исходя из наивысшей частоты в спектре входного сигнала е (t), руководствуясь правилами, изложенными в § 2.17. Таким образом, входящие в уравнение (8.112) величины е (t). Рис. 8.53. Структурная схема моделирования уравнения (8.112) "вых (О и 1д (О ДОЛЖНЫ быть заме-нены соответственно, на е (тТ), / - "вых (тТ) и (д (тГ). Для перехода от дифференци- ального уравнения (8.112) к эквивалентному ему разностному уравнению воспользуемся соответствием dusuxit) «вых («Л-«вых ("Т -л Интеаратор (8.113) Рис. 8.54. Структурная схема моделирования уравнения (8.112) которое имеет смысл при условии достаточной малости Т. Тогда уравнение (8.112) принимает вид -r ["вых (тТ)-иАтТ-Т)\ = /[е(тТ)-ив,(тТ)\- --г- «вых imT) "вых("1Т) = Rf\e{tnT)~UB,{tnT)] + вых(тГ-7), (8.114) Z+T где т = RC\ а„ = Т/ (х + Т); = х/ (х + Т). Уравнению (8.114) соответствует схема, представленная на рис. 8.55, а, а общая аналоговая схема (см. рис. 8.54) после перехода к дискретной обработке принимает вид, показанный на рис. 8.55, б. Из сопоставления аналоговой и дискретной схем следует,что звено с задержкой Т в кольце обратной связи для разностного уравнения аналогично интегратору для дифференциального уравнения. Характеристики и свойства подобных устройств обсуждаются в § 12.8. Моделирование рассматриваемой цепи целесообразно в тех случаях, когда аналитическое решение затруднительно или даже невозможно. Подобная ситуация имеет место, например, при детектировании относительно коротких радиоимпульсов, когда в пределах длительности импульса укладывается всего лишь несколько периодов высокочастотного заполнения и условие (8.56) оказывается невыполнимым (см. § 8.9). Пусть, например, средняя частота заполнения /о - 1/Го = 465 кГц (Промежуточная частота в обычно.м радиоприемнике), а длительность им- е(тТ)-иьь\-(тТ) иъыхСгпР) Рис. 8.55. Дискретный вариант схемы, представленной на рис. 8.54 1 Постоянная времени т= i?C, как правило, значительно больше периода высокочастотного заполнения входного сигнала е (i), так что спектр выходного напряжения «вых (О много же спектра сигнала е (t). Следовательно, при значении Г, обеспечивающем сохранение информации, которая содержится как в е (t), так и в / [е (t)], наверняка выполняется (8.113). [0] [1] [2] [3] [4] [5] [6] [7] [8] [9] [10] [11] [12] [13] [14] [15] [16] [17] [18] [19] [20] [21] [22] [23] [24] [25] [26] [27] [28] [29] [30] [31] [32] [33] [34] [35] [36] [37] [38] [39] [40] [41] [42] [43] [44] [45] [46] [47] [48] [49] [50] [51] [52] [53] [54] [55] [56] [57] [58] [59] [60] [61] [62] [63] [64] [65] [66] [67] [68] [69] [70] [71] [72] [73] [74] [75] [76] [77] [78] [79] [80] [81] [82] [83] [84] [85] [86] [87] [ 88 ] [89] [90] [91] [92] [93] [94] [95] [96] [97] [98] [99] [100] [101] [102] [103] [104] [105] [106] [107] [108] [109] [110] [111] [112] [113] [114] [115] [116] [117] [118] [119] [120] [121] [122] [123] [124] [125] [126] [127] [128] [129] [130] [131] [132] [133] [134] [135] [136] [137] [138] [139] [140] [141] [142] [143] [144] [145] [146] [147] [148] [149] [150] [151] [152] [153] [154] [155] [156] [157] [158] [159] [160] [161] [162] [163] [164] [165] [166] [167] [168] [169] 0.0014 |