mgupi.org.ru

ТАУ
ИС1 спец.2101 и 2108
3 курс
Филимонов А.Б.

Синтез модального регулятора и Синтез наблюдателя

Типовая курсовая работа

Программа курса. Весенний семестр 2006. Летняя сессия

Всё то же самое, но 6 семестр. Листинг двух m-файлов с каментами:
Первая лаба. Синтез модального регулятора

Код:

% Синтез модального регулятора
clc; clear;
n=2;
A0=[0 0;0 -1]
B0=[1;3]
C0=[]
D0=[]
% Конструирование ss-модели объекта
ObSS=ss(A0,B0,C0,D0)
rootG=-0.5;
LamdaG=rootG*ones(1,n);
% Построение матричного усилителя
K=acker(A0,B0,LamdaG);
A=A0-B0*K;
% Проверка спектра
eig(A)
% Конструирование ss-модели САР
B=B0
C=eye(2);
D=zeros(2,1);
SRSS=ss(A,B,C,D);
% Моделирование процесса регулирования
figure(1); clf;
x0=[1;0]
initial(SRSS,'r',x0)
ylabel('x');
xlabel('t');
legend('x1','x2');

Вторая лаба. Синтез наблюдателя

Код:

% Синтез наблюдателя
clc;clear;
n=2;
A0=[0 0;0 -0.1];
B0=[1;3];
C0=[2 1];
D0=[0];
% Конструирование ss-модели объекта
ObjSS=ss(A0,B0,C0,D0);
rootE=-1;
LamdaE=rootE*ones(1,n);
L=(acker(A0',C0',LamdaE))';
AE=A0-L*C0;
eig(AE)
% Конструирование модели оценивателя
BE=[L B0];
CE=eye(2);
DE=zeros(2,2);
EstSS=ss(AE,BE,CE,DE);
% Подключение оценивателя к объекту
sys1=[ObjSS;1];
sys=EstSS*sys1;
% Моделирование входного сигнала
tau=5
tF=36
tD=0.02
[U,T]=gensig('square',tau,tF,tD)
figure(1);
clf;
subplot(2,1,1);
plot(T,0,'k',T,U,'b');
axis([0 tF -0.5 1.5]);
ylabel('u');
% Моделирование процесса оценивания
xEst0=[0;0];
xObj0=[1;1]
x0=[xEst0;xObj0];
[YE,T,X]=lsim(sys,U,T,x0);
nEst=[1:n];
nObj=[(n+1):2*n];
XObj=X(:,nObj);
figure(2);clf;
plot(T,XObj(:,1),'b',T,XObj(:,2),'r',T,YE(:,1),':b',T,YE(:,2),':r');
grid;
ylabel('x,YE');
xlabel('t');
legend('x1','x2','yE1','yE2');

Типовая курсовая работа, для тех, кто сам ничего не нашёл.

Крайний срок сдачи курсового - 31.05.06 в 10:00 в Укудемии, в 230 ауд.

Программа курса
Весенний семестр 2006

на экзамене можно будет лекциями пользоваться? Rolling Eyes

CITROEN
Будем... Как обычно...

Осталось узнать, когда можно будет ещё защитить курсовую %)

что значит когда еще?
Разве не завтра?

вот бы завтра прокатила эта халява. Да поможет нам БОГ!

аминь...

Добавлено спустя 39 секунд:

надо покурить перед сном не забыть...

Харе молиться))) че никто не сделал чтоль?)

сделали, только вот прокатит ли, да и вообще никто же ни в зуб ногой о чем там речь...

Добавлено спустя 24 секунды:

а вообще чего не спим-то?
завтра кстати к 10, а не к 9...

Вот поэтому то и не спим что к 10=)

как вас тут много... тема плавно перетекает из тау в никуда %)

Код:

а вообще чего не спим-то?
завтра кстати к 10, а не к 9...

УРА!
А че спать надо курсовик доделывать.
А вообще, все здадим... рано или поздно Smile