下载文档原格式
《随机信号分析》实验报告二班级_______学号______姓名_______实验二高斯噪声的产生和性能测试１．实验目的（1）掌握加入高斯噪声的随机混合信号的分析方法。
（2）研究随机过程的均值、相关函数、协方差函数和方差。
⒉实验原理（1）利用随机过程的积分统计特性，给出随机过程的均值、相关函数、协方差函数和方差。
（2）随机信号均值、方差、相关函数的计算公式，以及相应的图形。
⒊实验报告要求（1）简述实验目的及实验原理。
（2）采用幅度为1，频率为25HZ的正弦信号为原信号，在其中加入均值为2，方差为0.04的高斯噪声得到混合随机信号X(t)。
试求随机过程的均值、相关函数、协方差函数和方差。
用MATLAB进行仿真，给出测试的随机过程的均值、相关函数、协方差函数和方差图形，与计算的结果作比较，并加以解释。
（3）分别给出原信号与混合信号的概率密度和概率分布曲线，并以图形形式分别给出原信号与混合信号均值、方差、相关函数的对比。
（4）读入任意一幅彩色图像，在该图像中加入均值为0，方差为0.01的高斯噪声，请给出加噪声前、后的图像。
（5）读入一副wav格式的音频文件，在该音频中加入均值为2，方差为0.04的高斯噪声，得到混合随机信号X(t)，请给出混合信号X(t)的均值、相关函数、协方差函数和方差，频谱及功率谱密度图形。
4、源程序及功能注释(逐句注释)（1）：clear all;clc;t=0:320;x=sin(2*pi*t*25);x1=wgn(1,321,0);z=x+x1;y=trapz(t,z);%y=int(z,x,0,t);subplot(3,2,1),plot(z);title('随机信号序列')meany=mean(z);subplot(3,2,3),plot(t,meany,'.');title('随机信号均值')vary=var(y); %方差subplot(3,2,4),plot(t,vary,'.');title('随机信号方差')cory=xcorr(z,'unbiased');%自相关函数subplot(3,2,2),plot(cory);title('随机信号自相关函数')covv=cov(y);subplot(3,2,5),plot(t,covv,'.');title('随机信号协方差')（2）：t=[0:0.0005:0.045];X1=sin(2*pi*25*t);%正弦subplot(3,4,1);plot(t,X1);gridtitle('正弦函数序列');X2=randn(1,length(t)); %产生均值为0，方差σ^2=1，标准差σ=1的正态分布的随机数或矩阵的函数高斯随机信号%X2=normrnd(2,0.04); %高斯随机序列均值,标准差subplot(3,4,2);plot(t,X2);title('高斯噪声序列');X=X1+X2; %混合随机信号X(t)subplot(3,4,3);plot(t,X);gridtitle('混合随机信号');meany1=mean(X1); %原信号的均值subplot(3,4,6),plot(t,meany1);title('原信号均值');vary1=var(X1); %原信号的方差subplot(3,4,7),plot(t,vary1);title('原信号方差');cory1=xcorr(X1,'unbiased'); %原信号的自相关函数subplot(3,4,8),plot(cory1);title('原信号自相关函数');meany=mean(X); %混合信号的均值subplot(3,4,10),plot(t,meany);title('混合信号均值');vary=var(X); %混合信号的方差subplot(3,4,11),plot(t,vary);title('混合信号方差')cory=xcorr(X,'unbiased'); %混合信号的自相关函数subplot(3,4,12),plot(cory);title('混合信号自相关函数')covy=cov(X1,X); %协方差subplot(3,4,4),plot(covy);title('协方差');[f1,xi]=ksdensity(X1); %原信号的概率密度subplot(3,4,5);plot(xi,f1);title('原信号的概率密度分布)');[f2,xi]=ksdensity(X); %混合信号的概率密度subplot(3,4,9);plot(xi,f2);title('混合信号概率密度分布');（3）：clcclear allclose allA = imread('dadian.jpg'); % 读入图像V=0.01;Noisy=imnoise(A,'gaussian',0,V);subplot(1,2,1),imshow(A),title('原图像');subplot(1,2,2),imshow(Noisy),title('加噪后图像'); （4）：clcclear allclose allt=0:320;A = wavread('alert.wav'); % 读入音频x = double(A);y=awgn(x,2,0.04);%x1 = double(z);%y=x+x1;subplot(2,3,1),plot(y);title('随机信号序列')meany=mean(y);subplot(2,3,2),plot(t,meany,'.');title('随机信号均值')vary=var(y); %方差subplot(2,3,3),plot(t,vary,'.');title('随机信号方差')cory=xcorr(y,'unbiased');%自相关函数subplot(2,3,4),plot(cory);title('随机信号自相关函数')fy=fft(y);ym=abs(fy);subplot(2,3,5),plot(ym);title('随机信号频谱图')fz=fft(cory);zm=abs(fz);subplot(2,3,6),plot(zm);title('随机信号功率谱密度图')5. 实验总结（手写）可给出实验过程中遇到的问题、解决方法、自己的收获、可否有改进办法等。
随机信号分析实验报告引言：随机信号是指信号在时间或空间上的其中一种特性是不确定的，不能准确地预测其未来行为的一类信号。
随机信号是一种具有随机性的信号，其值在一段时间内可能是不确定的，但是可以通过概率论和统计学的方法来描述和分析。
实验目的：通过实验，学习了解随机信号的基本概念和特性，学习了解和掌握常见的随机信号分析方法。
实验原理：随机信号可以分为离散随机信号和连续随机信号。
离散随机信号是信号在离散时间点上，在该时间点上具有一定的随机性；而连续随机信号是信号在连续时间上具有随机性。
常见的随机信号分析方法包括概率密度函数、功率谱密度函数等。
实验器材：计算机、MATLAB软件、随机信号产生器、示波器、电缆、电阻等。
实验步骤：1.配置实验仪器：将随机信号产生器和示波器与计算机连接。
2.生成随机信号：调节随机信号产生器的参数，产生所需的随机信号。
3.采集数据：使用示波器采集随机信号的样本数据，并将数据导入MATLAB软件。
4.绘制直方图：使用MATLAB软件绘制样本数据的直方图，并计算概率密度函数。
5.计算统计特性：计算随机信号的均值、方差等统计特性。
6.绘制功率谱密度函数：使用MATLAB软件绘制随机信号的功率谱密度函数。
实验结果和讨论：我们采集了一段长度为N的随机信号样本数据，并进行了相应的分析。
通过绘制直方图和计算概率密度函数，我们可以看出随机信号的概率分布情况。
通过计算统计特性，我们可以得到随机信号的均值、方差等重要参数。
通过绘制功率谱密度函数，我们可以分析随机信号的频谱特性。
结论：本实验通过对随机信号的分析，加深了对随机信号的理解。
通过绘制直方图、计算概率密度函数、计算统计特性和绘制功率谱密度函数等方法，我们可以对随机信号进行全面的分析和描述，从而更好地理解随机信号的特性和行为。
2.王五,赵六.随机信号分析方法.物理学报,2024,30(2):120-130.。
实验2 随机过程的计算机模拟一、实验目的1、给定功率谱（相关函数）和概率分布，通过计算机模拟分析产生相应的随机过程；2、通过该随即过程的实际功率谱（相关函数）和概率分布验证该实验的有效性；3、学会运用Matlab 函数对随机过程进行模拟。
二、实验原理1、标准正态分布随机序列的产生方法：利用随机变量函数变换的方法。
设r1，r2为两个相互独立的（0，1）均匀分布的随机数，如果要产生服从均值为m,方差为正态分布的随机数x，则可以按如下变换关系产生：2、正态随机矢量的模拟：设有一 N 维正态随机矢量，其概率密度为为协方差矩阵，且是正定的。
3、具有有理谱的正态随机序列的模拟根据随机过程通过线性系统的理论，白噪声通过线性系统后，输出是正态的，且输出功率谱只与系统的传递函数有关。
利用这一性质，我们可以产生正态随机过程。
如上图所示，输入W(n)为白噪声，假定功率谱密度为G (z) = 1 W ，通过离散线性系统后，输出X (n)是正态随机序列，由于要求模拟的随机序列具有有理谱，则G (z) X 可表示为：其中，G (z) X+ 表示有理谱部分，即所有的零极点在单位圆之内，G (z) X? 表示非有理谱部分，即所有零极点在单位圆之外。
4、满足一定相关函数的平稳正态随机过程的模拟，当已知平稳随机过程的相关函数而要确定该随机过程的模拟算法。
很显然，只要我们设计一个合适的滤波器，使得该白噪声通过滤波器后，输出的功率谱满足上述相关函数的傅里叶变换，就可以模拟得到该随机过程。
三、实验内容1、产生两组相互独立的（0，1）均匀分布的随机数（随机数个数：500）程序及图形如下:clear;x=randn(1,500);y= randn(1,500);subplot(2,1,1);plot(x);title('第二组');subplot(2,1,2);plot(y);title('第一组')2、按照实验原理中的方法产生一组均值为1，方差为1 的正态分布的随机序列（序列长度：500）程序及图形如下：clear;y=1+sqrt(1)*randn(1,500);plot(y);title(‘正态分布，均值方差都为1’)3、画出功率谱密度为G(w)=1/(1.25+cosw) 的功率谱图（一个周期内），采用均匀采样方法，采样点数为500程序及图形如下：clear;w=rand(1,500);M=1.25+cos(w);N=1;G=N./M;plot(G);title('均匀采样功率频谱');5、模拟产生一个功率谱为G(w)=1/(1.25+cosw) 的正态随机序列程序及图形如下:clear;w=randn(1,500);M=1.25+cos(w);N=1;G=N./M;plot(G);title('均匀采样功率频谱');4、实验中所遇到问题及解决方法问题1、对Matlab软件很生疏、编程也不熟悉。
随机信号实验报告学院通信工程学院专业信息工程班级 1401051班制作人文杰制作人晓鹏一、 摘要根据实验的要求与具体容，我们组做了一下分工，XXX 完成了本次的第一组实验即基于MATLAB 的信号通过线性系统与非线性系统的特性分析，具体容有（高斯白噪声n ，输入信号x ，通过线性与非线性系统的信号a,b,y1,y2的均值，均方值，方差，自相关函数，概率密度，功率谱密度以及频谱并把它们用波形表示出来），XXX 和XXX 两人合力完成了基于QUARTUS II 的2ASK 信号的产生及测试实验具体容有（XXX 负责M 序列发生器以及分频器，XXX 负责载波的产生以及开关函数和管脚配置），最后的示波器调试以及观察2ASK 信号的FFT 变换波形由我们组所有人一起完成的。
二、实验原理及要求实验一、信号通过线性系统与非线性系统的特性分析1、实验原理① 随机过程的均值（数学期望）：均值表示集合平均值或数学期望值。
基于随机过程的各态历经性，可用时间间隔T 的幅值平均值表示，即：均值表达了信号变化的中心趋势，或称之为直流分量。
② 随机过程的均方值：信号x(t)的均方值，或称为平均功率，其表达式为：均方值表达了信号的强度，其正平方根值，又称为有效值，也是信号的平均能量的一种表达。
③ 随机信号的方差： 信号x(t)的方差定义为：描述了信号的静态量，方差反映了信号绕均值的波动程度。
在已知均值和均方值的前提下，方差就很容易求得了。
④随机信号的自相关函数信号的相关性是指客观事物变化量之间的相依关系。
对于平稳随机过程X(t)和Y(t)在两个不同时刻t和t+τ的起伏值的关联程度，可以用相关函数表示。
在离散情况下，信号x(n)和y(n)的相关函数定义为：τ,t=0,1，2,……N-1。
⑤随机过程的频谱：信号频谱分析是采用傅立叶变换将时域信号x(t)从另一个角度来了解信号的特征。
时域信号x(t)的傅氏变换为：⑥随机过程的功率谱密度：随机信号的功率普密度是随机信号的各个样本在单位频带的频谱分量消耗在一欧姆电阻上的平均功率的统计均值，是从频域描述随机信号的平均统计参量，表示X(t)的平均功率在频域上的分布。
实验一 随机序列的产生及数字特征估计一、实验目的1、学习和掌握随机数的产生方法；2、实现随机序列的数字特征估计。
二、实验原理1. 随机数的产生随机数指的是各种不同分布随机变量的抽样序列（样本值序列）。
进行随机信号仿真分析时，需要模拟产生各种分布的随机数。
在计算机仿真时，通常利用数学方法产生随机数，这种随机数称为伪随机数。
伪随机数是按照一定的计算公式产生的，这个公式称为随机数发生器。
伪随机数本质上不是随机的，而且存在周期性，但是如果计算公式选择适当，所产生的数据看似随机的，与真正的随机数具有相近的统计特性，可以作为随机数使用。
(0，1)均匀分布随机数是最最基本、最简单的随机数。
(0，1)均匀分布指的是在[0，1]区间上的均匀分布，即U(0，1)。
实际应用中有许多现成的随机数发生器可以用于产生(0，1)均匀分布随机数，通常采用的方法为线性同余法，公式如下：Ny x N ky Mod y y n n n n /))((110===-， （1.1）序列{}n x 为产生的(0，1)均匀分布随机数。
下面给出了上式的3组常用参数： (1) 7101057k 10⨯≈==，周期，N ；(2) （IBM 随机数发生器）8163110532k 2⨯≈+==，周期，N ； (3) （ran0）95311027k 12⨯≈=-=，周期，N ；由均匀分布随机数，可以利用反函数构造出任意分布的随机数。
定理1.1 若随机变量X 具有连续分布函数F X (x)，而R 为(0，1)均匀分布随机变量，则有)(1R F X x -= （1.2）由这一定理可知，分布函数为F X (x)的随机数可以由(0，1)均匀分布随机数按上式进行变换得到。
2. MATLAB 中产生随机序列的函数(1) (0，1)均匀分布的随机序列 函数：rand用法：x = rand(m,n)功能：产生m ×n 的均匀分布随机数矩阵。
(2) 正态分布的随机序列 函数：randn用法：x = randn(m,n)功能：产生m ×n 的标准正态分布随机数矩阵。
一、实验目的1. 了解功率谱估计的基本原理和方法。
2. 掌握使用MATLAB进行功率谱估计的步骤。
3. 通过实验验证不同功率谱估计方法的性能。
二、实验原理功率谱估计是信号处理中的一个重要分支，它能够揭示信号在频域中的特性。
功率谱估计的基本原理是将信号通过傅里叶变换转换为频域信号，然后对频域信号进行功率计算，得到功率谱。
功率谱反映了信号在不同频率上的能量分布，对于信号分析、系统设计等具有重要意义。
常用的功率谱估计方法有周期图法、Welch方法、Bartlett方法等。
本实验主要研究周期图法、Welch方法和Bartlett方法的性能。
三、实验内容1. 数据采集：使用MATLAB自带的信号生成函数生成一段模拟信号，作为实验数据。
2. 周期图法：计算信号的功率谱，并与理论功率谱进行比较。
3. Welch方法：计算信号的功率谱，并与周期图法的结果进行比较。
4. Bartlett方法：计算信号的功率谱，并与Welch方法和周期图法的结果进行比较。
四、实验步骤1. 生成模拟信号：使用MATLAB自带的信号生成函数生成一段模拟信号，如正弦波、方波等。
2. 周期图法：a. 计算信号的自相关函数；b. 对自相关函数进行傅里叶变换，得到功率谱；c. 将计算得到的功率谱与理论功率谱进行比较。
3. Welch方法：a. 将信号分割成多个短时段；b. 对每个短时段进行自相关函数计算；c. 对所有短时段的自相关函数进行加权平均；d. 对加权平均后的自相关函数进行傅里叶变换，得到功率谱；e. 将计算得到的功率谱与周期图法的结果进行比较。
4. Bartlett方法：a. 将信号分割成多个短时段；b. 对每个短时段进行自相关函数计算；c. 对每个短时段的自相关函数进行 Bartlett 平滑；d. 对平滑后的自相关函数进行傅里叶变换，得到功率谱；e. 将计算得到的功率谱与Welch方法和周期图法的结果进行比较。
五、实验结果与分析1. 周期图法：通过实验发现，周期图法计算得到的功率谱在低频段存在较大的噪声，而在高频段噪声较小。
随机信号分析实验报告实验一：平稳随机过程的数字特征实验二：平稳随机过程的谱分析实验三：随机信号通过线性系统的分析实验四：平稳时间序列模型预测班级：姓名：学号：一、实验目的1、加深理解平稳随机过程数字特征的概念2、掌握平稳随机序列期望、自相关序列的求解3、分析平稳随机过程数字特征的特点二、实验原理平稳随机过程数字特征求解的相关原理三、实验过程function y = experimentnumber = 49; %学号49I = 8; %幅值为8u = 1/number;Ex = I*0.5 + (-I)*0.5;N = 64;C0 = 1; %计数p(1) = exp(-u);for m = 2:Nk = 1:m/2;p(m) = exp(-u*m) + sum((u*m).^(2*k)./factorial(2*k)*exp(-u*m));2222()[()()]{()()}{()()}X R m E X n X n m I P X n X n m I I P X n X n m I =+=+=-+=-E[X(n)]= I P{X(n)=+I}+(-I)P{X(n)=-I}=0⨯⨯0m >当时，/2220(){()()}(2)!m k m k m P X n X n m I e P k λλ⎢⎥⎣⎦-=+===∑222()(1)(21)X R m I P I P I P =--=-2()()X X X C m R m m =-me I m n X n X E m R λ22)]()([)(-=+=end;pp = [fliplr(p) C0 p];Rx = (2*pp - 1)*I^2;m = -N:N;Kx = Rx - Ex*Ex;rx = Kx/25;subplot(211), plot(m,Rx); axis([-N N 0 I*I]); title('自相关序列');subplot(212), plot(m,rx); axis([-N N 0 1]); title('自相关序数');四、实验结果及分析自相关序列的特点分析：m>0时Rx(m)随着m的增大而减小，m<0时Rx(m)随着m的增大而增大。
大连理工大学实验预习报告学院（系）：信息与通信工程学院 专业： 电子信息工程 班级： 电子1401 姓 名： ****** 学号： ****** 实验时间： 2016.11.4 实验室： c221 指导教师： 郭 成 安实验II ：系统对随机信号响应的统计特性分析、功率谱分析及应用实验一、 实验目的和要求掌握直接法估计随机信号功率谱的原理和实现方法；掌握间接法估计随机信号功率谱的原理和实现方法；掌握系统对随机信号响应的统计特性分析及仿真实现方法。
熟悉MATLAB 信号处理软件包的使用。
二、 实验原理和内容（一）实验原理：1. 直接法估计随机信号功率谱原理直接法又称为周期图法，它是把随机信号 x(n)的N 点观察数据xN(n)视为一能量有限信号， 直接取 xN(n)的傅里叶变换，得到 XN(ej ω)，然后取其模值的平方，并除以 N ，作为对 x(n)真实 的功率谱 P(ej ω)的估计。
工程上，常使用离散 Fourier变换（DFT ，编程上使用其快速算法 FFT ），即 PX(k)=2|)k (|1N X N进行计算。
2. 间接法估计随机信号功率谱间接法的理论基础是 Wiener-Khintchine 定理，具体的实现方法是先由 xN(n)估计出自相关函数(m)rˆ，然后对(m)r ˆ求傅里叶变换得到 xN(n)的功率谱，记之为 XN(ej ω)，并以此作为对真实功率谱 P(ej ω)的估计。
工程上，常使用离散 Fourier 变换（DFT ，编程上使用其快速算法FFT ），即122)(ˆ)(+--=∑=M km j M M m X e m r k P π，1||-≤N M ，进行计算。
因为由这种方法求出的功率谱是通过自相关函数间接得到的，所以又称为间接法或 Blackman-Tuckey(BT)法，该方法是 FFT 出现之前 常用的谱估计方法。
3. 时域中系统对随机信号响应的统计特性分析及仿真根据系统卷积性质，计算系统输出信号的统计特性。
有如下性质：∑=n X Y n h m )(m )()()()(k h j h k j m R m R j kx Y ++=∑∑4. 频域中系统对随机信号响应的统计特性分析及仿真根据卷积定理，输入、输出信号功率谱的关系为 RY(ej ω) = RX(ej ω)|H(ej ω)|2。
在计算系统输 出信号功率谱时，如果在时域时计算困难，可以按照上式在频域计算。
（二）实验内容：1． 直接法估计随机信号功率谱(1) 生成1024点数据的随机信号)()2cos(5)2cos(2)n (2211n N f f X ++++=ϕπϕπ其中f1=30Hz ，f2=100Hz ，21ϕϕ，为在[0,2п]内的均匀分布的随机变量，N(n)是数学期望为0，方差为1 的高斯白噪声。
(2) 用周期图法计算的功率谱，并绘图。
(3) 用MATLAB 函数periodogram 重新计算的功率谱，并与(2)做比较。
2．间接法估计随机信号功率谱(1) 计算以上的自相关函数。
(2) 通过计算自相关函数的Fourier 变换，求的功率谱并绘图。
(3) 利用MATLAB 函数psd 、pwelch 重新计算的功率谱，并与(2)做比较。
3．系统对随机信号响应的统计特性分析及仿真(1) 生成含500点数据的高斯分布白噪声随机信号。
(2) 设计一个带通系统，其上、下截止频率分别为4KHz 和3KHz 。
(3) 计算通过以上带通滤波器的自相关函数和功率谱密度。
三、 实验步骤(1) 生成1024点随机信号X （n ），高斯白噪声N （n ）。
(2) 周期法绘制其功率谱。
(3) 利用MATLAB 函数periodogram 计算功率谱，与（2）比较。
(4) 计算其自相关函数。
(5) 计算自相关函数的傅里叶变换。
(6) 求其功率谱。
(7) 利用MATLAB 函数psd 、pwelch 重新计算的功率谱，并与(2)做比较。
(8) 生成含500点数据的高斯分布白噪声随机信号。
(9) 设计一个带通系统，其上、下截止频率分别为4KHz 和3KHz 。
(10) 计算通过以上带通滤波器的自相关函数和功率谱密度。
大连理工大学实验报告学院（系）： 信息与通信工程学院 专业： 电 子 信 息 工 程 班级： 1401 姓 名： ****** 学号： ****** 实验时间： 2016/11/04 实验室： C221 指导教师： 郭 成 安实验II ：系统对随机信号响应的统计特性分析、功率谱分析及应用实验一、 实验目的和要求掌握直接法估计随机信号功率谱的原理和实现方法；掌握间接法估计随机信号功率谱的原理和实现方法；掌握系统对随机信号响应的统计特性分析及仿真实现方法。
熟悉MATLAB 信号处理软件包的使用。
二、 实验原理和内容（一）实验原理：1. 直接法估计随机信号功率谱原理直接法又称为周期图法，它是把随机信号 x(n)的N 点观察数据xN(n)视为一能量有限信号， 直接取 xN(n)的傅里叶变换，得到 XN(ej ω)，然后取其模值的平方，并除以 N ，作为对 x(n)真实 的功率谱 P(ej ω)的估计。
工程上，常使用离散 Fourier变换（DFT ，编程上使用其快速算法 FFT ），即 PX(k)=2|)k (|1N X N进行计算。
2. 间接法估计随机信号功率谱间接法的理论基础是 Wiener-Khintchine 定理，具体的实现方法是先由 xN(n)估计出自相关函数(m)rˆ，然后对(m)r ˆ求傅里叶变换得到 xN(n)的功率谱，记之为 XN(ej ω)，并以此作为对真实功率谱 P(ej ω)的估计。
工程上，常使用离散 Fourier 变换（DFT ，编程上使用其快速算法FFT ），即122)(ˆ)(+--=∑=M km j M M m X e m r k P π，1||-≤N M ，进行计算。
因为由这种方法求出的功率谱是通过自相关函数间接得到的，所以又称为间接法或 Blackman-Tuckey(BT)法，该方法是 FFT 出现之前 常用的谱估计方法。
3. 时域中系统对随机信号响应的统计特性分析及仿真根据系统卷积性质，计算系统输出信号的统计特性。
有如下性质：∑=n X Y n h m )(m )()()()(k h j h k j m R m R j kx Y ++=∑∑5. 频域中系统对随机信号响应的统计特性分析及仿真根据卷积定理，输入、输出信号功率谱的关系为 RY(ej ω) = RX(ej ω)|H(ej ω)|2。
在计算系统输 出信号功率谱时，如果在时域时计算困难，可以按照上式在频域计算。
（二）实验内容：2． 直接法估计随机信号功率谱(1) 生成1024点数据的随机信号)()2cos(5)2cos(2)n (2211n N f f X ++++=ϕπϕπ其中f1=30Hz ，f2=100Hz ，21ϕϕ，为在[0,2п]内的均匀分布的随机变量，N(n)是数学期望为0，方差为1 的高斯白噪声。
(2) 用周期图法计算的功率谱，并绘图。
(3) 用MATLAB 函数periodogram 重新计算的功率谱，并与(2)做比较。
2．间接法估计随机信号功率谱(1) 计算以上的自相关函数。
(2) 通过计算自相关函数的Fourier 变换，求的功率谱并绘图。
(3) 利用MATLAB 函数psd 、pwelch 重新计算的功率谱，并与(2)做比较。
3．系统对随机信号响应的统计特性分析及仿真(1) 生成含500点数据的高斯分布白噪声随机信号。
(2) 设计一个带通系统，其上、下截止频率分别为4KHz 和3KHz 。
(3) 计算通过以上带通滤波器的自相关函数和功率谱密度。
三、 实验步骤(1) 生成1024点随机信号X （n ），高斯白噪声N （n ）。
(2) 周期法绘制其功率谱。
(3) 利用MATLAB 函数periodogram 计算功率谱，与（2）比较。
(4) 计算其自相关函数。
(5) 计算自相关函数的傅里叶变换。
(6) 求其功率谱。
(7) 利用MATLAB 函数psd 、pwelch 重新计算的功率谱，并与(2)做比较。
(8) 生成含500点数据的高斯分布白噪声随机信号。
(9) 设计一个带通系统，其上、下截止频率分别为4KHz 和3KHz 。
(10) 计算通过以上带通滤波器的自相关函数和功率谱密度。
四、 主要仪器设备微型计算机、Matlab 开发环境（本报告采用MATLAB2016a ）五、编程代码与操作方法1． 实验内容（1）：直接法估计随机信号功率谱(1) 生成1024点数据的随机信号)()2cos(5)2cos(2)n (2211n N f f X ++++=ϕπϕπ其中f1=30Hz ，f2=100Hz ，21ϕϕ，为在[0,2п]内的均匀分布的随机变量，N(n)是数学期望为0，方差为1 的高斯白噪声。
(2) 用周期图法计算的功率谱，并绘图。
(3) 用MATLAB 函数periodogram 重新计算的功率谱，并与(2)做比较。
代码如下：clc;clear; %清屏及清除缓存N=1024;f1=30;f2=100;fs=1000;t=(0:N-1)/fs;%第一问Nn=randn(1,N); %高斯白噪声fai=random('unif',0,1,1,2)*2*pi; %均匀随机变量xn=2*cos(2*pi*f1*t+fai(1))+5*cos(2*pi*f2*t+fai(2))+Nn;Xk=fft(xn);Sk=(abs(Xk).^2)/N; %随机信号分析书第二版上式2.6-28 %第二问f=(0:N/2-1)*fs/N;w=10*log10(Sk(1:N/2));figure,plot(f,w);title('周期法估算随机信号功率谱');xlabel('f/Hz'),ylabel('Sx(f)(dB/Hz)');%第三问Sk=periodogram(xn);figure,plot(f,w);title('periodogram 函数算随机信号功率谱');xlabel('f/Hz'),ylabel('Sx(f)(dB/Hz)');figure,subplot(2,1,1),plot(f,w);title('直接法估算随机信号功率谱');xlabel('f/Hz'),ylabel('Sx(f)(dB/Hz)');subplot(2,1,2),plot(f,w,'r'); %延迟与原随机信号同一张图中比较title('periodogram函数算随机信号功率谱')xlabel('f/Hz'),ylabel('Sx(f)(dB/Hz)');实验内容(2)：间接法估计随机信号功率谱(1) 计算以上的自相关函数。
