大家看看这是怎么回事啊？
下午做了个小程序如下：load('IR014_1.mat','y')
N=1200;
n=0:1199;
f0=120;
fs=12000;
f=linspace(-fs/2,fs/2,N);
for i=1:1200
x(i)=y(2*i);
end
figure;
plot(n,x);xlim();xlabel('点数');ylabel('幅值');title('原信号');grid off;
X=fft(x,N);
Y=fftshift(X);
figure;
plot(n,X);xlim();xlabel('点数');ylabel('幅值');title('原信号FFT变换');grid off;
figure;
plot(f,abs(Y));xlim();xlabel('频率/fz');ylabel('幅值');title('原信号FFT谱');grid off;
输出的FFT谱如下图：
在此程序后面添加了
w0=2*pi*f0/fs;%归一化频率
m=0.99; %input('please input a number : m=');
B1=2*(1-m); %带宽计算
b1=;
a1=;
impz(b1,a1); %单位冲激响应
=freqz(b1,a1);
z=filter(b1,a1,x);
figure;
plot(n,z);xlim();xlabel('点数');ylabel('幅值');title('原信号经低通滤波器后的输出');grid off;
这段低通滤波器程序后，上面的FFT谱输出的怎么是这样的啊？
请问大家这是怎么回事啊？低通滤波器的程序对吗？要做一个低通1000HZ的低通滤波器程序该怎么写啊？求帮助！
页:
[1]