一段matlab低通滤波器程序，求改编成C语言。

这个我刚好做过一个滤波器，事实上对时域信号做FFT，截取一定点数再做逆FFT相当于理想滤波。设计滤波器代码如下：

f1=100;f2=200;%待滤波正弦信号频率

fs=2000;%采样频率

m=(03f1)/(fs/2);%定义过度带宽

M=round(8/m);%定义窗函数的长度

N=M-1;%定义滤波器的阶数

b=fir1(N,f2/fs);%使用fir1函数设计滤波器

%输入的参数分别是滤波器的阶数和截止频率

figure(1)

[h,f]=freqz(b,1,512);%滤波器的幅频特性图

%[H,W]=freqz(B,A,N)当N是一个整数时函数返回N点的频率向量和幅频响应向量

plot(ffs/(2pi),20log10(abs(h)))%参数分别是频率与幅值

xlabel('频率/赫兹');ylabel('增益/分贝');title('滤波器的增益响应');

figure(2)

subplot(211)

t=0:1/fs:05;%定义时间范围和步长

s=sin(2pif1t)+sin(2pif2t);%滤波前信号

plot(t,s);%滤波前的信号图像

xlabel('时间/秒');ylabel('幅度');title('信号滤波前时域图');

subplot(212)

Fs=fft(s,512);%将信号变换到频域

AFs=abs(Fs);%信号频域图的幅值

f=(0:255)fs/512;%频率采样

plot(f,AFs(1:256));%滤波前的信号频域图

xlabel('频率/赫兹');ylabel('幅度');title('信号滤波前频域图');

figure(3)

sf=filter(b,1,s);%使用filter函数对信号进行滤波

%参数分别为滤波器系统函数的分子和分母多项式系数向量和待滤波信号输入

subplot(211)

plot(t,sf)%滤波后的信号图像

xlabel('时间/秒');ylabel('幅度');title('信号滤波后时域图');

axis([02 05 -2 2]);%限定图像坐标范围

subplot(212)

Fsf=fft(sf,512);%滤波后的信号频域图

AFsf=abs(Fsf);%信号频域图的幅值

f=(0:255)fs/512;%频率采样

plot(f,AFsf(1:256))%滤波后的信号频域图

xlabel('频率/赫兹');ylabel('幅度');title('信号滤波后频域图');

wp=02pi;ws=03pi;

wd=ws-wp; %计算过渡带宽度

N=ceil(8pi/wd); %根据汉宁窗计算所需h(n)长度去大于等于x的最小整数

wn=(02+03)pi/2; %计算理想理想低通滤波器通带截止频率

b=fir1(N,wn/pi,hanning(N+1)); %调用fir1计算低通FIR数字滤波器

freqz(b,025,512)

long lvbo_test(long NEW_DATA,unsigned char k)

{

long result; //声明一个变量作为滤波缓存

if( NEW_DATA < OLD_DATA )

{

result = OLD_DATA - NEW_DATA;

result = k;

result += 128; //+128是为了四舍五入

result /= 256;

result = OLD_DATA - result;

}

else if( NEW_DATA > OLD_DATA )

{

result = NEW_DATA - OLD_DATA;

result = k;

result += 128; //+128是为了四舍五入

result /= 256;

result = OLD_DATA + result;

}

else

{

result = OLD_DATA;

}

OLD_DATA = result;

return result;

}

以上就是关于一段matlab低通滤波器程序，求改编成C语言。全部的内容，包括:一段matlab低通滤波器程序，求改编成C语言。、求汉宁窗设计FIR低通滤波器 程序、二阶低通滤波器的C程序求一个，急用！！！等相关内容解答，如果想了解更多相关内容，可以关注我们，你们的支持是我们更新的动力！