طراحی و کاربرد عملی پالایه‌های رقمی در پردازش سیگنال‌های هواشناسی

روش های پالایش رقمی برای جداسازی اطلاعات مجزا از سیگنال های پیچیده هواشناسی بسیار موثر هستند. هدف از مطالعه حاضر طراحی پالایه‌های مناسب برای اعمال بر روی سیگنال های هواشناسی دما و بارش به منظور هموارسازی، جداسازی مولفه های بسامدی با اهمیت و حذف نوفه های بسامد بالای بی اهمیت است. طراحی این پالایه ها، هم در حوزه زمان و هم در حوزه بسامد امکان پذیر است که در این تحقیق طراحی و اعمال پالایه در حوزه زمان انجام شده است. نتایج اهمیت انتخاب پالایه مناسب در تحلیل محتوای طیفی موجود در سیگنال هواشناسی را نشان می دهد. با پالایش سیگنال خام دما و بارش ماهانه طی دوره آماری 1979 تا 2021 در ایستگاه های تهران و تبریز به مولفه های فصلی، سالانه و بین سالی، ارتباط خطی بین این مولفه ها به‌طور مجزا توسط ضریب همبستگی پیرسون مورد بررسی قرار گرفته و شباهت ها و تفاوت ها با هم مقایسه شده اند. همچنین پالایه همینگ نسبت به پالایه میانگین متحرک ساده، عملکرد بهتری را در تضعیف نوسان های گیبس در بخش های جانبی پاسخ بسامدی پالایه نشان می دهد. البته انتخاب تعداد وزن بیشتر در طراحی پالایه می تواند در بهبود عملکرد پالایه اثرگذار باشد. این مسئله در زمانی که تعداد مولفه های با اهمیت بیشتری در طیف بسامدیِ سیگنال هواشناسی وجود دارند، بایستی بیشتر مورد توجه قرار گیرد.

On the design and implementation of digital filters to process meteorological signals

Separation of different frequency bands in the complex and combined signals related to meteorological variables and also climatic indices requires the use of digital filtering methods. In this way, the information on different frequency bands can be organized and used. Given that these signals generally exhibit complex and nonlinear behavior, the use of mathematical filtering methods to identify their stochastic and periodic components leads to a better understanding of their behavior and helps modeling them as well. Therefore, the use of digital filters in order to recognize regular variabilities and facilitate statistical forecasting is one of the main goals in this field.The design and implementation of these filters are possible both in time and frequency domains. In the frequency space, this process is performed based on the Fourier transform of the signals on the basis of the Fast Fourier Transform Algorithm (FFT), in which the variances of the desired signal can be extracted based on spectral analysis in different frequencies. By employing different types of nonrecursive and recursive digital filters, which they can be implemented as lowpass, highpass, bandpass, and bandstop, the related signal inthe  time domain for each state can be constructed, and the corresponding spectrum can be studied. The isolated spectrum can be related to the effect of a special phenomenon that influences the main signal. In addition, it is possible to remove high frequency components from the original signal, which include noises and may not contain important information. Moreover, the process of optimal smoothing the original signal can also be carried out.In this study, different digital filters have been designed and then applied to meteorological data such as monthly surface temperature and precipitation. Two synoptic stations over Iran are selected and the related discrete monthly signals are constructed for 504 months during 19792021. Then, the moving average (MA) filter is used as a main filter, because it is the most common filter in digital signal processing (DSP), and also it is the easiest digital filter to understand and use. In spite of its simplicity, the moving average filter is optimal for a common task such as reducing random noise while retaining a sharp step response. This makes it the premier filter for time domain encoded signals. The filtering process in this study is conducted to denoise the original signals, and also to examine seasonal, annual, and interannual components of the original signals. Since employed filters are digital, they must be applied to the initial discrete signal in the form of convolution with the finite impulse response (FIR) of the filter in the time domain, or they can be applied in the form of multiplication in the frequency domain based on discrete Fourier transform and then using of the inverse Fourier transform to recover the desired signal.The results of this study show the importance of using digital filters in analyzing the spectral contents of meteorological signals. Furthermore, the Hamming filter, which is defined based on the cosine truncation and windowing, shows better performance in attenuating Gibbs oscillations in the lateral sidelobes of the filter frequency response than the simple moving average (MA) filter. In addition, the correlation analysis is carried out separately to indicate the linear relationships between different frequency components of the signals. The higher correlations are observed in annual frequency bands of the temperature and precipitation signals for the selected stations. It shows the effect of external climate forcing on both temperature and precipitation that is stemmed from the earth’s motion around the sun during a year. Obviously, choosing more weights in the design of a filter can improve the filtering performance, but it should be avoided to use more weights than necessary.