Canmv k230 案例4.2——实时FIR 滤波器原理（二）(初版)

Canmv k230 案例4.2——实时FIR 滤波器原理

基本原理

理论待补充

Canmv k230 代码

import array
import math
from ulab import numpy as np
import utime
import time
from machine import UART
from machine import FPIOA
fpioa = FPIOA()
fpioa.set_function(5, fpioa.UART2_TXD)  
fpioa.set_function(6, fpioa.UART2_RXD)  
uart = UART(UART.UART2, baudrate=115200, bits=UART.EIGHTBITS, parity=UART.PARITY_NONE, stop=UART.STOPBITS_ONE)
PI = 3.14159265358979323846264338327950288419716939937510

rx = []   
def input_data():   
    for i in range(10000):
        data0 = 100 * math.cos(2 * PI * i / 64)
        data1  =100 * math.cos(400 * 2 * PI * i / 64)
        data2  = 0 * math.cos(3 * 2 * PI * i / 64)
        data3  = 0 * math.cos(4 * 2 * PI * i / 64)
        data4  = 0 * math.cos(5 * 2 * PI * i / 64)
        rx.append(data0 + data1 + data2 + data3 + data4)   
input_data()                             
data = np.array(rx,dtype=np.float) 
yy=[]
b=[0.0400,0.1810,0.3448,0.3448,0.1810,0.0400]
x=[0,0,0,0,0]      
y=[0,0,0,0,0]      

st=time.time_ns()
for i in range(9900):
    x=rx[0+i:5+i]
    for j in range(5):
        y[j]=b[j]*x[j]

    yy.append((sum(y)))
    uart.write("{0},{1}\n".format(rx[i],sum(y)))

VOFA+的输出结果

串口工具VOFA+的输出波形如下图所示

待续。。。