1 / 23

2 раздел. Математические методы обработки результатов параметрических исследований

Базовый курс «Основы автоматизированного проектирования» Семестр : 9 Название части курса: «Параметрический анализ систем». 2 раздел. Математические методы обработки результатов параметрических исследований. 2.3. Статистический анализ.

boris
Download Presentation

2 раздел. Математические методы обработки результатов параметрических исследований

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Базовый курс «Основы автоматизированного проектирования»Семестр : 9Название части курса: «Параметрический анализ систем».

  2. 2 раздел. Математические методы обработки результатов параметрических исследований

  3. 2.3. Статистический анализ Первый класс задач параметрического анализа в статистической постановке: анализ в условиях параметрической неопределенности. Задача № 1. Исследование влияния внешних параметров : на вход системы поступает стационарный центрированный случайный процесс x(t), имеющий спектральную плотность Sx(ω). Требуется определить статистические характеристики выходного процесса (для задач автоматического регулирования — регулируемой величины).

  4. Этап № 1. Формирование случайного процесса с заданной спектральной плотностью Sx(ω). Для этой цели используется формирующий фильтр (см. формулу (8) из «Статистические методы в параметрическом анализе»), на вход которого подается белый шум с характеристиками (7) (та же ссылка). На выходе фильтра получается желаемый процесс со спектральной плотностью (6) (та же ссылка). Для практического моделирования на ЭВМ используется цифровой формирующий фильтр. Sвых(z)=|Φ(z) |2 Sвх(z)

  5. Задание № 1 Сформировать случайный процесс с автокорреляционной функцией в виде экспоненты, т.е. имеющий спектральную плотность Gx(f)=4a/(a2+4π2f2) Очевидные преобразования : Gx(f)=4a/(a2+ω2)=(2√(a))2/(a2-(jω)2)=[2√(a)/(a-jω)] [2√(a)/(a+jω)=Φ(jω)Φ(-jω) Отсюда передаточная функция аналогового фильтра: Φ(s)=2√(a)/(a+s), где s — переменная Лапласа.

  6. Алгоритм • Сформировать аналоговый фильтр • Найти дискретный аналог формирующего фильтра • Извлечь информацию о полиномах числителя и знаменателя • Сгенерировать входной процесс в виде нормального дискретного белого шума со следующими параметрами: среднее – 0, дисперсия σ2=1, среднеквадратичное отклонение – σ=1. • Входной процесс пропустить через формирующий цифровой фильтр Параметры a=10 и Ts=0.01 Создать программный код- см.следующий слайд

  7. Программный код a=10; Fils=tf([2*sqrt(a)],[1 a]) Ts=0.01; Filz=c2d(Fils,Ts) % дискретный аналог формирующего фильтра [ndcell,ddcell]=tfdata(Filz) nd=cell2mat(ndcell) dd=cell2mat(ddcell) bzf=nd; azf =dd X0 = randn(1,1000); X = filter(bzf,azf, X0) Отчет: созданный в виде m-файла программный код с комментариями к операторам

  8. Этап № 2. Определение статистических характеристик процесса. Часть 1. Определение статистических характеристик полученного с помощью формирующего фильтра процесса X, который является входным для основной системы. Задание № 2 Определить характеристики процесса X с помощью метода Берга с функцией обращения [Pxx,f] = pburg(x,p,nfft,fs) [Pxx,w] = pburg(x,p,nfft) pburg(x,p) Создать программный код- см.следующий слайд

  9. Программный код fs=1/Ts; pburg(X,1,[],fs) % график 1 [Pxx,f1]=pburg(X,1,[],fs) figure pburg(X,1) [Pxx,w] = pburg(X,1) Отчет: созданный в виде m-файла программный код с комментариями к операторам и графики спектральной плотности

  10. Задание № 3 Выполнить расчет спектральной плотности с помощью пакета Simulink. Создать модель, соответствующую рисунку следующего слайда.

  11. Настройки параметров: Для формирования случайного процесса используются два блока: - блок входного белого шума White Noise; - блок цифрового фильтра All-pole Filter с заданной передаточной функцией Параметры блока приведены на следующем рисунке Отчет: График сравнительного анализа расчетов спектральной плотности с помощью блоков пакета Simulink и выводы о точности методов

  12. Задание № 4 Провести расчет характеристик случайных процессов: среднее (Mean), дисперсия (Variance), стандартное отклонение (Standard Deviation) , а также построение гистограммы (Histogram). Использовать пакет Simulink Создать модель, соответствующую рисунку следующего слайда. Параметры Подать сумму единичного сигнала (блок DSP Constant) и белого шума со средним 0 и дисперсией 0.05.

  13. Отчет: Показатели дисплеев с характеристиками, график с векторного осциллографа - данные моделирования суммарного процесса, график гистограммы по команде hist(yout)

  14. Задание № 5 Провести статистические исследования влияния внешних параметров на основную модель model1. Использовать отладочные модели.модифицированную модель model1 сохранить под именем model1stat1.mdl – см. следующий слайд Моделирование произвести при параметрах шума noisemean=0.2; noisevariance=0.1; Отчет: Показатели дисплеев с характеристиками, график моделирования процесса

  15. Оценка времени переходного процесса осуществляется с помощью видоизмененного блока Конец Задачи № 1.

  16. Задача № 2 Исследованиевлияния внутренних параметров : дана непрерывная система управления со случайными (неопределенными) параметрами. Требуется определить статистические характеристики выходного процесса (для задач автоматического регулирования — регулируемой величины). Задание № 6 Исследовать влияние вариаций первой постоянной времени объекта управления, субмодель которого представлена теперь в общем виде на следующем слайде: Базовые параметры модели не изменились: Tpf1=0.3; Tpf2=0.1; Tpf3=0.03;

  17. Задание № 7 Создать программный код с именем moddanstat2.m для исследования влияния случайной составляющей в первой постоянной времени Tpf1 на характеристики переходного процесса. Отчет: график значений Tpf1 со случайной составляющей и времени переходного процесса tpp для 20 экспериментов Текст кода приведен на следующем слайде

  18. noisemean=0; noisevariance=0; Tpf1=0.3; Tpf2=0.1; Tpf3=0.03; Tpf1b=Tpf1; Nrand=20; Tpf1var=Tpf1*0.3; Tpf1rand=0+sqrt(Tpf1var)*rand(Nrand,1); for i=1:Nrand Tpf1=Tpf1b+Tpf1rand(i) sim('model1stat2',tmod) [n,m]=size(tpptek1); tppstat(i)=tpptek1(n,1); end; tppmean=mean(tppstat) tppstd=std(tppstat) plot(Tpf1b+Tpf1rand,'r') hold on plot(tppstat,'g') hold off R = corrcoef(Tpf1b+Tpf1rand,tppstat)

  19. Задача № 3. Решение первого класса задач в случае задания интервала неопределенности. Задание № 8 Задачу исследования критерия качества системы устойчивость свести к исследованию ограниченного числа полиномов наихудшего вида Система имеет характеристическое уравнение s3+a2s2+a1s+a0=0 границы 4≤a0≤20; 10≤a1≤150 12≤a2≤30

  20. Программный код для исследования наихудших полиномов alfm=[4 10 12] betm=[20 150 30] pol1=[1 alfm(3) betm(2) betm(1)] pol2=[1 betm(3) alfm(2) alfm(1)] pol3=[1 betm(3) betm(2) alfm(1)] pol4=[1 alfm(3) alfm(2) betm(1)] kor1=roots(pol1) kor2=roots(pol2) kor3=roots(pol3) kor4=roots(pol4) Отчет: значения корней для 4 наихудших полиномов и их анализ на устойчивость.

More Related