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

#define Inputlength 20//輸入信號(hào)的個(gè)數(shù)

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

int main() {

int j,k;

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

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

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];

/*********輸入函數(shù)***********/

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

{

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

}

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

/***********輸出循環(huán)**************/

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

{

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

}

while (1);

return 0;

}

/**************濾波器函數(shù)****************/

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++)//對(duì)20個(gè)輸入函數(shù)逐一處理

{

y[i] = 0.0;

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

if (i >= j) {

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

}//分子系數(shù)與輸入函數(shù)的乘積和

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

if (i > l) {

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

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

}

return;

}

總結(jié)

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

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。