Octave Encyclopedia

demo-12c

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

course/basic/demo-12c.m

코드 인덱스로 돌아가기

카테고리

Course Basic

학습 소스 그룹

코드 길이

63

lines

작성자

won sunggyu

2025-04-22

패키지

none

pkg load 기준

전체 코드

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

# filename: demo-12c.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;
omr = sqrt(K/M);
# ze = C/2/sqrt(K*M);  # 0.05
xref = 1;

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

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

db = 0.01;
Bn = 5;
nn = floor(Bn/db);

ze_arr = 0.1:0.1:1.0;  # 감쇠비
bb = 0:db:db*(nn-1);  # 무차원 주파수축

Xp = zeros(length(ze_arr), length(bb));  # 무차원 변위응답
Xd = zeros(length(ze_arr), length(bb));  # 무차원 변위응답 (dB)

for i=1:length(ze_arr)
    ze = ze_arr(i);
    Xp(i, :) = 1 ./ (1 - bb.^2 + 2*1j*ze*bb);
    Xd(i, :) = 20 * log10(abs(Xp(i, :)) / xref);
end

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

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

# 데이터를 미리 준비하고 그리기
for i=1:length(ze_arr)
    plot(ax1(1), bb, abs(Xp(i, :)))
    plot(ax1(2), bb, Xd(i, :))
end

코드 해설

목적

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

입력

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

출력

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

실행 흐름

  1. 초기화
  2. 데이터 준비
  3. 데이터 연산
  4. M a + C v + K x = f
  5. 그래프 그리기
  6. 그래프
  7. 데이터를 미리 준비하고 그리기

핵심 함수/주제

lengthXpabsax1plotsqrtXdzeros

실습 과제

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

학습 팁

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

같은 카테고리의 다른 코드