mirror of
				https://github.com/f4exb/sdrangel.git
				synced 2025-10-29 03:50:46 -04:00 
			
		
		
		
	
		
			
				
	
	
		
			86 lines
		
	
	
		
			1.8 KiB
		
	
	
	
		
			C++
		
	
	
	
	
	
			
		
		
	
	
			86 lines
		
	
	
		
			1.8 KiB
		
	
	
	
		
			C++
		
	
	
	
	
	
| /*
 | |
|  * Filters from Fldigi.
 | |
| */
 | |
| 
 | |
| #ifndef	_FFTFILT_H
 | |
| #define	_FFTFILT_H
 | |
| 
 | |
| #include <complex>
 | |
| #include "gfft.h"
 | |
| 
 | |
| //----------------------------------------------------------------------
 | |
| 
 | |
| class fftfilt {
 | |
| enum {NONE, BLACKMAN, HAMMING, HANNING};
 | |
| 
 | |
| public:
 | |
| 	typedef std::complex<float> cmplx;
 | |
| 
 | |
| 	fftfilt(float f1, float f2, int len);
 | |
| 	fftfilt(float f2, int len);
 | |
| 	~fftfilt();
 | |
| // f1 < f2 ==> bandpass
 | |
| // f1 > f2 ==> band reject
 | |
| 	void create_filter(float f1, float f2);
 | |
| 	void create_dsb_filter(float f2);
 | |
|     void create_asym_filter(float fopp, float fin); //!< two different filters for in band and opposite band
 | |
| 
 | |
| 	int noFilt(const cmplx& in, cmplx **out);
 | |
| 	int runFilt(const cmplx& in, cmplx **out);
 | |
| 	int runSSB(const cmplx& in, cmplx **out, bool usb, bool getDC = true);
 | |
| 	int runDSB(const cmplx& in, cmplx **out);
 | |
| 	int runAsym(const cmplx & in, cmplx **out, bool usb); //!< Asymmetrical fitering can be used for vestigial sideband
 | |
| 
 | |
| protected:
 | |
| 	int flen;
 | |
| 	int flen2;
 | |
| 	g_fft<float> *fft;
 | |
| 	g_fft<float> *ift;
 | |
| 	cmplx *filter;
 | |
|     cmplx *filterOpp;
 | |
| 	cmplx *data;
 | |
| 	cmplx *ovlbuf;
 | |
| 	cmplx *output;
 | |
| 	int inptr;
 | |
| 	int pass;
 | |
| 	int window;
 | |
| 
 | |
| 	inline float fsinc(float fc, int i, int len) {
 | |
| 		return (i == len/2) ? 2.0 * fc:
 | |
| 				sin(2 * M_PI * fc * (i - len/2)) / (M_PI * (i - len/2));
 | |
| 	}
 | |
| 
 | |
| 	inline float _blackman(int i, int len) {
 | |
| 		return (0.42 -
 | |
| 				 0.50 * cos(2.0 * M_PI * i / len) +
 | |
| 				 0.08 * cos(4.0 * M_PI * i / len));
 | |
| 	}
 | |
| 
 | |
| 	void init_filter();
 | |
| 	void init_dsb_filter();
 | |
| };
 | |
| 
 | |
| 
 | |
| 
 | |
| /* Sliding FFT filter from Fldigi */
 | |
| class sfft {
 | |
| #define K1 0.99999
 | |
| public:
 | |
| 	typedef std::complex<float> cmplx;
 | |
| 	sfft(int len);
 | |
| 	~sfft();
 | |
| 	void run(const cmplx& input);
 | |
| 	void fetch(float *result);
 | |
| private:
 | |
| 	int fftlen;
 | |
| 	int first;
 | |
| 	int last;
 | |
| 	int ptr;
 | |
| 	struct vrot_bins_pair;
 | |
| 	vrot_bins_pair *vrot_bins;
 | |
| 	cmplx *delay;
 | |
| 	float k2;
 | |
| };
 | |
| 
 | |
| #endif
 |