使用C++Matlab中的lp2lp函数教程详解.docx

第

使用C++Matlab中的lp2lp函数教程详解

目录1.matlab的lp2lp函数的作用2.matlab的lp2lp函数的使用方法3.C++实现3.1complex.h文件3.2lp2lp.h文件4.测试结果4.1测试文件4.2测试3阶的情况4.3测试9阶的情况

1.matlab的lp2lp函数的作用

去归一化H(s)的分母

2.matlab的lp2lp函数的使用方法

[z,p,k]=buttap(3);

disp(零点：+z);

disp(极点：+p);

disp(增益：+k);

[Bap,Aap]=zp2tf(z,p,k);%由零极点和增益确定归一化Han(s)系数

disp(Bap=+Bap);

disp(Aap=+Aap);

[Bbs,Abs]=lp2lp(Bap,Aap,86.178823974858318);%低通到低通计算去归一化Ha(s),最后一个参数就是去归一化的截止频率

disp(Bbs=+Bbs);

disp(Abs=+Abs);

3.C++实现

3.1complex.h文件

#pragmaonce

#includeiostream

typedefstructComplex

doublereal;//实数

doubleimg;//虚数

Complex()

real=0.0;

img=0.0;

Complex(doubler,doublei)

real=r;

img=i;

}Complex;

/*复数乘法*/

intcomplex_mul(Complex*input_1,Complex*input_2,Complex*output)

if(input_1==NULL||input_2==NULL||output==NULL)

std::coutcomplex_mulerror!std::endl;

return-1;

output-real=input_1-real*input_2-real-input_1-img*input_2-

output-img=input_1-real*input_2-img+input_1-img*input_2-real;

return0;

3.2lp2lp.h文件

实现方法很简单，将H(s)的分母的系数乘以pow(wc,这一项的指数)即可

#pragmaonce

#includeiostream

#includevector

#includealgorithm

#includecomplex.h

usingnamespacestd;

vectorpairComplex*,intlp2lp(vectorpairComplex*,inttf,doublewc)

vectorpairComplex*,intresult;

if(tf.size()=0||wc=0.001)

returnresult;

result.resize(tf.size());

for(inti=0;itf.size();i++)

doublecoeff=pow(wc,tf[i].second);

Complex*c=(Complex*)malloc(sizeof(Complex));

c-real=coeff*tf[i].first-real;

c-img=coeff*tf[i].first-

pairComplex*,intp(c,tf[i].second);

result[i]=p;

returnresult;

4.测试结果

4.1测试文件

#includestdio.h

#includestdlib.h

#includemath.h

#includevector

#includebuttap.h

#includezp2tf.h

#includelp2lp.h

usingnamespacestd;

#definepi((double)3.141592653589793)

intmain()

vectorComplex*poles=buttap(3);

vectorpairComplex*,int