Code Detail

demo-12a

주파수 응답 함수 그려보기 (1자유도)

course/basic/demo-12a.m

목록으로

코드를 복사해 Octave에서 바로 실행할 수 있습니다.

카테고리

Course Basic

course-basic

코드 길이

lines

작성자

won sunggyu

2025-04-22

패키지

none

pkg load

전체 코드

59 lines

# filename: demo-12a.m
# writer: won sunggyu
# date: 2025-04-22
# language: octave
# description: 주파수 응답 함수 그려보기 (1자유도)

#------------------------------------------------------------------------------
# 초기화
#------------------------------------------------------------------------------

run("startup.m");
printf(fmt("{mfilename}\n", "#FF5733"));

#------------------------------------------------------------------------------
# 데이터 준비
#------------------------------------------------------------------------------

M = 1;
C = 10;
K = 10000;
F = 1;
xref = 1e-6;

# disp(sqrt(K/M)/(2*pi))

#------------------------------------------------------------------------------
# 데이터 연산
#------------------------------------------------------------------------------

# M a + C v + K x = f
# X = F / (K - omega^2 * M + 1j*omega * C)

df = 0.1;
Fn = 50;
nn = floor(Fn/df);

ff = 0:df:df*(nn-1);
X = zeros(size(ff));

for i=1:nn
    om = 2*pi*ff(i);
    X(i) = F / (K - om^2*M + 1j*om*C);
end

Xd = 20 * log10(abs(X) / xref);

#------------------------------------------------------------------------------
# 그래프 그리기
#------------------------------------------------------------------------------

# 그래프
figured("Size", [960, 960], "Move", [-1280, 0], "Name", mfilename);
ax1 = subplots(2, 1);

plot(ax1(1), ff, abs(X))
plot(ax1(2), ff, Xd)

# filename: demo-12a.m
# writer: won sunggyu
# date: 2025-04-22
# language: octave
# description: 주파수 응답 함수 그려보기 (1자유도)

#------------------------------------------------------------------------------
# 초기화
#------------------------------------------------------------------------------

run("startup.m");
printf(fmt("{mfilename}\n", "#FF5733"));

#------------------------------------------------------------------------------
# 데이터 준비
#------------------------------------------------------------------------------

M = 1;
C = 10;
K = 10000;
F = 1;
xref = 1e-6;

# disp(sqrt(K/M)/(2*pi))

#------------------------------------------------------------------------------
# 데이터 연산
#------------------------------------------------------------------------------

# M a + C v + K x = f
# X = F / (K - omega^2 * M + 1j*omega * C)

df = 0.1;
Fn = 50;
nn = floor(Fn/df);

ff = 0:df:df*(nn-1);
X = zeros(size(ff));

for i=1:nn
    om = 2*pi*ff(i);
    X(i) = F / (K - om^2*M + 1j*om*C);
end

Xd = 20 * log10(abs(X) / xref);


#------------------------------------------------------------------------------
# 그래프 그리기
#------------------------------------------------------------------------------

# 그래프
figured("Size", [960, 960], "Move", [-1280, 0], "Name", mfilename);
ax1 = subplots(2, 1);

plot(ax1(1), ff, abs(X))
plot(ax1(2), ff, Xd)

코드 해설

목적

주파수 응답 함수 그려보기 (1자유도)

입력

스크립트 상단에서 정의한 파라미터/입력 데이터를 사용합니다.

출력

그래프/figure 출력
콘솔 텍스트 출력

실행 흐름

핵심 함수

실습 과제

질량/감쇠/강성 또는 전달함수 계수를 바꿔 응답 변화를 확인해보세요.
축 범위와 라벨을 바꿔 그래프 해석성이 어떻게 달라지는지 확인해보세요.
핵심 함수 abs의 인자를 한 가지 바꿔 결과 변화를 기록해보세요.

학습 팁

그래프 비교 시 축 범위(XLim/YLim)와 단위를 먼저 고정하면 해석 오류를 줄일 수 있습니다.