![窄带随机信号性能分析_第1页](http://file2.renrendoc.com/fileroot_temp3/2021-5/3/ae3caf5c-436c-4ff4-a769-a41fa8634764/ae3caf5c-436c-4ff4-a769-a41fa86347641.gif)
![窄带随机信号性能分析_第2页](http://file2.renrendoc.com/fileroot_temp3/2021-5/3/ae3caf5c-436c-4ff4-a769-a41fa8634764/ae3caf5c-436c-4ff4-a769-a41fa86347642.gif)
![窄带随机信号性能分析_第3页](http://file2.renrendoc.com/fileroot_temp3/2021-5/3/ae3caf5c-436c-4ff4-a769-a41fa8634764/ae3caf5c-436c-4ff4-a769-a41fa86347643.gif)
![窄带随机信号性能分析_第4页](http://file2.renrendoc.com/fileroot_temp3/2021-5/3/ae3caf5c-436c-4ff4-a769-a41fa8634764/ae3caf5c-436c-4ff4-a769-a41fa86347644.gif)
![窄带随机信号性能分析_第5页](http://file2.renrendoc.com/fileroot_temp3/2021-5/3/ae3caf5c-436c-4ff4-a769-a41fa8634764/ae3caf5c-436c-4ff4-a769-a41fa86347645.gif)
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
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(幅度/
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 山西省朔州市怀仁县大地学校2024-2025学年高一生物下学期期末考试试题
- 海南省万宁市北京师范大学万宁附属中学2024-2025学年高一数学下学期开学考试试题无答案
- 山东省济南市章丘市第四中学2024-2025学年高二英语下学期第五次质量检测试题含解析
- 云南省昭通市昭阳区2024-2025学年高一化学上学期期末联考试题
- 区矿业市场清理整治工作方案-范文
- 自动磨刀转让合同范本
- 2022年理论力学课后习题答案-理论力学课后习题答案详解
- 饭店合作关系合同范本
- 青岛房屋赠与合同范本
- 山东省滨州市名校2024-2025学年初三5月模拟考英语试题含答案
- 民宿管家服务培训课件
- 2024年中国南水北调集团水网智慧科技有限公司招聘笔试参考题库附带答案详解
- 航天产品壳体数字化生产线关键技术
- 拍摄写真前期策划方案
- 履带抛丸机安全培训
- 课题研究报告和论文区别
- 合同廉洁承诺书
- 竞选自管会演讲稿
- 酒店投资项目融资计划书
- 失眠患者的护理查房课件
- 中职生就业指导(第3版)课件全套 项目1-8 讲究礼仪- 运用法律武器维护自身利益
评论
0/150
提交评论