asum.c

来自「g729音频编解码」· C语言 代码 · 共 22 行

#include <mex.h>
void mexFunction(int nlhs,mxArray *plhs[],int nrhs,const mxArray *prhs[])
{
    double *xb,*xn,*h;
    int i,j,len;
    //一个输出
    plhs[0]=mxCreateDoubleMatrix(1,40,mxREAL);
    xb=mxGetPr(plhs[0]);
    //两个输入
    xn=mxGetPr(prhs[0]);
    h=mxGetPr(prhs[1]);
    //xn的长度
    len=40;
    for (i=0;i<40;i++)
    {
        //xb(i)=xn(i:end)*h(1:end-i+1)';
        for (j=i;j<len;j++)
        {
            xb[i]=xb[i]+xn[j]*h[j-i];
        }
    }
}