목차
1.주제선정 - 우리별 1, 2호 인공위성
2. 시스템 모델링
3. 상태행렬의 가제어성(Controllability) 판별
4. 상태행렬의 가관측성(Observability) 판별
5. 설계 제한 조건
6. Reguation Problem
7. Proper poles
8. Realize control law, estimator, and compensator
[m-file]
[참고문헌]
2. 시스템 모델링
3. 상태행렬의 가제어성(Controllability) 판별
4. 상태행렬의 가관측성(Observability) 판별
5. 설계 제한 조건
6. Reguation Problem
7. Proper poles
8. Realize control law, estimator, and compensator
[m-file]
[참고문헌]
본문내용
잡았다. 즉 정착시간을 40초이내, 30%의 오버슈트를 고려한다. 그리고 출력에서 센서의 외란에 의한 영향을 1% 이내로 하였다.
6. Reguation Problem
-프로젝트의 설계에서 reference input을 도입하지 않았기 때문에 regulation problem이다.
7. Proper poles
-SRL설계에서 원하는 폴의 위치를 매트랩의 lqr(), lqe()함수를 이용하여 Try & Error 방식으로 위의 설계제한조건에 맞는 값을 찾으려고 하였다. 값이 크게 되면 정착시간이 빨라지지만 그만큼 오버슈트가 크게 생기고 작으면 그 반대의 상황이 발생한다. 그리고 q값이 크게 되면 시스템의 pole이 플랜트와 제어기의 pole에 의해 결정되어 estimator가 없는 것과도 같이 작용하지만 센서의 외란에 의한 영향이 그만큼 커지게 된다. 반면 q값이 작으면 센서의 외란에 의한 영향은 작아지지만 estimator의 폴의 영향이 커져서 전체 시스템이 반응이 느려지는 상황이 발생한다. 이러한 제한 조건들을 고려하여 =1, q=200으로 최종 결정하였다.
8. Realize control law, estimator, and compensator
[Control law]
R= 1;
rho = 1;
Q=rho* H' * H;
K=lqr(F, G, Q, R)
[Estimator]
qe=200;
gam=qe*G;
q1=gam'*gam;
rv=1;
L=lqe(F, gam, H, q1, rv)
[compensator]
ac=F-G*K-L*H ;bc=L;cc=K;dc=0;
[cnum, cden]=ss2tf(ac, bc, cc, dc)
[Simulink]
[Output1]
-위의 결과 그래프를 보면 overshoot가 30%이고 정착시간이 40초임을 알 수 있다. 그리고 외란에 의한 영향은 0.5% 이내로 위에서 설정한 제한조건을 만족시키는 것을 알 수 있다.
[Output2]
-위의 결과는 ramp 입력에 의한 결과로 16초만 되어도 입력을 잘 추종하는 것을 알 수 있다.
[m-file]
F=[0 1 0 0 ;
-0.091/0.4565 -0.0036/0.4565 0.091/0.4565 0.0036/0.4565 ;
0 0 0 1 ;
0.091/4.565 0.0036/4.565 -0.091/4.565 -0.0036/4.565];
G=[0 ; 0 ; 0 ; 1/4.565];
H=[1 0 0 0];
J=0;
[num,den]=ss2tf(F,G,H,J) %전달함수
C=[G F*G F^2*G F^3*G] % Controllability matrix
c=det(C) % Controllable 판별
O=[H ; H*F ; H*F^2 ; H*F^3] % Observability matrix
o=det(O) % Observable 판별
clf;
% satellite system matrices
j1=4.827;
j2=j1/10;
k=0.091;
b=0.0036;
F =[ 0 1.0000 0 0;
-k/j2 -b/j2 k/j2 b/j2;
0 0 0 1.0000 ;
k/j1 b/j1 -k/j1 -b/j1];
G =[0;
0;
0;
1/j1];
H =[1 0 0 0];
J =[0];
%Contorl law
R= 1;
rho = 1;
Q=rho* H' * H;
K=lqr(F, G, Q, R);
%Estimator
qe=200;
gam=qe*G;
q1=gam'*gam;
rv=1;
L=lqe(F, gam, H, q1, rv);
%Plant transper function
[gnum, gden]=ss2tf(F, G, H, J);
%Compensator transper function
ac=F-G*K-L*H ;bc=L;cc=K;dc=0;
[cnum, cden]=ss2tf(ac, bc, cc, dc);
[참고문헌]
1. Feedback Control of Dynamic Systems / Gene F. Franklin 외 / Prentice Hall / 2006
2. 우리별 위성의 자세제어 시스템 검증 방법 / 이현우 외(한국과학기술원 인공위성연구센터) / ICASE MAGAZINE Vol. 6, No. 5, 2000
3. 우리별 1, 2호의 기계구조 개발 및 환경시험 결과 / 배정석 외(한국과학기술원 인공위성 연구센터) / J. Astron. space sci. 13(2)S94-S104(1996)
6. Reguation Problem
-프로젝트의 설계에서 reference input을 도입하지 않았기 때문에 regulation problem이다.
7. Proper poles
-SRL설계에서 원하는 폴의 위치를 매트랩의 lqr(), lqe()함수를 이용하여 Try & Error 방식으로 위의 설계제한조건에 맞는 값을 찾으려고 하였다. 값이 크게 되면 정착시간이 빨라지지만 그만큼 오버슈트가 크게 생기고 작으면 그 반대의 상황이 발생한다. 그리고 q값이 크게 되면 시스템의 pole이 플랜트와 제어기의 pole에 의해 결정되어 estimator가 없는 것과도 같이 작용하지만 센서의 외란에 의한 영향이 그만큼 커지게 된다. 반면 q값이 작으면 센서의 외란에 의한 영향은 작아지지만 estimator의 폴의 영향이 커져서 전체 시스템이 반응이 느려지는 상황이 발생한다. 이러한 제한 조건들을 고려하여 =1, q=200으로 최종 결정하였다.
8. Realize control law, estimator, and compensator
[Control law]
R= 1;
rho = 1;
Q=rho* H' * H;
K=lqr(F, G, Q, R)
[Estimator]
qe=200;
gam=qe*G;
q1=gam'*gam;
rv=1;
L=lqe(F, gam, H, q1, rv)
[compensator]
ac=F-G*K-L*H ;bc=L;cc=K;dc=0;
[cnum, cden]=ss2tf(ac, bc, cc, dc)
[Simulink]
[Output1]
-위의 결과 그래프를 보면 overshoot가 30%이고 정착시간이 40초임을 알 수 있다. 그리고 외란에 의한 영향은 0.5% 이내로 위에서 설정한 제한조건을 만족시키는 것을 알 수 있다.
[Output2]
-위의 결과는 ramp 입력에 의한 결과로 16초만 되어도 입력을 잘 추종하는 것을 알 수 있다.
[m-file]
F=[0 1 0 0 ;
-0.091/0.4565 -0.0036/0.4565 0.091/0.4565 0.0036/0.4565 ;
0 0 0 1 ;
0.091/4.565 0.0036/4.565 -0.091/4.565 -0.0036/4.565];
G=[0 ; 0 ; 0 ; 1/4.565];
H=[1 0 0 0];
J=0;
[num,den]=ss2tf(F,G,H,J) %전달함수
C=[G F*G F^2*G F^3*G] % Controllability matrix
c=det(C) % Controllable 판별
O=[H ; H*F ; H*F^2 ; H*F^3] % Observability matrix
o=det(O) % Observable 판별
clf;
% satellite system matrices
j1=4.827;
j2=j1/10;
k=0.091;
b=0.0036;
F =[ 0 1.0000 0 0;
-k/j2 -b/j2 k/j2 b/j2;
0 0 0 1.0000 ;
k/j1 b/j1 -k/j1 -b/j1];
G =[0;
0;
0;
1/j1];
H =[1 0 0 0];
J =[0];
%Contorl law
R= 1;
rho = 1;
Q=rho* H' * H;
K=lqr(F, G, Q, R);
%Estimator
qe=200;
gam=qe*G;
q1=gam'*gam;
rv=1;
L=lqe(F, gam, H, q1, rv);
%Plant transper function
[gnum, gden]=ss2tf(F, G, H, J);
%Compensator transper function
ac=F-G*K-L*H ;bc=L;cc=K;dc=0;
[cnum, cden]=ss2tf(ac, bc, cc, dc);
[참고문헌]
1. Feedback Control of Dynamic Systems / Gene F. Franklin 외 / Prentice Hall / 2006
2. 우리별 위성의 자세제어 시스템 검증 방법 / 이현우 외(한국과학기술원 인공위성연구센터) / ICASE MAGAZINE Vol. 6, No. 5, 2000
3. 우리별 1, 2호의 기계구조 개발 및 환경시험 결과 / 배정석 외(한국과학기술원 인공위성 연구센터) / J. Astron. space sci. 13(2)S94-S104(1996)
추천자료
공업제어계측제어기능사
[홈 오토메이션]인터넷을 통한 무선 제어
[공학법제] 국내개발품의 특허 사례(김치냉장고의 숙성제어 방법)
P(Proportinal 비례) I(Integral 적분) D(Differential 미분) 제어
공정제어화공 실험-다단 액위 제어
전기기기 설계 - matlab을 이용한 DC 모터 속도제어 설계 - 계자저항 변화에 따른 속도제어 /...
[실험] 직교 로봇의 동작 원리와 구조 : 2축 로봇 제어 - 실험목적, 로봇의 정의, ROBOT의 구...
무인항공기(UAV)의 시스템 구성, 항공기 자세 측정 방법 및 자세제어
제어 시스템 설계 Term project - 온도 제어시스템에 대한 제어기 설계
[재료공학실험] 비례제어기를 갖는 폐루프 제어 - 2차 지연요소를 제어대상으로 하여 비례제...
[컴퓨터 프로그래밍] C++개요, C++의 기본 형식, C++의 변수, 연산자, 제어문, 배열과 포인터
제어공학실험Ⅱ <비례제어기를 갖는 폐루프제어> - 개루프 전달함수, 폐루프 전달함수
궤환이 시스템에 미치는 영향 (피드백,외란, 궤환 제어시스템)
AVR을 이용하여 서보모터 제어하기 (서보모터 제어방법,AVR,ATmega128,소스코드, 회로도,PWM,...
소개글