Kenji Arai
CW Decoder (Morse code decoder) 1st release version. Only run on Nucleo-F446RE mbed board.
Dependencies: Array_Matrix F446_AD_DA ST7565_SPI_LCD TextLCD UIT_FFT_Real
Fork of
F446_MySoundMachine
by 
不韋 呂
Base on F446_MySoundMachine program created by 不韋 呂-san.
Thanks to 不韋 呂-san making fundamental part such as FFT and ADC high speed interrupt driven program.
I just combined LCD and show CW code.
Diff: SignalProcessing/BilinearDesignLH.cpp
- Revision:
- 0:fa74b1130cc3
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/SignalProcessing/BilinearDesignLH.cpp Sun Jan 29 09:11:30 2017 +0000
@@ -0,0 +1,76 @@
+//------------------------------------------------------------------------
+// Design of Butterworth LPF and HPF using bilinear transform
+//
+// 2016/03/31, Copyright (c) 2016 MIKAMI, Naoki
+//------------------------------------------------------------------------
+
+#include "BilinearDesignLH.hpp"
+
+namespace Mikami
+{
+ // Execute design
+ // input
+ // fc: Cutoff frequency
+ // pb: Passband (LPF or HPF)
+ // output
+ // c : Coefficients for cascade structure
+ // g : Gain factor for cascade structure
+ void BilinearDesign::Execute(float fc, Type pb, Coefs c[], float& g)
+ {
+ Butterworth();
+ Bilinear(fc);
+ ToCascade(pb);
+ GetGain(pb);
+ GetCoefs(c, g);
+ }
+
+ // Get poles for Butterworth characteristics
+ void BilinearDesign::Butterworth()
+ {
+ float pi_2order = PI_/(2.0f*ORDER_);
+ for (int j=0; j<ORDER_/2; j++) // Pole with imaginary part >= 0
+ {
+ float theta = (2.0f*j + 1.0f)*pi_2order;
+ sP_[j] = Complex(-cosf(theta), sinf(theta));
+ }
+ }
+
+ // Bilinear transform
+ // fc: Cutoff frequency
+ void BilinearDesign::Bilinear(float fc)
+ {
+ float wc = tanf(fc*PI_FS_);
+ for (int k=0; k<ORDER_/2; k++)
+ zP_[k] = (1.0f + wc*sP_[k])/(1.0f - wc*sP_[k]);
+ }
+
+ // Convert to coefficients for cascade structure
+ void BilinearDesign::ToCascade(Type pb)
+ {
+ for (int j=0; j<ORDER_/2; j++)
+ {
+ ck_[j].a1 = 2.0f*real(zP_[j]); // a1m
+ ck_[j].a2 = -norm(zP_[j]); // a2m
+ ck_[j].b1 = (pb == LPF) ? 2.0f : -2.0f; // b1m
+ ck_[j].b2 = 1.0f; // b2m
+ }
+ }
+
+ // Calculate gain factor
+ void BilinearDesign::GetGain(Type pb)
+ {
+ float u = (pb == LPF) ? 1.0f : -1.0f;
+ float g0 = 1.0f;
+ for (int k=0; k<ORDER_/2; k++)
+ g0 = g0*(1.0f - (ck_[k].a1 + ck_[k].a2*u)*u)/
+ (1.0f + (ck_[k].b1 + ck_[k].b2*u)*u);
+ gain_ = g0;
+ }
+
+ // Get coefficients
+ void BilinearDesign::GetCoefs(Coefs c[], float& gain)
+ {
+ for (int k=0; k<ORDER_/2; k++) c[k] = ck_[k];
+ gain = gain_;
+ }
+}