matlab数信号处理周期图程序.docVIP

1、假设一平稳随机信号为，其中是均值为0，方差为1的白噪声，数据长度为1024。 （1）、产生符合要求的和； （2）、给出信号x(n)的理想功率谱； （3）、编写周期图谱估计函数，估计数据长度N=1024及256时信号功率谱，分析估计效果。 （4）、编写Bartlett平均周期图函数，估计当数据长度N=1024及256时，分段数L分别为2和8时信号的功率谱，分析估计效果。 解:思路 在matlab中提供的有 randn（m.n）函数，其为均值为零，方差为1的函数，所以w(n)可以通过随机序列randn(1,N)来产生，x(n)可以通过对 w(n)滤波产生，也可以直接由递推式迭代产生。 由于线性系统的输出功率谱等于输入功率谱乘以传递函数模的平方，X(n)可以看做w(n) 通过一线性系统的输出， H(z)=1/(1-0.8z) 所以x(n)的理想功率谱 直接产生 Matlab 程序： clear;close all; F=500; %采样率 N=1024; %观测数据 subplot(2,1,1); w=sqrt(1)+randn(1,N); plot(w); xlabel（观察次数）；ylabel(功率db)；title(白噪声w的分布情况); subplot(2,1,2); x=[w(1) zeros(1,N-1)]; %初始化x