Файл: Элементы математического моделирования в программных средах MATLAB 5 и Scilab (Андриевский Фрадков).pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 05.04.2024
Просмотров: 425
Скачиваний: 1
Рис. 3.14. Континуализация моделей.
убедиться, что такая система имеет колебательную переходную характеристику. Получим непрерывную модель этой системы. Выполним операторы:
s l = t f ( 1 . 5 , [1 0.5],Т0) csl=d2c(sl)
В результате появляется сообщение "System order was increased to handle real negative poles": "Порядок системы был увеличен, чтобы обработать отрицательные вещественные полюсы". Непрерывная система csl имеет второй порядок и описывается передаточной функцией WJs) = s*+ 2ns + 4.14^ • ,10 • Переходные характеристики исходной дискретной системы и ее непрерывной модели в моменты tk = кТ0)к = 0,1,2,... совпадают (рис. 3.14, г).
З а м е ч а н и е . Стоит обратить внимание на то, что в некоторых версиях тулбокса CONTROL SYSTEMS функция d2c для систем, заданных в НПК-форме (в виде совокупности нулей, полюсов и коэффициента передачи), работает не-
верно. |
Например, выполнив оператор s2=zpk(sl), получим |
||
ту же дискретную систему, что и исходная, W(z) |
= z |
q ^, |
|
однако для нее оператор cs2=d2c(s2) дает ошибочный |
ответ |
||
WH(s) |
= s ^ q ^. В версии 4.2 тулбокса эта ошибка |
устранена. |
93
3.5.Детерминированные и стохастические модели
Модели систем, о которых мы говорили до сих пор, были детерминированными (определенными), т.е. задание входного воздействия определяло выход системы однозначно. Однако на практике так бывает редко: описанию реальных систем обычно присуща неопределенность. Например, для статической модели неопределенность можно учесть, записывая вместо (3.1) соотношение
y(t) = F(u(t)) + <p(t), |
(3.44) |
где ip(t) - погрешность, приведенная к выходу системы. Причины неопределенности разнообразны:
-погрешности и помехи измерений входов и выходов системы (естественные погрешности);
-неточность самой модели системы, учитываемая путем искусственного введения в модель погрешности;
-неполнота информации о параметрах системы и т.д. Среди различных способов уточнения и формализации не-
определенности наибольшее распространение получил стохастический (вероятностный) подход, при котором неопределенные величины считаются случайными. Развитый понятийный и вычислительный аппараты теории вероятностей и математической статистики позволяют дать конкретные рекомендации по выбору структуры системы и оценке ее параметров. Классификация стохастических моделей систем и методов их исследования представлена в табл. 3.3.
Выводы и рекомендации основаны на эффекте усреднения: случайные отклонения результатов измерения некоторой величины от ее ожидаемого значения при суммировании взаимно уничтожаются и среднее арифметическое большого числа измерений оказывается близким к ожидаемому значению. Математические формулировки этого эффекта даются зако-
ном |
больших |
чисел |
и центральной |
предельной |
теоремой. За- |
кон |
больших |
чисел |
гласит, что |
если |
~ случайные |
величины с математическим ожиданием (среднее значение)
= а и дисперсией М(£, |
- |
а)2 |
= <72, то |
|
+ |
- |
+ |
0 |
(3.45) |
при достаточно больших N. Это говорит о принципиальной |
||||
возможности сколь угодно точной оценки |
по измерениям. |
94
Т а б л и ц а З . З . Стохастические модели систем
|
Статические |
Дискретные по Т |
Динамические |
|
||
|
|
|
Непрерывные по Т |
|||
|
Дискретные |
Непрерывные |
Дискретные |
Непрерывные |
Дискретные |
Непрерывные |
Математи- |
по С/, У |
по l/,Y |
по U,Y |
по UyY |
по U, Y |
по U,Y |
Схема незави- |
Регрессион- |
Марковские |
Стохастичес- |
Системы массо- |
Стохастичес- |
|
ческий ап- |
симых испыта- |
ные модели |
цепи; стохас- |
кие разност- |
вого обслужи- |
кие диффе- |
парат опи- |
ний [101, 102] |
[101, 102] |
тические ав- |
ные уравне- |
вания; полумар- |
ренциальные |
сания |
|
|
томаты [88] |
ния [75] |
ковские про- |
уравнения |
|
|
|
|
|
цессы [15, 94] |
[53] |
Методы |
Статистические |
Регрессион- |
Оценка пере- |
Статистическое |
Теория массо- |
Теория устой- |
оценки па- |
оценки вероят- |
ный анализ |
ходных веро- |
оценивание со- |
вого обслу- |
чивости [75] |
раметров |
ности; дисперси- |
[1, 86, 78, 110, 119] |
ятностей; ста- |
стояний и пара- |
живания; ими- |
|
и анализа |
онный анализ |
|
тистическое |
метров; анализ |
тационное мо- |
|
|
[46, 110] |
|
моделирова- |
стох астической |
делирование |
|
|
|
|
ние [15, 27] |
устойчивости [75, 86] [103] |
|
|
Методы |
Стохастическое |
Планирование |
Динамическое |
Динамическое |
Перебор; мето- |
Оптимальное |
синтеза |
программирова- |
эксперимента; |
программиро- |
программиро- |
ды оптималь- |
и адаптивное |
|
ние [80] |
стохастическое |
вание [88] |
вание [75, 95] |
ного управле- |
управление |
|
|
программиро- |
|
|
ния [103] |
[75, 95] |
|
|
вание [78, 110] |
|
|
|
|
Области |
Задачи выбора |
Обработка |
Компьютеры |
Импульсные |
Системы обслужи- |
САУ; механи- |
примене- |
из конечного |
результатов |
|
и цифровые |
вания (вычисл. |
ческие, электронные, |
ния |
числа вариантов |
измерений |
|
САУ |
системы, про- |
тепловые и другие |
|
(испытания, |
и испытаний |
|
|
изводственные |
процессы |
|
управление) |
|
|
|
системы и т.д.) |
|
Центральная предельная теорема, уточняя (3.45), утверждает, что
|
^ K i |
+ ... + Ы - в |
» ^ , |
|
(3.46) |
где |
стандартная |
= О, М£2 = |
1) нормально |
распреде- |
|
ленная случайная величина. (Если |
величина |
а2 неизвестна, |
|||
то следует заменить в (3.46) о на оценку s= yjjjzГ |
]Cili(&"~02i |
||||
где £ = |
jj'YliLi^i- При этом величина f будет |
распределена |
уже не нормально, а по закону Стьюдента с числом степней свбоды N — 1, который, к счастью, при 7V > 20 практически неотличим от нормального).
|
Функция распределения |
нормальной случайной |
величины |
|
£ |
хорошо |
известна и детально затабулирована. |
Имеются |
|
в |
книгах |
по статистике и |
таблицы функций распределения |
других часто встречающихся случайных величин, например стьюдентовой. Однако нет необходимости добывать статистические таблицы, если под рукой есть компьютер с установленной системой MATLAB. В составе тулбокса STATISTICS есть модули, вычисляющие функции распределения и другие характеристики более двадцати распространенных типов случайных величин. Узнать перечень всех функций тулбокса можно с помощью команды help tooIbox\stats или встроенной системы помощи через меню Help, окно Help Window. Проиллюстрируем применение MATLAB на простой статистической задаче.
Пусть требуется найти, при каком числе измерений N по-
грешность оценки математического ожидания |
независимых |
и одинаково распределенных величин |
окажется |
меньше чем 0.02, с вероятностью не меньше чем 0.99. Пусть среднеквадратическая погрешность каждого измерения из-
вестна: а = 0.5, а в качестве оценки берется среднее |
арифме- |
|
тическое величин |
• • • |
|
Допустим, что N > 30, так что можно пользоваться цен- |
||
тральной предельной |
теоремой. Из (3.46) следует, что иско- |
|
мое N должно удовлетворять соотношению Р |
< 0.02^ > |
|
>0.99, или эквивалентному соотношению Р j|£| < |
j > |
|
> 0.99, где £ - нормальная случайная величина со |
средним |
96
О и дисперсией 1. Вероятность, стоящую в левой части последнего соотношения, легко определить с помощью функции norminv, вычисляющей по величине р значение х) удовлетворяющее равенству Р {£ < х) = р. В силу симметрии нормального распределения Р{|f| < х} = 2р — 1. Поэтому значение р, удовлетворяющее заданному условию 2р— 1 = 0.99, равно 0.995. Выполнение оператора norminv(0.995) дает результат
ans = 2.578, откуда y/N |
> |
или N> ( 0 5 ц (ji>578)2 « |
« 4154. Поскольку 4154 > |
30, применение центральной пре- |
дельной теоремы законно (в противном случае вместо функции norminv следовало бы воспользоваться функцией tinv, обращающей функцию распределения Стьюдента). В заключение приведем полезное равенство Р{|£| < 1.96} = 0.05, позволяющее в прикидочных расчетах обходиться и без таблиц, и без компьютера. Интересно, что для решения задачи нам не понадобилась точная формулировка центральной предельной теоремы. Разумеется, формулировкам (3.45), (3.46) можно придать более строгий вид и это легко достижимо с помощью понятий вероятностной сходимости. Однако при попытке проверить условия этих строгих утверждений могут возникнуть трудности. В частности, в законе больших чисел и центральной предельной теореме требуется независимость отдельных измерений (реализаций) случайной вели-
чины и конечность ее дисперсии. Если |
эти условия |
наруша- |
||
ются, то могут нарушаться и выводы. |
|
Например, |
если все |
|
измерения совпадают: ^ = ... = |
то, |
хотя все остальные |
условия выполняются, об усреднении не может быть и речи. Другой пример: закон больших чисел несправедлив, если случайные величины распределены по закону Коши
(с плотностью распределения р(х) = 1 /тг(1 + я2)), не обладающему конечным математическим ожиданием и дисперсией. А ведь такой закон встречается в жизни! Например, по Коши распределена интегральная освещенность точек прямолинейного берега равномерно вращающимся прожектором, находящимся в море (на корабле) и включающимся в случайные моменты времени.
Но еще большие трудности вызывает проверка обоснованности самого употребления термина "случайный". Что такое случайная величина, случайное событие и т.д.? Часто говорят, что событие А случайно, если в результате экспе-
4 Б. Р. Андриевский и др. |
97 |
|
римента оно может наступить (с вероятностью р) или не наступить (с вероятностью 1— р). Все, однако, не так просто. Сама вероятность события может быть связана с результатами экспериментов лишь через частоту наступления события в некотором ряде (серии) экспериментов: vN = NA/N, где NA - число экспериментов, в которых событие наступило; N - общее число экспериментов. Если числа и^ при достаточно большом N приближаются к некоторому постоянному числу
РА- |
|
^ « Р л , |
(3.47) |
то событие А можно назвать случайным, а число р - его вероятностью. При этом частоты, наблюдавшиеся в различных сериях экспериментов, должны быть близки между собой
(это свойство называется статистической |
устойчивостью, |
или однородностью). Сказанное относится |
и к понятию слу- |
чайной величины, поскольку величина f является случайной, если случайными являются события {а < £ < 6} для любых чисел а, 6. Частоты наступления таких событий в длинных сериях экспериментов должны группироваться около некоторых постоянных значений.
Итак, для применимости стохастического подхода должны выполняться следующие требования:
1)массовость проводимых экспериментов, т.е. достаточно большое их число;
2)повторяемость условий экспериментов, оправдывающая сравнение результатов различных экспериментов;
3)статистическая устойчивость.
Стохастический подход заведомо нельзя применять к единичным экспериментам: бессмысленны выражения типа "вероятность того, что завтра будет дождь", "с вероятностью 0.8 « З е н и т » выиграет кубок" и т.п. Но даже если массовость и повторяемость экспериментов имеются, статистической устойчивости может и не быть, а проверить это - непростое дело. Известные оценки допустимого отклонения частоты от вероятности основаны на центральной предельной теореме1 или неравенстве Чебышева и требуют дополнительных гипотез о независимости или слабой зависимости измерений.
|
1 Легко |
заметить, что |
(3.47) есть частный |
случай (3.45), когда берет- |
|
ся |
= 1, |
если событие |
А наступило в i-м |
эксперименте, |
и £, = 0 - в |
противном случае. При этом Л/& = рд, |
- рА)2 = рл(1 - |
р^) . |
98