陷波器的DSP实现.doc

陷波器的DSP实现 陷波器传递函数 其中，， 为陷波器中心频率点，为陷波器带宽，为陷波器作用程度(取0) 采用双线性变换把传递函数离散化，带入 有 其中 ，，， ， 带入参数： wc= 4.398; T=0.0001; k1=40; ， 得 a0=0.991280695424675 a1=-1.982561199111843 a2=0.991280695424675 b1=-1.982561199111843 b2=0.982561390849351 转化为差分方程为 即 编程： float x_k; float x_k1=0; float x_k2=0; float y_k; float y_k1=0; float y_k2=0; float a0=0.991280695424675; float a1=-1.982561199111843; float a2=0.991280695424675; float b1=-1.982561199111843; float b2=0.982561390849351; x_k=t.DT; y_k=a0 * x_k+a1 * x_k1+a2 * x_k2-b1 * y_k1-b2*y_k2; x_k2=x_k1; x_k1=x_k; y_k2=y_k1; y_k1=y_k; EPwm1Regs.CMPA.half.CMPA = t.DT; EPwm3Regs.CMPA.half.CMPA = y_k;//+3750; // EPwm4Regs.CMPA.half.CMPA = x_k+3750; 仿真波形 实验波形