窄带信号包络及相位检波.docxVIP

随机实验窄带信号及包络和相位检波分析孙博彪2012年11月23日窄带信号及包络和相位检波分析I：摘要当窄带系统(接收机)的输入噪声(如热噪声)的功率谱分布在足够宽的频带(相对于接收机带宽)上时，系统的输出即为窄带过程。对于一个窄带信号，通过包络检波器之后，在检波器的输出端可以得到包络A（t）。当窄带信号通过一个相位检波器之后，可以得到有关相位的信息。论文通过用matlab软件对窄带信号的包络和相位，以及窄带信号的数字特征、概率密度、功率谱密度等进行了画图分析，进一步研究了窄带信号的包络和相位的特性。II：实验背景与目的通过实验掌握窄带随机信号的特点，关键在于包络和相位检波分析。分析并了解了解窄带信号的特性，包括均值、均方值、方差、相关函数、概率密度、频谱及功率谱密度等。熟悉运用常用的信号处理仿真软件平台：matlab软件。III：窄带信号的产生与验证一、窄带随机过程一个实平稳随机过程X(t)，若它的功率谱密度具有下述性质：中心频率为ωc，带宽为△ω=2ω0，当△ωωc时，就可认为满足窄带条件。若随机过程的功率谱满足该条件则称为窄带随机过程。若带通滤波器的传输函数满足该条件则称为窄带滤波器。随机过程通过窄带滤波器传输之后变成窄带随机过程。图1 为典型窄带随机过程的功率谱密度图。若用一示波器来观测次波形，则可看到，它接近于一个正弦波，但此正弦波的幅度和相位都在缓慢地随机变化，图2所示为窄带随机过程的一个样本函数。图1 典型窄带随机过程的功率谱密度图图2 窄带随机过程的一个样本函数二、窄带随机过程的数学表示1、用包络和相位的变化表示由窄带条件可知,窄带过程是功率谱限制在ωc附近的很窄范围内的一个随机过程,从示波器观察(或由理论上可以推知):这个过程中的一个样本函数(一个实现)的波形是一个频率为?c且幅度和相位都做缓慢变化的余弦波。写成包络函数和随机相位函数的形式：X(t)=A(t)*cos[ωct+ Φ(t)]其中:A(t)称作X(t)的包络函数; Φ(t)称作X(t)的随机相位函数。包络随时间做缓慢变化，看起来比较直观，相位的变化，则看不出来。2、莱斯（Rice）表示式任何一个实平稳随机过程X(t)都可以表示为：X(t)=Ac(t) cosωct-AS(t) sinωct其中同相分量：Ac(t)= X(t) cosφt= X(t) cosωct＋sinωct=LP[X(t) *2cosωct]正交分量：AS(t) = X(t)sinφt= cosωct— X(t) sinωct= LP[-X(t) *2sinωct]（LP[A]表示取A的低频部分）。Ac(t)和AS(t)都是实随机过程，均值为0，方差等于X(t)的方差。三、窄带随机过程仿真建模要求1、根据要求得到X(t)的表达式：x= (l+a) .*cos (2*pi*4000*t+2) +noisy/10;其中：noisy为高斯白噪声，由wgn函数生成，a=cos (2*pi*l000*t)，均值：Ex=mean (x)，方差：Dx=var (x)，计算可得：X(t)的均值为0.0019，X(t)的方差为0.7590。如图4所示，其中蓝色线为X(t)一个样本的时域波形，红色点连成的线为X(t)的均值，绿色点连成的线为X(t)的方差。图4 窄带随机信号时域波形2、求X(t)的概率密度，方法是将最大最小区间分成14等份，然后分别计算各个区间的个数，如图02中柱形条所示，利用曲线拟合，得到合适的概率密度函数。为了得到光滑的曲线，利用了多项式拟合，经过测试，9次拟合曲线比较符合要求，获得的曲线如图5中曲线所示：图5 X(t)的概率分布密度函数3、对X(t)进行频谱分析，在Matlab中，利用fft函数可以很方便得求得X(t)的频谱，然后用abs和angle函数求得幅值和相位，画出图像如图6-1所示，图6-2为无噪声叠加时X(t)的频谱图：图6 -1X(t)的频谱图（有噪声叠加）图6 -2X(t)的频谱图（无噪声叠加）4、求X(t)的自相关函数，用xcorr函数求出自相关序列，得到X(t)自相关函数的时域波形，如图7-1所示，图7-2为无噪声叠加时X(t)的自相关函数图。图7 -1X(t)自相关函数的时域波形（有噪声叠加）图7-2 X(t)自相关函数的时域波形（无噪声叠加）5、对X(t)自相关函数进行fft变换，得到X(t)的功率谱密度，如图8-1所示。图8-2为无噪声叠加时X(t)的功率谱密度图：图8 -1X(t)的功率谱密度（有噪声叠加）图8 -2X(t)的功率谱密度（无噪声叠加）6、建立滤波器，建立一个巴特沃思滤波器，对产生的x(t)进行检测。滤波器的幅度谱和相位谱所示：图9 地通滤波器的幅度谱和相位谱7、求Ac(t)的统计特性，Ac(t)