Octave Encyclopedia

main_03_pjm

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

ex-recv/03/03_pjm/main_03_pjm.m

코드 인덱스로 돌아가기

카테고리

Submission Archive

학습 소스 그룹

코드 길이

60

lines

작성자

-

날짜 정보 없음

패키지

none

pkg load 기준

전체 코드

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

######  시간/주파수 베이스 Power Sum 구하기 및 hanning window에 의한 비교
######  2025. 04. 07. pjm, 5th, ECO
clc; clear; close all;

dt = 1e-4;
N = 2e4;
t = [1:N]*dt;
df = 1/t(end);


#### Org.
acc = Rand_Harm_Sigs( t , 5 );            ## acc( 1,N ary )
[ Acc_mir, f ] = two_side_fft( acc, t );  ## Acc_mir( 1,N ary )  / ## f( 1,N ary )
Auto_PS_mir = conj(Acc_mir) .*Acc_mir;    ## Auto_PS_mir( 1,N ary )

PS_time = PowerSum( acc,dt );
PS_freq = PowerSum( Acc_mir,1/df );
ratio_PS_org = PS_freq/PS_time;

#### hanning wind w/o coef.
acc_hann = acc .* hanning(N).';
[ Acc_mir_hann, f ] = two_side_fft( acc_hann , t );   ## Acc_mir_hann( 1,N ary )  / ## f( 1,N ary )
Auto_PS_mir_hann = conj(Acc_mir_hann) .*Acc_mir_hann; ## Auto_PS_mir_hann( 1,N ary )

%PS_time_hann = PowerSum( acc_hann,dt );
PS_freq_hann = PowerSum( Acc_mir_hann,1/df );
ratio_PS_hann_wo = PS_freq_hann/PS_time;


#### hanning wind w. coef.
Acc_mir_hann_ecf = sqrt(8/3)*Acc_mir_hann;   ## Acc_mir_hann( 1,N ary )  / ## f( 1,N ary )
Auto_PS_mir_hann_ecf = conj(Acc_mir_hann_ecf) .*Acc_mir_hann_ecf; ## Auto_PS_mir_hann( 1,N ary )

%PS_time_hann = PowerSum( acc_hann,dt );
PS_freq_hann_ecf = PowerSum( Acc_mir_hann_ecf,1/df );
ratio_PS_hann_ecf = PS_freq_hann_ecf/PS_time;

#### Power Sum 결과 비교  ( Sum of PowerSpectrum / raw_data )
#### raw_data가 랜덤이므로 오차율 변동할 수 있음
%ratio_PS_org      = 1.0000
%ratio_PS_hann_wo  = 0.3751
%ratio_PS_hann_ecf = 1.0004


figure(1,'position',[300 200 800 600],'defaultaxesfontsize',14)
hold on; grid on; grid minor;
plot(t,acc,';no wind;')
ylabel('Acc. [m/s^2]'); xlabel('Time [sec]')
hold off

figure(2,'position',[300 200 800 600],'defaultaxesfontsize',14)
hold on; grid on; grid minor;
plot(f,Auto_PS_mir,';no wind;','linewidth',3)
plot(f,Auto_PS_mir_hann,';hanning w/o coef.;','linewidth',2)
plot(f,Auto_PS_mir_hann_ecf,';hanning w. coef.;','linewidth',1)
ylabel('Mag. [(m/s^2).^2]'); xlabel('Freq. [Hz]')
hold off

코드 해설

목적

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

입력

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

출력

  • 그래프/figure 출력

실행 흐름

  1. ##### 시간/주파수 베이스 Power Sum 구하기 및 hanning window에 의한 비교
  2. ### Power Sum 결과 비교 ( Sum of PowerSpectrum / raw_data )

핵심 함수/주제

PowerSumplotconjAcc_mir_hannAuto_PS_mir_hannfiguretwo_side_fftxlabel

실습 과제

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

학습 팁

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

같은 카테고리의 다른 코드