[논문 리뷰] 3rd-order Spectral Representation Method: Part II -- Ergodic Multi-variate random processes with fast Fourier transform

Lohit Vandanapu, Michael D. Shields|arXiv (Cornell University)|2019. 11. 22.

Scientific Research and Discoveries참고 문헌 10인용 수 4

한 줄 요약

이 논문은 주어진 상호 전력 스펙트럼 밀도 및 상호 비스펙트럼 밀도를 바탕으로, 빠른 푸리에 변환(FFT)을 사용하여 비가역성 다변량 스토케스틱 과정을 효율적으로 시뮬레이션하기 위한 3차 스펙트럼 표현 방법을 제안한다. 이 방법은 시뮬레이션된 과정이 기대값에서 목표로 하는 2차 및 3차 통계적 성질을 정확히 만족시키며, 다중 색인 주파수 이산화를 통해 비가역성을 확보하고, FFT 구현을 통해 상당한 계산 비용 절감을 달성한다.

ABSTRACT

The second in a two-part series, this paper extends the 3rd-order Spectral Representation Method for simulation of ergodic multi-variate stochastic processes according to a prescribed cross power spectral density and cross bispectral density. The 2nd and 3rd order ensemble properties of the simulated stochastic vector processes are shown to satisfy the target cross correlation properties in expectation. A multi-indexed frequency discretization is introduced to ensure ergodicity of the sample functions. This is first shown for uni-variate processes and then the simulation formula for multi-variate processes is provided. Ensemble properties and ergodicity of the sample functions are proven. Additionally, it is shown that the simulations can be implemented efficiently with the Fast Fourier Transform, which greatly reduces computational effort. An example involving the simulation of turbulent wind velocity fluctuations is presented to further highlight the features and applications of the algorithm.

연구 동기 및 목표

• 3차 스펙트럼 표현 방법을 다변량, 비가역성 스토케스틱 과정으로 확장하기.
• 시뮬레이션된 과정이 기대값에서 목표로 하는 2차 및 3차 집합 통계를 정확히 만족시키기.
• 다중 색인 주파수 이산화 방식을 통해 표본 함수의 비가역성을 강제하기.
• 빠른 푸리에 변환(FFT) 알고리즘의 통합을 통해 효율적인 계산을 가능하게 하기.
• 난류 풍속 속도 변동 시뮬레이션과 같은 실제 적용 사례를 통해 방법의 유용성을 입증하기.

제안 방법

• 다변량 과정의 시뮬레이션된 표본 함수에서 비가역성을 확보하기 위해 다중 색인 주파수 이산화를 사용한다.
• 주어진 상호 전력 스펙트럼 밀도 및 상호 비스펙트럼 밀도에서 유도된 무작위 위상과 진폭을 갖는 삼각함수 성분을 조합하여 스토케스틱 과정을 구성한다.
• 3차 누적량을 사용하여 스펙트럼 표현을 수립함으로써, 다변량 과정 내의 비정규성 및 비선형 의존성을 정확히 모델링할 수 있다.
• 시간 도메인 실현치의 계산을 가속화하기 위해 빠른 푸리에 변환(FFT)을 적용하여 계산 비용을 크게 감소시킨다.
• 시뮬레이션된 과정의 2차 및 3차 집합 모멘트가 목표로 하는 상호상관 및 비스펙트럼 성질을 정확히 만족시킴을 보장한다.
• 스펙트럼 밀도 함수를 역 FFT 연산과 적절한 위상 무작위화를 통해 시간 도메인 실현치로 매핑하는 수치 알고리즘을 유도한다.

실험 결과

연구 질문

• RQ1 3차 스펙트럼 표현 방법은 어떻게 주어진 상호스펙트럼 및 비스펙트럼 밀도를 갖는 다변량 비가역성 스토케스틱 과정을 시뮬레이션할 수 있도록 확장될 수 있는가?
• RQ2 어떤 주파수 이산화 전략이 시뮬레이션된 과정의 표본 함수에서 비가역성을 보장하는가?
• RQ3 통계적 정밀도를 훼손하지 않으면서 이러한 과정을 시뮬레이션하는 데 드는 계산 비용은 어떻게 감소시킬 수 있는가?
• RQ4 시뮬레이션된 과정이 기대값에서 목표로 하는 2차 및 3차 통계적 성질을 어느 정도 정확히 재현하는가?
• RQ5 실제 응용에 적합하게 빠른 푸리에 변환(FFT)을 효율적으로 구현할 수 있는가?

주요 결과

• 시뮬레이션된 다변량 과정은 기대값에서 목표로 하는 2차 상호상관 성질을 정확히 일치시킨다.
• 3차 상호비스펙트럼 밀도는 시뮬레이션된 과정에서 유지되어 비정규성 및 비선형 의존성을 정확히 표현함을 보장한다.
• 다중 색인 주파수 이산화의 사용은 수학적으로 증명된 바와 같이 표본 함수의 비가역성을 성공적으로 강제한다.
• 빠른 푸리에 변환의 통합으로 계산 복잡도가 감소하여 장기간에 걸친 고해상도 과정의 효율적 시뮬레이션을 가능하게 한다.
• 이 방법은 난류 풍속 속도 변동의 시뮬레이션에 성공적으로 적용되어 실용적 적용 분야에서의 관련성과 정확성을 입증하였다.
• 이론적 증명을 통해 제안된 프레임워크 하에서 집합 통계 및 시간 평균 통계 성질이 모두 목표 값으로 수렴하는 것으로 확인되었다.