Rand_Harm_Sigs

function [ sig ] = Rand_Harm_Sigs( time, m )
######  m개의 주파수 성분 조합으로 이루어진 시간 신호 생성
######  2025. 04. 07. pjm, 5th, ECO

####  Output
####  sig ( 1D array )

##  Inputs
##  time ( 1D array )
##  m ( scalar ) : Number of Frequency Components

  dt = time(2)-time(1);
  Fs = 1/dt;
  Fn = Fs/2;

  t_rsh = reshape(time,1,[]);
  sig = zeros( size(t_rsh) );

  Amplitude = Frquency = Phase = Frquency_tmp = zeros(m,1);
  Frquency_tmp = Fn*rand(m,1);
  Amplitude = rand(m,1);
  Frquency = sort(Frquency_tmp);
  Phs_Sinu = 360*rand(m,1);

  fprintf('\n  m  Amplitudes  Frequency  Phase\n')
  for ii = 1:m
    fprintf( '% 3d  % 10.2f  % 9.2f  % 5.f\n',ii,Amplitude(ii),Frquency(ii),Phase(ii))
  end

  sig = sum( Amplitude .*cosd( 360*Frquency.*t_rsh + Phase ), 1);
end