编程实现任意确定信号的频谱分析算法

(N2,Rp2,Wn2); %切贝雪夫滤波器系统函数[h2,w2]=freqz(num2,den2); %计算频谱响应 figure(4); %第四个图形 subplot(3,1,1); %在窗口中生成 3行 1列共 3个子图,当前激活第 1个子图 plot(abs(h2));grid; %画出幅值图 xlabel('\omega/\pi'); %在 xabel 中显示归一化?/? ylabel(' 振幅');%定义图形纵坐标为''振幅' title(' 契比雪夫Ⅰ型带通滤波器的幅频响应');%定义标题 f3=filter(num2,den2,x); %滤波 9 y3=fft(f3,8000); %作 fft变换 subplot(3,1,2); %在窗口中生成 3行 1列共 3个子图,当前激活第 2个子图 plot(f3); %做原始语音信号的时域图形 title(' 带通滤波后的信号');%定义标题 xlabel(' 时间');%定义图形横坐标为'时间’ ylabel(' 幅值');%定义图形纵坐标为''幅值' subplot(3,1,3); %在窗口中生成 3行 1列共 3个子图,当前激活第 2个子图 plot(abs(y3(1:4000))); %画出频谱图像 title(' 带通滤波后的信号频谱')%定义标题 xlabel(' 频率');%定义图形横坐标为'频率’ ylabel(' 幅值');%定义图形纵坐标为''幅值' wavwrite(f3,' 带通.wav'); %写出过滤后带通信号 x3=wavread(' 带通.wav'); %读出过滤后带通信号 sound(x3); %播放带通后信号执行如下: 4.3. 对滤波后的音频信号再滤出三个能量最集中的频簇 3 、在低、中、高三个频带中,各滤出三个能量最集中的频簇, 显示滤波后信号的时域和频域曲线。程序代码如下: