📄 pred_lt_3.c
字号:
#include <mex.h>
#include "matrix.h"
void mexFunction(int nlhs,mxArray *plhs[],int nrhs,const mxArray *prhs[])
{
double *exc,*b30,*T,*fraction;
int k,frac;
double *u1,*v;
int i,j;
exc=mxGetPr(prhs[0]);
T=mxGetPr(prhs[1]);
fraction=mxGetPr(prhs[2]);
b30=mxGetPr(prhs[3]);
//创建输出
plhs[0]=mxCreateDoubleMatrix(1,40,mxREAL);
v=mxGetPr(plhs[0]);
u1=(double *)mxCalloc(200,sizeof(double));
//主程序
k=(int) T[0];
frac=(int) fraction[0];
frac=-frac;
if (frac==-1)
{
k=k+1;
frac=2;
}
for (i=0;i<200;i++)
{
u1[i]=exc[i];
}
for (i=0;i<40;i++)
{
for (j=0;j<10;j++)
{
v[i]=v[i]+u1[(i-k+160-j)]*b30[(frac+3*j)]+u1[i-k+160+j]*b30[3-frac+3*j];
u1[159+i]=v[i];
}
}
mxFree(u1);
}
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -