随机信号分析与处理实验报告
随机信号分析实验报告
《随机信号分析》实验报告二班级_______学号______姓名_______实验二高斯噪声的产生和性能测试1.实验目的(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.。
《随机信号分析与处理》实验报告完整版(GUI)内附完整函数代码
《随机信号分析与处理》实验报告指导教师:班级:学号:姓名:实验一 熟悉MA TLAB 的随机信号处理相关命令一、实验目的1、熟悉GUI 格式的编程及使用。
2、掌握随机信号的简单分析方法3、熟悉语音信号的播放、波形显示、均值等的分析方法及其编程 二、实验原理 1、语音的录入与打开在MATLAB 中,[y,fs,bits]=wavread('Blip',[N1 N2]);用于读取语音,采样值放在向量y 中,fs 表示采样频率(Hz),bits 表示采样位数。
[N1 N2]表示读取从N1点到N2点的值。
2,均匀分布白噪声在matlab 中,有x=rand (a ,b )产生均匀白噪声序列的函数,通过与语言信号的叠加来分析其特性。
3、均值随机变量X 的均值也称为数学期望,它定义为对于离散型随机变量,假定随机变量X 有N 个可能取值,各个取值的概率为则均值定义为上式表明,离散型随机变量的均值等于随机变量的取值乘以取值的概率之和,如果取值是等概率的,那么均值就是取值的算术平均值,如果取值不是等概率的,那么均值就是概率加权和,所以,均值也称为统计平均值。
4、方差定义为随机过程的方差。
方差通常也记为D 【X (t )】 ,随机过程的方差也是时间 t 的函数, 由方差的定义可以看出,方差是非负函数。
5、自相关函数设任意两个时刻1t ,2t ,定义为随机过程X (t )的自相关函数,简称为相关函数。
自相关函数可正,可负,其绝对值越大表示相关性越强。
6.哈明(hamming)窗(10.100)121212121212(,)[()()](,,,)X R t t E X t X t x x f x x t t dx dx +∞+∞-∞-∞==⎰⎰(10.101)B = 1.3Δf,A = -43dB,D= -6dB/oct.哈明窗本质上和汉宁窗是一样的,只是系数不同。
哈明窗比汉宁窗消除旁瓣的效果好一些而且主瓣稍窄,但是旁瓣衰减较慢是不利的方面。
随机信号分析实验报告
随机信号分析实验报告实验一:平稳随机过程的数字特征实验二:平稳随机过程的谱分析实验三:随机信号通过线性系统的分析实验四:平稳时间序列模型预测班级:姓名:学号:一、实验目的1、加深理解平稳随机过程数字特征的概念2、掌握平稳随机序列期望、自相关序列的求解3、分析平稳随机过程数字特征的特点二、实验原理平稳随机过程数字特征求解的相关原理三、实验过程function y = experiment number = 49; %学号49 I = 8; %幅值为8 u = 1/number;Ex = I*0.5 + (-I)*0.5; N = 64; C0 = 1; %计数 p(1) = exp(-u);for m = 2:N k = 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 >当时,/222(){()()}(2)!m k mk 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 XC 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的增大而增大。
-随机信号分析实验报告
-随机信号分析实验报告H a r b i n I n s t i t u t e o f T e c h n o l o g y实验报告课程名称:随机信号分析院系:电⼦与信息⼯程学院班级:姓名:学号:指导教师:实验时间:实验⼀、各种分布随机数的产⽣(⼀)实验原理1.均匀分布随机数的产⽣原理产⽣伪随机数的⼀种实⽤⽅法是同余法,它利⽤同余运算递推产⽣伪随机数序列。
最简单的⽅法是加同余法)(mod 1M c y y n n +=+My x n n 11++= 为了保证产⽣的伪随机数能在[0,1]内均匀分布,需要M 为正整数,此外常数c 和初值y0亦为正整数。
加同余法虽然简单,但产⽣的伪随机数效果不好。
另⼀种同余法为乘同余法,它需要两次乘法才能产⽣⼀个[0,1]上均匀分布的随机数)(mod 1M ay y n n =+ My x n n 11++= 式中,a 为正整数。
⽤加法和乘法完成递推运算的称为混合同余法,即 )(mod 1M c ay y n n +=+ M y x n n 11++=⽤混合同余法产⽣的伪随机数具有较好的特性,⼀些程序库中都有成熟的程序供选择。
常⽤的计算语⾔如Basic 、C 和Matlab 都有产⽣均匀分布随机数的函数可以调⽤,只是⽤各种编程语⾔对应的函数产⽣的均匀分布随机数的范围不同,有的函数可能还需要提供种⼦或初始化。
Matlab 提供的函数rand()可以产⽣⼀个在[0,1]区间分布的随机数,rand(2,4)则可以产⽣⼀个在[0,1]区间分布的随机数矩阵,矩阵为2⾏4列。
Matlab 提供的另⼀个产⽣随机数的函数是random('unif',a,b,N,M),unif 表⽰均匀分布,a 和b 是均匀分布区间的上下界,N 和M 分别是矩阵的⾏和列。
2.随机变量的仿真根据随机变量函数变换的原理,如果能将两个分布之间的函数关系⽤显式表达,那么就可以利⽤⼀种分布的随机变量通过变换得到另⼀种分布的随机变量。
随机信号分析实验报告
实验一 随机噪声的产生与性能测试一、实验内容1.产生满足均匀分布、高斯分布、指数分布、瑞利分布的随机数,长度为N=1024,并计算这些数的均值、方差、自相关函数、概率密度函数、概率分布函数、功率谱密度,画出时域、频域特性曲线; 2.编程分别确定当五个均匀分布过程和5个指数分布分别叠加时,结果是否是高斯分布; 3.采用幅度为2, 频率为25Hz 的正弦信号为原信号,在其中加入均值为2 , 方差为0.04 的高斯噪声得到混合随机信号()X t ,编程求 0()()tY t X d ττ=⎰的均值、相关函数、协方差函数和方差,并与计算结果进行比较分析。
二、实验步骤 1.程序N=1024; fs=1000; n=0:N —1;signal=chi2rnd (2,1,N); %rand(1,N)均匀分布 ,randn(1,N )高斯分布,exprnd(2,1,N )指数分布,raylrnd (2,1,N)瑞利分布,chi2rnd(2,1,N )卡方分布 signal_mean=mean(signal ); signal_var=var (signal );signal_corr=xcorr(signal,signal ,'unbiased ’); signal_density=unifpdf(signal ,0,1); signal_power=fft(signal_corr); %[s,w]=periodogram (signal); [k1,n1]=ksdensity(signal);[k2,n2]=ksdensity (signal,’function ’,'cdf ’); figure ;hist(signal);title (’频数直方图’); figure ;plot (signal);title(’均匀分布随机信号曲线’); f=n *fs/N ; %频率序列 figure;plot(abs (signal_power)); title('功率幅频’); figure;plot(angle (signal_power)); title ('功率相频'); figure;plot (1:2047,signal_corr); title ('自相关函数’); figure;plot(n1,k1);title('概率密度’);figure;plot(n2,k2);title('分布函数’);结果(1)均匀分布(2)高斯分布(3)指数分布(4)瑞利分布(5)卡方分布2.程序N=1024;signal_1=rand(1,N);signal_2=rand(1,N);signal_3=rand(1,N);signal_4=rand(1,N);signal_5=rand(1,N);signal=signal_1+signal_2+signal_3+signal_4+signal_5; [k1,n1]=ksdensity(signal);figure(1)subplot(1,2,1);hist(signal);title('叠加均匀分布随机数直方图');subplot(1,2,2);plot(n1,k1);title(’叠加均匀分布的概率密度');结果指数分布叠加均匀分布叠加结果:五个均匀分布过程和五个指数分布分别叠加时,结果是高斯分布。
随机信号分析实验报告(基于MATLAB语言)
随机信号分析实验报告——基于MATLAB语言姓名:_班级:_学号:专业:目录实验一随机序列的产生及数字特征估计 (2)实验目的 (2)实验原理 (2)实验内容及实验结果 (3)实验小结 (6)实验二随机过程的模拟与数字特征 (7)实验目的 (7)实验原理 (7)实验内容及实验结果 (8)实验小结 (11)实验三随机过程通过线性系统的分析 (12)实验目的 (12)实验原理 (12)实验内容及实验结果 (13)实验小结 (17)实验四窄带随机过程的产生及其性能测试 (18)实验目的 (18)实验原理 (18)实验内容及实验结果 (18)实验小结 (23)实验总结 (23)实验一随机序列的产生及数字特征估计实验目的1.学习和掌握随机数的产生方法。
2.实现随机序列的数字特征估计。
实验原理1.随机数的产生随机数指的是各种不同分布随机变量的抽样序列(样本值序列)。
进行随机信号仿真分析时,需要模拟产生各种分布的随机数。
在计算机仿真时,通常利用数学方法产生随机数,这种随机数称为伪随机数。
伪随机数是按照一定的计算公式产生的,这个公式称为随机数发生器。
伪随机数本质上不是随机的,而且存在周期性,但是如果计算公式选择适当,所产生的数据看似随机的,与真正的随机数具有相近的统计特性,可以作为随机数使用。
(0,1)均匀分布随机数是最最基本、最简单的随机数。
(0,1)均匀分布指的是在[0,1]区间上的均匀分布, U(0,1)。
即实际应用中有许多现成的随机数发生器可以用于产生(0,1)均匀分布随机数,通常采用的方法为线性同余法,公式如下:,序列为产生的(0,1)均匀分布随机数。
定理1.1若随机变量X 具有连续分布函数,而R 为(0,1)均匀分布随机变量,则有2.MATLAB中产生随机序列的函数(1)(0,1)均匀分布的随机序列函数:rand用法:x = rand(m,n)功能:产生m×n 的均匀分布随机数矩阵。
(2)正态分布的随机序列函数:randn用法:x = randn(m,n)功能:产生m×n 的标准正态分布随机数矩阵。
随机信号分析实验报告
一、实验名称微弱信号的检测提取及分析方法二、实验目的1.了解随机信号分析理论如何在实践中应用2.了解随机信号自身的特性,包括均值、方差、相关函数、频谱及功率谱密度等3.掌握随机信号的检测及分析方法三、实验原理1.随机信号的分析方法在信号与系统中,我们把信号分为确知信号和随机信号。
其中随机信号无确定的变化规律,需要用统计特新进行分析。
这里我们引入随机过程的概念,所谓随机过程就是随机变量的集合,每个随机变量都是随机过程的一个取样序列。
随机过程的统计特性一般采用随机过程的分布函数和概率密度来描述,他们能够对随机过程作完整的描述。
但由于在实践中难以求得,在工程技术中,一般采用描述随机过程的主要平均统计特性的几个函数,包括均值、方差、相关函数、频谱及功率谱密度等来描述它们。
本实验中算法都是一种估算法,条件是N要足够大。
2.微弱随机信号的检测及提取方法因为噪声总会影响信号检测的结果,所以信号检测是信号处理的重要内容之一,低信噪比下的信号检测是目前检测领域的热点,而强噪声背景下的微弱信号提取又是信号检测的难点。
噪声主要来自于检测系统本身的电子电路和系统外空间高频电磁场干扰等,通常从以下两种不同途径来解决①降低系统的噪声,使被测信号功率大于噪声功率。
②采用相关接受技术,可以保证在信号功率小于噪声功率的情况下,人能检测出信号。
对微弱信号的检测与提取有很多方法,常用的方法有:自相关检测法、多重自相法、双谱估计理论及算法、时域方法、小波算法等。
对微弱信号检测与提取有很多方法,本实验采用多重自相关法。
多重自相关法是在传统自相关检测法的基础上,对信号的自相关函数再多次做自相关。
即令:式中,是和的叠加;是和的叠加。
对比两式,尽管两者信号的幅度和相位不同,但频率却没有变化。
信号经过相关运算后增加了信噪比,但其改变程度是有限的,因而限制了检测微弱信号的能力。
多重相关法将当作x(t),重复自相关函数检测方法步骤,自相关的次数越多,信噪比提高的越多,因此可检测出强噪声中的微弱信号。
随机信号处理实验报告
随机信号处理实验报告院系名称学生姓名学号指导教师目录一、实验要求: (3)二、实验原理: (3)2.1 随机信号的分析方法 (3)2.2 随机过程的频谱 (3)2.3 随机过程的相关函数和功率谱 (4)(1)随机信号的相关函数: (4)(2)随机信号的功率谱 (4)三、实验步骤与分析 (5)3.1实验方案 (5)3.2实验步骤与分析 (5)任务一:(s1 变量)求噪声下正弦信号的振幅和频率 (5)任务二:(s1 变量)求噪声下正弦信号的相位 (8)任务三:(s1 变量)求信号自相关函数和功率谱 (11)任务四:(s变量)求噪声下信号的振幅和频率 (14)任务五:(s变量)求信号的自相关函数和功率谱 (17)3.3实验结果与误差分析 (19)(1)实验结果 (19)(2)结果验证 (19)(3)误差分析 (21)四、实验总结和感悟 (22)1、实验总结 (22)2、实验感悟 (23)五、附低通滤波器的Matlab程序 (23)一、实验要求:(学号末尾3,7)两个数据文件,第一个文件数据中只包含一个正弦波,通过MA TLAB 仿真计算信号频谱和功率谱来估计该信号的幅度,功率,频率和相位?对第二个文件数据估计其中正弦波的幅度,功率和频率?写出报告,包含理论分析,仿真程序及说明,误差精度分析等。
第一文件调用格式load FileDat01_1 s1,数据在变量s1中;第二文件调用格式load FileDat01_2 s ,数据在变量s 中。
二、实验原理:2.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的增大而增大。
北京理工大学随机信号分析实验报告
北京理工大学随机信号分析实验报告本科实验报告实验名称:随机信号分析实验实验一随机序列的产生及数字特征估计一、实验目的1、学习和掌握随机数的产生方法。
2、实现随机序列的数字特征估计。
二、实验原理1、随机数的产生随机数指的是各种不同分布随机变量的抽样序列(样本值序列)。
进行随机信号仿真分析时,需要模拟产生各种分布的随机数。
在计算机仿真时,通常利用数学方法产生随机数,这种随机数称为伪随机数。
伪随机数是按照一定的计算公式产生的,这个公式称为随机数发生器。
伪随机数本质上不是随机的,而且存在周期性,但是如果计算公式选择适当,所产生的数据看似随机的,与真正的随机数具有相近的统计特性,可以作为随机数使用。
(0,1)均匀分布随机数是最最基本、最简单的随机数。
(0,1)均匀分布指的是在[0,1]区间上的均匀分布,即 U(0,1)。
实际应用中有许多现成的随机数发生器可以用于产生(0,1)均匀分布随机数,通常采用的方法为线性同余法,公式如下:)(m od ,110N ky y y n n -=Ny x n n /=序列{}nx 为产生的(0,1)均匀分布随机数。
下面给出了上式的3组常用参数: 1、10N 10,k 7==,周期7510≈⨯;2、(IBM 随机数发生器)3116N 2,k 23,==+周期8510≈⨯;3、(ran0)315N 21,k 7,=-=周期9210≈⨯;由均匀分布随机数,可以利用反函数构造出任意分布的随机数。
定理 1.1 若随机变量 X 具有连续分布函数F X (x),而R 为(0,1)均匀分布随机变量,则有)(1R F X x -=由这一定理可知,分布函数为F X (x)的随机数可以由(0,1)均匀分布随机数按上式进行变换得到。
2、MATLAB 中产生随机序列的函数(1)(0,1)均匀分布的随机序列函数:rand用法:x = rand(m,n)功能:产生m×n 的均匀分布随机数矩阵。
南京理工大学随机信号处理实验报告
题目:雷达线性调频信号的脉冲压缩处理线性调频脉冲信号,时宽10us,带宽40MHz,对该信号进行匹配滤波后,即脉压处理,处理增益为多少,脉压后的脉冲宽度为多少?用图说明脉压后的脉冲宽度,内差点看4dB带宽,以该带宽说明距离分辨率与带宽的对应关系。
1.程序为:T=10e-6;B=112e6;Rmin=8500;Rmax=11500;R=[9000,10000,10020];RCS=[1 1 1 ];C=3e8;K=B/T;Rwid=Rmax-Rmin;Twid=2*Rwid/C;Fs=10*B;Ts=1/Fs;Nwid=ceil(Twid/Ts);t=linspace(2*Rmin/C,2*Rmax/C,Nwid); M=length(R);td=ones(M,1)*t-2*R'/C*ones(1,Nwid);Srt1=RCS*(exp(1i*pi*K*td.^2).*(abs(td)<T/2));Srt=Srt1;Nchirp=ceil(T/Ts);Nfft=2^nextpow2(Nwid+Nwid-1); Srw=fft(Srt,Nfft);Srw1=fft(Srt1,Nfft);t0=linspace(-T/2,T/2,Nchirp);St=exp(1i*pi*K*t0.^2);Sw=fft(St,Nfft);Sot=fftshift(ifft(Srw.*conj(Sw)));Sot1=fftshift(ifft(Srw1.*conj(Sw)));N0=Nfft/2-Nchirp/2;Z=abs(Sot(N0:N0+Nwid-1));Z=Z/max(Z);Z=20*log10(Z+1e-6);figuresubplot(211)plot(t*1e6,real(Srt));axis tight;xlabel('us');ylabel('幅度')title(['线性信号压缩前']);subplot(212)plot(t*C/2,Z)xlabel('Range in meters');ylabel('幅度 ')title(['线性信号压缩后']);选取0.9*10^4HZ 的一个脉冲进行放大分析(调整Y 轴与X 轴的范围)58606264666870727476us幅度线性调频信号压缩前0.850.90.9511.05 1.1 1.15x 104-150-100-5050Range in meters 幅度 线性调频信号压缩后选取主瓣调整:大致可以看出压缩后的带宽为0.1hz理论上分析处理增益为:D=10*10e -6*112*10e6=1120D=112/B1=1120.B1=0.1HZ2.分辩率。
大连理工大学随机信号实验报告完整
大连理工大学实验预习报告学院(系):信息与通信工程学院专业:电子信息工程班级:姓名:学号:组:___实验时间:2015.12.14 实验室:C221 实验台:指导教师:实验I:随机信号的产生、相关分析及其应用实验实验1 均匀分布随机数的产生,统计特性分析及计算机仿真一、实验目的和要求掌握均匀分布随机信号的基本产生方法二、实验原理和内容较简单的伪随机序列产生方法是采用数论中基于数环理论的线性同余法(乘同余法、混合同余法),其迭代公式的一般形式为f(x) = (r·x+ b) Mod M,其离散形式为s(n + 1) = [r·s(n) + b] Mod M。
其中,s(n)为n 时刻的随机数种子,r 为扩展因子,b 为固定扰动项,M 为循环模,Mod M 表示对M 取模。
为保证s(n)的周期为M,r 的取值应满足r = 4k + 1,M 2p,k 与p的选取应满足:r < M,r(M-1) + 1< 231-1。
通常公式中参数常用取值为s(0) =12357,r = 2045,b = 1,M =1048576。
三、实验步骤1. 编程实现产生10000个在(0, 1)区间均匀分布随机数。
2. 计算生成随机数的1~4阶矩,最大值,最小值,频度直方图。
实验2 高斯分布随机数的产生,统计特性分析及计算机仿真一、实验目的和要求掌握高斯白噪声的基本产生方法二、实验原理和内容1.变换法2.较简单的高斯白噪声产生方法是基于概率论中的中心极限定理。
即无穷多个同分布随机变量之和构成随机变量服从高斯分布。
方便起见,可以使用多个均匀分布随机变量之和近似高斯分布随机变量。
三、实验步骤1.编程实现产生10000 个N(3, 4) 高斯随机数。
2.计算生成随机数的1~4 阶矩,最大值,最小值,频度直方图。
实验3 随机信号相关函数计算、相关分析及计算机仿真一、实验目的和要求掌握随机信号相关函数计算、相关分析及实现二、实验原理和内容根据自相关和互相关的定义,自相关,互相关计算随机信号的自相关和互相关。
北理工随机信号分析实验报告
本科实验报告实验名称:随机信号分析实验实验一 随机序列的产生及数字特征估计一、实验目的1、学习和掌握随机数的产生方法。
2、实现随机序列的数字特征估计。
二、实验原理1、随机数的产生随机数指的是各种不同分布随机变量的抽样序列(样本值序列)。
进行随机信号仿真分析时,需要模拟产生各种分布的随机数。
在计算机仿真时,通常利用数学方法产生随机数,这种随机数称为伪随机数。
伪随机数是按照一定的计算公式产生的,这个公式称为随机数发生器。
伪随机数本质上不是随机的,而且存在周期性,但是如果计算公式选择适当,所产生的数据看似随机的,与真正的随机数具有相近的统计特性,可以作为随机数使用。
(0,1)均匀分布随机数是最最基本、最简单的随机数。
(0,1)均匀分布指的是在[0,1]区间上的均匀分布,即 U(0,1)。
实际应用中有许多现成的随机数发生器可以用于产生(0,1)均匀分布随机数,通常采用的方法为线性同余法,公式如下:)(m od ,110N ky y y n n -=N y x n n /=序列{}n x 为产生的(0,1)均匀分布随机数。
下面给出了上式的3组常用参数: 1、10N 10,k 7==,周期7510≈⨯;2、(IBM 随机数发生器)3116N 2,k 23,==+周期8510≈⨯; 3、(ran0)315N 21,k 7,=-=周期9210≈⨯;由均匀分布随机数,可以利用反函数构造出任意分布的随机数。
定理 1.1 若随机变量 X 具有连续分布函数F X (x),而R 为(0,1)均匀分布随机变量,则有)(1R F X x -=由这一定理可知,分布函数为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.熟悉几种常用的时间序列。
二、实验内容1.白噪声通过线性系统后的统计特性分析。
(1)白噪声通过低通系统后的统计特性变化:对比输入输出的波形,自相关函数,功率谱密度,功率,互相关函数等;(2)白噪声通过不同带宽的低通系统后的概率密度;(3)窄带随机过程的产生与特性分析。
(调制,滤波)2.典型时间序列模型分析。
(1)模拟产生AR,ARMA模型序列,画出波形,并估计其均值,方差,自相关函数,功率谱密度;*(2)模拟产生指定功率密度的正态随机序列。
三、实验设备Matlab软件四、实验步骤以及实验结果分析1.白噪声通过线性系统后的统计特性分析。
>>l=(0:length(a2)-1)*200/length(a>>l=(0:length(a2)-1)*200/length(a2.典型时间序列的模拟分析模拟产生AR,ARMA模型序列:五、实验收获(本次实验的感受,对你的哪方面技能或知识有提高。
)本次实验我们收获很多,不仅理解了白噪声通过线性系统后统计特性的变化规律,同时也熟悉了如何使用matlab求信号的波形,自相关函数,功率谱密度,功率,互相关函数等等的统计特性。
深刻地理解到了线性系统对白噪声的影响。
除此之外,我们也深入地了解了AR 和ARMA模型序列。
最重要的是让我们加深了对课本知识的理解。
总之,本次实验我们受益匪浅。
随机信号分析实验报告范文
随机信号分析实验报告范文HaarrbbiinnIInnttiittuutteeooffTTeecchhnnoollooggyy实验报告告课程名称:院系:电子与信息工程学院班级:姓名:学号:指导教师:实验时间:实验一、各种分布随机数得产生(一)实验原理1、、均匀分布随机数得产生原理产生伪随机数得一种实用方法就是同余法,它利用同余运算递推产生伪随机数序列.最简单得方法就是加同余法为了保证产生得伪随机数能在[0,1]内均匀分布,需要M为正整数,此外常数c与初值y0亦为正整数。
加同余法虽然简单,但产生得伪随机数效果不好。
另一种同余法为乘同余法,它需要两次乘法才能产生一个[0,1]上均匀分布得随机数ﻩﻩﻩ式中,a为正整数。
用加法与乘法完成递推运算得称为混合同余法,即ﻩﻩﻩ用混合同余法产生得伪随机数具有较好得特性,一些程序库中都有成熟得程序供选择。
常用得计算语言如Baic、C与Matlab都有产生均匀分布随机数得函数可以调用,只就是用各种编程语言对应得函数产生得均匀分布随机数得范围不同,有得函数可能还需要提供种子或初始化。
Matlab提供得函数rand()可以产生一个在[0,1]区间分布得随机数,rand(2,4)则可以产生一个在[0,1]区间分布得随机数矩阵,矩阵为2行4列。
Matlab提供得另一个产生随机数得函数就是random(’unif’,a,b,N,M),unif表示均匀分布,a与b就是均匀分布区间得上下界,N与M分别就是矩阵得行与列。
2、、随机变量得仿真根据随机变量函数变换得原理,如果能将两个分布之间得函数关系用显式表达,那么就可以利用一种分布得随机变量通过变换得到另一种分布得随机变量。
若X就是分布函数为F(某)得随机变量,且分布函数F(某)为严格单调升函数,令Y=F(某),则Y必为在[0,1]上均匀分布得随机变量.反之,若Y就是在[0,1]上均匀分布得随机变量,那么即就是分布函数为F某(某)得随机变量。
随机信号分析报告实验:随机过程通过线性系统地分析报告
实验三 随机过程通过线性系统的分析实验目的1. 理解和分析白噪声通过线性系统后输出的特性。
2. 学习和掌握随机过程通过线性系统后的特性,验证随机过程的正态化问题。
实验原理1.白噪声通过线性系统设连续线性系统的传递函数为)(ωH 或)(s H ,输入白噪声的功率谱密度为2)(0N S X =ω,那么系统输出的功率谱密度为2)()(02N H S Y ⋅=ωω (3.1) 输出自相关函数为⎰∞∞-=ωωπτωτd e H N R j Y 20)(4)( (3.2)输出相关系数为)0()()(Y Y Y R R ττγ=(3.3) 输出相关时间为⎰∞=00)(ττγτd Y (3.4)输出平均功率为[]⎰∞=202)(2)(ωωπd H N t Y E (3.5)上述式子表明,若输入端是具有均匀谱的白噪声,则输出端随机信号的功率谱主要由系统的幅频特性)(ωH 决定,不再是常数。
2.等效噪声带宽在实际中,常常用一个理想系统等效代替实际系统的)(ωH ,因此引入了等效噪声带宽的概念,他被定义为理想系统的带宽。
等效的原则是,理想系统与实际系统在同一白噪声的激励下,两个系统的输出平均功率相等,理想系统的增益等于实际系统的最大增益。
实际系统的等效噪声带宽为⎰∞=∆022max)()(1ωωωωd H H e (3.6)或⎰∞∞--=∆j j e ds s H s H H j )()()(212maxωω (3.7)3.线性系统输出端随机过程的概率分布 (1)正态随机过程通过线性系统若线性系统输入为正态过程,则该系统输出仍为正态过程。
(2)随机过程的正态化随机过程的正态化指的是,非正态随机过程通过线性系统后变换为正态过程。
任意分布的白噪声通过线性系统后输出是服从正态分布的;宽带噪声通过窄带系统,输出近似服从正态分布。
实验内容设白噪声通过图3.1所示的RC 电路,分析输出的统计特性。
图3.1 RC 电路(1)试推导系统输出的功率谱密度、相关函数、相关时间和系统的等效噪声带宽。
随机信号分析与处理实验报告
随机信号分析与处理实验题目:对音频信号的随机处理班级:0312412姓名:肖文洲学号:031241217指导老师:钱楷时间:2014年11月25日实验目的:1、学会利用MATLAB模拟产生各类随机序列。
2、熟悉和掌握随机信号数字特征估计的基本方法。
3、熟悉掌握MATLAB的函数及函数调用、使用方法。
4、学会在MATLAB中创建GUI文件。
实验内容:1、选用任意一个音频信号作为实验对象,进行各种操作并画出信号和波形。
2、操作类型:(1)、概率密度;(2)、希尔伯特变换;(3)、误差函数;(4)、randn;(5)、原始信号频谱;(6)、axis;(7)、原始信号;(8)、normpdf;(9)、unifpdf;(10)、unifcdf;(11)、raylpdf;(12)、raylcdf;(13)、exppdf;(14)、截取声音信号的频谱;(15)、expcdf;(16)、periodogram;(17)、weibrnd;(18)、rand;(19)、自相关函数;(20)、截取信号的均方值。
实验步骤:1、打开MATLAB软件,然后输入guide创建一个GUI文件。
2、在已经创建好的GUI文件里面穿件所需要的.fig面板(以学号姓名格式命名)。
入下图所示:图为已经创建好的.fig面板3、右击“概率密度”,查看回调,然后点击“callback”.在相应的位置输入程序。
然后点击运行,出现下图:4、依次对后续操作方式进行类似的操作。
5、当完成所有按键的“callback”后,出现的均为上图。
实验程序:function varargout = xiaowenzhou(varargin)% XIAOWENZHOU M-file for xiaowenzhou.fig% XIAOWENZHOU, by itself, creates a new XIAOWENZHOU or raises the existing% singleton*.%% H = XIAOWENZHOU returns the handle to a new XIAOWENZHOU or the handle to% the existing singleton*.%% XIAOWENZHOU('CALLBACK',hObject,eventData,handles,...) calls the local% function named CALLBACK in XIAOWENZHOU.M with the given input arguments.%% XIAOWENZHOU('Property','Value',...) creates a new XIAOWENZHOU or raises the% existing singleton*. Starting from the left, property value pairs are% applied to the GUI before xiaowenzhou_OpeningFunction gets called. An% unrecognized property name or invalid value makes property application% stop. All inputs are passed to xiaowenzhou_OpeningFcn via varargin.%% *See GUI Options on GUIDE's Tools menu. Choose "GUI allows only one% instance to run (singleton)".%% See also: GUIDE, GUIDATA, GUIHANDLES% Edit the above text to modify the response to help xiaowenzhou% Last Modified by GUIDE v2.5 02-Dec-2014 23:14:41% Begin initialization code - DO NOT EDITgui_Singleton = 1;gui_State = struct('gui_Name', mfilename, ...'gui_Singleton', gui_Singleton, ...'gui_OpeningFcn', @xiaowenzhou_OpeningFcn, ...'gui_OutputFcn', @xiaowenzhou_OutputFcn, ...'gui_LayoutFcn', [] , ...'gui_Callback', []);if nargin && ischar(varargin{1})gui_State.gui_Callback = str2func(varargin{1});endif nargout[varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:});elsegui_mainfcn(gui_State, varargin{:});end% End initialization code - DO NOT EDIT% --- Executes just before xiaowenzhou is made visible.function xiaowenzhou_OpeningFcn(hObject, eventdata, handles, varargin)% This function has no output args, see OutputFcn.% hObject handle to figure% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)% varargin command line arguments to xiaowenzhou (see VARARGIN)% Choose default command line output for xiaowenzhouhandles.output = hObject;% Update handles structureguidata(hObject, handles);% UIWAIT makes xiaowenzhou wait for user response (see UIRESUME)% uiwait(handles.figure1);% --- Outputs from this function are returned to the command line. function varargout = xiaowenzhou_OutputFcn(hObject, eventdata, handles) % varargout cell array for returning output args (see VARARGOUT);% hObject handle to figure% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)% Get default command line output from handles structurevarargout{1} = handles.output;% --- Executes on button press in pushbutton1.function pushbutton1_Callback(hObject, eventdata, handles)[y,Fs,bits]=wavread('Íõ·Æ.wav');y=y(1:5000);N=length(y);i=1:N;[f,i]=ksdensity(y);plot(i,f);grid;xlabel('x');ylabel('f(x)');axis();title('¸ÅÂÊÃܶÈ');% hObject handle to pushbutton1 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)% --- Executes on button press in pushbutton2.function pushbutton2_Callback(hObject, eventdata, handles)x=wavread('Íõ·Æ.wav');x=x(20000:40000);y=hilbert(x);y=real(y);plot(x);% hObject handle to pushbutton2 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA)% --- Executes on button press in pushbutton3.function pushbutton3_Callback(hObject, eventdata, handles)x=wavread('Íõ·Æ.wav');x=randn(500,1);plot(x);% hObject handle to pushbutton3 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA)% --- Executes on button press in pushbutton4.function pushbutton4_Callback(hObject, eventdata, handles)x=wavread('Íõ·Æ.wav');y=erf(x);plot(y);% hObject handle to pushbutton4 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA)% --- Executes on button press in pushbutton5.function pushbutton5_Callback(hObject, eventdata, handles)x=wavread('Íõ·Æ.wav');plot(x);axis([0 5000 -0.01 0.01]);% hObject handle to pushbutton5 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA)% --- Executes on button press in pushbutton6.x=wavread('Íõ·Æ.wav');x=x(20000:40000);plot(x);% hObject handle to pushbutton6 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA)% --- Executes on button press in pushbutton7.function pushbutton7_Callback(hObject, eventdata, handles)x=wavread('Íõ·Æ.wav');x=-6:0.01:7;y=normpdf(x,1,2);plot(y);% hObject handle to pushbutton7 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA)% --- Executes on button press in pushbutton8.function pushbutton8_Callback(hObject, eventdata, handles)x=wavread('Íõ·Æ.wav');x=0:0.1:5;y=unifpdf(x,1,30);plot(y);% hObject handle to pushbutton8 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA)% --- Executes on button press in pushbutton9.function pushbutton9_Callback(hObject, eventdata, handles)x=wavread('Íõ·Æ.wav');x=0:0.1:5;y=unifcdf(x,1,5);plot(y);% hObject handle to pushbutton9 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA)% --- Executes on button press in pushbutton10.x=wavread('Íõ·Æ.wav');x=0:0.1:3;y=raylpdf(x,2);plot(y);% hObject handle to pushbutton10 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA)% --- Executes on button press in pushbutton11.function pushbutton11_Callback(hObject, eventdata, handles)x=wavread('Íõ·Æ.wav');x=0:0.1:3;y=raylcdf(x,10);plot(y);% hObject handle to pushbutton11 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA)% --- Executes on button press in pushbutton12.function pushbutton12_Callback(hObject, eventdata, handles)x=wavread('Íõ·Æ.wav');x=0:0.01:5;y=exppdf(x,1);plot(y);% hObject handle to pushbutton12 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA)% --- Executes on button press in pushbutton13.function pushbutton13_Callback(hObject, eventdata, handles)x=wavread('Íõ·Æ.wav');x=0:0.01:5;y=expcdf(x,1);plot(y);% hObject handle to pushbutton13 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA)% --- Executes on button press in pushbutton14.[y,Fs,bits]=wavread('Íõ·Æ.wav');y1=y(1:1000);t=0:1/Fs:1;y1=periodogram(y1,[],1000,Fs);plot(y1);% hObject handle to pushbutton14 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA)% --- Executes on button press in pushbutton15.function pushbutton15_Callback(hObject, eventdata, handles)x=wavread('Íõ·Æ.wav');x=0:0.01:5;x=weibrnd(1,1.5,100,1);plot(x);% hObject handle to pushbutton15 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA)% --- Executes on button press in pushbutton16.function pushbutton16_Callback(hObject, eventdata, handles)x=wavread('Íõ·Æ.wav');x=rand(200,1);plot(x);% hObject handle to pushbutton16 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA)% --- Executes on button press in pushbutton17.function pushbutton17_Callback(hObject, eventdata, handles)[x,Fs,bits]=wavread ('Íõ·Æ.wav');x=x (:,1);X=fft (x,4096);magX=abs (X);angX=angle (X);plot (X); title ('Ô-ʼÐźÅƵÆ×');% hObject handle to pushbutton17 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA)% --- Executes on button press in pushbutton18.function pushbutton18_Callback(hObject, eventdata, handles)[y,Fs,bits]=wavread('Íõ·Æ.wav');y=y(500:1000);h=[ones(1,20) zeros(1,20)];y2=conv(h,y);stem(y2,'.');grid;title('½ØÈ¡ÉùÒôÐźŵľí»ý');% hObject handle to pushbutton18 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA)% --- Executes on button press in pushbutton19.function pushbutton19_Callback(hObject, eventdata, handles)[y,Fs,bits]=wavread('Íõ·Æ.wav');y=y(1:5000);R=xcorr(y);plot(R);grid;title('×ÔÏà¹Øº¯Êý');% hObject handle to pushbutton19 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA)% --- Executes on button press in pushbutton20.function pushbutton20_Callback(hObject, eventdata, handles)[y,Fs,bits]=wavread('Íõ·Æ.wav');y=y(1:1000);n=length(y);x=randn(50,n);square=zeros(1,50);for i=1:50for j=1:1000square(i)=square(i)+x(i,j).^2;endsquare(i)=square(i)/1000;endRMS=sum(square)/30;plot(square);grid;title('½ØÈ¡ÉùÒôÐźŵľù·½Öµ');% hObject handle to pushbutton20 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)% --- Executes during object creation, after setting all properties. function pushbutton1_CreateFcn(hObject, eventdata, handles)% hObject handle to pushbutton1 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles empty - handles not created until after all CreateFcns called 实验结果:(1)概率密度(2)希尔伯特变换(3)randn(4)误差函数(5)axis(6)原始信号(7)normpdfd(8)unifpdf(9)unifcdf(10)raylpdf(11)raylcdf(12)exppdf(13)截取声音信号的卷积(14)expcdf(15)periodogram(16)weibrnd(17)rand(18)原始信号频谱(19)自相关函数(20)截取信号的均方值实验总结:随机信号分析与处理是研究随机信号的特点及其处理方法的专业基础课,是目标检测、估计、滤波等信号处理理论的基础。
哈尔滨工业大学(威海)随机信号分析实验一报告
《随机信号分析》实验报告班级: 1302502学号:姓名:《随机信号分析》实验报告实验一一、实验目的:熟悉并练习使用随机信号Matlab的函数二、实验内容:1、熟悉并练习使用下列Matlab 的函数,给出各个函数的功能说明和内部参数的意义,并给出至少一个使用例子和运行结果:1)rand() 11)unifpdf()2)randn() 12)unifcdf()3)normrnd() 13)raylpdf()4)mean() 14)raylcdf()5)var() 15)exppdf()6)xcorr() 16)expcdf()7)periodogram() 17)chol()8)fft()18)ksdensity()9)normpdf() 19)hist()10)normcdf() 20)int()用法、功能、程序如下:1)randn(m,n)功能:返回一个从标准正态分布中得到的伪随机标量。
>> r = randn(5) %由标准正态分布随机数组成的5×5 矩阵。
r =-1.0689 -0.7549 0.3192 0.6277 -1.2141-0.8095 1.3703 0.3129 1.0933 -1.1135-2.9443 -1.7115 -0.8649 1.1093 -0.00681.4384 -0.1022 -0.0301 -0.8637 1.53260.3252 -0.2414 -0.1649 0.0774 -0.76972)rand(m,n)功能:返回一个从开区间(0,1) 上的标准均匀分布得到的伪随机标量。
r = rand(5) %生成一个由介于0 和1 之间的均匀分布的随机数组成的5×5 矩阵>>r =0.5469 0.9572 0.9157 0.8491 0.39220.9575 0.4854 0.7922 0.9340 0.65550.9649 0.8003 0.9595 0.6787 0.17120.1576 0.1419 0.6557 0.7577 0.70600.9706 0.4218 0.0357 0.7431 0.03183)normrnd(mu,sigma,m,n)功能:以均值μ和标准差σ为参数的正态分布随机数mxn>> normrnd(0,1,3,4) %生成均值μ=0,σ=1的3x4正态分布随机数ans =0.2761 0.3919 -0.7411 0.0125-0.2612 -1.2507 -0.5078 -3.02920.4434 -0.9480 -0.3206 -0.45704)mean(A,dim)功能:数组的平均值mean(A,dim) dim=1,返回列平均数,默认为1dim=2,返回列平均数dim>2,返回AA = [0 1 1; 2 3 2; 1 3 2; 4 2 2] %M = mean(A) 沿A 的大小不等于1 的第一个数组维度返回均值。
四川大学随机信号分析实验报告
随机信号分析基础实验报告课程随机信号分析基础实验题目随机信号通过线性系统学生姓名笔墨东韵专业电子信息科学与技术一、实验目的1.理解白噪声通过线性系统后统计特性的变化规律。
2.熟悉几种常用的时间序列。
二、实验内容1.白噪声通过线性系统后的统计特性分析。
(1)白噪声通过低通系统后的统计特性变化:对比输入输出的波形,自相关函数,功率谱密度,功率,互相关函数等;(2)白噪声通过不同带宽的低通系统后的概率密度;(3)窄带随机过程的产生与特性分析。
(调制,滤波)2.典型时间序列模型分析。
(1)模拟产生AR,ARMA模型序列,画出波形,并估计其均值,方差,自相关函数,功率谱密度;*(2)模拟产生指定功率密度的正态随机序列。
三、实验设备Matlab软件四、实验步骤以及实验结果分析1.白噪声通过线性系统后的统计特性分析。
>>l=(0:length(a2)-1)*200/length(a>>l=(0:length(a2)-1)*200/length(a2.典型时间序列的模拟分析模拟产生AR,ARMA模型序列:五、实验收获(本次实验的感受,对你的哪方面技能或知识有提高。
)本次实验我们收获很多,不仅理解了白噪声通过线性系统后统计特性的变化规律,同时也熟悉了如何使用matlab求信号的波形,自相关函数,功率谱密度,功率,互相关函数等等的统计特性。
深刻地理解到了线性系统对白噪声的影响。
除此之外,我们也深入地了解了AR 和ARMA模型序列。
最重要的是让我们加深了对课本知识的理解。
总之,本次实验我们受益匪浅。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
随机信号分析与处理实验题目:对音频信号的随机处理班级:0312412姓名:肖文洲学号:031241217指导老师:钱楷时间:2014年11月25日实验目的:1、学会利用MATLAB模拟产生各类随机序列。
2、熟悉和掌握随机信号数字特征估计的基本方法。
3、熟悉掌握MATLAB的函数及函数调用、使用方法。
4、学会在MATLAB中创建GUI文件。
实验内容:1、选用任意一个音频信号作为实验对象,进行各种操作并画出信号和波形。
2、操作类型:(1)、概率密度;(2)、希尔伯特变换;(3)、误差函数;(4)、randn;(5)、原始信号频谱;(6)、axis;(7)、原始信号;(8)、normpdf;(9)、unifpdf;(10)、unifcdf;(11)、raylpdf;(12)、raylcdf;(13)、exppdf;(14)、截取声音信号的频谱;(15)、expcdf;(16)、periodogram;(17)、weibrnd;(18)、rand;(19)、自相关函数;(20)、截取信号的均方值。
实验步骤:1、打开MATLAB软件,然后输入guide创建一个GUI文件。
2、在已经创建好的GUI文件里面穿件所需要的.fig面板(以学号姓名格式命名)。
入下图所示:图为已经创建好的.fig面板3、右击“概率密度”,查看回调,然后点击“callback”.在相应的位置输入程序。
然后点击运行,出现下图:4、依次对后续操作方式进行类似的操作。
5、当完成所有按键的“callback”后,出现的均为上图。
实验程序:function varargout = xiaowenzhou(varargin)% XIAOWENZHOU M-file for xiaowenzhou.fig% XIAOWENZHOU, by itself, creates a new XIAOWENZHOU or raises the existing% singleton*.%% H = XIAOWENZHOU returns the handle to a new XIAOWENZHOU or the handle to% the existing singleton*.%% XIAOWENZHOU('CALLBACK',hObject,eventData,handles,...) calls the local% function named CALLBACK in XIAOWENZHOU.M with the given input arguments.%% XIAOWENZHOU('Property','Value',...) creates a new XIAOWENZHOU or raises the% existing singleton*. Starting from the left, property value pairs are% applied to the GUI before xiaowenzhou_OpeningFunction gets called. An% unrecognized property name or invalid value makes property application% stop. All inputs are passed to xiaowenzhou_OpeningFcn via varargin.%% *See GUI Options on GUIDE's Tools menu. Choose "GUI allows only one% instance to run (singleton)".%% See also: GUIDE, GUIDATA, GUIHANDLES% Edit the above text to modify the response to help xiaowenzhou% Last Modified by GUIDE v2.5 02-Dec-2014 23:14:41% Begin initialization code - DO NOT EDITgui_Singleton = 1;gui_State = struct('gui_Name', mfilename, ...'gui_Singleton', gui_Singleton, ...'gui_OpeningFcn', @xiaowenzhou_OpeningFcn, ...'gui_OutputFcn', @xiaowenzhou_OutputFcn, ...'gui_LayoutFcn', [] , ...'gui_Callback', []);if nargin && ischar(varargin{1})gui_State.gui_Callback = str2func(varargin{1});endif nargout[varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:});elsegui_mainfcn(gui_State, varargin{:});end% End initialization code - DO NOT EDIT% --- Executes just before xiaowenzhou is made visible.function xiaowenzhou_OpeningFcn(hObject, eventdata, handles, varargin)% This function has no output args, see OutputFcn.% hObject handle to figure% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)% varargin command line arguments to xiaowenzhou (see VARARGIN)% Choose default command line output for xiaowenzhouhandles.output = hObject;% Update handles structureguidata(hObject, handles);% UIWAIT makes xiaowenzhou wait for user response (see UIRESUME)% uiwait(handles.figure1);% --- Outputs from this function are returned to the command line. function varargout = xiaowenzhou_OutputFcn(hObject, eventdata, handles) % varargout cell array for returning output args (see VARARGOUT);% hObject handle to figure% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)% Get default command line output from handles structurevarargout{1} = handles.output;% --- Executes on button press in pushbutton1.function pushbutton1_Callback(hObject, eventdata, handles)[y,Fs,bits]=wavread('Íõ·Æ.wav');y=y(1:5000);N=length(y);i=1:N;[f,i]=ksdensity(y);plot(i,f);grid;xlabel('x');ylabel('f(x)');axis();title('¸ÅÂÊÃܶÈ');% hObject handle to pushbutton1 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)% --- Executes on button press in pushbutton2.function pushbutton2_Callback(hObject, eventdata, handles)x=wavread('Íõ·Æ.wav');x=x(20000:40000);y=hilbert(x);y=real(y);plot(x);% hObject handle to pushbutton2 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA)% --- Executes on button press in pushbutton3.function pushbutton3_Callback(hObject, eventdata, handles)x=wavread('Íõ·Æ.wav');x=randn(500,1);plot(x);% hObject handle to pushbutton3 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA)% --- Executes on button press in pushbutton4.function pushbutton4_Callback(hObject, eventdata, handles)x=wavread('Íõ·Æ.wav');y=erf(x);plot(y);% hObject handle to pushbutton4 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA)% --- Executes on button press in pushbutton5.function pushbutton5_Callback(hObject, eventdata, handles)x=wavread('Íõ·Æ.wav');plot(x);axis([0 5000 -0.01 0.01]);% hObject handle to pushbutton5 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA)% --- Executes on button press in pushbutton6.x=wavread('Íõ·Æ.wav');x=x(20000:40000);plot(x);% hObject handle to pushbutton6 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA)% --- Executes on button press in pushbutton7.function pushbutton7_Callback(hObject, eventdata, handles)x=wavread('Íõ·Æ.wav');x=-6:0.01:7;y=normpdf(x,1,2);plot(y);% hObject handle to pushbutton7 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA)% --- Executes on button press in pushbutton8.function pushbutton8_Callback(hObject, eventdata, handles)x=wavread('Íõ·Æ.wav');x=0:0.1:5;y=unifpdf(x,1,30);plot(y);% hObject handle to pushbutton8 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA)% --- Executes on button press in pushbutton9.function pushbutton9_Callback(hObject, eventdata, handles)x=wavread('Íõ·Æ.wav');x=0:0.1:5;y=unifcdf(x,1,5);plot(y);% hObject handle to pushbutton9 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA)% --- Executes on button press in pushbutton10.x=wavread('Íõ·Æ.wav');x=0:0.1:3;y=raylpdf(x,2);plot(y);% hObject handle to pushbutton10 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA)% --- Executes on button press in pushbutton11.function pushbutton11_Callback(hObject, eventdata, handles)x=wavread('Íõ·Æ.wav');x=0:0.1:3;y=raylcdf(x,10);plot(y);% hObject handle to pushbutton11 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA)% --- Executes on button press in pushbutton12.function pushbutton12_Callback(hObject, eventdata, handles)x=wavread('Íõ·Æ.wav');x=0:0.01:5;y=exppdf(x,1);plot(y);% hObject handle to pushbutton12 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA)% --- Executes on button press in pushbutton13.function pushbutton13_Callback(hObject, eventdata, handles)x=wavread('Íõ·Æ.wav');x=0:0.01:5;y=expcdf(x,1);plot(y);% hObject handle to pushbutton13 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA)% --- Executes on button press in pushbutton14.[y,Fs,bits]=wavread('Íõ·Æ.wav');y1=y(1:1000);t=0:1/Fs:1;y1=periodogram(y1,[],1000,Fs);plot(y1);% hObject handle to pushbutton14 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA)% --- Executes on button press in pushbutton15.function pushbutton15_Callback(hObject, eventdata, handles)x=wavread('Íõ·Æ.wav');x=0:0.01:5;x=weibrnd(1,1.5,100,1);plot(x);% hObject handle to pushbutton15 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA)% --- Executes on button press in pushbutton16.function pushbutton16_Callback(hObject, eventdata, handles)x=wavread('Íõ·Æ.wav');x=rand(200,1);plot(x);% hObject handle to pushbutton16 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA)% --- Executes on button press in pushbutton17.function pushbutton17_Callback(hObject, eventdata, handles)[x,Fs,bits]=wavread ('Íõ·Æ.wav');x=x (:,1);X=fft (x,4096);magX=abs (X);angX=angle (X);plot (X); title ('Ô-ʼÐźÅƵÆ×');% hObject handle to pushbutton17 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA)% --- Executes on button press in pushbutton18.function pushbutton18_Callback(hObject, eventdata, handles)[y,Fs,bits]=wavread('Íõ·Æ.wav');y=y(500:1000);h=[ones(1,20) zeros(1,20)];y2=conv(h,y);stem(y2,'.');grid;title('½ØÈ¡ÉùÒôÐźŵľí»ý');% hObject handle to pushbutton18 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA)% --- Executes on button press in pushbutton19.function pushbutton19_Callback(hObject, eventdata, handles)[y,Fs,bits]=wavread('Íõ·Æ.wav');y=y(1:5000);R=xcorr(y);plot(R);grid;title('×ÔÏà¹Øº¯Êý');% hObject handle to pushbutton19 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA)% --- Executes on button press in pushbutton20.function pushbutton20_Callback(hObject, eventdata, handles)[y,Fs,bits]=wavread('Íõ·Æ.wav');y=y(1:1000);n=length(y);x=randn(50,n);square=zeros(1,50);for i=1:50for j=1:1000square(i)=square(i)+x(i,j).^2;endsquare(i)=square(i)/1000;endRMS=sum(square)/30;plot(square);grid;title('½ØÈ¡ÉùÒôÐźŵľù·½Öµ');% hObject handle to pushbutton20 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)% --- Executes during object creation, after setting all properties. function pushbutton1_CreateFcn(hObject, eventdata, handles)% hObject handle to pushbutton1 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles empty - handles not created until after all CreateFcns called 实验结果:(1)概率密度(2)希尔伯特变换(3)randn(4)误差函数(5)axis(6)原始信号(7)normpdfd(8)unifpdf(9)unifcdf(10)raylpdf(11)raylcdf(12)exppdf(13)截取声音信号的卷积(14)expcdf(15)periodogram(16)weibrnd(17)rand(18)原始信号频谱(19)自相关函数(20)截取信号的均方值实验总结:随机信号分析与处理是研究随机信号的特点及其处理方法的专业基础课,是目标检测、估计、滤波等信号处理理论的基础。