- 1、本文档共32页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
MATLAB语言实现对模拟信号的处理
题目要求:? ?1.录制一段自己的语音信号,并对录制的信号进行采样;? ?2.画出采样后的语音信号的时域波形和频谱图;? ?3.给定滤波器的性能指标,采用窗函数法和双线性变换法设计滤波器,? ?并划出滤波器的频域响应;? ?4.用该滤波器对采集的信号进行滤波,画出滤波后信号的时域波形和频谱,? ?并对滤波前后的信号进行对比,分析信号的变化;? ?5.回放语音信号;? ?6.设计一个信号处理系统界面。? ?报告要求:??根据自己录制的语音信号选择合适的滤波器,并按照题目要求完成报告;
程序1:用MATLAB对原始语音信号进行分析,画出它的时域波形和频谱
fs=22050;? ?? ?? ?? ?? ?? ?%语音信号采样频率为22050x1=wavread(C:\teng.wav); %读取语音信号的数据,赋给变量x1sound(x1,22050);? ?? ?? ???%播放语音信号y1=fft(x1,1024);? ?? ?? ???%对信号做1024点FFT变换f=fs*(0:511)/1024;figure(1)plot(x1)? ?? ?? ?? ?? ?? ? %做原始语音信号的时域图形title(原始语音信号);xlabel(time n);ylabel(fuzhi n);figure(2)freqz(x1)? ?? ?? ?? ?? ?? ?%绘制原始语音信号的频率响应图title(频率响应图)figure(3)subplot(2,1,1);plot(abs(y1(1:512)))? ?? ? %做原始语音信号的FFT频谱图title(原始语音信号FFT频谱)subplot(2,1,2);plot(f,abs(y1(1:512)));title(原始语音信号频谱)xlabel(Hz);ylabel(fuzhi);
程序2:给原始的语音信号加上一个高频余弦噪声,频率为5kHz。画出加噪后的语音信号时域和频谱图,与原始信号对比,可以很明显的看出区别。
fs=22050; x1=wavread(C:\teng.wav);f=fs*(0:511)/1024;t=0:1/22050:(size(x1)-1)/22050; %将所加噪声信号的点数调整到与原始信号相同Au=0.03; d=[Au*cos(2*pi*5000*t)]; %噪声为5kHz的余弦信号x2=x1+d;sound(x2,22050); %播放加噪声后的语音信号y2=fft(x2,1024);figure(1)plot(t,x2)title(加噪后的信号);xlabel(time n);ylabel(fuzhi n);figure(2)subplot(2,1,1);plot(f,abs(y1(1:512)));title(原始语音信号频谱);xlabel(Hz);ylabel(fuzhi);subplot(2,1,2);plot(f,abs(y2(1:512)));title(加噪后的信号频谱);xlabel(Hz);ylabel(fuzhi);
程序3:双线性变换法设计Butterworth滤波器
fs=22050;x1=wavread(h:\课程设计2\shuzi.wav);t=0:1/22050:(size(x1)-1)/22050;Au=0.03;d=[Au*cos(2*pi*5000*t)];x2=x1+d;wp=0.25*pi;ws=0.3*pi;Rp=1;Rs=15;Fs=22050;Ts=1/Fs;wp1=2/Ts*tan(wp/2);? ?? ?? ?? ?? ???%将模拟指标转换成数字指标ws1=2/Ts*tan(ws/2); [N,Wn]=buttord(wp1,ws1,Rp,Rs,s);??%选择滤波器的最小阶数[Z,P,K]=buttap(N);? ?? ?? ?? ?? ?? ?%创建butterworth模拟滤波器[Bap,Aap]=zp2tf(Z,P,K);[b,a]=lp2lp(Bap,Aap,Wn);? ?[bz,az]=bilinear(b,a,Fs);? ?? ?? ???%用双线性变换法实现模拟滤波器到数字滤波器的转换[H,W]=freqz(bz,az);? ?? ?? ?? ?? ???%绘制频率响应曲线figure(1)plot(W*Fs/(2*pi),abs(H))gridxlabel(频率/Hz)ylabel(频率响应幅度)title(Butterworth)f1=filter(bz,az,x2);
您可能关注的文档
- IMD动画.ppt
- Inertnet域名IP搜索程序实现.doc
- inside job中文字幕.doc
- Ikebana日本插花艺术.ppt
- Imagewar入门教程.doc
- INTECONT PLUS申克.doc
- Integrating design and construction through vi英文文献.doc
- INTERNET技术与应用课程论文.doc
- Internet网络流量的自相似性分析.docx
- iPad评测报告:挑战笔记本霸主地位.doc
- 新的一年工作展望.docx
- 医生年终个人工作的述职报告(3篇).docx
- 2023年消防设施操作员之消防设备中级技能考前冲刺练习题附答案详解.docx
- 2022-2023年环境影响评价工程师之环评技术导则与标准通关练习题包括详细解答.docx
- 2023年中级注册安全工程师之安全生产管理考前冲刺检测卷和答案.docx
- 2023年中级银行从业资格之中级银行管理考前冲刺测试卷提供答案解析.docx
- 2023年公共营养师之二级营养师通关模拟考试试卷附带答案.docx
- 证券分析师之发布证券研究报告业务考前冲刺模拟题库.docx
- 2022-2023年二级建造师之二建建设工程法规及相关知识综合提升测试卷附答案.docx
- 2023年二级建造师之二建机电工程实务通关模拟考试试卷提供答案解析.docx
文档评论(0)