窄带随机信号性能分析

1、窄带随机信号性能分析一摘要窄带信号在通信系统中有着重要的意义，信号处理技术及通信网络系统与计算机分析技术的相互融合，都要求我们对研究分析随机信号经过系统的响应有一个深入的了解。本实验包括四部分：窄带信号及包络和相位检波分析，窄带随机信号的仿真与分析，希尔伯特变换在单边带系统中的应用，随机信号的dsb分析。主要涉及窄带滤波器的设计，高斯窄带信号包络的均值，均方值和方差的测定，相位概率密度函数的测定等。通过实验了解窄带信号在信号处理领域的应用。 复杂的实际通信系统可以通过抽象与仿真来研究它的特性。本实验通过matlab中的仿真出理想信号，并对其进行分析与测量。2 实验特点与原理1. 窄带信号及包络

2、和相位检波分析 一般无线电接收机中，通常都有高频或中频放大器，它们的通频带往往远小于中心频率，既有 这种线性系统通称为窄带线性系统。在通信、雷达等许多电子系统中，都常常用一个宽带平稳随机过程来激励一个窄带滤波器，这是在滤波器输出端得到的便是一个窄带随机过程。若用示波器观测此波形，则可看到，它接近一个正弦波，但此正弦波的幅度和相位都在缓慢的随机变化。我们可以证明，任何一个实窄带随机过程x(t)都可以表示为： 式中， 是固定值，对于窄带随机过程来说，一般取窄带滤波器的中心频率或载波频率。在实际应用中，常常需要检测出包络和的信息。若将窄带随机过程x(t)送入包络检波器，则在检波器的输出端可得到包络；

3、若将窄带随机过程x(t)送入一个相位检波器，便可检测出相位信息。如图10所示： 图10 窄带信号及包络和相位检波器图10中，在相位检波器之前加入一个理想限幅器，其作用是消除包络起伏对相位检波器的影响。2. 窄带随机信号的仿真与分析 （1）.窄带随机过程一个实平稳随机过程x(t)，若它的功率谱密度具有下述性质：中心频率为c，带宽为=20，当c 时，就可认为满足窄带条件。若随机过程的功率谱满足该条件则称为窄带随机过程。若带通滤波器的传输函数满足该条件则称为窄带滤波器。随机过程通过窄带滤波器传输之后变成窄带随机过程。图1 为典型窄带随机过程的功率谱密度图。若用一示波器来观测次波形，则可看到，它接近于

4、一个正弦波，但此正弦波的幅度和相位都在缓慢地随机变化，图2所示为窄带随机过程的一个样本函数。图1 典型窄带随机过程的功率谱密度图图2 窄带随机过程的一个样本函数(2).窄带随机过程的数学表示1).用包络和相位的变化表示由窄带条件可知,窄带过程是功率谱限制在c附近的很窄范围内的一个随机过程,从示波器观察(或由理论上可以推知):这个过程中的一个样本函数(一个实现)的波形是一个频率为c且幅度和相位都做缓慢变化的余弦波。写成包络函数和随机相位函数的形式：x(t)=a(t)*cosct+ (t)其中:a(t)称作x(t)的包络函数; (t)称作x(t)的随机相位函数。包络随时间做缓慢变化，看起来比较直观

5、，相位的变化，则看不出来。2).莱斯（rice）表示式任何一个实平稳随机过程x(t)都可以表示为：x(t)=ac(t) cosct-as(t) sinct其中同相分量：ac(t)= x(t) cost= x(t) cosct sinct=lpx(t) *2cosct正交分量：as(t) = x(t)sint= cosct x(t) sinct= lp-x(t) *2sinct（lpa表示取a的低频部分）。ac(t)和as(t)都是实随机过程，均值为0，方差等于x(t)的方差。3. 希尔伯特变换子单边带系统中的应用 在单边带幅度调制中，可以保留上边带，也可以保留下边带。信号单边带调制可以提高信道

6、的利用率。信号单边调制（ssb）有上边带（usb）和下边带（lsb）两种，一般利用希尔伯特变换来实现。1).利用希尔伯特实现单边带调制的原理框图如下所示：低通滤波器 希尔伯特变换xt(t)a点单边输出 图1 利用希尔伯特变换实现单边带调制框图 其中输入信号x(t)： x(t)=s(t)+n(t)。s(t)为频率为1khz、幅值为1v的正弦波信号。载波为4 khz、幅值为1v的正弦波信号。n(t)为高斯噪声。2).希尔伯特变换器的时域特性h(t)为对上式进行傅里叶变化，可得希尔伯特变换器的频率特性h(jw)为：由以上可知，希尔伯特变换器的幅度响应为h(jw)=1,相位响应为（w）=-sgn(w)

7、，因此，希尔伯特变换器是一个全通系统，称为90度相移器。3).如果调制信号的频谱为x(jw)，则对usb(t)及lsb(t)的时域表达式两边进行傅里叶变换可得下式：yusb(jw)= 利用上面xh(jw)与x(jw)的关系，将xh(jw)用x(jw)替换得： =4.随机信号dsb分析1）.dsb信号的模型在am信号中，载波分量并不携带信息，信息完全由边带传送。如果将载波抑制，只需在将直流去掉，即可输出抑制载波双边带信号，简称双边带信号（dsb）。dsb调制器模型如图1所示。图1 dsb调制器模型 其中，设正弦载波为式中，为载波幅度；为载波角频率；为初始相位（假定为0）。调制过程是一个频谱搬移的

8、过程，它是将低频信号的频谱搬移到载频位置。而解调是将位于载频的信号频谱再搬回来，并且不失真地恢复出原始基带信号。双边带解调通常采用相干解调的方式，它使用一个同步解调器，即由相乘器和低通滤波器组成。在解调过程中，输入信号和噪声可以分别单独解调。相干解调的原理框图如图2所示：图2 相干解调器的数学模型信号传输信道为高斯白噪声信道，其功率为。2).dsb信号调制过程分析假定调制信号的平均值为0，与载波相乘，即可形成dsb信号，其时域表达式为式中，的平均值为0。dsb的频谱为dsb信号的包络不再与调制信号的变化规律一致，因而不能采用简单的包络检波来恢复调制信号， 需采用相干解调(同步检波)。另外，在调

9、制信号的过零点处，高频载波相位有180的突变。除了不再含有载频分量离散谱外，dsb信号的频谱与am信号的频谱完全相同，仍由上下对称的两个边带组成。所以dsb信号的带宽与am信号的带宽相同，也为基带信号带宽的两倍， 即式中，为调制信号的最高频率。 3).高斯白噪声信道特性分析在实际信号传输过程中，通信系统不可避免的会遇到噪声，例如自然界中的各种电磁波噪声和设备本身产生的热噪声、散粒噪声等，它们很难被预测。而且大部分噪声为随机的高斯白噪声，所以在设计时引入噪声，才能够真正模拟实际中信号传输所遇到的问题，进而思考怎样才能在接受端更好地恢复基带信号。信道加性噪声主要取决于起伏噪声，而起伏噪声又可视为高

10、斯白噪声，因此我在此环节将对双边带信号添加高斯白噪声来观察噪声对解调的影响情况。为了具体而全面地了解噪声的影响问题，我将分别引入大噪声（信噪比为20db）与小噪声（信噪比为2db）作用于双边带信号，再分别对它们进行解调，观察解调后的信号受到了怎样的影响。在此过程中，我用函数来添加噪声，此函数功能为向信号中添加噪声功率为其方差的高斯白噪声。正弦波通过加性高斯白噪声信道后的信号为故其有用信号功率为噪声功率为信噪比满足公式则可得到公式我们可以通过这个公式方便的设置高斯白噪声的方差。4).dsb解调过程分析所谓相干解调是为了从接收的已调信号中，不失真地恢复原调制信号，要求本地载波和接收信号的载波保证同

11、频同相。相干解调的一般数学模型如图所示。 图5 dsb相干解调模型设图四的输入为dsb信号乘法器输出为通过低通滤波器后当常数时，解调输出信号为3 实验的设计思想与实现1. 窄带信号及包络和相位检波分析程序实现clear allwp1=2*pi*1*103;ws1=4*pi*1*103;rp=1;rs=30;wp1=1;ws1=ws1/wp1;n1,wc1=buttord(wp1,ws1,rp,rs,s);z1,p1,k1=buttap(n1);b1,a1=zp2tf(z1,p1,k1);w1=0:pi/512:2*pi;h1,w1=freqs(b1,a1,w1);plot(w1,20*log1

12、0(abs(h1),k);grid;xlabel(lambda);ylabel(a(lambda)/db);title(低通滤波器幅度谱)grid on;plot(w1,angle(h1);%在 02 区间内作相位谱title(低通滤波器相位谱);grid on;fs=20000; %设定采样频率n=100; %取的样本点数n=0:n-1; %取样点序列t=n/fs; %获得以1/16000为时间间隔采样序列noisy=wgn(1,n,0); %产生高斯白噪声a=cos(2*pi*1000*t); %获取a(t)的采样点x=(1+a).*cos(2*pi*4000*t+2)+noisy/10;

13、 %获取x(t)的采样点figure(1); %画图subplot(2,1,1);plot(n,x);%以t为横坐标画出x(t)的时域图型axis(0 100 -2.5 2.5);xlabel(采样点);ylabel(x(t)/v);title(图1-窄带随机信号波形);%窄带系统检测xc=filter(b1,a1,x);%对信号进行滤波ac=xc.*2.*cos(2*pi*4000*t);as=xc.*2.*sin(2*pi*4000*t);y=ac.*cos(2*pi*4000*t)-as.*sin(2*pi*4000*t);%作输出信号时域波形figure(10)subplot(2,1,

14、1);n=0:n-1;plot(n,y);%作信号时域波形axis(0 100 -3 3);%axis tight; xlabel(采样点 );ylabel(幅值) title(y时域谱);grid on;%进行 fft 变换并做频谱图yy=fft(y,length(y);%对滤波后信号进行 fft 变换longy=length(yy);%求傅里叶变换后的序列长度labelx=(0:longy-1)*16000/longy;magn3=abs(yy);% 求滤波后信号幅值subplot(2,1,2);plot(labelx,magn3);%做频谱图axis tight; xlabel(频率(h

15、z); ylabel(幅值 ); title(y频谱图);grid on;%求输出信号的自相关函数c3,lags3=xcorr(y,coeff);%求出自相关序列figure(11);subplot(2,1,1);plot(lags3/fs,c3);%在时域内画自相关函数xlabel(); ylabel(rx(); axis tight;title(y的自相关函数);grid on;%11.作窄带随即过程y(t)的相位图xiangwei=angle(y);figure(9);plot(n,xiangwei);%以t为横坐标画出y(t)的时域图型axis(0 50 0 4);title(图9-y

16、(t)的相位);grid on;%10.做相位得一维概率函数figure(6);hist(xiangwei);title(相位的一维概率密度);grid on;each1=linspace(min(xiangwei),max(xiangwei),10);%将最大最小区间分成14等份,然后分别计算各个区间的个数nr1=hist(xiangwei,each1);%计算各个区间的个数nr1=nr1/length(xiangwei);%计算各个区间的个数归一化subplot(2,1,2);%多项式拟合bar(each1,nr1);hold on;plot(each1,nr1,g)hold on;p1=

17、polyfit(each1,nr1,5);eachi1=0:0.1:3.5;nri1=polyval(p1,eachi1);plot(eachi1,nri1,r)axis tight; title(相位概率密度分布);xlabel(xiangwei);ylabel(p1(x);grid on;%9.窄带随即过程y(t)的包络a=abs(y);figure(7);plot(n,a);axis(0 100 0 5);xlabel(时间t/s);ylabel(包络a(t)/v);title(图7-y(t)的包络);grid on;%10.作包络的一维概率密度图figure(8);hist(a);ti

18、tle(包络的一维概率密度);grid on;each=linspace(min(a),max(a),14);%将最大最小区间分成14等份,然后分别计算各个区间的个数nr=hist(a,each);%计算各个区间的个数nr=nr/length(a);%计算各个区间的个数归一化subplot(2,1,2);%多项式拟合bar(each,nr);hold on;plot(each,nr,g)hold on;p=polyfit(each,nr,6);eachi=0:0.1:3.5;nri=polyval(p,eachi);plot(eachi,nri,r)axis tight; title(包络概率

19、密度分布);xlabel(a);ylabel(p(x)2. 窄带随机信号的仿真与分析1） .程序设计clcfs=16000; %设定采样频率n=1300;n=0:n-1; %取的样本点数t=n/fs; %获得以1/16000为时间间隔采样序列noisy=wgn(1,n,0); %产生高斯白噪声w1=2*pi*1000;%低频信号角频率w2=2*pi*4000;%载波角频率a=cos(2*pi*1000*t); %获取a(t)的采样点x=(1+a).*cos(2*pi*4000*t+2)+noisy/10; %获取x(t)的采样点%以t为横坐标画出x(t)的时域图型figure(1); subp

20、lot(2,1,1); plot(n,x); axis(0 140 -3 3);xlabel(采样点);ylabel(x(t)/v);title(窄带随机信号波形);grid on;%求x(t)的统计特性 并画出来disp(x(t)的均值为); ex=mean(x); disp(ex);%求x(t)均值hold on; plot(n,ex,r.); disp(x(t)的方差为);dx=var(x); disp(dx);%求x(t)方差hold on; plot(n,dx,g.);%画出x(t)的概率分布函数each=linspace(min(x),max(x),14); %将最大最小区间分成1

21、4等份,然后分别计算各个区间的个数nr=hist(x,each); %计算各个区间的个数nr=nr/length(x); %计算各个区间的个数归一化subplot(2,1,2); p=polyfit(each,nr,9); %画出概率分布直方图bar(each,nr); %多项式拟合hold on; plot(each,nr,g)eachi=-2:0.1:2;nri=polyval(p,eachi);plot(eachi,nri,r)axis tight;title(x(t)概率密度分布);xlabel(x(t);ylabel(p(x);grid on;%对x(t)进行频谱分析fx=fft(x

22、,n); %对x(t)进行fft变换，在016000区间内得到2n-1个频率值magn=abs(fx); %求x(t)幅值xangle=angle(fx); %求x(t)相位labelang=(0:length(x)-1)*16000/length(x); %在016000区间内求横坐标刻度figure(2); plot(labelang,magn*10); %在016000区间内做频谱和相位图axis(0 16000 -0.5 600); xlabel(频率/hz);ylabel(幅值);title(x(t)频谱图);grid on;%求x(t)的自相关函数c,lags=xcorr(x,co

23、eff); %求出自相关序列figure(3); subplot(2,1,1); plot(lags/fs,c); %在时域内画自相关函数axis tight; xlabel(t);ylabel(rx(t);title(x(t)的自相关函数);grid on;%求x(t)的功率谱密度long=length(c); sx=fft(c,long); labelx=(0:long-1)*2*pi;plot_magn=10*log10(abs(sx);subplot(2,1,2); plot(labelx,plot_magn); %画功率谱密度axis tight;xlabel(w);ylabel(s

24、x(w);title(x(t)的功率谱密度);grid on; %窄带系统检测z1=2.*cos(2*pi*4000*t); z2=-2.*sin(2*pi*4000*t); ac=z1.*x; %滤波后生成ac(t)as=z2.*x; %滤波后生成as(t)y=ac.*cos(2*pi*4000*t)-as.*sin(2*pi*4000*t);%滤波器设计f_p=1000;f_s=1600;r_p=1;r_s=35; %设定滤波器参数; 通、阻带截止频率,通、阻带衰减ws=2*f_s/fs;wp=2*f_p/fs; %频率归一化n,wn=buttord(wp,ws,r_p,r_s); %采用

25、巴特沃思滤波器b,a=butter(n,wn); %求得滤波器传输函数的多项式系数figure(4);h,w=freqz(b,a); %求得滤波器传输函数的幅频特性subplot(2,1,1); plot(w*fs/(2*pi),abs(h); %在02pi区间内作幅度谱title(低通滤波器幅度谱); grid on;subplot(2,1,2); plot(w*fs/(2*pi),angle(h); %在02pi区间内作相位谱title(低通滤波器相位谱); grid on;%求ac(t)滤波后的统计特性mc=filter(b,a,ac); %上支路通过滤波器 ac(t)disp(ac(t

26、)的均值);eh=mean(mc) %求ac(t)的均值disp(ac(t)的均方值是);e2h=mc*mc/n %求ac(t)的均方值disp(ac(t)的方差);dh=var(mc) %求ac(t)的方差%画ac(t)的时域波形figure(6); subplot(2,1,1); n=0:n-1; plot(n,mc); axis(0 300 -1 1);xlabel(采样点);ylabel(幅值);title(ac(t)的时域波形);grid on;%画ac(t)的频谱图yc=fft(mc,length(mc); %对ac(t)进行fft变换longc=length(yc); %求傅里叶

27、变换后的序列长度labelx=(0:longc-1)*16000/longc;magnl=abs(yc); %求ac(t)的幅值subplot(2,1,2); plot(labelx,magnl); %画ac(t)的频谱图axis tight; xlabel(频率(hz); ylabel(幅值); title(ac(t)频谱图); grid on;%求ac(t)的自相关函数c1,lags1=xcorr(mc,coeff); %求出ac(t)的自相关序列figure(7); subplot(2,1,1); plot(lags1/fs,c1); %在时域内画ac(t)的自相关函数xlabel(t)

28、;ylabel(rx(t);axis tight;title(ac(t)的自相关函数);grid on;%求ac(t)的双边功率谱sac=fft(c1,length(c1); %对ac(t)的自相关函数进行傅里叶变换magnc=abs(sac); %求ac(t)的双边功率谱幅值long=length(sac); %求傅里叶变换后的序列长度labelc=(0:long-1)*16000/long;subplot(2,1,2); plot(labelc,10*log10(magnc); %画ac(t)的自相关函数频谱 即为ac(t)的双边功率谱xlabel(频率(hz);ylabel(功率谱(db

29、w);axis tight;title(ac(t)的双边功率谱);grid on;%求得as(t)的统计特性ms=filter(b,a,as); %对下支路信号进行滤波得as(t)disp(as(t)的均值); eh=mean(ms) %求as(t)的均值disp(as(t)的均方值是); e2h=ms*ms/n %求as(t)的均方值disp(as(t)的方差); dh=var(ms) %求as(t)的方差 %作as(t)的时域波形figure(8);subplot(2,1,1); n=0:n-1;plot(n,ms); %画出as(t)的时域波形axis(0 300 -0.5 2); xl

30、abel(采样点);ylabel(幅值);title(as(t)的时域波形);grid on;%对as(t)进行fft变换并做频谱图ys=fft(ms,length(ms); %对as(t)进行fft变换longs=length(ys); %求傅里叶变换后的序列长度labelx=(0:longs-1)*16000/longs;magn2=abs(ys); %求as(t)的幅值subplot(2,1,2); plot(labelx,magn2); %画出as(t)的频谱图axis tight;xlabel(频率(hz);ylabel(幅值);title(as(t)的频谱图);grid on;%求

31、as(t)的自相关函数c2,lags2=xcorr(ms,coeff); %求出as(t)的自相关序列figure(9);subplot(2,1,1);plot(lags2/fs,c2); %画出as(t)自相关函数的时域波形xlabel(t);ylabel(rx(t);axis tight;title(as(t)的的自相关函数);grid on;%求as(t)的双边功率谱sas=fft(c2,length(c2); %对as(t)的自相关函数进行傅里叶变换magnc=abs(sac); %求as(t)的双边功率谱幅值long=length(sas); %求傅里叶变换后的序列长度labels=

32、(0:long-1)*16000/long;subplot(2,1,2); plot(labelc,10*log10(magnc); %画as(t)的自相关函数频谱xlabel(频率(hz);ylabel(功率谱(dbw);axis tight;title(as(t)的双边功率谱); % 求y(t)的统计特性disp(输出信号y(t)的均值);eh=mean(y) %求输出信号y(t)的均值 disp(输出信号y(t)的均方值);e2h=y*y/n %求输出信号y(t)的均方值disp(输出信号y(t)的方差);dh=var(y) %求输出信号y(t)的方差 %作输出信号y(t)的时域波形fi

33、gure(10); subplot(2,1,1);n=0:n-1;plot(n,y); axis(0 150 -2 2);xlabel(采样点);ylabel(幅值);title(y(t)的时域波形);grid on;%进行fft变换并做频谱图yy=fft(y,length(y); %对相加后的信号进行fft变换longy=length(yy); %y(t)傅里叶变换后的序列长度labelx=(0:longy-1)*16000/longy;magn3=abs(yy); %求y(t)的幅值subplot(2,1,2); plot(labelx,magn3); %做y(t)的频谱图axis tig

34、ht;xlabel(频率(hz);ylabel(幅值);title(y(t)的频谱图);grid on;%求输出信号y(t)的自相关函数c3,lags3=xcorr(y,coeff); %求出y(t)的自相关序列figure(11); subplot(2,1,1); plot(lags3/fs,c3); %画y(t)自相关函数的时域波形xlabel(t);ylabel(rx(t);axis tight;title(y(t)的的自相关函数);grid on;%求输出信号y(t)的双边功率谱sy=fft(c3,length(c3); %对y(t)的自相关函数进行傅里叶变换magny=abs(sy)

35、; %求y(t)双边功率谱幅值long=length(sy);labely=(0:long-1)*16000/long; subplot(2,1,2); plot(labely,10*log10(magny); %*画y(t)的功率谱密度xlabel(频率(hz);ylabel(功率谱(dbw);axis tight;title(y(t)的双边功率谱);grid on;f1=2800;f3=5200;%通带截止频率上下限fsl=2500;fsh=5500;%阻带截止频率上下限rp=0.1;rs=30;%通带边衰减db值和阻带边衰减db值wp1=2*pi*f1/fs;wp3=2*pi*f3/fs

36、;wsl=2*pi*fsl/fs;wsh=2*pi*fsh/fs;wp=wp1 wp3;ws=wsl wsh;n,wn=cheb1ord(ws/pi,wp/pi,rp,rs);%切比雪夫1型带通滤波器bz1,az1=cheby1(n,rp,wp/pi);%切比雪夫1型带通滤波器h,w=freqz(bz1,az1,256,fs);figure(12);subplot(2,1,1);plot(w,abs(h);%绘制带通滤波器的幅频曲线xlabel(频率/hz);ylabel(h/db);title(带通滤波器的幅频曲线);grid on;subplot(2,1,2);plot(w,angle(h

37、);%绘制带通滤波器的相频曲线xlabel(频率/hz);ylabel(w/db);title(带通滤波器的相频曲线);grid on;figure(13);y=filter(bz1,az1,x);%信号通过带通滤波器y=fft(y,length(y);%信号通过带通滤波器的频谱labelk=(-(length(y)-1)/2:(length(y)-1)/2)*fs/length(y);plot(labelk,abs(y),b-);%信号通过带通滤波器的幅频曲线xlabel(频率/hz);ylabel(y(h);title(x(t)通过带通滤波器后的幅频特性);grid on;窄带随机过程仿真

38、建模要求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中柱形条所示，利用曲线

39、拟合， 得到合适的概率密度函数。 为了得到光滑的曲线，利用了多项式拟合，经过测试，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)的自相关函

40、数图。图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)为x(t) *2cosct通过低通滤波器的信号，ac(t)的均值eh = -0.4075 4（带有直流分量），

41、ac(t)的均方值是e2h =0.2458ac(t)的方差dh = 0.0798ac(t)的波形如图10、图11所示：图10 -1ac(t)的时域波形图和频谱图（有噪声叠加）图10-2 ac(t)的时域波形图和频谱图（无噪声叠加）图11-1 ac(t)的自相关函数的时域波形图和ac(t)的功率谱密度（有噪声叠加）图11 -2ac(t)的自相关函数的时域波形图和ac(t)的功率谱密度（无噪声叠加）8、求as(t)的统计特性，as(t)为x(t) *2cosct通过低通滤波器的信号，as(t)的均值eh =0.8972（带有直流分量），as(t)的均方值是e2h = 1.1565as(t)的方差d

42、h = 0.3518 as(t)的波形如图13、图14所示：图13-1 as(t)的时域波形图和频谱图（有噪声叠加）图13 -2as(t)的时域波形图和频谱图（无噪声叠加）图14-1 as(t)的自相关函数的时域波形图和as(t)的功率谱密度（有噪声叠加）图14 -2as(t)的自相关函数的时域波形图和as(t)的功率谱密度（无噪声叠加）9、求出y(t)的统计特性，y (t)=ac(t) cosct-as(t) sinct，其统计特性如下输出信号y(t)的均值eh = -4.4011e-004s输出信号y(t)的均方值e2h = 3.0280输出信号y(t)的方差dh = 3.0303y(t)

43、的仿真图形如图15、图16所示。图15-1 y(t)的时域波形图和频谱图（有噪声叠加）图15-2 y(t)的时域波形图和频谱图（无噪声叠加）图16 -1y(t)的自相关函数的时域波形图和y(t)的功率谱密度（有噪声叠加）图16 -2y(t)的自相关函数的时域波形图和y(t)的功率谱密度（无噪声叠加）调幅调相波的相干解调在实际应用中，常常需要检测出包络和的信息。若将窄带随机过程x(t)送入包络检波器，则在检波器的输出端可得到包络；若将窄带随机过程x(t)送入一个相位检波器，便可检测出相位信息。如图10所示： 窄带信号及包络和相位检波器图中，在相位检波器之前加入一个理想限幅器，其作用是消除包络起伏

44、对相位检波器的影响。带通滤波器：此带通滤波器的作用是滤除已调信号频带以外的噪声。带通滤波器的参数如下：通带截止频率上下限：2800hz,5000hz阻带截止频率上下线：2500hz,5200hz阻带衰减：30db通带衰减：0.1db该带通滤波器的幅频和相频特性曲线如图9所示：信号通过带通滤波器后的幅频特性：3. 希尔伯特变换子单边带系统中的应用3.1实现的程序框图如下图所示 输 入画出a点的是与信号和频域信号并与未通过时进行比较设置参数（采样频率、频率点数、观察时间长度）利用希尔伯特变换产生ssb调制设置基带信号、载波信号的频率进行傅里叶变换对其频谱分析产生高斯白噪声，基带信号加上噪声并频谱分

45、析 结 束 基带信号加上噪声后通过低通滤波器3.2利用matlab具体实现的代码如下：1）参数设定clc clear allfs=15000; %采样频率ts=1/fs;%采样周期t=0:ts:0.01;%时间序列df=0.2;%采样分辨率m=2048;%频率点数fc = 4000; %载波频率lt=length(t); %时间序列长度l=2*min(at); r=2*max(abs(at);2）产生高斯白噪声n(t)并进行频谱分析nt = wgn(1,length(t),0.1);n_1=nt/max(abs(nt);%噪声figure(1);subplot(211);plot(t,n_1)

46、;title(高斯白噪声nt信号);xlabel(t/s);ylabel(幅度/v);grid on;pausen=0:m-1; %t=n/fs; %时间序列y0=fft(n_1,m);mag0=(abs(y0);f=n*fs/(1000*m);subplot(212);plot(f,mag0);title(高斯白噪声频谱分析);xlabel(f/khz);ylabel(幅度/v);axis(0 10 0 20);grid on;pause运行结果如下：3）产生基带信号s(t)并进行频谱分析st=sin(1000*2*pi*t);subplot(211);plot(t,st);title(初始

47、信号st=sin(1000*2*pi*t); xlabel(t/s);ylabel(幅度/v);grid on;pausey1=fft(st,m);mag1=(abs(y1);f=n*fs/(1000*m);subplot(212);plot(f,mag1);title(初始信号频谱分析);xlabel(f/khz);ylabel(幅度/v);grid on;axis(0 10 0 100);pause运行结果如下：4）调制信号（s(t)+n(t)）进行频谱分析xt=st+n_1;subplot(211);plot(t,xt);title(调制信号xt=st+nt（初始信号+噪声）); xla

48、bel(t/s);ylabel(幅度/v);grid on;pausey3=fft(xt,m);mag3=(abs(y3);f=n*fs/(1000*m);subplot(212);plot(f,mag3);title(调制信号频谱分析);xlabel(f/khz);ylabel(幅度/v);axis(0 10 0 100);grid on;pause运行结果如下：5）调制信号通过滤波器哟后a点的信号分析wp=2*2200/fs; %通带边界频率ws=2*2800/fs; %阻带边界频率rp=1; %通带最大衰减度as=30; %阻带最小衰减度v,wc=buttord(wp,ws,rp,as)

49、;b,a=butter(v,wc);h,w=freqz(b,a);at=filter(b,a,xt); %经过低通滤波器的a点信号figure(2)subplot(311);plot(w,abs(h);title(低通滤波器信号); xlabel(t/s);ylabel(幅度/v);grid on;pausey3=fft(at,m);mag3=(abs(y3);f=n*fs/(1000*m);subplot(312);plot(t,at);title(经过滤波器后的调制信号)xlabel(t/s);ylabel(幅度/v);grid on;subplot(313);plot(f,mag3);t

50、itle(调制信号经过低通滤波器后频谱分析); xlabel(f/khz);ylabel(幅度/v);grid on;axis(0 10 0 100);pause运行结果如下：6）信号经过希尔伯特变换产生ssb调制figure(3);subplot(3,2,1);plot(t,at);title(经过滤波器后的调制信号)xlabel(t/s);ylabel(幅度/v);grid on;pause %按任意键可看到调制信号的曲线c1=cos(2*pi*fc*t);c2=sin(2*pi*fc*t);subplot(3,2,3);u1=at(1:lt).*c1(1:lt)+imag(hilbert(at(1:lt).*c2(1:lt);plot(t,u1);title(下边带调制信号);xlabel(t/s);ylabel(幅度/