功率谱分析(Matlab)_第1页
功率谱分析(Matlab)_第2页
功率谱分析(Matlab)_第3页
功率谱分析(Matlab)_第4页
功率谱分析(Matlab)_第5页
已阅读5页,还剩5页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

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. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论