




已阅读5页,还剩8页未读, 继续免费阅读
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
实验五 谱分析 一 实验目的 1 研究不同类型的窗函数,研究一些不同的方法来测试窗的性能; 2 专注于有关窄带信号的几个不同的情形。 二 实验原理 信号是无限长的,而在进行信号处理时只能采用有限长信号,所以需要将 窄带信号“截断”。在信号处理中,“截断”被看成是用一个有限长的“窗口”看无限长的信倍号,或者从分析的角度是无限长的信号x(t)乘以有限长的窗函数w(t),由傅立叶变换性质可知:x(t)w(t)=1/2X(jw)W(jw) 如果x(t)是频带有限信号,而w(t)是频带无限函数,截断后的信号也必是频带无限信号,从而产生所谓的频谱泄露。频谱泄露是不可避免的,但是尽量减小,因此设计了不同的窗函数满足不同的要求。从能量的角度,频谱泄露也是能量泄露,因为加窗后,是原来的信号集中在宅频带内的能量分散到无限的频带范围。三 实验内容、程序及结果用MATLAB编程绘制各种窗函数的形状及其幅度响应。w1=boxcar(30);subplot(421)stem(w1)title(juxingchang)xlabel(t);ylabel(w1);w2=hanning(30);subplot(422)stem(w2)title(hanningchuang)xlabel(t);ylabel(w2);w3=hamming(30);subplot(423)stem(w3)title(hanmingchuang);xlabel(t);ylabel(w3);w4=bartlett(30);subplot(424)stem(w4)title(batelitechuang)xlabel(t);ylabel(w4);w5=blackman(30);subplot(425)stem(w5)title(balaikemanchuang)xlabel(t);ylabel(w5);w6=triang(30);subplot(426)stem(w6)title(triang)xlabel(t);ylabel(w6);w7=kaiser(30,60);subplot(427)stem(w7)title(kaiser)xlabel(t);ylabel(w7);w8=chebwin(30,60);subplot(428)stem(w8)title(chebwin)xlabel(t);ylabel(w8);2,各种窗函数的幅频响应 w1=boxcar(30);X,W=dtft(w1,300)subplot(421)plot(W/2/pi,abs(X);title(juxingchuang);xlabel(w);ylabel(|W(jw)|);w2=hanning(30);X,W=dtft(w2,300)subplot(422)plot(W/2/pi,abs(X);title(hanningchuang);xlabel(w);ylabel(|W(jw)|);w3=hamming(30);X,W=dtft(w3,300)subplot(423)plot(W/2/pi,abs(X);title(hanmingchuang);xlabel(w);ylabel(|W(jw)|);w4=bartlett(30);X,W=dtft(w4,300)subplot(424)plot(W/2/pi,abs(X);title(batelite);xlabel(w);ylabel(|W(jw)|);w5=blackman(30);X,W=dtft(w5,300)subplot(425)plot(W/2/pi,abs(X);title(bulaikemanchuang);xlabel(w);ylabel(|W(jw)|);w6=triang(30);X,W=dtft(w6,300)subplot(426)plot(W/2/pi,abs(X);title(triang);xlabel(w);ylabel(|W(jw)|);w7=kaiser(30,60);X,W=dtft(w7,300)subplot(427)plot(W/2/pi,abs(X);title(kaiser);xlabel(w);ylabel(|W(jw)|);w8=chebwin(30,60);X,W=dtft(w8,300)subplot(428)plot(W/2/pi,abs(X);title(chebwin);xlabel(w);ylabel(|W(jw)|);3矩形窗的幅频响应w1=boxcar(10);X,W=dtft(w1,300)subplot(221)plot(W/2/pi,abs(X);title(juxingchuang);xlabel(w);ylabel(|W(jw)|);w2=boxcar(20);X,W=dtft(w2,300)subplot(222)plot(W/2/pi,abs(X);title(juxingchuang);xlabel(w);ylabel(|W(jw)|);w3=boxcar(50);X,W=dtft(w3,300)subplot(223)plot(W/2/pi,abs(X);title(juxingchuang);xlabel(w);ylabel(|W(jw)|);w4=boxcar(100);X,W=dtft(w4,300)subplot(224)plot(W/2/pi,abs(X);title(juxingchuang);xlabel(w);ylabel(|W(jw)|);4.1信号周期的0.9矩形窗和汉宁窗f=25/16;T=2.56;fs=8;N=0.9*T*fs;n=0:(N-1);x=0.75+3.4*cos(2*pi*f*n/fs)+2.7*cos(4*pi*f*n/fs)+1.5*sin(3.5*pi*f*n/fs)+2.5*sin(7*pi*f*n/fs);w1=boxcar(N);y1=x.*w1;Y1,W=dtft(y1,1000);subplot(221),plot(W/2/pi,abs(Y1);title(信号周期0.9的矩形窗幅频);xlabel(f);ylabel(|Y1|);subplot(222),plot(W/2/pi,angle(Y1);title(信号周期0.9的矩形窗相频);xlabel(f);ylabel(Y1);w2=hanning(N);y2=x.*w2;Y2,W=dtft(y2,1000);subplot(223),plot(W/2/pi,abs(Y2);title(信号周期0.9的汉宁窗幅频);xlabel(f);ylabel(|Y2|);subplot(224),plot(W/2/pi,angle(Y2);title(信号周期0.9的汉宁窗相频);xlabel(f);ylabel(Y2);汉明窗和巴特利特窗f=25/16;T=2.56;fs=8;N=0.9*T*fs;n=0:(N-1);x=0.75+3.4*cos(2*pi*f*n/fs)+2.7*cos(4*pi*f*n/fs)+1.5*sin(3.5*pi*f*n/fs)+2.5*sin(7*pi*f*n/fs);w1=hamming(N);y1=x.*w1;Y1,W=dtft(y1,1000);subplot(221),plot(W/2/pi,abs(Y1);title(信号周期0.9的汉明窗幅频);xlabel(f);ylabel(|Y1|);subplot(222),plot(W/2/pi,angle(Y1);title(信号周期0.9的汉明窗相频);xlabel(f);ylabel(Y1);w2=bartlett(N);y2=x.*w2;Y2,W=dtft(y2,1000);subplot(223),plot(W/2/pi,abs(Y2);title(信号周期0.9的巴特利特窗幅频);xlabel(f);ylabel(|Y2|);subplot(224),plot(W/2/pi,angle(Y2);title(信号周期0.9的巴特利特窗相频);xlabel(f);ylabel(Y2);布莱克曼窗和Triang窗f=25/16;T=2.56;fs=8;N=0.9*T*fs;n=0:(N-1);x=0.75+3.4*cos(2*pi*f*n/fs)+2.7*cos(4*pi*f*n/fs)+1.5*sin(3.5*pi*f*n/fs)+2.5*sin(7*pi*f*n/fs);w1=blackman(N);y1=x.*w1;Y1,W=dtft(y1,1000);subplot(221),plot(W/2/pi,abs(Y1);title(信号周期0.9的布莱克曼窗幅频);xlabel(f);ylabel(|Y1|);subplot(222),plot(W/2/pi,angle(Y1);title(信号周期0.9的布莱克曼窗相频);xlabel(f);ylabel(Y1);w2=triang(N);y2=x.*w2;Y2,W=dtft(y2,1000);subplot(223),plot(W/2/pi,abs(Y2);title(信号周期0.9的Triang窗幅频);xlabel(f);ylabel(|Y2|);subplot(224),plot(W/2/pi,angle(Y2);title(信号周期0.9的Triang窗相频);xlabel(f);ylabel(Y2);Kaiser窗和切比雪夫窗f=25/16;T=2.56;fs=8;N=0.9*T*fs;n=0:(N-1);x=0.75+3.4*cos(2*pi*f*n/fs)+2.7*cos(4*pi*f*n/fs)+1.5*sin(3.5*pi*f*n/fs)+2.5*sin(7*pi*f*n/fs);w1=kaiser(N);y1=x.*w1;Y1,W=dtft(y1,1000);subplot(221),plot(W/2/pi,abs(Y1);title(信号周期0.9的Kaiser窗幅频);xlabel(f);ylabel(|Y1|);subplot(222),plot(W/2/pi,angle(Y1);title(信号周期0.9的Kaiser窗相频);xlabel(f);ylabel(Y1);w2=chebwin(N);y2=x.*w2;Y2,W=dtft(y2,1000);subplot(223),plot(W/2/pi,abs(Y2);title(信号周期0.9的切比雪夫窗幅频);xlabel(f);ylabel(|Y2|);subplot(224),plot(W/2/pi,angle(Y2);title(信号周期0.9的切比雪夫窗相频);xlabel(f);ylabel(Y2);4.2信号周期的1.1矩形窗和汉宁窗f=25/16;T=2.56;fs=8;N=1.1*T*fs;n=0:(N-1);x=0.75+3.4*cos(2*pi*f*n/fs)+2.7*cos(4*pi*f*n/fs)+1.5*sin(3.5*pi*f*n/fs)+2.5*sin(7*pi*f*n/fs);w1=boxcar(N);y1=x.*w1;Y1,W=dtft(y1,1000);subplot(221),plot(W/2/pi,abs(Y1);title(信号周期1.1的矩形窗幅频);xlabel(f);ylabel(|Y1|);subplot(222),plot(W/2/pi,angle(Y1);title(信号周期1.1的矩形窗相频);xlabel(f);ylabel(Y1);w2=hanning(N);y2=x.*w2;Y2,W=dtft(y2,1000);subplot(223),plot(W/2/pi,abs(Y2);title(信号周期1.1的汉宁窗幅频);xlabel(f);ylabel(|Y2|);subplot(224),plot(W/2/pi,angle(Y2);title(信号周期1.1的汉宁窗相频);xlabel(f);ylabel(Y2);汉明窗和巴特利特窗f=25/16;T=2.56;fs=8;N=1.1*T*fs;n=0:(N-1);x=0.75+3.4*cos(2*pi*f*n/fs)+2.7*cos(4*pi*f*n/fs)+1.5*sin(3.5*pi*f*n/fs)+2.5*sin(7*pi*f*n/fs);w1=hamming(N);y1=x.*w1;Y1,W=dtft(y1,1000);subplot(221),plot(W/2/pi,abs(Y1);title(信号周期1.1的汉明窗幅频);xlabel(f);ylabel(|Y1|);subplot(222),plot(W/2/pi,angle(Y1);title(信号周期1.1的汉明窗相频);xlabel(f);ylabel(Y1);w2=bartlett(N);y2=x.*w2;Y2,W=dtft(y2,1000);subplot(223),plot(W/2/pi,abs(Y2);title(信号周期1.1的巴特利特窗幅频);xlabel(f);ylabel(|Y2|);subplot(224),plot(W/2/pi,angle(Y2);title(信号周期1.1的巴特利特窗相频);xlabel(f);ylabel(Y2);布莱克曼窗和Triang窗f=25/16;T=2.56;fs=8;N=1.1*T*fs;n=0:(N-1);x=0.75+3.4*cos(2*pi*f*n/fs)+2.7*cos(4*pi*f*n/fs)+1.5*sin(3.5*pi*f*n/fs)+2.5*sin(7*pi*f*n/fs);w1=blackman(N);y1=x.*w1;Y1,W=dtft(y1,1000);subplot(221),plot(W/2/pi,abs(Y1);title(信号周期1.1的布莱克曼窗幅频);xlabel(f);ylabel(|Y1|);subplot(222),plot(W/2/pi,angle(Y1);title(信号周期1.1的布莱克曼窗相频);xlabel(f);ylabel(Y1);w2=triang(N);y2=x.*w2;Y2,W=dtft(y2,1000);subplot(223),plot(W/2/pi,abs(Y2);title(信号周期1.1的Triang窗幅频);xlabel(f);ylabel(|Y2|);subplo
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 通讯行业会计的工作计划
- 建立互助养老服务网络的计划
- 水务产业发展战略计划
- 促进学校与家长的紧密合作的学期工作计划
- 2025年合成橡胶型胶粘剂项目建议书
- 强化供需监测机制保障市场秩序
- 2025年人造岗石树脂项目建议书
- 2025年硅粉系列合作协议书
- 国有土地使用权租赁合同
- 第二季度部门工作计划表
- 快手申诉文本
- 现金调拨系统操作手册教学课件
- 学校物业管理宿舍管理方案995
- 荔枝依旧年年红
- SMT贴片线项目可行性研究报告
- 小学语文必备文学常识常考100题汇总(含答案)
- 第五版-FMEA-新版FMEA【第五版】
- 英语人教版高中必修三(2019新编)第一单元教案
- 口腔修复学第三章牙体缺损的修复
- 建设部环卫劳动定额
- 美容心理学概述(美容心理学课件)
评论
0/150
提交评论