小弟正在做毕业论文，是关于语音信号滤波。现在想编写一个MATLAB程序，滤除高于5000Hz频率，改怎么编写？
在网上已经找到了相关程序，只要改改就行了
[y,fs,bits]=wavread('chushi.wav');
sound(y,fs)
n=length(y)
y_p=fft(y,n);
f=fs*(0:n/2-1)/n;figure(1)
subplot(2,1,1);
plot(y);
title('原始语音信号采样后的时域波形');
xlabel('时间轴')
ylabel('幅值A')
subplot(2,1,2);
plot(f,abs(y_p(1:n/2)));
title('原始语音信号采样后的频谱图');
xlabel('频率Hz');
ylabel('频率幅值');
L=length(y)
noise=0.1*randn(L,2);
y_z=y+noise;
sound(y_z,fs)
n=length(y);
y_zp=fft(y_z,n);
f=fs*(0:n/2-1)/n;figure(2)
subplot(2,1,1);plot(y_z);
title('加噪语音信号时域波形');
xlabel('时间轴')
ylabel('幅值A')
subplot(2,1,2);
plot(f,abs(y_zp(1:n/2)));
title('加噪语音信号频谱图');
xlabel('频率Hz');
ylabel('频率幅值');
Ft=5000;
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=[1 0];dev=[p s];
[n21,wn21,beta,ftype]=kaiserord(fpts,mag,dev);
b21=fir1(n21,wn21,Kaiser(n21+1,beta));
[h,w]=freqz(b21,1);plot(w/pi,abs(h));
title('FIR低通滤波器');x=fftfilt(b21,y_z);
X=fft(x,n);figure(4);
subplot(2,2,1);plot(f,abs(y_zp(1:n/2)));
title('滤波前信号的频谱');
subplot(2,2,2);plot(f,abs(X(1:n/2)));
title('滤波后信号的频谱');
subplot(2,2,3);plot(y_z);
title('滤波前信号的时域波形')
subplot(2,2,4);plot(x);
title('滤波后信号的时域波形')
sound(x,fs,bits)
① 中国知网也好、万方数据也好都有例子！
② 并且大部分的院校都有免费的接口！
③ 如果真没有免费的接口，那就百度知道悬赏求助下载吧！
④ 要是要外文的论文准备翻译的话，最好的办法就是【谷歌学术】
⑤ 需要什么语言的论文直接就用相应的语言搜索！100% 能找到类似的！
⑥ 至于翻译，可以直接谷歌翻一下，弄完在自己缕缕就可以了！
⑦ 要是计算机类的代码什么的到CSDN或者51CTO下载！
【友情提示】==================论文写作方法===========================
{首先就不要有马上毕业，最后一次花点钱就得了的想法}
① 其实，原创的论文网上没有免费为你代谢的！谁愿意花时间给你写这个呢？难道你在空闲的时间原以为别人提供这种毫无意义的服务么？所以：还不如自己写。主要是网上的不可靠因素太多，万一碰到骗人的，就不上算了。
② 写作论文的简单方法，首先大概确定自己的选题【这个很可能老师已经给你确定了】，然后在网上查找几份类似的文章。
③ 通读一些相关资料，对这方面的内容有个大概的了解！看看别人都从哪些方面写这个东西！
④ 参照你们学校的论文的格式，列出提纲，接着要将提纲给你们老师看看，再修改。等老师同意你这个提纲之后，你就可以补充内容！
⑤ 也可以把这几份论文综合一下，从每篇论文上复制一部分，组成一篇新的文章！然后把按自己的语言把每一部分换下句式或词，经过换词不换意的办法处理后,网上就查不到了！
⑥ 最后，到万方等数据库进行检测【这里便宜啊，每一万字才1块钱】，将扫红部分进行再次修改！
⑦ 祝你顺利完成论文!
【WARNING】=========================================================
[Ⅰ] 如果确认找人代笔，交押金的要慎重考虑！
[Ⅱ] 淘宝交易的话，一定看好他的打款时间，有的设定为【3天】，到期之后人家自然收到钱！
[Ⅲ] 希望用我的回答能让童鞋们多个心眼！
。。我，，
组哟的，
