get_signal

# m개의 주파수 성분을 포함한 신호 생성 함수
# 작성자 : 신찬영
# 작성날짜 : 25.04.07.

function [t, signal, fs] = get_signal(fs, duration, m)
%   Input :
%           fs       - 샘플링 주파수 (Hz)
%           duration - 신호 길이 (s)
%           m        - sin 함수 개수

%   Output :
%           Signal   - 중첩된 m개의 sin신호 / (1/fs : 1/fs : duration)
%           fs       - 샘플링 주파수 (Hz) / (-fs/2 : duration/fs : fs/2-duration/fs)
%           t        - 신호 길이 (s) / (1/fs : 1/fs : duration)


    t = 1/fs:1/fs:duration;  % 시간 벡터
    signal = zeros(size(t));  % 신호 초기화
    freqs = round(rand(1, m) * (fs/2));  % (0, fs/2) 사이에서 m개의 주파수

    for i = 1:m
        A = rand();  % 진폭 (0~1 사이 랜덤)
        phi = 2*pi*rand();  % 위상 (0~2pi 사이 랜덤)
        signal = signal + A * sin(2*pi*freqs(i)*t + phi); %   Output (1/fs : 1/fs : duration)
    end

end