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

UptoLike

ВУЗ: 

ПГУ | Пенза

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

Рубрика: 

w=tf(2,[10,1]) % Передаточная функция системы
tm=1000;
t=1:tm; %
Задание времени наблюдения
n=length(t); %
Вычисление длины вектора времени
u=randn(1,tm); % Формирование входного сигнала
y=lsim(w,u,t); %
Формирование выходного сигнала
plot(t,u,t,y)
tau=-tm+1:1:tm-1;
mu=mean(u) %
Вычисление среднего значения входного сигнала
my=mean(y)
% Вычисление среднего значения выходного сигнала
du=std(u)%
Вычисление среднеквадратичного отклонения входного сигнала
dy=std(y) %
Вычисление среднеквадратичного отклонения выходного сигнала
ruu=xcorr(u,u,'biased'); %
Вычисление корреляционной функции входного сигнала
ryu=xcorr(y,u,'biased'); %
Вычисление корреляционной функции выходного
сигнала
ryy=xcorr(y,y,'biased'); %
Вычисление взимной корреляционной функции
subplot (3,1,1)
plot(tau,ruu),grid
subplot(3,1,2)
plot(tau,ryu),grid
subplot(3,1,3)
plot(tau,ryy),grid
pause
subplot(2,1,2)
[Su,f]=psd(u,n,1); %
Вычисление спектральной плотности входного сигнала
[Sy,f]=psd(y,n,1); %
Вычисление спектральной плотности выходного сигнала
subplot(2,1,1)
plot(f,Su),grid
subplot(2,1,2)
plot(f,Sy),grid
m
u
=-0,0376:
m
y
= -0,0736:
σ
u
= 1,0191:
σ
y
= 0,3901/ 
      w=tf(2,[10,1]) % Передаточная функция системы
      tm=1000;
      t=1:tm; % Задание времени наблюдения
      n=length(t); % Вычисление длины вектора времени
      u=randn(1,tm); % Формирование входного сигнала
      y=lsim(w,u,t); % Формирование выходного сигнала
      plot(t,u,t,y)
      tau=-tm+1:1:tm-1;
      mu=mean(u) % Вычисление среднего значения входного сигнала
      my=mean(y) % Вычисление среднего значения выходного сигнала
      du=std(u)% Вычисление среднеквадратичного отклонения входного сигнала
      dy=std(y) % Вычисление среднеквадратичного отклонения выходного сигнала
      ruu=xcorr(u,u,'biased'); % Вычисление корреляционной функции входного сигнала
      ryu=xcorr(y,u,'biased'); % Вычисление корреляционной функции выходного
сигнала
      ryy=xcorr(y,y,'biased'); % Вычисление взимной корреляционной функции
      subplot (3,1,1)
      plot(tau,ruu),grid
      subplot(3,1,2)
      plot(tau,ryu),grid
      subplot(3,1,3)
      plot(tau,ryy),grid
      pause
      subplot(2,1,2)
      [Su,f]=psd(u,n,1); % Вычисление спектральной плотности входного сигнала
      [Sy,f]=psd(y,n,1); % Вычисление спектральной плотности выходного сигнала
      subplot(2,1,1)
      plot(f,Su),grid
      subplot(2,1,2)
      plot(f,Sy),grid

            mu=-0,0376: my= -0,0736: σu= 1,0191: σy= 0,3901/

Страницы