当前位置:
文档之家› MATLAB仿真实现经典谱估计(采用周期图法)
MATLAB仿真实现经典谱估计(采用周期图法)
(3)了解周期图法(直接法)估计功率谱的缺陷。
三、实验原理
把随机信号 的N点观察数据 视为一能量有限信号,直接取得 傅里叶变换,得 ,然后再取其幅值的平方,并除以N,作为对 真实的功率谱 的估计。即为:
四、涉及实验的相关情况介绍(包含使用软件或实验设备等情况)
一台安装MATLAB软件的电脑
五、实验记录(程序、相关的图形、相关数据记录及分析)
axis([0 0.3 -2 2]);
grid on;wiFra bibliotekdow = boxcar( length( xn) ) ;%矩形窗
nfft = 512;
[Pxx f]= periodogram( xn,window,nfft,Fs) ;%直接法
subplot(312)
plot( f,10* log10( Pxx) ) ;
%采用直接法(周期图法)估计功率谱;
clear
Fs = 1000;%采样频率
n = 0:1 /Fs: .3;%产生含有噪声的序列
xn = cos(200*pi*n)+0.1*randn(size(n));
subplot(311);%输出随机信号xn;
plot(n,xn);xlabel('时间');ylabel('幅度');title('输入信号x(n)');
plot( f,10*log10( Pxx) ) ;
title('直接法经典谱估计,1024点');
xlabel('频率(Hz)');
ylabel('功率谱密度');
grid on;
六、实验总结
从上图我们可以得到这样的结论:在增加数据长度N时,就会使互不相关的点数增加,提高谱曲线的分辨力,但是加剧谱曲线 的起伏。经典功率谱估计不是一致估计,这是周期图法(直接法)的一个严重的缺点。
title('直接法经典谱估计,512点');
xlabel('频率(Hz)');
ylabel('功率谱密度');
grid on;
window = boxcar( length( xn) );
nfft = 1024;
[Pxx,f]= periodogram( xn,window,nfft,Fs) ;
subplot(313) ;
数字信号处理课程实验报告
实验指导教师:黄启宏
实验名称
MATLAB仿真实现经典谱估计(采用周期图法)
专业、班级
电子与通信工程
姓名
张帅
实验地点
仿古楼301
实验日期
2013.11.17
一、实验内容
采用周期图法(直接法)实现经典谱估计。
二、实验目的
(1)掌握周期图法(直接法)估计出功率谱的步骤和方法;
(2)在实验的过程中找到影响经典谱估计的因素;