matlab 现代功率谱估计,matlab经典、现代功率谱估计

上机作业：

1、假设一平稳随机信号为()()()0.81x n x n w n =?+，其中是均值为0，方差为1的白噪声，数据长度为1024。

(1)、产生符合要求的)(n w 和)(n x ；

(2)、给出信号)(n x 的理想功率谱；

(3)、编写周期图谱估计函数，估计数据长度N=1024及256时信号功率谱，分析估计效果。

(4)、编写Bartlett 平均周期图函数，估计当数据长度N=1024及256时，分段数L 分别为2和8时信号的功率谱，分析估计效果。

一、解题思路

w(n)可以通过随机序列randn(1,N)来产生，x(n)可以通过对w(n)滤波产生(由递推式可

得系统的传递函数)，也可以直接由递推式迭代产生。

由于线性系统的输出功率谱等于输入功率谱乘以传递函数模的平方，X(n)可以看做w(n)通过一线性系统的输出，H(z)=1/(1-0.8z)。所以x(n)的理想功率谱P(e jw )=σw 2|H(e jw )|2。

周期图方法：直接对观测数据做FFT 变换，变换的结果取模的平方再除以数据长度，作为估计的功率谱。256个观测点时可以对原观测数据以4为间隔提取得到。

Bartlett 法：将L 组独立的观测数据分别求周期图，再将L 个周期图求平均作为信号的功率谱估计。L 组数据可以通过对原观测数据以L 为间隔提取得到。二、MATLAB 实现程序及注解 clc;

clear;close all;

Fs=500; %采样率

N=1024; %观测数据

w=sqrt(1)+randn(1,N); %0均值,方差为1的白噪声,长度1024

x=[w(1) zeros(1,N-1)]; %初始化x(n),长度1024,x(1)=w(1)

for i=2:N

x(i)=0.8*x(i-1)+w(i); %迭代产生观测数据x(n)

end

%% 理想功率谱

[h,w1]=freqz(x);

figure,plot(w1*500/(2*pi),10*log10(abs(h).^2));grid on;

title('理想功率谱');

xlabel('频率'); ylabel('功率db');

%% 周期图法

%1024个观测点

Pxx=abs(fft(x)).^2/N; %周期图公式

Pxx=10*log10(Pxx(index+1)); %化为db

figure;plot(k,Pxx);grid on;

title('周期图1024点');

matlab 现代功率谱估计,matlab经典、现代功率谱估计相关推荐

spss22.0统计分析从入门到精通_数据分析最全资料：SPSS/MATLAB/SQL/SAS/EXCEL经典教材+视频教程，快速入门！...
可以说当今社会数据就是力量.数据分析能力已经成为各个行业必备的技能,如果能熟练掌握数据分析技能,不管是留学申请还是以后的就业,都会是一个很大的加分项. 但是目前国内很少有高校专门开设数据分析课程,想要 ...
svm算法 matlab,SVM算法matlab程序
经典SVM算法matlab程序\svm\binomial.m 经典SVM算法matlab程序\svm\centrefig.m 经典SVM算法matlab程序\svm\cmap.mat 经典SVM算法m ...
matlab 摄像头拍照,matlab摄像头拍照
图像采集: 该单元是指道路上安装的摄像头在检测到有车辆通过的同时进行拍照并...在 MATLAB 中我们可以调用 im2gray 函数对图像进行灰度化处理. 4.2.2图像的...... 2 车辆牌 ...
【 MATLAB 】使用 MATLAB 得到高密度谱（补零得到DFT）和高分辨率谱（获得更多的数据得到DFT）的方式对比（附MATLAB脚本）
上篇博文分析了同一有限长序列在不同的N下的DFT之间的不同: MATLAB ]使用 MATLAB 作图讨论有限长序列的 N 点 DFT(强烈推荐)(含MATLAB脚本) 那篇博文中,我们通过补零的方式 ...
【 MATLAB 】使用 MATLAB 实现模拟信号的近似及其连续傅里叶变换
严格来说,利用MATLAB是不可能用来分析模拟信号的.然而,如果有足够小的时间增量在足够细的栅格上对采样而产生一种平滑的图,并有足够大的时间来展示所有的模式,那么就能对模拟信号作近似分析.令是栅格间隔 ...
【 MATLAB 】使用 MATLAB 求由差分方程表示的滤波器的响应的两种方法
例题: 一个3阶低通滤波器由下面差分方程描述: y(n) = 0.0181 x(n) + 0.0543 x(n-1) + 0.0543 x(n-2) + 0.0181 x(n-3) + 1.76 y( ...
matlab中fdyn,Matlab的用法总结
1. 对序列进行洗牌 randperm() randperm()产生随机的序列 %if filepaths 是一个5*1的结构体,then cshuffle = randperm(length(fil ...
matlab 数学库,matlab数学函数库
(n) 求 n 的阶乘如何用 matlab 配方没有发现 matlab 有这一命令,不过我们可以调用 maple 的命令,调用方法如下: 首先加载 maple 中的 student 函数库,加载. ...
matlab字符衔接,matlab字符串连接(多个字符串)的经验，matlab字符串
matlab字符串连接(多个字符串)的经验,matlab字符串参考 matlab中字符串连接的3种方法文章,实验遇到的一些小白问题,如有问题希望大神们给予指导. Matlab中,两个或多个字符串 ...
matlab模块 python,Matlab 和Python结合使用
今天介绍的内容比较基础,希望能达到抛砖引玉的效果. Matlab 和 Python 结合使用无非分为两种情形: Matlab中使用Python包: Python中使用Matlab. 1. Matlab ...

matlab 现代功率谱估计,matlab经典、现代功率谱估计

matlab 现代功率谱估计,matlab经典、现代功率谱估计相关推荐

最新文章

热门文章