본문 바로가기
반응형

Matlab32

Matlab 실습 8. 감쇠비에 따른 단위계단 응답 문제 : 2차시스템에서 고유주파수가 2로 일정할때, 감쇠비의 변화에 따른 응답을 그래프로 그려보자. 코드 : clc; clear wn=2; z1=0.6; z2=1.0; z3=1.4; num=wn^2; d1=[1 2*z1*wn wn^2]; d2=[1 2*z2*wn wn^2]; d3=[1 2*z3*wn wn^2]; t=0:0.01:10; y1=step(num, d1, t); % step(분자, 분모, 시간) y2=step(num, d2, t); y3=step(num, d3, t); plot(t, y1, t, y2,'--', t, y3, ':'), grid xlabel('t(sec)'), ylabel('output y(t)'), title('Step Response') legend('Underdamped','.. 2022. 12. 24.
Matlab 실습 7.Bode 선도 문제 : \( G(s)=1/s(s+2) \)에 대한 Bode 선도를 그려보자 코드 : clc; clear w=logspace(-4,3,700); % log sacale num=1; den=[1 2 0]; sysp=tf(num,den); bode(sysp), grid 그래프 - x축은 logspace를 사용한다. 첫번째 요소는 \( 10^-4 \)를 의미하고, 두번째요소는 \( 10^3 \)을 말한다. 세번째요소는 700개의 점으로 구간을 나눈다는 의미이다. 2022. 12. 24.
Matlab 실습 6. if, else, switch를 이용한 짝수 홀수 판별 문제 : if, else, switch를 이용하여 짝수, 홀수를 판별해보자. 코드 : % 짝수 홀수 판별 clear all, close all x=input('number='); if x ~= abs(x) % 복소수 판별 error('the input x is not real') elseif x~=fix(x) % 실수 판별, 소수점 버림 error('the input x is not interger') elseif x 2022. 12. 24.
Matlab 실습 5. Bode 선도, 대역폭 문제 : \( T(s)=300/(s^2+20s+300) \)에서 Bode선도와 -3dB이 되는 주파수(대역폭)을 구해보자. 코드 : clear all, close all w=logspace(-1,3,400); nt=300; dt=[1 20 300]; syst=tf(nt,dt); [mag,phase]=bode(syst,w); magdB=20*log10(mag); index=1; for i=1:length(w)-1 if magdB(i)>-3 index=index+1; else index=index; end end BW=w(index), bode(syst,w), grid 결과 : BW = 20.6914 그래프 - Bode 선도는 x좌표가 log scale이다. - Bdoe 선도는 크기와 위상이 나타난다. - .. 2022. 12. 24.
Matlab 실습 4.Step function(2) 문제 : \( T(s)=300/(s^2+20s+300) \)에서 계단응답의 최종값과 정착시간을 계산하고, 최종값의 \(\pm2%\)인 상하한을 그려보자. 코드 : clc; clear all, close all, t=[0:0.01:2]; nt=300; dt=[1 20 300]; syst=tf(nt,dt); y=step(syst,t); finalt=length(t); yss=y(finalt); % 계단응답의 최종값 for i=1:finalt if y(i)>1.02*yss, ts=t(i); elseif y(i) 2022. 12. 24.
Matlab 실습 3. 구간별로 다른 함수의 그래프 문제 : \( y=-1 (x 2022. 12. 24.
반응형