Файл: Элементы математического моделирования в программных средах 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