s= tf ('s') ; Kp =2 ; G = 1 /(s^2 +3*s +2 ); T= feedback(G* Kp , 1 ,-1 ) T = 2 ------------- s^2 + 3 s + 4 연속시간 전달 함수입니다 . 모델 속성 e= (1-T) e = s^2 + 3 s + 2 ------------- s^2 + 3 s + 4 step(T) grid on 점검문제 5. 1 다음 공정을 비례 제어기 kp =2 로 제어해보자 . 주어진 G(S) 값이 있고 스텝응답을 그리고 정상상태 오차를 구하라 . (ess 구하기 ) 입력이 스텝 응답 일때 기준 ess 는 1 / 2 가 됩니다 .점검문제 5.2 다음 공정을 미분 제어기 kd =2 로 제어해보자 . 주어진 G(S) 값이 있고 스텝응답을 그리고 정상상태 오차를 구하라 . s= tf ('s') ; KD =2*s ; G = 1 /(s^2 +3*s +2 ); T= feedback(G*KD , 1 ,-1 ) T = 2 s ------------- s^2 + 5 s + 2 연속시간 전달 함수입니다 . 모델 속성 e= (1-T) e = s^2 + 3 s + 2 ------------- s^2 + 5 s + 2 step(T) grid on (ess 구하기 ) 입력이 스텝응답일때 기준으로 ess =1 이 됩니다 .점검문제 5.3 다음 공정을 적분 제어기 kI =2 로 제어해보자 . 주어진 G(S) 값이 있고 스텝응답을 그리고 정상상태 오차를 구하라 . s= tf ('s') ; KI =2/s ; G = 1 /(s^2 +3*s +2 ); T= feedback(G*KI , 1 ,-1 ) T = 2 --------------------- s^3 + 3 s^2 + 2 s + 2 연속시간 전달 함수입니다 . 모델 속성 e= (1-T) e = s^3 + 3 s^2 + 2 s --------------------- s^3 + 3 s^2 + 2 s + 2 step(T) grid on (ess 구하기 ) 입력이 스텝응답 일때 ess= 0 이 됩니다 .점검문제 5.4 다음 공정을 비례 미분 제어기로 제어해보자 . 주어진 G(s) 값이 있고 스텝응답을 그리고 정상상태 오차를 구하시오 . 오버슈트를 5% 이내로 만들기 위한 게인값을 구해보시오 . s= tf ('s') ; KD =1*s ; KP = 5 ; C= KD +KP ; G = 1 /(s^2 +3*s +2 ); T= feedback(G*C , 1 ,-1 ) T = s + 5 ------------- s^2 + 4 s + 7 연속시간 전달 함수입니다 . 모델 속성 e = (1 -T) e = s^2 + 3 s + 2 ------------- s^2 + 4 s + 7 연속시간 전달 함수입니다 . 모델 속성 step(T) 입력이 스텝 응답 일때 ess= 2 / 7점검문제 5.5 다음 공정을 비례 적분 제어기로 제어해보자 . 주어진 G(s) 값이 있고 스텝응답을 그리고 정상상태 오차를 구하시오 . 오버슈트를 5% 이내로 만들기 위한 게인값을 구해보시오 . s= tf ('s') ; KI =1; KP = 1 ; C= KP + KI/s ; G = 1 /(s^2 +3*s +2 ); T= feedback(G*C , 1 ,-1 ) T = s + 1 --------------------- s^3 + 3 s^2 + 3 s + 1 연속시간 전달 함수입니다 . 모델 속성 step(T) ess= 0점검문제 5.6 다음 공정을 비례 적분 미분 제어기로 제어해보자 . 주어진 G(s) 값이 있고 스텝응답을 그리고 정상상태 오차를 구하시오 . 오버슈트를 5% 이내로 만들기 위한 게인값을 구해보시오 . s= tf ('s') ; KP = 1 KI =1.23; KD = 1 ; C= KP + KD*s + KI/s ; G = 1 /(s^2 +3*s +2 ); T= feedback(G*C , 1 ,-1 ) KP = 1 T = s^2 + s + 1.23 ------------------------ s^3 + 4 s^2 + 3 s + 1.23 연속시간 전달 함수입니다 . 모델 속성 e = 1-T e = s^3 + 3 s^2 + 2 s ------------------------ s^3 + 4 s^2 + 3 s + 1.23 연속시간 전달 함수입니다 . 모델 속성 ess = 0연습문제 5.1 다음 이차 공정의 안정성을 확인하고 오버슈트가 5% 이내로 오도록 k 값을 선정하여 보시오 . routh array 계산 결과 k 0 일때 안정 된다 . s = tf (' s '); k= 0.5 ; G = 1/( s *(s+1)); T = feedback (k* G , 1, -1) T = 0.5 ------------- s^2 + s + 0.5 연속시간 전달 함수입니다. 모델 속성 step ( T )연습문제 5.2 다음 삼차 공정의 안정성을 확인하고 안정화하기 위해 k 값을 선정하여 보시오 . s = tf (' s ') ; k = 10 ; G = (s+1) / (s+2) ; D = (k) / ( s *(s^2 +5* s +6) ); T = feedback ( G * D , 1 ,-1 ) T = 10 s + 10 -------------------------------- s^4 + 7 s^3 + 16 s^2 + 22 s + 10 연속시간 전달 함수입니다. 모델 속성 step ( T )연습문제 5.3 다음 이차 공정의 안정성을 확인하고 안정화 하기 위해 PD 제어기를 선정하여 근을 확인해 보시오 . s = tf (' s ') ; K= 2; kp = 3; kd = 2.5; G =1/((s+1) * (s-1) ) ; D = kp+kd * s ; T = feedback ( G * D , 1, -1) T = 2.5 s + 3 --------------- s^2 + 2.5 s + 2 연속시간 전달 함수입니다. 모델 속성 step ( T )연습문제 5.4 공정이 다음과 같이 비최소위상 이다 . 먼저 시스템의 안정성을 조사하고 , 만약 불안정하다면 왜 불안정 한지를 예를 들어 설명하고 , 시스템이 안정화 하기 위해 제어기를 선정한 뒤 설계해 보시오 . s = tf (' s ') ; G = 1 /(s-1) ; T = feedback ( G , 1, -1) T = 1 - s 연속시간 전달 함수입니다. 모델 속성 ( 풀이 0 이기에 neutral 상태 ) s= tf ('s') ; kp = 100; kd =0.1; KI= 1; G=1/(s-1); D = kp+kd * s+KI /s; T = feedback(G*D,1,-1) step(T) T = 0.1 s^2 + 100 s + 1 ------------------- 1.1 s^2 + 99 s + 1 연속시간 전달 함수입니다 . 모델 속성{nameOfApplication=Show}
점검문제 1.2 A= [ 1 2 3 ; 2 3 2; 3 2 1] A = 1 2 3 2 3 2 3 2 1 save data A load data A clear load data A clear load data A who 사용자의 변수 : A a= A(3,:) a = 3 2 1점검문제 1.5 help plot plot - 2 차원 선 플롯 X 값에 대한 Y 데이터의 2 차원 선 플롯을 생성합니다 . 벡터 및 행렬 데이터 plot(X,Y) plot( X,Y,LineSpec ) plot(X1,Y1,..., Xn,Yn ) plot(X1,Y1,LineSpec1,..., Xn,Yn,LineSpecn ) plot(Y) plot( Y,LineSpec ) 테이블 데이터 plot( tbl,xvar,yvar ) plot( tbl,yvar ) 추가 옵션 plot(ax,___) plot(___, Name,Value ) p = plot(___)점검문제 1.7 a = [ 1 2 3 ] ; disp (a) 1 2 3 disp ( ' MATLAB is so easy') MATLAB is so easy점검문제 1.8 A= [1 2 ; 3 4]; B= [5 6 ; 7 8]; F= A'*B + (B'*A*B)^(-1) F = -40.5000 87.5000 95.2500 -5.5000점검문제 1.9 temp F = -40.5000 87.5000 95.2500 -5.5000점검문제 1.11 m=3 ; n=2; for i = 1:m, for j = 1:n , H ( i,j ) = 1/(1+j-1); end end H = 1.0000 0.5000 1.0000 0.5000 1.0000 0.5000점검문제 1.12 t = linspace (0,20,20); k=1 ; tmax = 20; while t (k) tmax , y (k) = cos ( pi * t (k)); k=k+1; end plot ( y ) xlabel (' time ( s )') ylabel (' y ( t )') grid점검문제 1.14 x=[-4: 0.1 :4]; k=1 ; a=2 ; for i =1:length(x) if x( i ) = -k y( i ) = 0; elseif x( i ) = k y( i ) = ( a/(2*k) ) * x( i ) ; else x( i ) k y( i ) = a ; end end plot ( x,y ) xlabel ( 'x' ) ylabel (' y' ) grid ;{nameOfApplication=Show}