Roofing filters are a type of bandpass filter conventionally used in HF radio receivers in the first IF stage to limit the frequency spectrum passed on to later stages in the receiver.
The goal in applying roofing filters to a price signal is to simultaneously attenuate high frequency noise and low frequency distortion to pass an oscillating signal with a nearly zero mean for analysis and/or further calculation.
In this study, there are three filter types to choose from:
-> Roofing Filter, which passes data through a 2 pole high pass filter, then through a Super Smoother filter.
-> Gaussian Roofing Filter, which passes data through a 2 pole Gaussian high pass filter, then through a 2 pole Gaussian low pass filter.
-> Butterworth Roofing Filter, which passes data through a 2 pole Butterworth high pass filter, then through a 2 pole Butterworth low pass filter.
Each filter type has different amplitude and delay characteristics, so play around with each type and see which response suits your needs best.
There is an option to normalize the scale of the output as well. The normalization process in this script is computed by comparing positive and negative outputs to the filter's moving RMS value.
The resulting oscillator can be fed through numerous conventional indicators including Oscillator, , , etc. to generate smoother, less distorted indicators for a clearer view of turning points.
Alternatively, it can also act as an indicator itself, as implied by the corresponding color scheme included in the script.
Although roofing filters are not conventionally used in the analysis of market data, applying such spectral analysis techniques may prove to be quite useful for the design of more efficient indicators and more reliable predictions.
-> Fixed bug in RMS calculation. Big thanks to @yatrader2 for bringing this issue to my attention!
-> Reoganized script structure.
-> Added an optional compression algorithm to the script. This can be used to reduce transients (short duration, high amplitude values) by a specified level of intensity.
There are 2 small bugs in the normalization feature.
1) sum should be sma because you need the average of the squares that is sum/len.
2) you should round (or floor) the hp+lp/2 to make sure you have an integer length
You can correct this by changing line 97 to: rms = sqrt(sma(pow(roof, 2), round((hp_per + lp_per)/2)))
Here is a general rms function for reuse:
I appreciate the feedback!
I'm wondering if line 139 should be:
rms = RMS(roof, round(sqrt(hp_per * lp_per)))
ie Should it use the geometric mean rather than the arithmetic since we are averaging wave lengths. I know it is pretty subtle, but since you care about details, I thought I'd raise it...