#include#define Width 3 //濾波器門寬

#define Inputlength 20//輸入信號的個數

void myfilter(int, float *, float *, int, float *, float *);

int main() {

int j,k;

float a[3] = { 1,-1.5,0.25 };//分母系數向量，可自行設置

float b[3] = { 8,-3.5,0.5 };//分子系數向量，可自行設置

float x[20] = { 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20 };//1-20的序列

float x1[20];

float y[20];

/*********輸入函數***********/

for (j = 0; j < 19; j++)

{

x1[j] = cos(0.2*3.14*x[j]);

}

myfilter(Width, a, b, Inputlength, x1, y);

/***********輸出循環**************/

for (k = 0; k < 20; k++)

{

printf("%f\n", y[k]);

}

while (1);

return 0;

}

/**************濾波器函數****************/

void myfilter(int width, float *a, float *b, int inputlength, float *x, float *y)

{

int i, j,l;

y[0] = b[0] * x[0];

for (i = 1; i < inputlength; i++)//對20個輸入函數逐一處理

{

y[i] = 0.0;

for (j = 0; j < width; j++)

if (i >= j) {

y[i] = y[i] + b[j] * x[i - j];

}//分子系數與輸入函數的乘積和

for (l = 0; l < width-1; l++)

if (i > l) {

y[i] = (y[i] - a[l + 1] * y[i-l-1])/a[0];

}//用之前的乘積和減去分母系數與過去輸出函數的乘積差

}

return;

}

總結

以上是生活随笔為你收集整理的matlab中fliter函数,用C语言实现MATLAB中Filter()函数的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。