matlab实现iir系数[b,a], 但差分方程实现有误！

savior120

#define IO_BUF SHORT
#define IIR_FLOATING DOUBLE

#define IIR_L_N  16   //低通滤波器阶数为16
#define IIR_H_N  6    //高通滤波器阶数为6

//高通滤波器的A系数
const IIR_FLOATING IIR_H_A[IIR_L_N]={1,-4.998,9.9919,-9.9878,4.9919,-0.99797};//

//高通滤波器的B系数
const IIR_FLOATING IIR_H_B[IIR_L_N]={0.9990,-4.9949,9.9899,-9.9899,4.9949,-0.9990};//




void iir(IO_BUF *BUF,INT N_BUF, INT N,IIR_FLOATING *coef,INT INC,IIR_FLOATING *T_BUF)
{
    INT i,j;
    IIR_FLOATING out;
    for(i=0;i<N_BUF;i=i+INC)
    {
        T_BUF[0]=(IIR_FLOATING)(*(BUF+i));
        out=0.0;
        for(j=0;j<(N*2-1);j++)
        {
            out+=T_BUF[j]*coef[j];
        }
        if(out>32767)
        {
            *(BUF+i)=32767;//float正数溢出
        }
        else if(out<-32768)
        {
            *(BUF+i)=-32768;//float负数溢出
        }
        else
        {
            *(BUF+i)=(IO_BUF)out;
        }
        //移动输入x1 x2...和y1 y2...
        for(j=N-1;j>0;j--)
        {
            T_BUF[j]=T_BUF[j-1];
        }
        for(j=(N*2-2);j>N;j--)
        {
            T_BUF[j]=T_BUF[j-1];
        }
        T_BUF[N]=out;
    }
}


void CIIRView::OnIirFilter()
{
    IIR_FLOATING coef[IIR_H_N*2];
    for(int i=0;i<IIR_H_N;i++)
    {
        coef=IIR_H_B;
    }
    for(i=IIR_H_N;i<IIR_H_N*2-1;i++)
    {
        coef=0-IIR_H_A[i-IIR_H_N+1];
    }
    IIR_FLOATING T_BUF[IIR_H_N*2];
    for(i=0;i<IIR_H_N*2;i++)
    {
        T_BUF=0.0;
    }

    iir(BUF,500, IIR_H_N,coef,1, T_BUF);
}


在iir实现过程中，老是出现(out>32767)和(out<-32768)，造成滤波实现不了，请问专家们，
上面实现是matlab求出的[b,a]系数问题，还是iir函数差分方程的实现问题？请指教，多谢了！

songzy41

楼主有否在MATLAB中仿真IIR差分方程的滤波过程？

savior120

songzy41  你好！
matlab里面模拟过，我上面iir函数在matlab里面用filter函数模拟了，
信号没出现什么问题，是对的，所以我确信a和b值的正确性。
但是，用C实现差分方程就出现问题，请高手指教，多谢了。

songzy41

用filter函数模拟还不够，要用差分方程去仿真模拟，当差分方程仿真模拟成功了，可报留差分方程仿真模拟开始时的中间过程，然而在C语言运行时把开始的中间数据也显示出来（与MATLAB差分模拟时用相同的输入），把中间数据进行对比，这样能较快地调试C语言程序，也容易发现问题。