随机信号分析实验报告二 2

《随机信号分析》实验报告二

班级：

学号：

姓名：

实验二高斯噪声的产生和性能测试

１．实验目的

（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、源程序及功能注释

（2）源程序：

clear all;

clc;

t=0:320; %t=0:320

x=sin(2*pi*t/25); %x=sin(2*p1*t/25)

x1=wgn(1,321,0); %产生一个一行32列的高斯白噪声矩阵，输出的噪声强度为0dbw

z=x+x1; %z=x+x1;

y=trapz(t,z); %y=int(z,x,o,t),返回到从0到t的定积分z，积分变量为x

subplot(2,3,1),plot(z); %将图像窗口分为6个画图区，在第一个区域画图，以z为纵坐标

title('随机信号序列') %命名图像

meany=mean(z); %返回y中表示的平均值

subplot(2,3,2),plot(t,meany,'.'); %在第二个区域画图

title('随机信号均值') %命名

vary=var(y); %方差

subplot(2,3,3),plot(t,vary,'.'); %在第三个区域作图

title('随机信号方差') %命名

cory=xcorr(z,'unbiased');%自相关函数

subplot(2,3,4),plot(cory); %自相关函数

title('随机信号自相关函数') %命名

covv=cov(y);%协方差

subplot(2,3,5),plot(t,covv,'.'); %在第五个区域作图

title('随机信号协方差') %命名

（3）源程序：

t=[0:0.0005:0.045]; %t=[0:0.0005:0.045]

X1=sin(2*pi*25*t); %正弦

subplot(3,4,1); %将图像窗口分为12个画图区，在第一个区域画图，以z为纵坐标plot(t,X1);grid %以t为横坐标，X1为纵坐标，在画图中设置网格线

title('正弦函数序列');%命名

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);grid

title('混合随机信号'); %命名

meany1=mean(X1); %原信号的均值subplot(3,4,5),plot(t,meany1); %在第六个区域作图title('原信号均值'); %命名

vary1=var(X1); %原信号的方差subplot(3,4,6),plot(t,vary1); %在第6个区域作图title('原信号方差'); %命名

cory1=xcorr(X1,'unbiased'); %原信号的自相关函数subplot(3,4,7),plot(cory1); %在第七个区域作图title('原信号自相关函数'); %命名

meany=mean(X); %混合信号的均值subplot(3,4,9),plot(t,meany); %在第九个区域作图title('混合信号均值'); %命名

vary=var(X); %混合信号的方差subplot(3,4,10),plot(t,vary); %在第十个区域作图

title('混合信号方差') %命名

cory=xcorr(X,'unbiased'); %混合信号的自相关函数subplot(3,4,11),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('混合信号概率密度分布') %图像命名

（4）源程序：

[I,M]=imread('大殿.jpg');

J=imnoise(I,'gaussian',0,0.01);

subplot(2,1,1),imshow(I,M),title('大殿1');

subplot(2,1,2),imshow(J,M),title('大殿2');

K=zeros(242,308);

i=1:100;

J=imnoise(I,'gaussian',0,0.01);

J1=im2double(J);

K=K+J1;

endK=K/100;

（5）源程序：

clc;%清除所有变量

clear all;%清屏

t=[0:0.005:0.5]; %t为0-0.045的由0开始每隔0.0005一个点的数组

A=wavread('F:\好好学习\随机信号\实验\实验2\alert.wav '); %读取音频信号

subplot(3,3,1),plot(A); %所画图像共有3行3列，该图是第1个，画出A的图像title('音频信号'); %显示图像名称“音频信号”

noise=normrnd(2,0.2); %生成均值为2，方差为0.04的高斯噪声

x=A+noise; %将音频与噪声混合

subplot(3,3,3),plot(x); %所画图像共有3行3列，该图是第3个，画出x的图像title('混合信号'); %显示图像名称“混合信号”

meany=mean(x); %计算混合信号的均值

subplot(3,3,4),plot(t,meany); %所画图像共有3行3列，该图是第4个，画出mean-t的图像

title('混合信号均值'); %显示图像名称“混合信号均值”

vary=var(x); %计算混合过程的方差

subplot(3,3,5),plot(t,vary); %所画图像共有3行3列，该图是第5个，画出vary-t 的图像

title('混合信号方差'); %显示图像名称“混合信号方差”

cory=xcorr(x,'unbiased');%计算混合信号的自相关函数

subplot(3,3,6),plot(cory); %所画图像共有3行3列，该图是第6个，画出cory 的图像

title('混合信号自相关函数'); %显示图像名称“混合信号自相关函数”

covv=cov(x);%计算混合信号的协方差

subplot(3,3,7),plot(t,covv,'.'); %所画图像共有3行3列，该图是第7个，画出covv-t 的图像

title('混合信号协方差'); %显示图像名称“混合信号协方差”

fx=fft(x);%求x的离散傅里叶变换

ym=abs(x);%求x的绝对值

subplot(3,3,8),plot(ym);%所画图像共有3行3列，该图是第8个

title('混合信号频谱图')%显示标题为：混合信号频谱图

fz=fft(cory);%求cory的离散傅里叶变换

zm=abs(fz);%求fz的绝对值

subplot(3,3,9),plot(zm);%所画图像共有2行3列，该图是第6个,画出zm的图像

title('混合信号功率谱密度图')%显示标题为：混合信号功率谱密度图

5 、实验结果（2）结果：

（3）结果：

（4）结果：6、实验总结

北理工随机信号分析实验报告

本科实验报告实验名称：随机信号分析实验

实验一 随机序列的产生及数字特征估计 一、实验目的 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、10 N 10,k 7==，周期7 510≈?； 2、（IBM 随机数发生器）31 16 N 2,k 23,==+周期8 510≈?； 3、（ran0）31 5 N 21,k 7,=-=周期9 210≈?； 由均匀分布随机数，可以利用反函数构造出任意分布的随机数。 定理 1.1 若随机变量 X 具有连续分布函数F X (x)，而R 为(0,1)均匀分布随机变量，则有 )(1R F X x -= 由这一定理可知，分布函数为F X (x)的随机数可以由(0,1)均匀分布随机数按上式进行变

随机信号分析实验报告

一、实验名称 微弱信号的检测提取及分析方法 二、实验目的 1．了解随机信号分析理论如何在实践中应用 2．了解随机信号自身的特性，包括均值、方差、相关函数、频谱及功率谱密度等 3．掌握随机信号的检测及分析方法 三、实验原理 1．随机信号的分析方法 在信号与系统中，我们把信号分为确知信号和随机信号。其中随机信号无确定的变化规律，需要用统计特新进行分析。这里我们引入随机过程的概念，所谓随机过程就是随机变量的集合，每个随机变量都是随机过程的一个取样序列。 随机过程的统计特性一般采用随机过程的分布函数和概率密度来描述，他们能够对随机过程作完整的描述。但由于在实践中难以求得，在工程技术中，一般采用描述随机过程的主要平均统计特性的几个函数，包括均值、方差、相关函数、频谱及功率谱密度等来描述它们。本实验中算法都是一种估算法，条件是N要足够大。 2．微弱随机信号的检测及提取方法 因为噪声总会影响信号检测的结果，所以信号检测是信号处理的重要内容之一，低信噪比下的信号检测是目前检测领域的热点，而强噪声背景下的微弱信号提取又是信号检测的难点。 噪声主要来自于检测系统本身的电子电路和系统外空间高频电磁场干扰等，通常从以下两种不同途径来解决 ①降低系统的噪声，使被测信号功率大于噪声功率。 ②采用相关接受技术，可以保证在信号功率小于噪声功率的情况下，人能检测出信号。 对微弱信号的检测与提取有很多方法，常用的方法有：自相关检测法、多重自相法、双谱估计理论及算法、时域方法、小波算法等。 对微弱信号检测与提取有很多方法，本实验采用多重自相关法。 多重自相关法是在传统自相关检测法的基础上，对信号的自相关函数再多次做自相关。即令： 式中，是和的叠加；是和的叠加。对比两式，尽管两者信号的幅度和相位不同，但频率却没有变化。信号经过相关运算后增加了信噪比，但其改变程度是有限的，因而限制了检测微弱信号的能力。多重相关法将 当作x（t），重复自相关函数检测方法步骤，自相关的次数越多，信噪比提高的越多，因此可检测出强噪声中的微弱信号。

信号与系统实验报告1

学生实验报告 （理工类） 课程名称：信号与线性系统专业班级：M11通信工程 学生学号：1121413017 学生姓名：王金龙 所属院部：龙蟠学院指导教师：杨娟

20 11 ——20 12 学年第 1 学期 金陵科技学院教务处制 实验报告书写要求 实验报告原则上要求学生手写，要求书写工整。若因课程特点需打印的，要遵照以下字体、字号、间距等的具体要求。纸张一律采用A4的纸张。 实验报告书写说明 实验报告中一至四项内容为必填项，包括实验目的和要求；实验仪器和设备；实验内容与过程；实验结果与分析。各院部可根据学科特点和实验具体要求增加项目。 填写注意事项 （1）细致观察，及时、准确、如实记录。 （2）准确说明，层次清晰。 （3）尽量采用专用术语来说明事物。 （4）外文、符号、公式要准确，应使用统一规定的名词和符号。 （5）应独立完成实验报告的书写，严禁抄袭、复印，一经发现，以零分论处。 实验报告批改说明 实验报告的批改要及时、认真、仔细，一律用红色笔批改。实验报告的批改成绩采用百分制，具体评分标准由各院部自行制定。 实验报告装订要求

实验批改完毕后，任课老师将每门课程的每个实验项目的实验报告以自然班为单位、按学号升序排列，装订成册，并附上一份该门课程的实验大纲。

实验项目名称：常用连续信号的表示 实验学时： 2学时 同组学生姓名： 无 实验地点： A207 实验日期： 11.12.6 实验成绩： 批改教师： 杨娟 批改时间： 一、实验目的和要求 熟悉MATLAB 软件；利用MATLAB 软件，绘制出常用的连续时间信号。 二、实验仪器和设备 586以上计算机，装有MATLAB7.0软件 三、实验过程 1. 绘制正弦信号)t Asin t (f 0?ω+=（），其中A=1，πω2=，6/π?=； 2. 绘制指数信号at Ae t (f =），其中A=1，0.4a -=； 3. 绘制矩形脉冲信号，脉冲宽度为2； 4. 绘制三角波脉冲信号，脉冲宽度为4；斜度为0.5； 5. 对上题三角波脉冲信号进行尺度变换，分别得出)2t (f ，)2t 2(f -； 6. 绘制抽样函数Sa （t ）,t 取值在-3π到+3π之间； 7. 绘制周期矩形脉冲信号，参数自定； 8. 绘制周期三角脉冲信号，参数自定。 四、实验结果与分析 1．制正弦信号)t Asin t (f 0?ω+=（），其中A=1，πω2=，6/π?= 实验代码： A=1;

随机信号通过线性和非线性系统后地特性分析报告 实验报告材料

实验三 随机信号通过线性和非线性系统后的特性分析 一、实验目的 1、了解随机信号的均值、均方值、方差、自相关函数、互相关函数、概率密度、频谱及功率谱特性。 2、研究随机信号通过线性系统和非线性系统后的均值、均方值、方差、自相关函数、互相关函数、概率密度、频谱及功率谱有何变化，分析随机信号通过线性系统和非线性系统后的特性 二、实验仪器与软件平台 1、 微计算机 2、 Matlab 软件平台 三、实验步骤 1、 根据本实验内容和要求查阅有关资料，设计并撰写相关程序流程。 2、 选择matlab 仿真软件平台。 3、 测试程序是否达到设计要求。 4、 分析实验结果是否与理论概念相符 四、实验内容 1、 随机信号通过线性系统和非线性系统后的特性分析 （1）实验原理 ①随机信号的分析方法 在信号系统中，可以把信号分成两大类：确定信号和随机信号。确定信号具有一定的变化规律，二随机信号无一定的变化规律，需要用统计特性进行分析。在这里引入了一个随机过程的概念。所谓随机过程，就是随机变量的集合，每个随机变量都是随机过程的一个采样序列。随机过程可以分为平稳的和非平稳的，遍历的和非遍历的。如果随机信号的统计特性不随时间的推移而变化。则随机过程是平稳的。如果一个平稳的随机过程的任意一个样本都具有相同的统计特性。则随机过程是遍历的。下面讨论的随机过程都认为是平稳的遍历的随机过程，因此，可以随机取随机过程的一个样本值来描述随机过程中的统计特性。 随机过程的统计特性一般采用主要的几个平均统计特性函数来描述，包括、均方值、方差、自相关系数、互相关系数、概率密度、频谱及功率谱密度等。 a.随机过程的均值 均值E[x(t)]表示集合平均值或数学期望值。基于过程的各态历经行，可用时间间隔T 内的幅值平均值表示，即 ∑-==1 /)()]([N t N t x t x E 均值表达了信号变化的中心趋势，或称之为直流分量。

信号系统实验报告

电子工程系 信号与系统课程实验报告 2011-----2012学年第一学期 专业: 电子信息工程技术班级: 学号 : 姓名: 指导教师: 实常用连续时间信号的实现

一、实验目的 （1）了解连续时间信号的特点； （2）掌握连续时间信号表示的向量法和符号法； （3）熟悉MATLAB Plot函数等的应用。 二、实验原理 1、信号的定义 信号是随时间变化的物理量。信号的本质是时间的函数。 2、信号的描述 1）时域法 时域法是将信号表示成时间的函数f(t)来对信号进行描述的方法。信号的时间特性指的是信号的波形出现的先后，持续时间的长短，随时间变化的快慢和大小，周期的长短等。 2）频域（变换域）法 频域法是通过正交变换，将信号表示成其他变量的函数来对信号进行描述的方法。一般常用的是傅立叶变换。信号的频域特性包括频带的宽窄、频谱的分布等。 信号的频域特性与时域特性之间有着密切的关系。 3、信号的分类 按照特性的不同，信号有着不同的分类方法。 （1）确定性信号：可以用一个确定的时间函数来表示的信号。 随机信号：不可以用一个确定的时间函数来表示，只能用统计特性加以描述的信号。 （2）连续信号：除若干不连续的时间点外，每个时间点在t上都有对应的数值信号。离散信号：只在某些不连续的点上有数值，其他时间点上信号没有定义的信号。 （3）周期信号：存在T，使得等式f(t+T)=f(t)对于任意时间t都成立的信号。非周期信号：不存在使得等式f(t+T)=f(t)对于任意时间t都成立的信号。 绝对的周期信号是不存在的，一般只要在很长时间内慢走周期性就可以了。 （4）能量信号：总能量有限的信号。 功率信号:平均功率有限切非零的信号。 （5）奇信号：满足等式f(t)=--f(--t)的信号。偶信号：满足等式f(t)=f(--t)的信号。 三、涉及的MATLAB函数 1、plot函数 功能：在X轴和Y轴方向都按线性比例绘制二维图形。 调用格式： Plot（x,y）:绘出相x对y的函数线性图。 Plot（x1,y1,x2,y2,…..）:会出多组x对y的线性曲线图。 2、ezplot函数 功能：绘制符号函数在一定范围内的二维图形。简易绘制函数曲线。 调用格式： Ezplot (fun):在[-2π，2π]区间内绘制函数。 Ezplot (fun，［min，max］)：在［min,max］区间内绘函数。 Ezplot (funx,funy):定义同一曲面的函数，默认的区间是[0, 2π]。】 3、sym函数 功能：定义信号为符号的变量。 调用格式：sym(fun):fun为所要定义的表达式。 4、subplot函数

北京理工大学信号与系统实验实验报告

实验1 信号的时域描述与运算 一、实验目的 1. 掌握信号的MATLAB表示及其可视化方法。 2. 掌握信号基本时域运算的MA TLAB实现方法。 3. 利用MA TLAB分析常用信号，加深对信号时域特性的理解。 二、实验原理与方法 1. 连续时间信号的MATLAB表示 连续时间信号指的是在连续时间范围内有定义的信号，即除了若干个不连续点外，在任何时刻信号都有定义。在MATLAB中连续时间信号可以用两种方法来表示，即向量表示法和符号对象表示法。 从严格意义上来说，MATLAB并不能处理连续时间信号，在MATLAB中连续时间信号是用等时间间隔采样后的采样值来近似表示的，当采样间隔足够小时，这些采样值就可以很好地近似表示出连续时间信号，这种表示方法称为向量表示法。表示一个连续时间信号需要使用两个向量，其中一个向量用于表示信号的时间范围，另一个向量表示连续时间信号在该时间范围内的采样值。例如一个正弦信号可以表示如下： >> t=0:0.01:10; >> x=sin(t); 利用plot(t,x)命令可以绘制上述信号的时域波形，如图1所示。 如果连续时间信号可以用表达式来描述，则还可以采用符号表达式來表示信号。例如对于上述正弦信号，可以用符号对象表示如下： >> x=sin(t); >> ezplot(X); 利用ezplot(x)命令可以绘制上述信号的时域波形 Time(seconds) 图1 利用向量表示连续时间信号

t 图 2 利用符号对象表示连续时间信号 sin(t) 2.连续时间信号的时域运算 对连续时间信号的运算包括两信号相加、相乘、微分、积分，以及位移、反转、尺度变换（尺度伸缩）等。 1）相加和相乘 信号相加和相乘指两信号对应时刻的值相加和相乘，对于两个采用向量表示的可以直接使用算术运算的运算符“+”和“*”来计算，此时要求表示两信号的向量时间范围和采样间隔相同。采用符号对象表示的两个信号，可以直接根据符号对象的运算规则运算。 2）微分和积分 对于向量表示法表示的连续时间信号，可以通过数值计算的方法计算信号的微分和积分。这里微分使用差分来近似求取的，由时间向量[N t t t ,,,21?]和采样值向量[N x x x ,,,21?]表示的连续时间信号，其微分可以通过下式求得 1,,2,1,|)('1-?=?-≈ +=N k t x x t x k k t t k 其中t ?表示采样间隔。MA TLAB 中用diff 函数来计算差分 k k x x -+1。 连续时间信号的定积分可以由MATLAB 的qud 函数实现，调用格式为 quad ('function_name',a,b) 其中，function_name 为被积函数名，a 、b 为积分区间。

信号与系统实验报告_1(常用信号的分类与观察)

实验一：信号的时域分析 一、实验目的 1.观察常用信号的波形特点及产生方法 2.学会使用示波器对常用波形参数的测量 二、实验仪器 1.信号与系统试验箱一台（型号ZH5004） 2.40MHz双踪示波器一台 3.DDS信号源一台 三、实验原理 对于一个系统特性的研究，其中重要的一个方面是研究它的输入输出关系，即在一特定的输入信号下，系统对应的输出响应信号。因而对信号的研究是对系统研究的出发点，是对系统特性观察的基本手段与方法。在本实验中，将对常用信号和特性进行分析、研究。 信号可以表示为一个或多个变量的函数，在这里仅对一维信号进行研究，自变量为时间。常用信号有：指数信号、正弦信号、指数衰减正弦信号、复指数信号、Sa(t)信号、钟形信号、脉冲信号等。 1、信号：指数信号可表示为f(t)=Ke at。对于不同的a取值，其波形表现为不同的形式，如下图所示： 图1―1 指数信号 2、信号：其表达式为f（t）=Ksin（ωt＋θ），其信号的参数：振幅K、角频率ω、与初始相位θ。其波形如下图所示：

图1－2 正弦信号 3、指数衰减正弦信号：其表达式为其波形如下图： 图1－3 指数衰减正弦信号 4、Sa(t)信号：其表达式为：。Sa(t)是一个偶函数，t= ±π,±2π,…,±nπ时，函数值为零。该函数在很多应用场合具有独特的运用。其信号如下图所示：

图1－4 Sa（t）信号 5、钟形信号（高斯函数）：其表达式为：其信号如下图所示： 图1－5 钟形信号 6、脉冲信号：其表达式为f(t)=u(t)-u(t-T)，其中u(t)为单位阶跃函数。其信号如下图所示： 7、方波信号：信号为周期为T，前T/2期间信号为正电平信号，后T/2期间信号为负电平信号，其信号如下图所示 U(t)

随机信号实验报告

随机信号分析 实验报告 目录 随机信号分析 (1) 实验报告 (1) 理想白噪声和带限白噪声的产生与测试 (2) 一、摘要 (2) 二、实验的背景与目的 (2) 背景： (2) 实验目的: (2) 三、实验原理 (3) 四、实验的设计与结果 (4) 实验设计: (4) 实验结果: (5) 五、实验结论 (12) 六、参考文献 (13) 七、附件 (13) 1

理想白噪声和带限白噪声的产生与测试一、摘要 本文通过利用MATLAB软件仿真来对理想白噪声和带限白噪声进行研究。理想白噪声通过低通滤波器和带通滤波器分别得到低通带限白噪声和帯通带限白噪声。在仿真的过程中我们利用MATLAB工具箱中自带的一些函数来对理想白噪声和带限白噪声的均值、均方值、方差、功率谱密度、自相关函数、频谱以及概率密度进行研究，对对它们进行比较分析并讨论其物理意义。 关键词：理想白噪声带限白噪声均值均方值方差功率谱密度自相关函数、频谱以及概率密度 二、实验的背景与目的 背景： 在词典中噪声有两种定义：定义1：干扰人们休息、学习和工作的声音，引起人的心理和生理变化。定义2：不同频率、不同强度无规则地组合在一起的声音。如电噪声、机械噪声，可引伸为任何不希望有的干扰。第一种定义是人们在日常生活中可以感知的，从感性上很容易理解。而第二种定义则相对抽象一些，大部分应用于机械工程当中。在这一学期的好几门课程中我们都从不同的方面接触到噪声，如何的利用噪声，把噪声的危害减到最小是一个很热门的话题。为了加深对噪声的认识与了解，为后面的学习与工作做准备，我们对噪声进行了一些研究与测试。 实验目的: 了解理想白噪声和带限白噪声的基本概念并能够区分它们，掌握用MATLAB 或c/c++软件仿真和分析理想白噪声和带限白噪声的方法，掌握理想白噪声和带限白噪声的性质。

北京理工大学信号与系统实验报告2 LTI系统的时域分析

实验2 LTI 系统的时域分析 （基础型实验） 一． 实验目的 1. 掌握利用MATLAB 对系统进行时域分析的方法。 2. 掌握连续时间系统零状态响应、冲击响应和阶跃响应的求解方法。 3. 掌握求解离散时间系统响应、单位抽样响应的方法。 4. 加深对卷积积分和卷积和的理解。掌握利用计算机进行卷积积分和卷积和计算的方法。 二． 实验原理与方法 1. 连续时间系统时域分析的MATLAB 实现 1） 连续时间系统的MA TLAB 表示 LTI 连续系统通常可以由系统微分方程描述，设描述系统的微分方程为： (N)(N 1)(M)(M 1)1010(t)(t)...(t)b (t)b (t)...b (t)N N M M a y a y a y x x x ----++=++ 则在MATLAB 中可以建立系统模型如下： 1010[b ,b ,...,b ];a [a ,a ,...,a ];sys tf(b,a); M M N N b --=== 其中，tf 是用于创建系统模型的函数，向量a 和b 的元素是以微分方程求导的降幂次序来排列的，如果有缺项，应用0补齐，例如由微分方程 2''(t)y'(t)3y(t)x(t)y ++= 描述的系统可以表示为： >> b=[1]; >> a=[2 1 3]; >> sys=tf(b,a); 而微分方程由 ''(t)y'(t)y(t)x''(t)x(t)y ++=- 描述的系统则要表示成 >> b=[1 0 -1]; >> a=[1 1 1]; >> sys=tf(b,a); 2） 连续时间系统的零状态响应 零状态响应指系统的初始状态为零，仅由初始信号所引起的响应。MATLAB 提供了一个用于求解零状态响应的函数lism ，其调用格式如下： lism （sys,x,t ）绘出输入信号及响应的波形，x 和t 表示输入信号数值向量及其时间向量。 y= lism （sys,x,t ）这种调用格式不绘出波形，而是返回响应的数值向量。 3） 连续时间系统的冲激响应与阶跃响应

随机信号分析实验报告(基于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)均匀分布随机数，通常采用的方法为线性同余法，公式如下： y0=1，y n=ky n(mod N) ? x n=y n N 序列{x n}为产生的(0,1)均匀分布随机数。 定理1.1若随机变量X 具有连续分布函数F x(x)，而R 为(0,1)均匀分布随机变量，则有 X=F x?1(R) 2.MATLAB中产生随机序列的函数 (1)(0,1)均匀分布的随机序列函数：rand 用法：x = rand(m,n) 功能：产生m×n 的均匀分布随机数矩阵。 (2)正态分布的随机序列 函数：randn 用法：x = randn(m,n) 功能：产生m×n 的标准正态分布随机数矩阵。 如果要产生服从N(μ,σ2)分布的随机序列，则可以由标准正态随机序列产生。 (3)其他分布的随机序列 分布函数分布函数 二项分布binornd 指数分布exprnd 泊松分布poissrnd 正态分布normrnd 离散均匀分布unidrnd 瑞利分布raylrnd 均匀分布unifrnd X2分布chi2rnd 3.随机序列的数字特征估计 对于遍历过程，可以通过随机序列的一条样本函数来获得该过程的统计特征。这里我们假定随机序列X(n)为遍历过程，样本函数为x(n)，其中n=0,1,2,……N-1。那么，

随机信号分析实验报告二 2

《随机信号分析》实验报告二 班级： 学号： 姓名：

实验二高斯噪声的产生和性能测试 １．实验目的 （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、源程序及功能注释 （2）源程序： clear all; clc; t=0:320; %t=0:320 x=sin(2*pi*t/25); %x=sin(2*p1*t/25) x1=wgn(1,321,0); %产生一个一行32列的高斯白噪声矩阵，输出的噪声强度为0dbw

信号与系统实验报告

实验三 常见信号的MATLAB 表示及运算 一、实验目的 1．熟悉常见信号的意义、特性及波形 2．学会使用MATLAB 表示信号的方法并绘制信号波形 3. 掌握使用MATLAB 进行信号基本运算的指令 4. 熟悉用MATLAB 实现卷积积分的方法 二、实验原理 根据MATLAB 的数值计算功能和符号运算功能，在MA TLAB 中，信号有两种表示方法，一种是用向量来表示，另一种则是用符号运算的方法。在采用适当的MA TLAB 语句表示出信号后，就可以利用MA TLAB 中的绘图命令绘制出直观的信号波形了。 1.连续时间信号 从严格意义上讲，MATLAB 并不能处理连续信号。在MATLAB 中，是用连续信号在等时间间隔点上的样值来近似表示的，当取样时间间隔足够小时，这些离散的样值就能较好地近似出连续信号。在MATLAB 中连续信号可用向量或符号运算功能来表示。 ⑴ 向量表示法 对于连续时间信号()f t ，可以用两个行向量f 和t 来表示，其中向量t 是用形如12::t t p t =的命令定义的时间范围向量，其中，1t 为信号起始时间，2t 为终止时间，p 为时间间隔。向量f 为连续信号()f t 在向量t 所定义的时间点上的样值。 ⑵ 符号运算表示法 如果一个信号或函数可以用符号表达式来表示，那么我们就可以用前面介绍的符号函数专用绘图命令ezplot()等函数来绘出信号的波形。 ⑶ 常见信号的MATLAB 表示 单位阶跃信号 单位阶跃信号的定义为：10()0 t u t t >?=? 0); %定义函数体，即函数所执行指令

信号与系统实验报告

信号与系统实验报告

信号与系统实验报告 姓名： 学号： 软件部分： 表示信号与系统的MATLAB 函数、工具箱 一、实验项目名称：表示信号、系统的MATLAB 函数、工具箱 二、实验目的与任务： 目的：1、加深对常用离散信号的理解； 2、熟悉表示信号的基本MATLAB 函数。 任务：基本MATLAB 函数产生离散信号；基本信号之间的简单运算；判 断信号周期。 三、实验原理： 利用MATLAB 强大的数值处理工具来实现信号的分析和处理，首先就是要学会应用MATLAB 函数来构成信号。 四、实验内容及步骤： 常见的基本信号可以简要归纳如下： 实验内容(一)、 编制程序产生上述5种信号（长度可输入确定），并绘出其图形。 其中5种信号分别为单位抽样序列、单位阶跃序列、正弦序列、指数序列和复正弦序列。 实验内容（二）、 在[0,31]出下列图像 1223[]sin( )cos() 4 4 []cos ( ) 4[]sin()cos() 48 n n x n n x n n n x n πππππ=== 五、项目需用仪器设备名称：计算机、MATLAB 软件。

六、所需主要元器件及耗材：无 七、实验程序及数据 函 数 程序图片 单位冲击函数x=zeros(1,10); x(1)=1; stem(x) 单位阶跃函数x=ones(1,30); plot(x)

正弦序列n=0:30-1; x=sin(2*pi*n/10); stem(x) x=cos(1/4*pi*n).*cos(1/4*pi*n) ; stem(x) 复正弦序列n=0:29; x=exp(j*5*n); stem(x) 指数序列n=0:10; x=2.^n; stem(x)

随机信号分析实验报告

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 +=+ M y x n n 1 1++= 为了保证产生的伪随机数能在[0,1]内均匀分布，需要M 为正整数，此外常数c 和初值y0亦为正整数。加同余法虽然简单，但产生的伪随机数效果不好。另一种同余法为乘同余法，它需要两次乘法才能产生一个[0,1]上均匀分布的随机数 )(mod 1M ay y n n =+ M y x n n 1 1++= 式中，a 为正整数。用加法和乘法完成递推运算的称为混合同余法，即 )(mod 1M c ay y n n +=+ M y x n n 1 1++= 用混合同余法产生的伪随机数具有较好的特性，一些程序库中都有成熟的程序供选择。 常用的计算语言如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.随机变量的仿真 根据随机变量函数变换的原理，如果能将两个分布之间的函数关系用显式表达，那么就可以利用一种分布的随机变量通过变换得到另一种分布的随机变量。 若X 是分布函数为F(x)的随机变量，且分布函数F(x)为严格单调升函数，令Y=F(X)，则Y 必为在[0,1]上均匀分布的随机变量。反之，若Y 是在[0,1]上 均匀分布的随机变量，那么)(1 Y F X X -= 即是分布函数为FX(x)的随机变量。式中F X -?1 ()为F X ()?的反函数。这样，欲求某个分布的随机变量，先产生在[0,1]区间上的均匀分布随机数，再经上式变 换，便可求得所需分布的随机数。 3.高斯分布随机数的仿真 广泛应用的有两种产生高斯随机数的方法，一种是变换法，一种是近似法。 如果X1,X2是两个互相独立的均匀分布随机数，那么下式给出的Y1,Y2

信号与系统实验报告3 (2)

信号与系统实验 实验三：信号的卷积 （第三次实验）

【实验目的】 1. 理解卷积的物理意义； 2. 掌握运用计算机进行卷积运算的原理和方法； 3. 熟悉卷积运算函数conv的应用； 【实验内容】 给定如下因果线性时不变系统： y[n]+0.71y[n-1]-0.46y[n-2]-0.62y[n-3=0.9x[n]-0.45x[n-1]+0.35x[n-2]+0.002x[n-3] (1)不用impz函数，使用filter命令，求出以上系统的单位冲激响应h[n]的前20个样本; 代码如下： clear all; N=[0:19]; num=[0.9 -0.45 0.35 0.002]; den=[1 0.71 -0.46 -0.62]; h=filter(num,den,N); stem(N,h); xlabel('ê±??Dòo?'); ylabel('??·ù'); title('μ￥??3??¤?ìó|'); grid; 图像如下：

(2)得到h[n]后，给定x[n],计算卷积输出y[n]；并用滤波器h[n]对输入x[n]滤波，求得y1[n]; 代码如下： clear all; N=[0:19]; num=[0.9 -0.45 0.35 0.002]; den=[1 0.71 -0.46 -0.62]; h=filter(num,den,N); x=[1 -2 3 -4 3 2 1]; y=conv(h,x); n=0:25; subplot(2,1,1); stem(n,y); xlabel('时间序号n');ylabel('振幅'); title('用卷积得到的输出');grid; x1=[x zeros(1,19)]; y1=filter(h,1,x1); subplot(2,1,2); stem(n,y1);

中科大信号与系统 实验报告2

信号与系统实验报告 一、实验目的 1. 熟悉连续时间系统的单位冲激响应、阶跃响应的意义及求解方法 2. 熟悉连续（离散）时间系统在任意信号激励下响应的求解方法 3. 熟悉应用MATLAB实现求解系统响应的方法 二、实验原理 1.连续时间系统求解各种响应 impulse( ) 函数 函数impulse( )将绘制出由向量a和b所表示的连续系统在指定时间范围内的单位冲激响应h(t)的时域波形图，并能求出指定时间范围内冲激响应的数值解。 以默认方式绘出由向量a和b所定义的连续系统的冲激响应的时域波形。 绘出由向量a和b所定义的连续系统在0 ~ t0时间范围内冲激响应的时域波形。 绘出由向量a和b所定义的连续系统在t1 ~ t2时间范围内，并且以时间间隔p均匀取样的冲激响应的时域波形。 只求出由向量a和b所定义的连续系统在t1 ~ t2时间范围内，并且以时间间隔p均匀取样的冲激响应的数值解，但不绘出其相应波形。 step( ) 函数 函数step( )将绘制出由向量a和b所表示的连续系统的阶跃响应，在指定的时间范围内的波形图，并且求出数值解。和impulse( )函数一样，step( )也有如下四种调用格式： step( b,a) step(b,a,t0) step(b,a,t1:p:t2) y=step(b,a,t1:p:t2) 上述调用格式的功能和impulse( )函数完全相同，所不同只是所绘制（求解）的是系统的阶跃响应g(t)，而不是冲激响应h(t)。 lsim( )函数 根据系统有无初始状态，lsim( )函数有如下两种调用格式： ①系统无初态时，调用lsim( )函数可求出系统的零状态响应，其格式如下：

随机信号处理模实验报告

随机信号分析与处理实验报告院系：信息工程学院 专业：电子信息科学与技术 姓名: 方静 学号：030941209 指导老师：廖红华

实验一 熟悉MATLAB 的随机信号处理相关命令 一、实验目的 1、利用Matlab 对随机熟悉各种随机信号函数的用法 2、掌握随机信号的简单分析方法 二、实验原理 1、语音的录入与打开 在MATLAB 中，wavread 函数用于读取语音信号，采样值放在向量y 中，s f 表示采样频率(Hz)，bits 表示 采样位数。[N1 N2]表示读取从N1点到N2点的值。 2、语音信号的频域分析 FFT 即为快速傅氏变换，是离散傅氏变换的快速算法，它是根据离散傅氏变换的奇、偶、虚、实等特性，对离散傅立叶变换的算法进行改进获得的。在Matlab 信号处理工具箱中，语音信号的频域分析就是对信号进行傅里叶变换后的分析。 4、方差 定义22)]}()({[t t m t X E X X -=）（δ 为随机过程的方差。方差通常也记为DX （t ） ，随机过程的方差也是时间 t 的函数, 由方差的 定义可以看出，方差是非负函数。 5、自相关与互相关 自相关和互相关分别表示的是两个时间序列之间和同一个时间序列在任意两个不同时刻的取值之间的相关程度，即互相关函数是描述随机信号x(t),y(t)在任意两个不同时刻t1，t2的取值之间的相关程度，自相关函数是描述随机信号x(t)在任意两个不同时刻t1，t2的取值之间的相关程度。 互相关函数给出了在频域内两个信号是否相关的一个判断指标，把两测点之间信号的互谱与各自的自谱联系了起来。它能用来确定输出信号有多大程度来自输入信号，对修正测量中接入噪声源而产生的误差非常有效. 事实上，在图象处理中，自相关和互相关函数的定义如下：设原函数是f(t)，则自相关函数定义为R(u)=f(t)*f(-t)，其中*表示卷积；设两个函数分别是f(t)和g(t)，则互相关函数定义为R(u)=f(t)*g(-t)，它反映的是两个函数在不同的相对位置上互相匹配的程度。 6. 短时过零率与短时能量 语音一般分为无声段，清音段和浊音段。由于语音信号是一个非平稳过程，不能用处理平稳信号的信号处理技术对其进行分析处理。但由于语音信号本身的特点，在10-30ms 的短时间范围内，其特性可以看作是一个准稳态过程，具有短时性，因此采用短时能量和过零率来对语音进行端点检测是可行的。 信号的短时能量定义为：设语音波形时域信号为x(t)，加窗分帧处理后得到第n 帧语音信号为xn(m),则定义的短时能量函数如下： ) ()()(x m n x m w m n +=，10-≤≤ N m ，,0)(),1(~0,1)(=-==n w N m m w m 为其他值，其中n=0,1T,2T……并且N 为帧长，T 为帧移长度。 短时过零率表示一帧语音中语音信号的波形穿过横轴的零电平的次数，他可以用来区分清音和浊音，因为语音信号中高音段有高的过零率，低音段有低的过零率，短时能量大的地方过零率小，短时能量小的地方过零率大。 过零率可以反映信号的频谱特性。当离散时间信号相邻两个样点的正负号相异时，我们称之为“过零”，即此时信号的时间波形穿过了零电平的横轴。统计单位时间内样点值改变符号的次数具可以得到平均过零

信号与系统实验报告(常用信号的分类与观察)

实验一:信号得时域分析 一、实验目得 1.观察常用信号得波形特点及产生方法 2.学会使用示波器对常用波形参数得测量 二、实验仪器 1.信号与系统试验箱一台（型号ＺH５0０4) 2.40MHｚ双踪示波器一台 3.ＤDS信号源一台 三、实验原理 对于一个系统特性得研究,其中重要得一个方面就是研究它得输入输出关系，即在一特定得输入信号下,系统对应得输出响应信号.因而对信号得研究就是对系统研究得出发点,就是对系统特性观察得基本手段与方法.在本实验中,将对常用信号与特性进行分析、研究。 信号可以表示为一个或多个变量得函数，在这里仅对一维信号进行研究,自变量为时间。常用信号有：指数信号、正弦信号、指数衰减正弦信号、复指数信号、Sa(ｔ)信号、钟形信号、脉冲信号等。 1、信号:指数信号可表示为ｆ(t）＝Ke at。对于不同得ａ取值，其波形表现为不同得形式,如下图所示: 图１―1 指数信号 2、信号:其表达式为f(t）=Ｋsin(ωt+θ),其信号得参数:振幅K、角频率ω、与初始相位θ。其波形如下图所示:

图1-2 正弦信号 3、指数衰减正弦信号:其表达式为其波形如下图: 图１-３指数衰减正弦信号 ４、Sａ(t）信号:其表达式为:。Sa（t）就是一个偶函数，t= ±π,±２π,…,±ｎπ时,函数值为零。该函数在很多应用场合具有独特得运用。其信号如下图所示: 图1－4 Ｓa(t)信号 5、钟形信号(高斯函数):其表达式为：其信号如下图所示：

图1-5 钟形信号 6、脉冲信号:其表达式为ｆ(t)=u(t)-ｕ（t—T),其中u（t)为单位阶跃函数。其信号如下图所示: f(ｔ） ? ……??…… 0 t 图1-6脉冲信号 7、方波信号:信号为周期为T,前T/2期间信号为正电平信号,后T/２期间信号为负电平信号，其信号如下图所示 U(t) ………… ?0?t 图1－７方波信号 四、实验内容及主要步骤 下列实验中信号产生器得工作模式为11 1、指数信号观察 通过信号选择键1，设置Ａ组输出为指数信号(此时信号输出指示灯为0００000）。用示波器测量“信号Ａ组”得输出信号。 输出波形为：

信号系统实验报告123(1)

实验一连续时间信号的时域基本运算 一、实验目的 （1）掌握连续时间信号时域运算的基本方法； （2）掌握相关函数的调用格式及作用； （3）掌握连续信号的基本运算。 二、实验原理 信号的基本运算包括信号的相加（减）和相乘(除).信号的时域变换包括信号的平移、翻转、倒相以及尺度变换。 （1）加减： f(t)=f1(t)±f2(t)（2）乘: f(t)=f1(t)×f2(t) （3）延时或平移：f(t)→(t-t0) t0>0时右移；t0<0时左移 （4）翻转： f(t)→f(-t) （5）尺度变换：f(t)→ f(at) |a|>1时尺度缩小；|a|<1时尺度放大；a<0时，尺度翻转。 （6）标量乘法：f(t)→af(t) （7）倒相： f(t)→-f(t) （8）微分： f(t)→df(t)/dt （9）积分： f(t)→∫t -∞f(t)d(t) 三、涉及的MATLAB函数及其运算 1、stepfun函数 功能：产生一个阶跃信号。 调用格式： Stepfun(t,t 0)其中，t是时间区间，在该区间内阶跃信号一定会产生；t 是信号 发生从0到1跳跃的时刻。 2、diff函数 调用格式： diff (f) : 求函数f对预设独立变数的一次微分值。 diff (f, ’t’) : 求函数f对独立变数t的一次微分值。 3、int函数 调用格式: Int（f）: 函数f对预设独立变数的积分值。 Int(f,’t’): 函数f对独立变数t的积分值。 4、heaviside函数 Heaviside（t）:产生没有移位的阶跃信号。 Heaviside（t-k）:产生向右平移K单位的阶跃信号。 四、实验内容与方法1、验证性实验 （1）移位 实现连续信号的移位，即f(t-t0),或者f(t+t0);常数t0>0。 MATLAB程序 clear all t=0:0.0001:2 y=sin(2*pi*(t)); y1=sin(2*pi*(t-0.2)); plot(t,y,'-',t,y1,'--'); ylabel('f(t)');xlabel('t');title('信号的移位') 运行结果：

随机信号分析与处理实验报告

随机信号分析与处理实 验 题目：对音频信号的随机处理 班级：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.