完善资料让更多小伙伴认识你,还能领取20积分哦, 立即完善>
[y,fs,nbits]=wavread('speech'); %窗函数高通滤波 n=length(y); %求出语音信号的长度 noise=0.01*randn(n,2); %随机函数产生噪声 s=y+noise; %语音信号加入噪声 S=fft(s); %傅里叶变换 Ft=8000; Fp=1000; Fs=1200; wp=2*Fp/Ft; ws=2*Fs/Ft; rp=1; rs=50; p=1-10.^(-rp/20); s=10.^(-rs/20); fpts=[wp ws]; mag=[0 1]; dev=[p s]; [n23,wn23,beta,ftype]=kaiserord(fpts,mag,dev);%kaiserord求阶数截止频率 b23=fir1(n23,wn23,Kaiser(n23+1,beta));%由firl设计滤波器 z23=fftfilt(b23,s); sound(z23); m23=fft(z23); %求滤波后的信号 figure; subplot(2,2,1); plot(abs(S),'g'); title('滤波前信号的频谱'); subplot(2,2,2); plot(abs(m23),'r'); title('滤波后信号的频谱'); subplot(2,2,3); plot(s); title('滤波前信号的波形'); subplot(2,2,4); plot(z23); title('滤波后信号的波形');[y,fs,nbits]=wavread('speech'); %窗函数高通滤波 n=length(y); %求出语音信号的长度 noise=0.01*randn(n,2); %随机函数产生噪声 s=y+noise; %语音信号加入噪声 S=fft(s); %傅里叶变换 Ft=8000; Fp=1000; Fs=1200; wp=2*Fp/Ft; ws=2*Fs/Ft; rp=1; rs=50; p=1-10.^(-rp/20); s=10.^(-rs/20); fpts=[wp ws]; mag=[0 1]; dev=[p s]; [n23,wn23,beta,ftype]=kaiserord(fpts,mag,dev);%kaiserord求阶数截止频率 b23=fir1(n23,wn23,Kaiser(n23+1,beta));%由firl设计滤波器 z23=fftfilt(b23,s); sound(z23); m23=fft(z23); %求滤波后的信号 figure; subplot(2,2,1); plot(abs(S),'g'); title('滤波前信号的频谱'); subplot(2,2,2); plot(abs(m23),'r'); title('滤波后信号的频谱'); subplot(2,2,3); plot(s); title('滤波前信号的波形'); subplot(2,2,4); plot(z23); title('滤波后信号的波形'); 为什么只输出一个波形?
|
|
相关推荐
1个回答
|
|
回帖奖励 +1 分积分
没用过 这个真心不会
|
|
|
|
你正在撰写答案
如果你是对答案或其他答案精选点评或询问,请使用“评论”功能。
1697 浏览 1 评论
200圆!求助大佬给一份VSG并网和离网模式的simulink仿真
1973 浏览 0 评论
MATLAB(3)--矩阵的引用(sub2ind、ind2sub、reshape函数使用)
2792 浏览 0 评论
3130 浏览 0 评论
4105 浏览 1 评论
小黑屋| 手机版| Archiver| 电子发烧友 ( 湘ICP备2023018690号 )
GMT+8, 2024-11-28 01:37 , Processed in 0.581204 second(s), Total 73, Slave 56 queries .
Powered by 电子发烧友网
© 2015 bbs.elecfans.com
关注我们的微信
下载发烧友APP
电子发烧友观察
版权所有 © 湖南华秋数字科技有限公司
电子发烧友 (电路图) 湘公网安备 43011202000918 号 电信与信息服务业务经营许可证:合字B2-20210191 工商网监 湘ICP备2023018690号