Идентификация объектов управления. Семенов А.Д - 138 стр.

UptoLike

ВУЗ: 

ПГУ | Пенза

Составители: 

Рубрика: 

uc=u-mu;
m=3;%
Задание размерности АРСС - модели
% Задание начальных условий в РМНК
P=1000*eye(2*m,2*m);
Q=zeros(2*m,1);
F=Q;
%
Метод РМНК
for i=1:n-m
F=[-yc(i+m-1:-1:i);uc(i+m-1:-1:i)];
ch=P*F;
zn=1+F'*P*F;
gm=ch/zn;
P=(eye(2*m)-gm*F')*P;
Q=Q+gm*(yc(m+i)-F'*Q);
kf(i,1:2*m)=Q'; %
Коэффициенты АРСС - модели
Tp(i)=F'*Q;
End
%
Ошибка идентификации и ее характеристики
e=yc(m+1:end)-Tp';
de=std(e);
plot(t(1:n-m),kf),grid
pause
sr=[yc(m+1:end),Tp'];
plot(sr),grid
pause
plot (e),grid
%
Вычисление передаточной функции объекта и его характеристик
nun=kf(end,m+1:end);
den=[1 kf(end,1:m)];
wod=tf(nun,den,dt) %
Вычисление дискретной передаточной функции объекта
won=d2c(wod) %
Вычисление передаточной функции объекта
sm=ss(won)%
Модель объекта в пространстве состояний
impulse(won,wo),grid % Функция веса
pause
step(won,wo),grid %
Переходная характеристика
pause
bode(won,wo),grid %
ЛАЧХ И ФЧХ
pause
nyquist(won,wo),grid %
АФЧХ
pause
wonz=zpk(won)
[z,p,k]=zpkdata(won,'v')%
Нули, полюса и коэффициент передачи объекта 
uc=u-mu;
m=3;% Задание размерности АРСС - модели
% Задание начальных условий в РМНК
P=1000*eye(2*m,2*m);
Q=zeros(2*m,1);
F=Q;
% Метод РМНК
for i=1:n-m
   F=[-yc(i+m-1:-1:i);uc(i+m-1:-1:i)];
   ch=P*F;
   zn=1+F'*P*F;
   gm=ch/zn;
   P=(eye(2*m)-gm*F')*P;
   Q=Q+gm*(yc(m+i)-F'*Q);
   kf(i,1:2*m)=Q'; % Коэффициенты АРСС - модели
   Tp(i)=F'*Q;
End
% Ошибка идентификации и ее характеристики
e=yc(m+1:end)-Tp';
de=std(e);
plot(t(1:n-m),kf),grid
pause
sr=[yc(m+1:end),Tp'];
plot(sr),grid
pause
plot (e),grid
% Вычисление передаточной функции объекта и его характеристик
nun=kf(end,m+1:end);
den=[1 kf(end,1:m)];
wod=tf(nun,den,dt) % Вычисление дискретной передаточной функции объекта
won=d2c(wod) % Вычисление передаточной функции объекта
sm=ss(won)% Модель объекта в пространстве состояний
impulse(won,wo),grid % Функция веса
pause
step(won,wo),grid % Переходная характеристика
pause
bode(won,wo),grid % ЛАЧХ И ФЧХ
pause
nyquist(won,wo),grid % АФЧХ
pause
wonz=zpk(won)
[z,p,k]=zpkdata(won,'v')% Нули, полюса и коэффициент передачи объекта

Страницы