main_03_mjy

% 주제: 신호 생성, 주파수 변환, 파워 계산을 합니다.
% 작성자 : 문제영
% 작성 날짜 : 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으로 수렴된다.