随机信号分析仿真
.
'. 随机信号分析原理大作业报告
专业：水声工程
姓名： xxx
学号：xxxxxxxxxx
题目要求：
给定一个白噪声信号，它的均值和方差自定。
1.设计一个线性滤波器，使该滤波器的输出为一个窄带信号。并给出
该窄带信号在不同的3个典型中心频率和带宽时的波形。
2.对该滤波器输出的上述窄带信号，用莱斯表示法对其进行建模，画
出)(t a和)(t b的波形。
3.计算上述3种窄带信号对应的瞬时频率和瞬时相位，并进行包络检
测。
1 窄带信号的生成 1.1 高斯白噪声的产生
若)(t N 为一个具有零均值的平稳随机过程，其功率谱密度均匀分布在（-）整个频率区间，即
02
1
)(N w S n =
(1) 其中，0N 为一个正实常数，则称)(t N 为白噪声。 白噪声的自相关函数为
)(2
1
)(0τδτN R N =
(2) 白噪声在任意两个相邻时刻（不管这两个时刻多么的近）的取值都是不相关的，这意味着白噪声过程随时间的起伏很快，过程的功率谱极宽。
这种形式定义的白噪声只是一种理想化的模型，实际上这种白噪声是不存在的，因为按照定义，白噪声的均方值为无限大，而物理上存在的随机过程，其均方值总是有限的。在实际工作中，当所研究的随机过程通过某一系统时，只要过程的功率谱密度在一个比系统宽大的多的频率范围内近似均匀分布，就可以把它作为白噪声来处理。
现产生一均值为0，方差为3高斯白噪声，如图1所示：
滤波器2：中心频率f0=700Hz，带通：650～750 Hz。滤波器2的幅度相位图以及滤波器2输出信号波形如下：
图6 滤波器2输出信号的时域波形
附件一滤波器1输出信号仿真程序clear all
close all
clc
%产生高斯白噪声
N=25000; %序列长度
my_var = 2;
noise = sqrt(my_var)*randn(1,N);%均值为0，方差为2 figure(1)
plot(noise)
title('均值为0方差为2的高斯白噪声')
grid on
fs = 25000;%采样频率
f0 = 1000;%中心频率
%滤波器
f_pass = [900 1100];
omega_pass = 2*f_pass/fs;
b = fir1(192,omega_pass);
figure(2)
freqz(b,1,1024)%滤波器幅度和相位图像
grid on
%噪声通过窄带滤波器
filter_outpu = filter(b,1,noise);
figure(3)
plot(filter_outpu)
title('窄带信号在时域的波形')
grid on
%做fft变换
Nfft = fs;
fft_x = fft(filter_outpu,Nfft);
ff = 0:fs/Nfft:fs-fs/Nfft;
figure(4)
plot(ff,20*log10(abs(fft_x)))%窄带信号的频谱
title('窄带信号的频谱')
xlabel('频率 Hz')
ylabel('幅度 dB')
grid on
%窄带信号在时域的波形
X_t = filter_outpu;
t = 0:1/fs:1-1/fs;
figure(5)
plot(t,X_t)
title('窄带信号在时域的波形')
xlabel('t / s')
grid on
%莱斯表示法
h_X = hilbert(X_t,Nfft) ;%希尔伯特变换
omega0 = 2*pi*f0;
A_t = X_t.*cos(omega0*t)+h_X.*sin(omega0*t);
B_t = -1*X_t.*sin(omega0*t)+h_X.*cos(omega0*t); figure(6)
subplot(2,1,1);
plot(t,A_t)
grid on
hold on
subplot(2,1,2);
plot(t,B_t)
grid on
%瞬时频率瞬时相位
theta_t = atan(h_X./X_t); xh1=unwrap(angle(h_X)); omega_t=fs*diff(xh1)/(2*pi); figure(7)
plot(omega_t);
title('瞬时频率')
omega_t = diff(theta_t); figure(8)
plot(t,theta_t)
title('瞬时相位')
grid on
%包络检测
am = abs(h_X);
figure(9)
plot(t,X_t,t,am,'r') %包络title('窄带信号的包络')
grid on
随机信号分析实验报告
一、实验名称 微弱信号的检测提取及分析方法 二、实验目的 1．了解随机信号分析理论如何在实践中应用 2．了解随机信号自身的特性，包括均值、方差、相关函数、频谱及功率谱密度等 3．掌握随机信号的检测及分析方法 三、实验原理 1．随机信号的分析方法 在信号与系统中，我们把信号分为确知信号和随机信号。其中随机信号无确定的变化规律，需要用统计特新进行分析。这里我们引入随机过程的概念，所谓随机过程就是随机变量的集合，每个随机变量都是随机过程的一个取样序列。 随机过程的统计特性一般采用随机过程的分布函数和概率密度来描述，他们能够对随机过程作完整的描述。但由于在实践中难以求得，在工程技术中，一般采用描述随机过程的主要平均统计特性的几个函数，包括均值、方差、相关函数、频谱及功率谱密度等来描述它们。本实验中算法都是一种估算法，条件是N要足够大。 2．微弱随机信号的检测及提取方法 因为噪声总会影响信号检测的结果，所以信号检测是信号处理的重要内容之一，低信噪比下的信号检测是目前检测领域的热点，而强噪声背景下的微弱信号提取又是信号检测的难点。 噪声主要来自于检测系统本身的电子电路和系统外空间高频电磁场干扰等，通常从以下两种不同途径来解决 ①降低系统的噪声，使被测信号功率大于噪声功率。 ②采用相关接受技术，可以保证在信号功率小于噪声功率的情况下，人能检测出信号。 对微弱信号的检测与提取有很多方法，常用的方法有：自相关检测法、多重自相法、双谱估计理论及算法、时域方法、小波算法等。 对微弱信号检测与提取有很多方法，本实验采用多重自相关法。 多重自相关法是在传统自相关检测法的基础上，对信号的自相关函数再多次做自相关。即令： 式中，是和的叠加；是和的叠加。对比两式，尽管两者信号的幅度和相位不同，但频率却没有变化。信号经过相关运算后增加了信噪比，但其改变程度是有限的，因而限制了检测微弱信号的能力。多重相关法将 当作x（t），重复自相关函数检测方法步骤，自相关的次数越多，信噪比提高的越多，因此可检测出强噪声中的微弱信号。
信号参考电源层的仿真分析
大多数layout工程师以及SI/硬件工程师都知道, 信号除了不能跨分割层布线之外，一般还不容许参考电源层布线的(当然，这里指的高速高频信号)，为什么不能参考电源层？究竟会带来多大影响？如果叠层空间限制的情况下可以容许哪 些信号参考电源？针对这些问题，本篇将结合ANSYS/Ansoft仿真软件进行理论及仿真方法介绍。 1 参考电源层的回流路径 首先，从信号回流路径的角度开始基本理论的回顾。一个简单四层PCB信 号通过过孔换层参考电源，其信号的回流路径如图1 示意: 图1 信号回流路径 由上图可见，当高速信号在信号线上传播时，在信号电流向前传播的过程当中，由于与参考平面之间存在容性耦合，所以当发生dV/dt时，就会有电流经耦合电容流向参考平面的现象，传输线正下方位置都会有瞬态电流流回到源端电路。如果信号的参考为电源平面，那么信号回流将首先流向电源层，然后再通过电源与地网络之间的Cpg流向地网络，最后再经地层流向源端电路，最终形成一个 完整的电流回路。我们都知道，控制好高速信号的回路阻抗非常关键，因为它直接影响到信号传输特性。 当信号参考电源层布线时，回流路径当中对信号影响最大的就是Cpg电源与地网络之间的容性通道。它可以是电源地网络上分布复杂的退耦电容，也可能包含电源地层平面之间的平板电容，构成非常复杂，在各个频点所表现的阻抗特性都不一样，难以量化与控制。所以不建议高速信号参考电源。 那么究竟有多大影响，下面通过仿真软件来帮忙我们看看具体信号传输差异的情况。
2，参考电源层的仿真分析 2.1 基础研究模型的建立 有了以上理论了解之后，接下来通过仿真技术协助研究，到底参考电源层会跟信号传输带来怎样的影响？ 为了说明问题，把模型简单化，这里利用板级仿真工具SIwave的自行建模功能（也可通过版图工具画一个类似PCB走线再导入）建立一个简单的10X10四层PCB, 叠层分布为SIG/GND/PWR/SIG，第二层全部为地，第三层电源平面为一小块不规则平面，如下图，并布置两根传输线，一根为表层走线，此案例中，它属于完全参考地层平面的微带线，一根为表层走线经过孔到底层走线的微带线，属于部分参考地层又部分参考电源层的走线。即建立了我们需要研究的参考电源的信号模型。如图2所示： 图2 简单的四层PCB模型 2.2 回流仿真分析 通过SIwave2014以上版本的AC CURRENTS 功能可以进行信号回流路径的仿真分析，只需要在两条传输线两端分别添加相应频率的信号源和负载，即可仿真得到信号源传输时，各个平面层上的电流分别情况。如图3所示，显示为地层的电流分布，跟前面理论分析结论非常一致。完全参考地层的传输线，回流路径主要集中在走线正下方，而参考电源层的信号回流会经电源地耦合到地层上，所以在电源与地层重叠的地方分布，不同频点的回流分布也不尽相同，这势必会影响信号传送质量，同时也可能对外界电路造成干扰。
五款信号完整性仿真工具介绍
现在的高速电路设计已经达到GHz的水平，高速PCB设计要求从三维设计理论出发对过孔、封装和布线进行综合设计来解决信号完整性问题。高速PCB设计要求中国工程师必须具备电磁场的理论基础，必须懂得利用麦克斯韦尔方程来分析PCB设计过程中遇到的电磁场问题。目前，Ansoft公司的仿真工具能够从三维场求解的角度出发，对PCB设计的信号完整性问题进行动态仿真。 (一)Ansoft公司的仿真工具 现在的高速电路设计已经达到GHz的水平，高速PCB设计要求从三维设计理论出发对过孔、封装和布线进行综合设计来解决信号完整性问题。高速PCB设计要求中国工程师必须具备电磁场的理论基础，必须懂得利用麦克斯韦尔方程来分析PCB设计过程中遇到的电磁场问题。目前，Ansoft公司的仿真工具能够从三维场求解的角度出发，对PCB设计的信号完整性问题进行动态仿真。 Ansoft的信号完整性工具采用一个仿真可解决全部设计问题： SIwave是一种创新的工具，它尤其适于解决现在高速PCB和复杂IC封装中普遍存在的电源输送和信号完整性问题。 该工具采用基于混合、全波及有限元技术的新颖方法，它允许工程师们特性化同步开关噪声、电源散射和地散射、谐振、反射以及引线条和电源/地平面之间的耦合。该工具采用一个仿真方案解决整个设计问题，缩短了设计时间。 它可分析复杂的线路设计，该设计由多重、任意形状的电源和接地层，以及任何数量的过孔和信号引线条构成。仿真结果采用先进的3D图形方式显示，它还可产生等效电路模型，使商业用户能够长期采用全波技术，而不必一定使用专有仿真器。 (二)SPECCTRAQuest Cadence的工具采用Sun的电源层分析模块： Cadence Design Systems的SpecctraQuest PCB信号完整性套件中的电源完整性模块据称能让工程师在高速PCB设计中更好地控制电源层分析和共模EMI。 该产品是由一份与Sun Microsystems公司签署的开发协议而来的，Sun最初研制该项技术是为了解决母板上的电源问题。 有了这种新模块，用户就可根据系统要求来算出电源层的目标阻抗；然后基于板上的器件考虑去耦合要求，Shah表示，向导程序能帮助用户确定其设计所要求的去耦合电容的数目和类型；选择一组去耦合电容并放置在板上之后，用户就可运行一个仿真程序，通过分析结果来发现问题所在。 SPECCTRAQuest是CADENCE公司提供的高速系统板级设计工具，通过它可以控制与PCB layout相应的限制条件。在SPECCTRAQuest菜单下集成了一下工具： （1）SigXplorer可以进行走线拓扑结构的编辑。可在工具中定义和控制延时、特性阻抗、驱动和负载的类型和数量、拓扑结构以及终端负载的类型等等。可在PCB详细设计前使用此工具，对互连线的不同情况进行仿真，把仿真结果存为拓扑结构模板，在后期详细设计中应用这些模板进行设计。 （2）DF/Signoise工具是信号仿真分析工具，可提供复杂的信号延时和信号畸变分析、IBIS 模型库的设置开发功能。SigNoise是SPECCTRAQUEST SI Expert和SQ Signal Explorer Expert进行分析仿真的仿真引擎，利用SigNoise可以进行反射、串扰、SSN、EMI、源同步及系统级的仿真。 （3）DF/EMC工具——EMC分析控制工具。 （4）DF/Thermax——热分析控制工具。 SPECCTRAQuest中的理想高速PCB设计流程： 由上所示，通过模型的验证、预布局布线的space分析、通过floorplan制定拓朴规则、由规
信号与系统仿真实验报告
信号与系统仿真实验报告1．实验目的 了解MATLAB的基本使用方法和编程技术，以及Simulink平台的建模与动态仿真方法，进一步加深对课程内容的理解。 2．实验项目 信号的分解与合成，观察Gibbs现象。 信号与系统的时域分析，即卷积分、卷积和的运算与仿真。 信号的频谱分析，观察信号的频谱波形。 系统函数的形式转换。 用Simulink平台对系统进行建模和动态仿真。 3．实验内容及结果 3.1以周期为T，脉冲宽度为2T1的周期性矩形脉冲为例研究Gibbs现象。 已知周期方波信号的相关参数为：x(t)=∑ak*exp(jkω),ω=2*π/T,a0=2*T1/T,ak=sin(kωT1)/kπ。画出x(t)的波形图（分别取m=1,3,7,19,79,T=4T1）,观察Gibbs现象。 m=1; T1=4; T=4*T1;k=-m:m; w0=2*pi/T; a0=2*T1/T; ak=sin(k*w0*T1)./(k*pi); ak(m+1)=a0; t=0:0.1:40; x=ak*exp(j*k'*w0*t); plot(t,real(x)); 3.2求卷积并画图 (1)已知：x1(t)=u(t-1)-u(t-2), x2(t)=u(t-2)-u(t-3)求：y(t)=x1(t)*x2(t)并画出其波形。 t1=1:0.01:2; f1=ones(size(t1)); f1(1)=0; f1(101)=0; t2=2:0.01:3; f2=ones(size(t2)); f2(1)=0; f2(101)=0; c=conv(f1,f2)/100;
t3=3:0.01:5; subplot(311); plot(t1,f1);axis([0 6 0 2]); subplot(312); plot(t2,f2);axis([0 6 0 2]); subplot(313); plot(t3,c);axis([0 6 0 2]); (2）已知某离散系统的输入和冲击响应分别为：x[n]=[1,4,3,5,1,2,3,5], h[n]=[4,2,4,0,4,2].求系 统的零状态响应，并绘制系统的响应图。 x=[1 4 3 5 1 2 3 5]; nx=-4:3; h=[4 2 4 0 4 2]; nh=-3:2; y=conv(x,h); ny1=nx(1)+nh(1); ny2=nx(length(nx))+nh(length(nh)); ny=[ny1:ny2]; subplot(311); stem(nx,x); axis([-5 4 0 6]); ylabel('输入') subplot(312); stem(nh,h); axis([-4 3 0 5]); ylabel('冲击效应') subplot(313); stem(ny,y); axis([-9 7 0 70]); ylabel('输出'); xlabel('n'); 3.3 求频谱并画图 (1) 门函数脉冲信号x1(t)=u(t+0.5)-u(t-0.5) N=128;T=1; t=linspace(-T,T,N); x=(t>=-0.5)-(t>=0.5); dt=t(2)-t(1); f=1/dt; X=fft(x); F=X(1:N/2+1); f=f*(0:N/2)/N; plot(f,F)
北理工随机信号分析实验报告
本科实验报告实验名称：随机信号分析实验
实验一 随机序列的产生及数字特征估计 一、实验目的 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)均匀分布随机数按上式进行变
信号分析与处理仿真实验
实验报告 实验名称MATLAB仿真实验 课程名称信号分析与处理 院系部: 专业班级：学生姓名：学号:同组人:实验台号:指导教师：成绩：实验日期:2015-11-29
实验一信号的产生与运算 1.单位阶跃信号 （1）源程序 t=-0.5:0.01:1.5; u=stepfun(t,0); u1=stepfun(t,0.5); figure(1) plot(t,u);axis([-0.5 1.5 -0.2 1.2]);title('单位阶跃信号波形'); figure(2) plot(t,u1);axis([-0.5 1.5 -0.2 1.2]);title('延迟单位阶跃信号波形'); （2）实验结果
2.单位冲激信号 （1）源程序 clear;clc; t=-1:0.001:1; for i=1:3 dt=1/(i^4); x=(1/dt)*((t>=-(1/2*dt))-(t>=(1/2*dt))); subplot(1,3,i); stairs(t,x); end （2）实验结果
3.抽样信号 （1）源程序 clear;clc; t=-20:0.01:20; x=sinc(t/pi); plot(t,x); title('抽样信号'); （2）实验结果
4.单位样值序列（1）源程序 clear;clc; n1=input('n1='); n2=('n2='); n=n1:n2; k=length(n); x1=zeros(1,k); x1(1,-n1+1)=1 subplot(1,2,1); stem(n,x1,'filled') （2）实验结果
雷达系统中杂波信号的建模与仿真
1.雷达系统中杂波信号的建模与仿真目的 雷达的基本工作原理是利用目标对雷达波的散射特性探测和识别目标。然而目标存在于周围的自然环境中，环境对雷达电磁波也会产生散射，从而对目标信号的检测产生干扰，这些干扰就称为雷达杂波。对雷达杂波的研究并通过相应的信号处理技术可以最大限度的压制杂波干扰，发挥雷达的工作性能。 雷达研制阶段的外场测试不仅耗费大量的人力、物力和财力，而且容易受大气状况影响，延长了研制周期。随着现代数字电子技术和仿真技术的发展，计算机仿真技术被广泛应用于包括雷达系统设计在内的科研生产的各个领域，在一定程度上可以替代外场测试，降低雷达研制的成本和周期。 长期以来，由于对杂波建模与仿真的应用己发展了多种杂波类型和多种建模与仿真方法。然而却缺少一个集合了各种典型杂波产生的成熟的软件包，雷达系统的研究人员在需要用到某一种杂波时，不得不亲自动手，从建立模型到计算机仿真，重复劳动，造成了大量的时间和人力的浪费。因此，建立一个雷达杂波库，就可以使得科研人员在用到杂波时无需重新编制程序，而直接从库中调用杂波生成模块，用来产生杂波数据或是用来构成雷达系统仿真模型，在节省时间和提高仿真效率上的效益是十分可观的。 从七十年代至今已经公布了很多杂波模型，其中有几类是公认的比较合适的模型。而且，杂波建模与仿真技术的发展己有三十多年的历史，己经有了一些比较成熟的理论和行之有效的方法，这就使得建立雷达杂波库具有可行性。 为了能够反映雷达信号处理机的真实性能，同时为改进信号处理方案提供理论依据，雷达杂波仿真模块输出的杂波模拟信号应该能够逼真的反映对象环境的散射环境。模拟杂波的一些重要散射特性影响着雷达对目标的检测和踉踪性能，比如模拟杂波的功率谱特性与雷达的动目标显示滤波器性能有关;模拟杂波的幅度起伏特性与雷达的恒虚警率检测处理性能有关。因此，杂波模拟方案的设计是雷达仿真设计中极其重要的内容，杂波模型的精确性、通用性和灵活性是衡量杂波产生模块的重要指标。 2.Simulink简介 Simulink是MATLAB最重要的组件之一，它提供一个动态系统建模、仿真和
DSB信号的仿真分析
《MATLAB课程设计》报告题目：基于MATLAB的DSB调制与解调分析专业班级: 通信1104班 学生姓名： 指导教师：
MATLAB课程设计任务书 学生姓名：专业班级： 指导教师：工作单位： 题目: 基于MATLAB的DSB调制与解调分析 设计内容和要求 DSB信号的仿真分析 调制信号：分别为300Hz正弦信号和矩形信号；载波频率：30kHz； 解调：同步解调； 要求：画出以下三种情况下调制信号、已调信号、解调信号的波形、频谱以及解调器输入输出信噪比的关系曲线； 1）调制信号幅度=×载波幅度；2）调制信号幅度=载波幅度； 3）调制信号幅度=×载波幅度； 时间安排 2013年12月25日：复习DSB的原理，初步构想设计的流程。 2013年12月26日至28日：程序编写及调试。 2013年12月29日：写报告。 指导教师签名：年月日
目录
摘要 调制在通信系统中有十分重要的作用。通过调制，不仅可以进行频谱搬移，把调制信号的频谱搬移到所希望的位置上，从而将调制信号转换成适合于传播的已调信号，而且它对系统的传输有效性和传输的可靠性有着很大的影响，调制方式往往决定了一个通信系统的性能。MATLAB软件广泛用于数字信号分析，系统识别，时序分析与建模，神经网络、动态仿真等方面有着广泛的应用。本课题利用MATLAB软件对DSB 调制解调系统进行模拟仿真，分别利用300HZ正弦波和矩形波，对30KHZ正弦波进行调制，观察调制信号、已调信号和解调信号的波形和频谱分布，并在解调时引入高斯白噪声，对解调前后信号进行信噪比的对比分析，估计DSB调制解调系统的性能。 Abstract Modulation in communication systems have an important role. Through the modulation, not only can move the spectrum, the modulated signal spectrum move to the desired position, which will convert into a modulated signal suitable for transmission of modulated signals, and that its transmission system, the effectiveness and reliability of transmission has a great impact, the modulation method is often decided on a communication system performance. MATLAB software is widely used in digital signal analysis, system identification, time series analysis and modeling, neural networks, dynamic simulation have a wide range of applications. This topic using MATLAB software DSB modulation and demodulation system simulation, use, respectively, 300HZ sine wave and rectangular wave, sine wave modulation of the 30KHZ observed modulated signal modulated signal and demodulate the signal waveform and spectrum distribution, and in the solution white Gaussian noise introduced when adjusted for demodulating the signal-noise ratio before and after the comparative analysis, it is estimated DSB modulation and demodulation performance of the system.
交通仿真实验报告
交通仿真实验报告 篇一：交通仿真实验报告 目录 1 上机性质与目的.................................. 2 2 上机内容....................................... 2 3 交叉口几何条件、信号配时和交通流数据描述.......... 3 3.1 交叉口几何数据................................ 3 3.2 交叉口信号配时系统............................ 3 3.3 交叉口交通流数据.............................. 4 4 交叉口交通仿真.................................. 4 4.1 交通仿真步骤.................................. 4 4.2 二维输出..................................... 13 4.3 3D输出...................................... 14 5 仿真结果分析................................... 15 6 实验总结和体会 (15) 实验上机名称：信号交叉口仿真 1 上机性质与目的 本实验属于计算机仿真实验，借助仿真系统模拟平面信号交叉口场景，学生将完成从道路条件设计到信号相位配置等一系列仿真实验。 实验目的： 1. 了解平面信号交叉口在城市交通中的地位； 2. 了解平面信号交叉口的主要形式、规模等基本情况； 3. 了解交叉口信号相位配时及对交叉口通行能力的影响；
随机信号分析实验报告二 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仿真实验报告
MATLA仿真实验报告 学院：计算机与信息学院 课程：—随机信号分析 姓名： 学号： 班级： 指导老师： 实验一
题目：编写一个产生均值为1,方差为4的高斯随机分布函数程序, 求最大值，最小值，均值和方差，并于理论值比较。 解：具体的文件如下，相应的绘图结果如下图所示 G仁random( 'Normal' ,0,4,1,1024); y=max(G1) x=mi n(G1) m=mea n(G1) d=var(G1) plot(G1);
实验二 题目：编写一个产生协方差函数为CC)=4e":的平稳高斯过程的程序，产生样本函数。估计所产生样本的时间自相关函数和功率谱密度，并求统计自相关函数和功率谱密度，最后将结果与理论值比较。 解：具体的文件如下，相应的绘图结果如下图所示。 N=10000; Ts=0.001; sigma=2; beta=2; a=exp(-beta*Ts); b=sigma*sqrt(1-a*a); w=normrnd(0,1,[1,N]); x=zeros(1,N); x(1)=sigma*w(1); for i=2:N x(i)=a*x(i-1)+b*w(i); end %polt(x); Rxx=xcorr(x0)/N; m=[-N+1:N-1]; Rxx0=(sigma A2)*exp(-beta*abs(m*Ts)); y=filter(b,a,x) plot(m*Ts,RxxO, 'b.' ,m*Ts,Rxx, 'r');
periodogram(y,[],N,1/Ts); 文件旧硯化）插入（1〕 ZMCD 克闻〔D ］窗口曲） Frequency (Hz) 50 100 150 200 250 300 350 400 450 500 NH---.HP)&UO 二 balj/ 」- □歹
五款信号完整性仿真分析工具
SI 五款信号完整性仿真工具介绍 （一）Ansoft公司的仿真工具 现在的高速电路设计已经达到GHz的水平，高速PCB设计要求从三维设计理论出发对过孔、封装和布线进行综合设计来解决信号完整性问题。高速PCB 设计要求中国工程师必须具备电磁场的理论基础，必须懂得利用麦克斯韦尔方程来分析PCB设计过程中遇到的电磁场问题。目前，An soft公司的仿真工具能够从三维场求解的角度出发，对PCB 设计的信号完整性问题进行动态仿真。 Ansoft 的信号完整性工具采用一个仿真可解决全部设计问题： Slwave是一种创新的工具，它尤其适于解决现在高速PCB和复杂IC封装中普遍存在的电源输送和信号完整性问题。 该工具采用基于混合、全波及有限元技术的新颖方法，它允许工程师们特性化同步开关噪声、电源散射和地散射、谐振、反射以及引线条和电源/地平面之间的耦合。该工具采用一个仿真方案解决整个设计问题，缩短了设计时间。 它可分析复杂的线路设计，该设计由多重、任意形状的电源和接地层，以及任何 数量的过孔和信号引线条构成。仿真结果采用先进的3D 图形方式显示，它还可产生等效电路模型，使商业用户能够长期采用全波技术，而不必一定使用专有仿 （二）SPECCTRAQuest Cade nee的工具采用Sun的电源层分析模块: Cade nee Design System 的SpeeetraQuest PCB信号完整性套件中的电源完整性模块据称能让工程师在高速PCB设计中更好地控制电源层分析和共模EMI 。 该产品是由一份与Sun Microsystems公司签署的开发协议而来的，Sun最初研制该项技术是为了解决母板上的电源问题。 有了这种新模块，用户就可根据系统要求来算出电源层的目标阻抗；然后基于板上的器件考虑去耦合要求，Shah表示，向导程序能帮助用户确定其设计所要求的去耦合电容的数目和类型；选择一组去耦合电容并放置在板上之后，用户就可运行一个仿真程序，通过分析结果来发现问题所在。 SPECCTRAQuest是CADENCE公司提供的高速系统板级设计工具，通过它可以控制与PCB layout相应的限制条件。在SPECCTRAQuest菜单下集成了一下工具： （1）SigXplorer 可以进行走线拓扑结构的编辑。可在工具中定义和控制延时、特性阻抗、驱动和负载的类型和数量、拓扑结构以及终端负载的类型等等。可在
MATLAB通信系统仿真实验报告1
MATLAB通信系统仿真实验报告
实验一、MATLAB的基本使用与数学运算 目的：学习MATLAB的基本操作，实现简单的数学运算程序。 内容： 1-1要求在闭区间[0，2π]上产生具有10个等间距采样点的一维数组。试用两种不同的指令实现。 运行代码：x=[0:2*pi/9:2*pi] 运行结果： 1-2用M文件建立大矩阵x x=[0.10.20.30.40.50.60.70.80.9 1.11.21.31.41.51.61.71.81.9 2.12.22.32.42.52.62.72.82.9 3.13.23.33.43.53.63.73.83.9] 代码：x=[0.10.20.30.40.50.60.70.80.9 1.11.21.31.41.51.61.71.81.9 2.12.22.32.42.52.62.72.82.9 3.13.23.33.43.53.63.73.83.9] m_mat 运行结果： 1-3已知A=[5，6;7，8]，B=[9，10;11,12]，试用MATLAB分别计算 A+B,A*B,A.*B,A^3，A.^3，A/B,A\B. 代码：A=[56;78]B=[910;1112]x1=A+B X2=A-B X3=A*B X4=A.*B X5=A^3 X6=A.^3X7=A/B X8=A\B
运行结果： 1-4任意建立矩阵A，然后找出在[10,20]区间的元素位置。 程序代码及运行结果： 代码：A=[1252221417;111024030;552315865]c=A>=10&A<=20运行结果： 1-5总结：实验过程中，因为对软件太过生疏遇到了些许困难，不过最后通过查书与同学交流都解决了。例如第二题中，将文件保存在了D盘，而导致频频出错，最后发现必须保存在MATLAB文件之下才可以。第四题中，逻辑语言运用到了ij，也出现问题，虽然自己纠正了问题，却也不明白错在哪了，在老师的讲解下知道位置定位上不能用ij而应该用具体的整数。总之第一节实验收获颇多。
随机信号实验报告
随机信号分析 实验报告 目录 随机信号分析 (1) 实验报告 (1) 理想白噪声和带限白噪声的产生与测试 (2) 一、摘要 (2) 二、实验的背景与目的 (2) 背景： (2) 实验目的: (2) 三、实验原理 (3) 四、实验的设计与结果 (4) 实验设计: (4) 实验结果: (5) 五、实验结论 (12) 六、参考文献 (13) 七、附件 (13) 1
理想白噪声和带限白噪声的产生与测试一、摘要 本文通过利用MATLAB软件仿真来对理想白噪声和带限白噪声进行研究。理想白噪声通过低通滤波器和带通滤波器分别得到低通带限白噪声和帯通带限白噪声。在仿真的过程中我们利用MATLAB工具箱中自带的一些函数来对理想白噪声和带限白噪声的均值、均方值、方差、功率谱密度、自相关函数、频谱以及概率密度进行研究，对对它们进行比较分析并讨论其物理意义。 关键词：理想白噪声带限白噪声均值均方值方差功率谱密度自相关函数、频谱以及概率密度 二、实验的背景与目的 背景： 在词典中噪声有两种定义：定义1：干扰人们休息、学习和工作的声音，引起人的心理和生理变化。定义2：不同频率、不同强度无规则地组合在一起的声音。如电噪声、机械噪声，可引伸为任何不希望有的干扰。第一种定义是人们在日常生活中可以感知的，从感性上很容易理解。而第二种定义则相对抽象一些，大部分应用于机械工程当中。在这一学期的好几门课程中我们都从不同的方面接触到噪声，如何的利用噪声，把噪声的危害减到最小是一个很热门的话题。为了加深对噪声的认识与了解，为后面的学习与工作做准备，我们对噪声进行了一些研究与测试。 实验目的: 了解理想白噪声和带限白噪声的基本概念并能够区分它们，掌握用MATLAB 或c/c++软件仿真和分析理想白噪声和带限白噪声的方法，掌握理想白噪声和带限白噪声的性质。
信号仿真实验报告
信号与系统仿真实验报告
实验一 (1)()t δ Function-M 文件 function [x,t]=dirac(t1,t2,t0) %y=dirac(t-t0),t1
分析：在新版的Matlab 函数库中有自带的阶跃函数，调用方法为f=heaviside(t)，这里为了方便画位移后0()t t ε-的图像，故自定义了一个阶跃函数。 (3)指数 ①a=1； >> f=sym('exp(t)'); >> ezplot(f,[-3,3]) >> xlabel('时间t') >> ylabel('函数f （x ）') ②a=-1； f=sym('exp((-1)*t)'); >> ezplot(f,[-3,3]) >> xlabel('时间t') >> ylabel('函数f （x ）') 图a ）a=1时的指数信号图像 图b ）a=-1时的指数函数图像 分析：y=sym （‘f （x ）’）是用了符号运算法 （4）(),5N R t N = >> t=-1:0.001:10; >> y=heaviside(t,0)-heaviside(t,5); >> plot(t,y) >> axis([0,10,-0.2,1.2]) 分析：采用两个跳变点不等的阶跃函数相减得到一个矩形函数的方法生成的门函数。
信号实验报告
大连理工大学 本科实验报告 课程名称：信号与系统实验 学院(系）：电子信息与电气工程学部专业: 通信工程 班级: 1401班 学号：20148３091 学生姓名：李睿 20１6年 5 月2１日 ?实验项目列表
?大连理工大学实验预习报告 学院（系）：电信专业:通信工程班级:1401班 姓名:李睿学号：20148309１组:5 _＿_ 实验时间：2016、5、6 实验室：创新园大厦c0221 实验台： 5 指导教师签字：成绩： 信号得频谱图 一、实验目得与要求 1、掌握周期信号得傅里叶级数展开 ２、掌握周期信号得有限项傅里叶级数逼近 3、掌握周期信号得频谱分析 4、掌握连续非周期信号得傅立叶变换 5、掌握傅立叶变换得性质 二、实验用得mａｔlａb命令与例子
１、ａ:b：ｃ：产生一个从ａ到 c，间隔为ｂ得等间隔数列例：5：1:11,产生一个从 5 到11,间隔为 1 得等间隔数列 ２、ｑuａrｅ（t，dutｙ)：周期性矩形脉冲信号(duty 表示占空比）调用形式： y=squaｒe(ｔ，dｕtｙ）例：产生一个周期为2π，幅值为±１得周期性方波。y＝sｑuarｅ(2＊pi＊30*t,７5）； plot（ｔ,y）,grｉd on aｘis（[—０、１,０、1，—1、5，1、5]） 3、plot（）：ｍaｔlab 中二维线画图函数ｐloｔ（x，ｙ，’颜色与标识’）：若 y 与ｘ为同维向量，则以ｘ为横坐标，y 为纵坐标绘制连线图. 若x 就是向量，y 就是行数或列数与ｘ长度相等得矩阵,则绘制多条不同色彩得连线图，x 被作为这些曲线得共同横坐标.若 x 与 y 为同型矩阵，则以ｘ，ｙ对应元素分别绘制曲线,曲线条数等于矩阵列数. 例:在０≤x≤2π区间内，绘制曲线 y=2ｅ-0、5xcoｓ(４πx）。 ｘ=0：2＊pi； y＝２＊eｘp(－０、5*x）、＊cos（４＊ｐｉ＊x); plot（x，y) ‘’：y 黄ｍ紫 c 青 r 红 g 绿 b 蓝ｗ白 k 黑—实线、点 <小于号 :点线ｏ圆s 正方形 -、点划线x 叉号 d 菱形- -虚线 +加号h 六角星 *星号 p 五角星 v 向下三角形 ^向上三角形〉大于号 4、grｉd ｏｎ:有网格 gｒｉｄ off:关掉格网下面就是加上命令grｉd on后画得图，有网格. 5、 aｘis（[a b c d］）:表明图线得x轴范围为a～by轴范围为c～ｄ例：ｐlｏt(x,ｙ）ａｘis(［０ 1 ２３]) ｇｒｉｄ oｎ ６、 lengｔｈ（a)：表示矩阵ａ得最大得长度比如lｅngth([1 2 3；4 5 ６]) 等于3，因为2行与3列中最大就是3。当a就是向量时,即表示向量得元素个数，因为向量总就是１×n或ｎ×1得，而n一定大于或等于1、所以得到得结果一定就是n. ７、 1、/tan(pi、*x)：表示点乘。点乘就是值对值得运算上面得式子中 X 可能就是一个向量或矩阵，PＩ后面得点就是一个ＰI 与一个向量相乘，得到得也就是一个向量;1 后面乘得自然也就是个向量所以要加点,也就就是对应不同得Ｘ，有不同得 Y 值. 8.fｉｇurｅ就是建立图形得意思. 系统自动从 1，2，3,4、、、来建立图形，数字代表第几幅图形,figure（1)，figurｅ（２）就就是第一第二副图得意思,在建立图形
随机信号分析实验报告(基于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。那么，
五款信号完整性仿真分析工具
SI五款信号完整性仿真工具介绍 (一)Ansoft公司的仿真工具 现在的高速电路设计已经达到GHz的水平，高速PCB设计要求从三维设计理论出发对过孔、封装和布线进行综合设计来解决信号完整性问题。高速PCB设计要求中国工程师必须具备电磁场的理论基础，必须懂得利用麦克斯韦尔方程来分析PCB设计过程中遇到的电磁场问题。目前，Ansoft公司的仿真工具能够从三维场求解的角度出发，对PCB设计的信号完整性问题进行动态仿真。 Ansoft的信号完整性工具采用一个仿真可解决全部设计问题： SIwave是一种创新的工具，它尤其适于解决现在高速PCB和复杂IC封装中普遍存在的电源输送和信号完整性问题。 该工具采用基于混合、全波及有限元技术的新颖方法，它允许工程师们特性化同步开关噪声、电源散射和地散射、谐振、反射以及引线条和电源/地平面之间的耦合。该工具采用一个仿真方案解决整个设计问题，缩短了设计时间。 它可分析复杂的线路设计，该设计由多重、任意形状的电源和接地层，以及任何数量的过孔和信号引线条构成。仿真结果采用先进的3D图形方式显示，它还可产生等效电路模型，使商业用户能够长期采用全波技术，而不必一定使用专有仿真器。 (二)SPECCTRAQuest Cadence的工具采用Sun的电源层分析模块： Cadence Design Systems的SpecctraQuest PCB信号完整性套件中的电源完整性模块据称能让工程师在高速PCB设计中更好地控制电源层分析和共模EMI。 该产品是由一份与Sun Microsystems公司签署的开发协议而来的，Sun最初研制该项技术是为了解决母板上的电源问题。 有了这种新模块，用户就可根据系统要求来算出电源层的目标阻抗；然后基于板上的器件考虑去耦合要求，Shah表示，向导程序能帮助用户确定其设计所要求的去耦合电容的数目和类型；选择一组去耦合电容并放置在板上之后，用户就可运行一个仿真程序，通过分析结果来发现问题所在。 SPECCTRAQuest是CADENCE公司提供的高速系统板级设计工具，通过它可以控制与PCB layout相应的限制条件。在SPECCTRAQuest菜单下集成了一下工具： （1）SigXplorer可以进行走线拓扑结构的编辑。可在工具中定义和控制延时、特性阻抗、驱动和负载的类型和数量、拓扑结构以及终端负载的类型等等。可在