试验一白噪声检验.docx
- 文档编号:23039475
- 上传时间:2023-04-30
- 格式:DOCX
- 页数:31
- 大小:647.56KB
试验一白噪声检验.docx
《试验一白噪声检验.docx》由会员分享,可在线阅读,更多相关《试验一白噪声检验.docx(31页珍藏版)》请在冰豆网上搜索。
试验一白噪声检验
白噪声测试
实验目的
⑴了解白噪声信号的特性,包括均值(数学期望)、均方值、方差、相关函数、概率密度、频谱及功率谱密度等。
⑵掌握白噪声信号的分析方法。
二、实验原理
所谓白噪声是指它的概率统计特性服从某种分布而它的功率谱密度又是均匀的。
确切的
说,白噪声只是一种理想化的模型,因为实际的噪声功率谱密度不可能具有无限宽的带宽,否
则它的平均功率将是无限大,是物理上不可实现的。
然而白噪声在数学处理上比较方便,所以它在通信及电子工程系统的分析中有十分重要的作用。
一般地说,只要噪声的功率谱密度的宽度远大于它所作用的系统的带宽,并且在系统的带内,它的功率谱密度基本上是常数,就可以作为白噪声处理了。
白噪声的功率谱密度为:
其中N。
为单边功率谱密度。
白噪声的自相关函数为:
No
白噪声的自相关函数是位于T=0处、强度为N0的冲击函数。
这表明白噪声在任何两个不
2
同的瞬间的取值是不相关的。
同时也意味着白噪声能随时间无限快的变化,因为它的带宽是无
限宽的。
下面我们给出几种分布的白噪声。
随机过程的几种分布
前人已证明,要产生一个服从某种分布的随机数,可以先求出其分布函数的反函数的解析式,再将一个在[0,1]区间内的均匀分布的随机数的值代入其中,就可以计算出服从某种分布的随机数。
下面我们就求解这些随机数。
[0,1]区间均匀分布随机信号的产生:
采用混合同余法产生[0,1]区间的均匀分布随机数。
混合同余法产生随机数的递推公式为:
yn1ayncn=0,1,2……
Xn普n=1,2,3……
M
由上式的出如下实用算法:
Xn1aXn
X0M
cc
[aXn]
MM
其中:
M
2k,
其中k为计算几种数字尾部的字长
a
4t1
,t为任意选定的正整数
y。
为任意非负整数
c,为奇数
Matlab语言中的rand()函数是服从[0,1]均匀分布的,所以在以后的实验中如果用到均匀分布的随机数,我们统一使用rand()函数
正态分布(高斯分布)随机信号的产生:
高斯分布的密度函数为:
f(x)--exp(+)弋22)
采用变换法产生正态分布随机数,若Ri、R2示[0,1]均匀分布随机数,则有正态分布随机数:
2lnRicos2R2
2lnRisin2R2
指数分布随机信号的产生:
指数分布的密度函数为:
三、实验内容与结果
1.产生五种概率分布的信号
Matlab程序:
%生成各种分布的随机数
x4=raylrnd(1,1,1024);%生成长度为1024的瑞利分布
%时域特性曲线:
figure;
subplot(3,2,1),plot(1:
1024,x1);grid
on;title('均匀分布');xlabel('时间(t)');ylabel('幅度');axis([01024-22]);
subplot(3,2,2),plot(1:
1024,x2);grid
on;title('正态分布');xlabel('时间(t)');ylabel('幅度');axis([01024-22]);
subplot(3,2,3),plot(1:
1024,x3);grid
on;title('指数分布');xlabel('时间(t)');ylabel('幅度');axis([01024-15]);
subplot(3,2,4),plot(1:
1024,x4);grid
on;title('瑞利分布');xlabel('时间(t)');ylabel('幅度');axis([01024-14]);
subplot(3,2,5),plot(1:
1024,x5);grid
on;title('卡方分布');xlabel('时间(t)');ylabel('幅度');axis([01024-15]);
瑞利分布
时间(t)卡方分布
2.均值:
均值E[x(t)]表示集合平均值或数学期望值。
基于随机过程的各态历经性,可用时间间隔t
内的幅值平均值表示:
N-1
EM)]=
r=0
均值表达了信号变化的中心趋势,或称之为直流分量
在MATLAB中,可以用
mean()函数来计算。
%求各种分布的均值
figure;
m1=mean(x1);m2=mean(x2);m3=mean(x3);m4=mean(x4);m5=mean(x5);
subplot(3,2,1),plot(1:
1024,m1);title(
'均匀分布均值');xlabel('时间(t)');ylabel('幅度');axis([01024-22]);
subplot(3,2,2),plot(1:
1024,m2);title(
'高斯分布均值');xlabel('时间(t)');ylabel('幅度');axis([01024-22]);
subplot(3,2,3),plot(1:
1024,m3);title(
'指数分布均值');xlabel('时间(t)');ylabel('幅度');axis([01024-22]);
subplot(3,2,4),plot(1:
1024,m4);title(
'瑞利分布均值');xlabel('时间(t)');ylabel('幅度');axis([01024-22]);
subplot(3,2,5),plot(1:
1024,m5);title(
'卡方分布均值');xlabel('时间(t)');ylabel('幅度');axis([01024-22]);
均匀分布均值
时间(t)指数分布均值
时间(t)卡方分布均值
瑞利分布均值
%10020030040050060070»8009001000
时间(t>
3.方差:
随机过程的方差函数描述了随机过程所有样本函数在t时刻的函数值相对于其数学期望
的偏离程度。
定义:
4L
t-Q
其中c(t)是随机过程的标准差。
当随即过程表征的是接收机输出端的噪声电压时,c2(t)表
示小号在单位电阻上的瞬时交流功率统计平均值,而c(t)表示噪声电压相对于电压统计平均值的交流分量
在MATLAB中,可以用std()函数计算出标准差c(t),再平方就可以得到方差
%求各种分布的方差
figure;
v1=var(x1);v2=var(x2);v3=var(x3);v4=var(x4);v5=var(x5);
subplot(3,2,1),plot(1:
1024,v1);grid
on;title('均匀分布方差');xlabel('时间(t)');ylabel('幅度');
subplot(3,2,2),plot(1:
1024,v2);grid
on;title('高斯分布方差');xlabel('时间(t)');ylabel('幅度');
subplot(3,2,3),plot(1:
1024,v3);grid
on;title('指数分布方差');xlabel('时间(t)');ylabel('幅度');
subplot(3,2,4),plot(1:
1024,v4);grid
on;title('瑞利分布方差');xlabel('时间(t)');ylabel('幅度');
subplot(3,2,5),plot(1:
1024,v5);grid
on;title('卡方分布方差');xlabel('时间(t)');ylabel('幅度');
时间Ct)指数分布方差
卡方分布方差
高斯分布方差
4•自相关:
信号的相关性是指客观事物变化量之间的相依关系。
对于平稳随机过程x(t)和y(t)在两个
不同时刻t和t+t的起伏值的关联程度,可以用相关函数表示。
在离散情况下,信号x(n)和y(n)
的相关函数定义为:
T=0t
随机信号的自相关函数表示波形自身不同时刻的相似程度。
与波形分析、频谱分析相比,它具有能够在强噪声干扰情况下准确地识别信号周期的特点。
%求各种分布的自相关函数
);[x_c2,lags]=xcorr(x2,200,'unbiased');[x_c3,lags]=xcorr(x3,200,
);[x_c5,lags]=xcorr(x5,200,'unbiased');
on;title('均匀分布自相关');
on;title('正态分布自相关');
on;title('指数分布自相关');
on;title('瑞利分布自相关');
on;title('卡方分布自相关');
均匀分布目相关
指数分布自相关
11
1k
i■
1fe
11
11
■iiii
h1ft
■11H
H1l>
1[|B
11・
11
ii
ii
i■
1p
1It
1■[
1H
ii■[
11
11l>
iiiii
V1•
■aii
91l>
i1l!
k1t1
■[IB
H4l>
■[|H
1iI
-200-1G0-100-50C50100150200
■111
1||■[
i<1I<
Iil1
fell''
111
1II||
H4b
■cii
l>41
11I
I11
11■
1IIH
14I'
i[ii
|ji|l>
i口I
」■■■■■■J.J1
iy■!
ft1t Illi --fa--i----—— |1||i IlHl! |1|! 1! lH 141 ■ciii fciI -200-150-100-50050100150200 瑞利分布自相关 111 111 11V ■11 1>1■ 111 111 11* I1|i iii1 p■|i Ii41 13|l HXl> 11■ 11■ 11l> I1I f1» iii I41 isii *4l> iiiii V4I h«Ii ■■■・■■■■t■孑■■■■■■/■■■■■■■』■■ 11■1 fell* ■il■II \\I\ 111 ■! il> ■II1 iii 〔41iiIiiii -200-150-10ft-50050100150200 5.概率密度函数: 一维分布函数为: 作(幻;6)=P{x(g)莖龙丄} 若Fx(xi;ti)对xi的一阶偏导存在,则一维概率密度为: 3xx 在MATLAB中,可以用ksdensity()函数来计算一维概率密度 %求各种分布的概率密度函数y1=unifpdf(x1,-1,1); y2=normpdf(x2,0,1); y3=exppdf(x3,1); y4=raylpdf(x4,1); y5=chi2pdf(x5,1);%各种分布的概率密度估计figure; [k1,n1]=ksdensity(x1); [k2,n2]=ksdensity(x2);[k3,n3]=ksdensity(x3);[k4,n4]=ksdensity(x4);[k5,n5]=ksdensity(x5); 揩数分布褫率密虞 时间 卡方分布概率密虞 吋间 时间 6.频谱: 信号频谱分析是采用傅立叶变换将时域信号x(t)变换为频域信号x(f),从另一个角度来了 解信号的特征。 时域信号x(t)的傅氏变换为: xCO=f 在MATLAB中,对信号进行快速傅立叶变换fft()就可以得到频谱函数。 %幅频特性曲线 x1=unifrnd(-1,1,1,1024);%生成长度为1024的均匀分布 x2=normrnd(0,1,1,1024);%生成长度为1024的正态分布 x4=raylrnd(1,1,1024);%生成长度为1024的瑞利分布 x5=chi2rnd(1,1,1024);%生成长度为1024的卡方分布 f1=fft(x1,1024); f2=fft(x2,1024); f3=fft(x3,1024); f4=fft(x4,1024); f5=fft(x5,1024); figure; subplot(3,2,1),plot(abs(f1)),axis([01023050]);grid幅值V') subplot(3,2,2),plot(abs(f2));axis([01023050]);grid subplot(3,2,3),plot(abs(f3)),axis([010230100]);grid subplot(3,2,4),plot(abs(f4)),axis([01023050]);grid on;title('均匀分布幅频特性');Xlabel('频率Hz');Ylabel( on;title('正态分布');Xlabel('频率Hz');Ylabel('幅值V') on;title('指数分布');Xlabel('频率Hz');Ylabel('幅值V') on;title('瑞利分布');Xlabel('频率Hz');Ylabel('幅值V') on;title('卡方分布');Xlabel('频率Hz');Ylabel('幅值V') subplot(3,2,5),plot(abs(f5)),axis([010230100]);grid 均匀外布幅瓠特性 频率Hn指数分布 卡方分布 正态分布 50 40 >30麵 @20 10 50 4Q >30 丽 彰20 10 0 010020030040050060070030C9001000 频率Hz 7.功率谱密度: 随机信号的功率谱密度是随机信号的各个样本在单位频带内的频谱分量消耗在一欧姆电 阻上的平均功率的统计均值,是从频域描述随机信号的平均统计参量,表示x(t)的平均功率在 频域上的分布。 它只反映随机信号的振幅信息,而没有反映相位信息。 在MATLAB中,可由下式得到功率谱密度: 错误! 未找到引用源 错误! 未找到引用源 %功率谱密度 figure; f1=fft(x1,1024); f2=fft(x2,1024); f3=fft(x3,1024); f4=fft(x4,1024); f5=fft(x5,1024); p1=mean(f1.*conj(f1))/1024;p2=mean(f2.*conj(f2))/1024;p3=mean(f3.*conj(f3))/1024; p4=mean(f4.*conj(f4))/1024;p5=mean(f5.*conj(f5))/1024; subplot(3,2,1),plot(1: 1024,abs(p1));grid on;title('均匀分布功率谱密度');xlabel('频率Hz');ylabel('幅值'); subplot(3,2,2),plot(1: 1024,abs(p2));grid on;title('正态分布功率谱密度');xlabel('频率Hz');ylabel('幅值'); subplot(3,2,3),plot(1: 1024,abs(p3));grid on;title('指数分布功率谱密度');xlabel('频率Hz');ylabel('幅值'); subplot(3,2,4),plot(1: 1024,abs(p4));grid on;title('瑞利分布功率谱密度');xlabel('频率Hz');ylabel('幅值'); subplot(3,2,5),plot(1: 1024,abs(p5));grid on;title('卡方分布功率谱密度');xlabel('频率Hz');ylabel('幅值'); 正态分布功率谱密度 2iiii Ims二s牛二二m彳 ■i log ■[IIII0L '020040060080010001200 频率缶 瑞制分布功率谱密度 jna I豐 8.取10240个点时的功率谱密度 %功率谱密度 x1=unifrnd(-1,1,1,10240);%生成长度为10240的均匀分布 x2=normrnd(0,1,1,10240);%生成长度为10240的正态分布 x3=exprnd(1,1,10240);%生成长度为10240的指数分布均值为零 x4=raylrnd(1,1,10240); %生成长度为10240的瑞利分布 x5=chi2rnd(1,1,10240); %生成长度为10240的卡方分布 figure; f1=fft(x1,10240); f2=fft(x2,10240); f3=fft(x3,10240); f4=fft(x4,10240); f5=fft(x5,10240); p1=mean(f1.*conj(f1))/10240; p4=mean(f4.*conj(f4))/10240; p2=mean(f2.*conj(f2))/10240; p5=mean(f5.*conj(f5))/10240; p3=mean(f3.*conj(f3))/10240; subplot(3,2,1),plot(1: 10240,abs(p1));grid ylabel('幅值'); subplot(3,2,2),plot(1: 10240,abs(p2));grid ylabel('幅值'); subplot(3,2,3),plot(1: 10240,abs(p3));grid ylabel('幅值'); subplot(3,2,4),plot(1: 10240,abs(p4));grid ylabel('幅值'); subplot(3,2,5),plot(1: 10240,abs(p5));grid on;title('N=10240均匀分布 on;title('N=10240正态分布 on;title('N=10240指数分布 on;title('N=10240瑞利分布 on;title('N=10240卡方分布 ylabel('幅值'); 功率谱密度');xlabel('频率Hz'); 功率谱密度');xlabel('频率Hz'); 功率谱密度');xlabel('频率Hz'); 功率谱密度');xlabel('频率Hz'); 功率谱密度');xlabel('频率Hz'); 2 d N=10240均匀俞布功率谱密度 111 ihl> ■a■ 1I* 111 ll・ 1 ■P・ 111 1PV 悝 111 0 ■IP ■■■■--——-f"™"■■■■■■■■-|-——-■■■■ 1pi 111 -1 1it> ■a■ 1[ i 0 2000 400060008010 10000 12000 频率吃 3 扌旨蹴分布功率谱密度 1R■ ■bk ■■・ 1 Ii I 2 111 1pp lab 1 ■ fa iik H ■pa I■h I il 1 11* 111 111 1 1 0 I 111 III ff・ ■ii 11* 111 i 1 ■ I I 1 0 2000 400060008000 10000 12000 频率业 4 10240卡方符布功率谱密度 1 ~Ir ■ii I99 iib 1ii ■i■ ■1i 111 ll>It ■a■ 2 111 1P1 ■ia ■F■9BW--0(,-BP1 111 IPP i 1pp iih 111 L ( ) 2000 4000600080CO 10000 *2000 频率H? N^IOZlOiE态分布功率谙密度 N=1024-DfS! 利分布功率谓密度 频率吃 9.随机信号叠加及高斯检验: %判断是否是高斯分布 s=zeros(1,1024); fori=0: 4 z=unifrnd(0,1,1,1024); s=s+z; end s1=zeros(1,1024); forj=0: 4 z1=exprnd(3,1,1024); on;title('均匀分布叠加'); on;title('均匀分布叠加') on;title('指数分布叠加'); on;title('指数分布叠加') '均匀分布正态检验'); '指数分布正态检验'); s1=s1+z1; end figure; [j1,
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 试验 噪声 检验