数字低通滤波器的指标如下:信号的取样频率是1000Hz，截止频率是1000Hz，截止频率是200Hz，滤波器的阶数是81，用矩形窗、汗宁窗、凯赛窗、三角窗分别设计低通滤波器。
程序如下：
wpass=0.4*pi; 通带截止频率为0.4π
w1=boxcar(81); 矩形窗时域函数
w2=hamming(81); 汉宁窗时域函数
w3=kaiser(81,7.865); 凯赛窗时域函数
w4=triang(81); 三角窗时域函数
n=1:1:81;
hd=sin(wpass*(n-41))./(pi*(n-41));
hd(41)=wpass/pi;
h1=hd.*rot90(w1);
h2=hd.*rot90(w2);
h3=hd.*rot90(w3);
h4=hd.*rot90(w4);
[a1,rad]=freqz(h1); 频率响应函数求其幅度和频率
[a2,rad]=freqz(h2);
[a3,rad]=freqz(h3);
[a4,rad]=freqz(h4);
subplot(2,2,1);
plot(rad,20*log10(abs(a1)));
xlabel('归一化频率');
ylabel('幅度(dB)');
title(' 利用矩形窗设计的滤波器');
grid on;
subplot(2,2,2);
plot(rad,20*log10(abs(a2)));
xlabel('归一化频率');
ylabel('幅度(dB)');
title('利用汗宁窗设计的滤波器');
grid on;
subplot(2,2,3);
plot(rad,20*log10(abs(a3)));
xlabel('归一化频率');
ylabel('幅度(dB)');
title('利用凯赛窗设计的滤波器');
grid on;
subplot(2,2,4);
plot(rad,20*log10(abs(a4)));
xlabel('归一化频率');
ylabel('幅度(dB)');
title('利用三角窗设计的滤波器');
grid on;
FFT(三个)
书上P140
第一个FFT
已知信号由频率15Hz、幅值为0.5的正xuan信号和频率为40Hz、幅值为2的正XUAN信号组成，数据取样频率为FS=100Hz，试分别绘制N=128点DFT的幅频图和N=1024点DFT幅频图
由题可知 信号可写为
x(t)=0.5sin(2πf1t)+2sin(2πf2t)
%MatLAB Example
clf
fs=100;
N=128;
n=0:N;
t=n/fs;
x=0.5*sin(2*pi*15*t)+2*sin(2*pi*40*t);
y=fft(x,N);
mag=abs(y);
f=(0:length(y)-1)'*fs/length(y);
subplot(221)
plot(f,mag);
xlabel(Frequency(Hz)');
ylabel(Magnitude');
title('N=128')
reid
subplot(222)
plot(f(1:N/2),mag(1:N/2));
xlabel('Frequency(Hz)');
ylabel('Magnitude');
title('N=128')
grid
fs=100;
N=1024
n=0:N-1;
t=n/fs;
x=0.5*sin(2*pi*15*t)+2*sin(2*pi*40*t);
y=fft(x,N);
mag=abs(y);
f=(0:length(y)-1)'*fs/length(y);
subplot(223)
plot(f,mag);
xlabel('Frequency(Hz)');
ylabel('Magnitude');
title('N=1024')
grod
subplot(224)
plt(f(1:N/2),mag(1:N/2));
xlabel('Frequency(Hz)');
ylabel('Magnitude');
title('N=1024')
grid
要改
