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

UptoLike

ВУЗ: 

ПГУ | Пенза

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

Рубрика: 

pr1=(lh1(n)-lh1(n-1))/dt % Первый (действительный) корень характеристического
уравнения
b1=(t(n)*lh1(n-1)-t(n-1)*lh1(n))/dt;
c1=exp(b1);%
Первая постоянная интегрирования
s1=c1*exp(pr1*t);
%
Вычисление второго слагаемого переходной характеристики
lh2=log(abs(ko-h'-s1));
plot(t,lh2),grid
%
Определение координат огибающей полулогарифмической переходной характеристики
z(1)=0;l=0;
for j=2:n
z(j)=lh2(j)-lh2(j-1);
if z(j)*z(j-1)<0 & z(j)<z(j-1)
l=l+1;
y(l)=lh2(j-1);
tt(l)=t(j-1);
end
end
%
Вычисление частоты и начальной фазы
m=7;
a11=sum(tt(1:m).*tt(1:m));
a12=sum(tt(1:m));
a21=a12;
a22=m;
j=1:m;
b1=sum((2*j-1).*tt(j))/4;
b2=sum(2*j-1)/4;
d=a11*a22-a12*a21;
d1=b1*a22-b2*a21;
d2=a11*b2-a21*b1;
x=d1/d;
z=d2/d;
w=2*pi*x; %
Частота
f=2*pi*z; % Начальная фаза
b1=sum(y(j).*tt(j));
b2=sum(y(j));
d=a11*a22-a12*a21;
d1=b1*a22-b2*a21;
d2=a11*b2-a21*b1;
x=d1/d;
z=d2/d;
p2r=x+i*w %
Корни (мнимые) характеристического уравнения
p3r=x-i*w
c2=exp(x); %
Вторая постоянная интегрирования
s2=c2*exp(x*t).*sin(w*t-f); 
pr1=(lh1(n)-lh1(n-1))/dt % Первый (действительный) корень характеристического
уравнения
b1=(t(n)*lh1(n-1)-t(n-1)*lh1(n))/dt;
c1=exp(b1);% Первая постоянная интегрирования
s1=c1*exp(pr1*t);
% Вычисление второго слагаемого переходной характеристики
lh2=log(abs(ko-h'-s1));
plot(t,lh2),grid
% Определение координат огибающей полулогарифмической переходной характеристики
z(1)=0;l=0;
for j=2:n
  z(j)=lh2(j)-lh2(j-1);
  if z(j)*z(j-1)<0 & z(j)

Страницы