Yoichi Nagashima / Mbed 2 deprecated CQ_nagasm_08_notch_test

add 60Hz/50Hz notch filter

Dependencies:   mbed

main.cpp@0:5ec74240d7f1, 2014-12-16 (annotated)

Committer:
nagasm
Date:
Tue Dec 16 08:28:20 2014 +0000
Revision:
0:5ec74240d7f1
Child:
1:8fa2f521009a
add 60Hz/50Hz notch filter

Who changed what in which revision?

UserRevisionLine numberNew contents of line
nagasm 0:5ec74240d7f1 1 #include "mbed.h"
nagasm 0:5ec74240d7f1 2 #include "sub.hpp"
nagasm 0:5ec74240d7f1 3 #include "FIR_LPF.hpp"
nagasm 0:5ec74240d7f1 4 #include "notch.hpp"
nagasm 0:5ec74240d7f1 5
nagasm 0:5ec74240d7f1 6 int main(){
nagasm 0:5ec74240d7f1 7 common_setup();
nagasm 0:5ec74240d7f1 8 xbee.baud(38400);
nagasm 0:5ec74240d7f1 9 xbee.attach(&rx_fifoset, xbee.RxIrq);
nagasm 0:5ec74240d7f1 10 timer_setup.attach_us(&timer_interrupt, 5); // 5usec
nagasm 0:5ec74240d7f1 11 coef_set(60); // Ham notch filter set
nagasm 0:5ec74240d7f1 12 while(1){
nagasm 0:5ec74240d7f1 13 tx_fifo_check();
nagasm 0:5ec74240d7f1 14 if(timer_value[1] > 19){ // 0.1msec sampling
nagasm 0:5ec74240d7f1 15 timer_value[1] = 0;
nagasm 0:5ec74240d7f1 16 float data1 = (float)gain * (analog_value2.read() - 0.5f); // A/D in (3)
nagasm 0:5ec74240d7f1 17 if (data1 < 0) data1 = -data1; // always detection ON
nagasm 0:5ec74240d7f1 18 if (fir_lpf != 0) data1 = FIR_calc1(data1); // FIR calc (1) call
nagasm 0:5ec74240d7f1 19 data1 = notch_filter1(data1);
nagasm 0:5ec74240d7f1 20 if(timer_value[2] > 2999){ // 15msec
nagasm 0:5ec74240d7f1 21 timer_value[2] = 0;
nagasm 0:5ec74240d7f1 22 if(max_count != 0) data1 = move_mean_calc1(data1);
nagasm 0:5ec74240d7f1 23 tx_message((uint16_t)((data1 + 1.0f) * 2047)<<4);
nagasm 0:5ec74240d7f1 24 }
nagasm 0:5ec74240d7f1 25 }
nagasm 0:5ec74240d7f1 26 if(timer_value[3] > 19){ // 0.1msec sampling
nagasm 0:5ec74240d7f1 27 timer_value[3] = 0;
nagasm 0:5ec74240d7f1 28 float data2 = (float)gain * (analog_value3.read() - 0.5f); // A/D in (4)
nagasm 0:5ec74240d7f1 29 if (data2 < 0) data2 = -data2; // always detection ON
nagasm 0:5ec74240d7f1 30 if (fir_lpf != 0) data2 = FIR_calc2(data2); // FIR calc (1) call
nagasm 0:5ec74240d7f1 31 data2 = notch_filter2(data2);
nagasm 0:5ec74240d7f1 32 if(timer_value[4] > 2999){ // 15msec
nagasm 0:5ec74240d7f1 33 timer_value[4] = 0;
nagasm 0:5ec74240d7f1 34 if(max_count != 0) data2 = move_mean_calc2(data2);
nagasm 0:5ec74240d7f1 35 tx_message(0x400000 + ((uint16_t)((data2 + 1.0f) * 2047)<<4));
nagasm 0:5ec74240d7f1 36 }
nagasm 0:5ec74240d7f1 37 }
nagasm 0:5ec74240d7f1 38 if(timer_value[0] > 199999){ // 1000msec
nagasm 0:5ec74240d7f1 39 timer_value[0] = 0;
nagasm 0:5ec74240d7f1 40 myled = !myled;
nagasm 0:5ec74240d7f1 41 }
nagasm 0:5ec74240d7f1 42 if(rx_fifo_check() == 1){
nagasm 0:5ec74240d7f1 43 int sum = 0;
nagasm 0:5ec74240d7f1 44 for (int i=0; i<6; i++) sum += conv_hex(raw_data[i])<<(4*(5-i));
nagasm 0:5ec74240d7f1 45 tx_message(sum); // Echo Back
nagasm 0:5ec74240d7f1 46 if(sum>>16 == 0x80){
nagasm 0:5ec74240d7f1 47 switch((sum & 0xff00)>>8){
nagasm 0:5ec74240d7f1 48 case 0x00:
nagasm 0:5ec74240d7f1 49 fir_lpf = sum & 0x01;
nagasm 0:5ec74240d7f1 50 break;
nagasm 0:5ec74240d7f1 51 case 0x01:
nagasm 0:5ec74240d7f1 52 gain = sum & 0x0f;
nagasm 0:5ec74240d7f1 53 break;
nagasm 0:5ec74240d7f1 54 case 0x02:
nagasm 0:5ec74240d7f1 55 max_count = sum & 0x7f;
nagasm 0:5ec74240d7f1 56 if (max_count>100) max_count = 100;
nagasm 0:5ec74240d7f1 57 sum_clear();
nagasm 0:5ec74240d7f1 58 break;
nagasm 0:5ec74240d7f1 59 case 0x03:
nagasm 0:5ec74240d7f1 60 coef_set(sum & 0x3f);
nagasm 0:5ec74240d7f1 61 break;
nagasm 0:5ec74240d7f1 62 }
nagasm 0:5ec74240d7f1 63 }
nagasm 0:5ec74240d7f1 64 }
nagasm 0:5ec74240d7f1 65 }
nagasm 0:5ec74240d7f1 66 }