savior120 发表于 2006-10-31 12:39

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

#define IO_BUF SHORT
#define IIR_FLOATING DOUBLE

#define IIR_L_N16   //低通滤波器阶数为16
#define IIR_H_N6    //高通滤波器阶数为6

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

//高通滤波器的B系数
const IIR_FLOATING IIR_H_B={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=(IIR_FLOATING)(*(BUF+i));
out=0.0;
for(j=0;j<(N*2-1);j++)
{
out+=T_BUF*coef;
}
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=T_BUF;
}
for(j=(N*2-2);j>N;j--)
{
T_BUF=T_BUF;
}
T_BUF=out;
}
}


void CIIRView::OnIirFilter()
{
IIR_FLOATING coef;
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;
}
IIR_FLOATING T_BUF;
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求出的系数问题,还是iir函数差分方程的实现问题?请指教,多谢了!
页: [1]
查看完整版本: matlab实现iir系数[b,a], 但差分方程实现有误!