Octave Encyclopedia

main_03_mjy

Signal Processing + Visualization 중심의 Octave 학습 예제

ex-recv/03/03_mjy/main_03_mjy.m

코드 인덱스로 돌아가기

카테고리

Submission Archive

학습 소스 그룹

코드 길이

70

lines

작성자

-

날짜 정보 없음

패키지

io, signal

pkg load 기준

전체 코드

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

% 주제: 신호 생성, 주파수 변환, 파워 계산을 합니다.
% 작성자 : 문제영
% 작성 날짜 : 20250402
clc;clear
pkg load io
pkg load signal
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 1. 주파수 성분 m개가 포함된 sinusoidal 신호를 출력하는 함수를 작성하라.

fs = 1000;
T = 1;
frequency = [50, 150, 300];  % 포함할 주파수 성분 (Hz)
amplitude = [1, 0.5, 0.8];  % 각 주파수 성분의 진폭
phases = [0, pi/4, pi/2];  % 각 주파수 성분의 위상 (라디안)

[x,t] = get_sinusoidal(amplitude,frequency,phases,fs,T);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 2. 함수의 출력이 가속도 물리량으로 가정하고 시간 그래프를 그려라.

width = 1300; height = 800;
fig1 = figure(1);
set(fig1,'Position',[200 100 width height])
plot_time(t, x,'b','시간 신호 그래프','Time [sec]','Acceleration [m/s^2]');

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 3. 시간 신호에 대해 auto-power spectrum을 two-side로 그려라.

fig2 = figure(2);
set(fig2,'Position',[200 100 width height])
plot_power_two_side(x, fs);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 4. 시간 도메인에서 power sum을 구하라. (함수 작성)
% 값 출력 포함

p_sum_time = calc_power_sum_time(x)

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 5. 주파수 도메인에서 power sum을 구하고 (함수 작성) 값을 비교하라.
% 값 출력 포함

p_sum_freq = calc_power_sum_freq(x, fs);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 6. 시간 신호에 hanning 윈도우를 적용하고 직전 과정을 반복하고 비교 결과를 주석에 남겨라. (1-5)

calc_hanning(amplitude,frequency,fs,T,phases);

% hanning 윈도우를 적용하기 전에는 스펙트럼이 퍼지는 현상이 있지만,
% 윈도우를 적용한 후에는 주파수 성분이 뚜렷하게 나타난다.

% 주파수 도메인에서 윈도우 적용하기 전에 누설 현상이 일어나는데 윈도우 적용을 하면 누설이 줄어든다.

% 시간 도메인에서는 윈도우를 적용하면 양 끝이 0으로 수렴된다.

코드 해설

목적

Signal Processing + Visualization 중심의 Octave 학습 예제

입력

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

출력

  • 그래프/figure 출력

실행 흐름

  1. 1. 주파수 성분 m개가 포함된 sinusoidal 신호를 출력하는 함수를 작성하라
  2. 2. 함수의 출력이 가속도 물리량으로 가정하고 시간 그래프를 그려라
  3. 값 출력 포함

핵심 함수/주제

figuresetcalc_hanningcalc_power_sum_freqcalc_power_sum_timeget_sinusoidalplot_power_two_sideplot_time

실습 과제

  • 샘플링 주파수나 입력 주파수를 바꿔 스펙트럼 변화를 비교해보세요.
  • 축 범위와 라벨을 바꿔 그래프 해석성이 어떻게 달라지는지 확인해보세요.
  • 핵심 함수 figure의 인자를 한 가지 바꿔 결과 변화를 기록해보세요.

학습 팁

  • 그래프 비교 시 축 범위(XLim/YLim)와 단위를 먼저 고정하면 해석 오류를 줄일 수 있습니다.
  • 입력 파일 경로가 현재 작업 디렉터리 기준인지 먼저 확인하세요.

같은 카테고리의 다른 코드