版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、功率谱分析 由题目内容,设采样频率fs=1000HZ,数据长度为256,模型阶数为14,f1=200,f2=300、250。 (1)用最大熵法进行谱估计运行程序后,观察图像f1和f2相差较小时,功率谱变化更剧烈;模型的阶数越高,图像中能够获得的信息就越多,但同时计算量也就越大;增加数据长度可以获得更多的信息,提高了谱分析的分辨率,这是因为AR模型的谱估计隐含着对数据和自相关函数的外推,其长度可能会超过给定长度,分辨率不受信源信号的限制。(2) 分别用Levinson递推法和Burg法进行功率谱分析 Levinson递推法运行程序后,观察图像,f1和f2相差较小时,功率谱变化更剧烈;模型的阶数越
2、高,图像中能够获得的信息就越多,但同时计算量也就越大;增加数据长度可以获得更多的信息,提高了谱分析的分辨率,但本题中信号为正弦信号加白噪声,故图像观察不明显。 Burg法运行程序后,观察图像,f1和f2相差较小时,功率谱变化更剧烈;模型的阶数越高,图像中能够获得的信息就越多,但同时计算量也就越大;增加数据长度可以获得更多的信息,提高了谱分析的分辨率。(3) 改变信号的相位、频率、信噪比,上述谱分析结果有何变化如果正弦信号的频率过大,超过fs/2,会产生频率混叠现象,输入f1=600HZ,会在400HZ处产生一个波峰;降低信噪比会导致谱分辨率下降;信号起始相位的变动可导致谱线的偏移和分裂(我的图
3、像观察不到)。最大熵法估计N=1024;Nfft=256;Fs=1000;n=0:N-1;t=n/Fs;x1=sin(2*pi*200*t);x2=sin(2*pi*300*t); %0.3xn=x1+awgn(x1,10)+x2+awgn(x2,10); Pxx1,f=pmem(xn,14,Nfft,Fs);subplot(4,1,1)plot(f,10*log10(Pxx1);xlabel('Frequency (Hz)');ylabel('Power Spectrum (dB)');title('MEM f2/fs=0.3,Nfft=256,Ode
4、r=14');grid N=1024;Nfft=256;Fs=1000;n=0:N-1;t=n/Fs;x1=sin(2*pi*200*t);x2=sin(2*pi*250*t); %0.25xn=x1+awgn(x1,10)+x2+awgn(x2,10); Pxx1,f=pmem(xn,14,Nfft,Fs);subplot(4,1,2)plot(f,10*log10(Pxx1);xlabel('Frequency (Hz)');ylabel('Power Spectrum (dB)');title('MEM f2/fs=0.25,Nfft=25
5、6,Oder=14');grid N=1024; Nfft=512; %修改数据长度512Fs=1000;n=0:N-1;t=n/Fs;x1=sin(2*pi*200*t);x2=sin(2*pi*300*t); %0.3xn=x1+awgn(x1,10)+x2+awgn(x2,10); Pxx1,f=pmem(xn,14,Nfft,Fs);subplot(4,1,3)plot(f,10*log10(Pxx1);xlabel('Frequency (Hz)');ylabel('Power Spectrum (dB)');title('MEM f2
6、/fs=0.3,Nfft=512,Oder=14');grid N=1024;Nfft=256;Fs=1000;n=0:N-1;t=n/Fs;x1=sin(2*pi*200*t);x2=sin(2*pi*300*t); %0.3xn=x1+awgn(x1,10)+x2+awgn(x2,10); Pxx1,f=pmem(xn,24,Nfft,Fs); %修改阶数为24subplot(4,1,4)plot(f,10*log10(Pxx1);xlabel('Frequency (Hz)');ylabel('Power Spectrum (dB)');title
7、('MEM f2/fs=0.3,Nfft=256,Oder=24');GridBurg法估计N=1024;Nfft=256;Fs=1000;n=0:N-1;t=n/Fs;x1=sin(2*pi*200*t);x2=sin(2*pi*300*t); %0.3xn=x1+awgn(x1,10)+x2+awgn(x2,10); Pxx1,f=pburg(xn,14,Nfft,Fs);subplot(4,1,1)plot(f,10*log10(Pxx1);xlabel('Frequency (Hz)');ylabel('Power Spectrum (dB)
8、39;);title('Burg f2/fs=300,Nfft=256, Oder=14');grid N=1024;Nfft=256;Fs=1000;n=0:N-1;t=n/Fs;x1=sin(2*pi*200*t);x2=sin(2*pi*250*t); %0.25xn=x1+awgn(x1,10)+x2+awgn(x2,10); Pxx1,f=pburg(xn,14,Nfft,Fs);subplot(4,1,2)plot(f,10*log10(Pxx1);xlabel('Frequency (Hz)');ylabel('Power Spectrum
9、 (dB)');title('Burg f2/fs=250,Nfft=256, Oder=14');grid N=1024; Nfft=512; %修改数据长度512Fs=1000;n=0:N-1;t=n/Fs;x1=sin(2*pi*200*t);x2=sin(2*pi*300*t); %0.3xn=x1+awgn(x1,10)+x2+awgn(x2,10); Pxx1,f=pburg(xn,14,Nfft,Fs);subplot(4,1,3)plot(f,10*log10(Pxx1);xlabel('Frequency (Hz)');ylabel(&
10、#39;Power Spectrum (dB)');title('Burg f2/fs=300,Nfft=512, Oder=14');grid N=1024;Nfft=256;Fs=1000;n=0:N-1;t=n/Fs;x1=sin(2*pi*200*t);x2=sin(2*pi*300*t); %0.3xn=x1+awgn(x1,10)+x2+awgn(x2,10); Pxx1,f=pburg(xn,24,Nfft,Fs); %修改阶数为24subplot(4,1,4)plot(f,10*log10(Pxx1);xlabel('Frequency (Hz)
11、');ylabel('Power Spectrum (dB)');title('Burg f2/fs=300,Nfft=256, Oder=24');gridLevinson递推法N=1024;Nfft=256;Fs=1000;n=0:N-1;t=n/Fs;x1=sin(2*pi*200*t);x2=sin(2*pi*300*t); %0.3xn=x1+awgn(x1,10)+x2+awgn(x2,10); Pxx1,f=pyulear(xn,14,Nfft,Fs);%Pxx1,f=Levinson(xn,14,Nfft,Fs);subplot(4,1,
12、1)plot(f,10*log10(Pxx1);xlabel('Frequency (Hz)');ylabel('Power Spectrum (dB)');title('Levinson Nfft=256,f2/fs=0.3,Oder=14');grid N=1024;Nfft=256;Fs=1000;n=0:N-1;t=n/Fs;x1=sin(2*pi*200*t);x2=sin(2*pi*250*t); %0.25xn=x1+awgn(x1,10)+x2+awgn(x2,10); Pxx1,f=pyulear(xn,14,Nfft,Fs);
13、subplot(4,1,2)plot(f,10*log10(Pxx1);xlabel('Frequency (Hz)');ylabel('Power Spectrum (dB)');title('Levinson Nfft=256,f2/fs=0.25,Oder=14');grid N=1024; Nfft=512; %修改数据长度512Fs=1000;n=0:N-1;t=n/Fs;x1=sin(2*pi*200*t);x2=sin(2*pi*300*t); %0.3xn=x1+awgn(x1,10)+x2+awgn(x2,10); Pxx1,f
14、=pyulear(xn,14,Nfft,Fs);subplot(4,1,3)plot(f,10*log10(Pxx1);xlabel('Frequency (Hz)');ylabel('Power Spectrum (dB)');title('Levinson Nfft=512,f2/fs=0.3,Oder=14');grid N=1024;Nfft=256;Fs=1000;n=0:N-1;t=n/Fs;x1=sin(2*pi*200*t);x2=sin(2*pi*300*t); %0.3xn=x1+awgn(x1,10)+x2+awgn(x2,
15、10); Pxx1,f=pyulear(xn,24,Nfft,Fs); %修改阶数为24subplot(4,1,4)plot(f,10*log10(Pxx1);xlabel('Frequency (Hz)');ylabel('Power Spectrum (dB)');title('Levinson Nfft=256,f2/fs=0.3,Oder=24');grid最大熵法改变信号的相位、频率、信噪比N=1024;Nfft=256;Fs=1000;n=0:N-1;t=n/Fs;x1=sin(2*pi*200*t);x2=sin(2*pi*300*
16、t); %0.3xn=x1+awgn(x1,10)+x2+awgn(x2,10); Pxx1,f=pmem(xn,14,Nfft,Fs);subplot(4,1,1)plot(f,10*log10(Pxx1);xlabel('Frequency (Hz)');ylabel('Power Spectrum (dB)');title('MEM f2/fs=0.3,Nfft=256,Oder=14');grid N=1024;Nfft=256;Fs=1000;n=0:N-1;t=n/Fs;x1=sin(2*pi*200*t+pi/6); %相位加了pi
17、/6x2=sin(2*pi*300*t); %0.3xn=x1+awgn(x1,10)+x2+awgn(x2,10); Pxx1,f=pmem(xn,14,Nfft,Fs);subplot(4,1,2)plot(f,10*log10(Pxx1);xlabel('Frequency (Hz)');ylabel('Power Spectrum (dB)');title('MEM f2/fs=0.3,Nfft=256,Oder=14,相位加pi/6');grid N=1024;Nfft=256;Fs=1000;n=0:N-1;t=n/Fs;x1=sin
18、(2*pi*200*t);x2=sin(2*pi*300*t); %0.3xn=x1+awgn(x1,5)+x2+awgn(x2,5); %性噪比改为5Pxx1,f=pmem(xn,14,Nfft,Fs);subplot(4,1,3)plot(f,10*log10(Pxx1);xlabel('Frequency (Hz)');ylabel('Power Spectrum (dB)');title('MEM f2/fs=0.3,Nfft=256,Oder=14,性噪比=5');grid N=1024;Nfft=256;Fs=1000;n=0:N-1
19、;t=n/Fs;x1=sin(2*pi*300*t);x2=sin(2*pi*400*t); %0.3xn=x1+awgn(x1,10)+x2+awgn(x2,10); Pxx1,f=pmem(xn,14,Nfft,Fs);subplot(4,1,4)plot(f,10*log10(Pxx1);xlabel('Frequency (Hz)');ylabel('Power Spectrum (dB)');title('MEM f1/fs=0.3,f2/fs=0.4,Nfft=256,Oder=14');gridBurg改变信号的相位、频率、信噪比
20、N=1024;Nfft=256;Fs=1000;n=0:N-1;t=n/Fs;x1=sin(2*pi*200*t);x2=sin(2*pi*300*t); %0.3xn=x1+awgn(x1,10)+x2+awgn(x2,10); Pxx1,f=pburg(xn,14,Nfft,Fs);subplot(4,1,1)plot(f,10*log10(Pxx1);xlabel('Frequency (Hz)');ylabel('Power Spectrum (dB)');title('Burg f2/fs=300,Nfft=256, Oder=14')
21、;grid N=1024;Nfft=256;Fs=1000;n=0:N-1;t=n/Fs;x1=sin(2*pi*200*t+pi/6); %相位加了pi/6x2=sin(2*pi*300*t); %0.3xn=x1+awgn(x1,10)+x2+awgn(x2,10); Pxx1,f=pburg(xn,14,Nfft,Fs);subplot(4,1,2)plot(f,10*log10(Pxx1);xlabel('Frequency (Hz)');ylabel('Power Spectrum (dB)');title('Burg f2/fs=300,Nf
22、ft=256, Oder=14,相位加pi/6');grid N=1024;Nfft=256;Fs=1000;n=0:N-1;t=n/Fs;x1=sin(2*pi*300*t);x2=sin(2*pi*400*t); %0.3xn=x1+awgn(x1,10)+x2+awgn(x2,10); Pxx1,f=pburg(xn,14,Nfft,Fs);subplot(4,1,3)plot(f,10*log10(Pxx1);xlabel('Frequency (Hz)');ylabel('Power Spectrum (dB)');title('Bur
23、g f1/fs=300,f2/fs=400,Nfft=256, Oder=14');grid N=1024;Nfft=256;Fs=1000;n=0:N-1;t=n/Fs;x1=sin(2*pi*200*t);x2=sin(2*pi*300*t); %0.3xn=x1+awgn(x1,5)+x2+awgn(x2,5); %性噪比改为5Pxx1,f=pburg(xn,14,Nfft,Fs);subplot(4,1,4)plot(f,10*log10(Pxx1);xlabel('Frequency (Hz)');ylabel('Power Spectrum (dB)
24、');title('Burg f2/fs=300,Nfft=256, Oder=14,性噪比=5');gridLevinson法改变信号的相位、频率、信噪比 N=1024;Nfft=256;Fs=1000;n=0:N-1;t=n/Fs;x1=sin(2*pi*200*t);x2=sin(2*pi*300*t); %0.3xn=x1+awgn(x1,10)+x2+awgn(x2,10); Pxx1,f=pyulear(xn,14,Nfft,Fs);%Pxx1,f=Levinson(xn,14,Nfft,Fs);subplot(4,1,1)plot(f,10*log10(P
25、xx1);xlabel('Frequency (Hz)');ylabel('Power Spectrum (dB)');title('Levinson Nfft=256,f2/fs=0.3,Oder=14');grid N=1024;Nfft=256;Fs=1000;n=0:N-1;t=n/Fs;x1=sin(2*pi*200*t+pi/6); %相位加了pi/6x2=sin(2*pi*300*t); %0.3xn=x1+awgn(x1,10)+x2+awgn(x2,10); Pxx1,f=pyulear(xn,14,Nfft,Fs);%Pxx1,f=Levinson(xn,14,Nfft,Fs);subplot(4,1,2)plot(f,10*log10(Pxx1);xl
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 自无证建房出售合同(2篇)
- 粘土冰淇淋课件
- 语文单据 课件
- 专题01 字音字形词语(考点串讲)-七年级语文上学期期末考点大串讲(统编版2024·五四学制)
- 第一讲 成长趣事(看图写话教学)-二年级语文上册(统编版)
- 第七讲 做手工(看图写话教学)-二年级语文上册(统编版)
- 西京学院《应用统计学》2021-2022学年第一学期期末试卷
- 西京学院《数字电子技术》2021-2022学年期末试卷
- 西京学院《机械工程测试技术》2021-2022学年第一学期期末试卷
- 冲突 作文 课件
- 机加工安全事故案例演示文稿
- 凯文杜兰特-英语介绍
- 剖宫产术后再次妊娠阴道分娩管理的专家共识
- 铁路工程挡土墙加固施工平安方案
- 最全的俄语教学课件
- 改进维持性血液透析患者贫血状况PDCA
- 再生资源回收利用体系建设项目方案
- 循证护理学(理论部分)智慧树知到答案章节测试2023年复旦大学
- 医院开展老年友善医疗机构建设工作总结
- 马克思主义基本原理概论智慧树知到答案章节测试2023年泰山学院
- 餐饮档口和门店消防安全培训
评论
0/150
提交评论