实验三 随机信号通过线性时不变系统
第5章 随机信号通过线性系统_20120509

5.2 平稳白噪声通过LTI系统
2019/2/9
22
2019/2/9
23
2019/2/9
24
2019/2/9
25
2019/2/9
26
2019/2/9
27
2019/2/9
28
2019/2/9
29
2019/2/9
30
2019/2/9
31
2019/2/9
32
2019/2/9
33
2019/2/9
L a1x1 (t ) a2 x2 (t ) a1L x1 (t ) a2 L x2 (t )
时不变性:L x(t ) y(t )
系统完全由算子L[ ]确定。
2019/2/9 4
线性时不变系统的其他表示形式
线性时不变系统的输出可看成是输入信号的冲击响 应(严格说是零状态响应),表示为:
2019/2/9
10
证明:(1)
E Y (t ) E X (t u) h(u)du mX h(u )du mX H ( j 0)
其中 H ( j 0) h(u)du (2)
E Y (t1 ) X (t2 ) E X (t1 u ) X (t2 )h(u )du
y(t ) x(t ) h(t )
x(t u)h(u)du
Y ( jw) X ( jw) H ( jw)
其中 y(t ) Y ( jw)
x(t ) X ( jw)
h(t ) H ( jw)
2019/2/9
5
线性系统的输出过程
当输入信号X(t,ξi)是随机过程时,系统的输出Y(t,ξi) 映射为Y(t,ξi)=L[X(t,ξi)]。由于ξi在X(t)对应的整个 样本空间取值,使Y(t,ξi)成为随机函数,构造出一个 新的随机过程。 若输入信号的均方值E|X(t)|2存在,则该类信号是稳定 的。对于稳定的LTI系统,当输入信号稳定时,其输出 信号也一定是稳定的。稳定系统满足 h(u ) du
实验三 随机过程通过线性系统

实验名称线性系统对随机过程的响应一、实验目的通过本仿真实验了解正态白色噪声随机过程通过线性系统后相关函数以及功率谱的变化;培养计算机编程能力。
二、实验平台MATLAB R2014a三、实验要求(1)运用正态分布随机数产生函数产生均值为m=0,根方差σ=1的白色正态分布序列{u(n)|n=1,2,…,2000},画出噪声u(n)的波形图。
(2)设离散时间线性系统的差分方程为x(n)=u(n)-0.36u(n-1)+0.85u(n-2)(n=3,4,…,2000).画出x(n)的波形图。
(3)随机过程x(n)的理论上的功率谱函数为在[0,π]范围内对w进行采样,采样间隔0.001π,计算S(i×0.001π) (i=1,2,…,1000);画出波形图。
(4)根据步骤二产生的数据序列x(n)计算相关函数的估计值与理论值1.1296、-0.666、0.85、0、0、0的差异。
(5)根据相关函数的估计值对随机过程的功率谱密度函数进行估计在[0,π]范围内对w进行采样,采样间隔0.001π,计算S(i×0.001π) (i=1,2,…,1000);画出波形图,比较其与理论上的功率谱密度函数S(w)的差异。
(6)依照实验1的方法统计数据x(n)在不同区间出现的概率,计算其理论概率,观察二者是否基本一致。
四、实验代码及结果A、运用正态分布随机数产生函数产生均值为m=0,根方差σ=1的白色正态分布序列{u(n)|n=1,2,…,2000},画出噪声u(n)的波形图。
代码实现:波形图:分析:运用正态分布随机数产生函数产生均值为0,根方差σ=1的白色噪声样本序列。
B、设离散时间线性系统的差分方程为x(n)=u(n)-0.36u(n-1)+0.85u(n-2)(n=3,4,…,2000).画出x(n)的波形图。
代码实现:波形图:分析:正态随机序列通过离散时间线性系统生成的仍是正态随机序列。
C、随机过程x(n)的理论上的功率谱函数为在[0,π]范围内对w进行采样,采样间隔0.001π,计算S(i×0.001π) (i=1,2,…,1000);画出波形图。
随机信号分析-3 随机信号通过线性系统

X(t ) X (t ) RX ( )
则输出 Y(t ) 的时间平均
1 Y(t ) lim T 2T 1 lim T 2T
T
T
T
Y (t )dt
0
[
T
h(u ) X (t u )du ]dt
1 T [ lim X(t u )dt ]h(u )du 0 T 2T T mX h(u )du mY
(1)线性性: L[ax1 (t ) bx2 (t )] aL[ x1 (t )] bL[ x2 (t )] L[] 称作算子 稳定系统:对于任意有界输入,响应也是有界的
4
连续时不变线性系统的分析方法 1. 时域分析
y(t ) x(t )h( )d x( )h(t )d x(t ) h(t )
RXY (t1 , t2 ) RX (t1, t2 )* h(t2 )
RYX (t1 , t2 ) RX (t1, t2 )* h(t1 )
9
时域分析法
证明:
RXY (t1 , t2 ) RX (t1, t2 )* h(t2 )
RYX (t1 , t2 ) RX (t1, t2 )* h(t1 )
14
时域分析法
结论1:若输入是
X(t )宽平稳的,则系统输出Y(t )也是宽平
稳的,且输入与输出联合宽平稳。
证明:若输入 X(t )为宽平稳随机过程,则有:
mX (t ) mX 常数
RX (t1 , t2 ) RX ( )
=t 2 t1
RX (0) E[ X 2 (t )]
h( ) E[X(t )]d
实验三线性系统

实验三 信号通过线性系统的特性分析一 实验目的1.掌握无失真传输的概念及无失真传输的线性系统满足的条件2.分析无失真传输的线性系统输入、输出频谱特性,给出系统的频谱特性 3.掌握系统幅频特性的测试及绘制方法二 实验原理通过频谱分析可以看出,在一般情况下线性系统的响应波形与激励波形是不同的,即:信号在通过线性系统传输的过程中产生了失真。
线性系统引起的信号失真是由两方面的因素造成的,一是系统对信号中各频率分量的幅度产生不同程度的衰减,使响应各频率分量的相对幅度产生变化,造成幅度失真;一是系统对各频率分量产生的相移不与频率成正比,使响应各频率分量在时间轴上的相对位置产生变化,造成相位失真。
线性系统的幅度失真与相位失真都不产生新的频率分量。
对于非线性系统,由于其非线性特性,对于传输信号产生非线性失真,非线性失真可能产生新的频率分量。
如果信号在传输过程中不失真,则响应)(t r 与激励)(t e 波形相同,只是幅度大小或出现的时间不同。
激励与响应的关系可表示为)()(0t t ke t r -= (3-15)为了实现信号无失真传输,线性系统应该满足什么条件?由式(3-1)得)()(t j ej kE j R ωωω-= (3-16)设)(t e 与)(t r 的傅里叶变换分别是)(ωj E 和)(ωj R ,则)()()(ωωωj E j H j R = (3-17)比较式(3-2)与式(3-3),在信号无失真传输时,系统函数应为)()()(t j j keej H j H ωωφωω-== (3-18)因此,为了实现任意信号通过线性系统不产生波形失真,该系统应满足以下两个理想条件,如图3-1。
⎩⎨⎧-==0)()(t kj H ωωφω图3.13 理想线性传输系统的系统函数的频率特性很显然,在传输有限频宽的信号时,上述的理想条件可以放宽,只要在信号占有频带范围内系统满足上述理想条件即可。
三 实验方法 实验电路采用如下电路图3.14 实验电路22211122221111)()()(C R j R C R j R C R j R j U j U j H ωωωωωω++++==若2211C R C R =,则212)(R R R j H +=ω,0)(=ωφ,该系统满足无失真传输的条件。
信号与系统实验之连续线性时不变系统的分析

信号与系统实验报告连续线性时不变系统的分析专业:电子信息工程(实验班)姓名:曾雄学号:14122222203班级:电实12-1BF目录一、实验原理与目的 (3)二、实验过程及结果测试 (3)三、思考题 (10)四、实验总结 (10)五、参考文献 (11)一、实验原理与目的深刻理解连续时间系统的系统函数在分析连续系统的时域特性、频域特性及稳定性中的重要作用及意义。
掌握利用MATLAB 分析连续系统的时域响应、频响特性和零极点的基本方法。
二、实验过程及结果测试1.描述某线性时不变系统的微分方程为: ''()3'()2()'()y t y t y t f t f t++=+ 且f(t)=t 2,y(0-)=1,y ’(0-)=1;试求系统的单位冲激响应、单位阶跃响应、全响应、零状态响应、零输入响应、自由响应和强迫响应。
编写相应MATLAB 程序,画出各波形图。
(1)单位冲激响应: 程序如下:%求单位冲激响应a=[1,3,2]; b=[1,2]; sys=tf(b,a); t=0:0.01:10; h=impulse(sys,t);%用画图函数plot( )画单位冲激响应的波形plot(h); %单位冲激响应曲线 xlabel('t'); ylabel('h');title('单位冲激响应h(t)') 程序运行所得波形如图一:200400600800100012000.10.20.30.40.50.60.70.80.91th单位冲激响应h(t )图一 单位冲激响应的波形(2)单位阶跃响应: 程序如下:%求单位阶跃响应a=[1,3,2]; b=[1,2]; sys=tf(b,a); t=0:0.01:10; G=step(sys,t);%用画图函数plot( )画单位阶跃响应的波形plot(G); %单位阶跃响应曲线 xlabel('t'); ylabel('g');title('单位阶跃响应g(t)') 程序运行所得波形如图二:2004006008001000120000.10.20.30.40.50.60.70.80.91tg单位阶跃响应g(t )图二 单位阶跃响应的波形 (3)零状态响应: 程序如下:%求零状态响应yzs=dsolve('D2y+3*Dy+2*y=2*t+2*t^2','y(0)=0,Dy(0)=0') %用符号画图函数ezplot( )画各种响应的波形 t=0:0.01:3;ezplot(yzs,t); %零状态响应曲线 axis([0,3,-1 5]);title('零状态响应曲线yzs'); ylabel('yzs');程序运行所得波形如图三:00.511.522.53-112345t零状态响应曲线yzsy z s图三 零状态响应的波形(4)零输入响应: 程序如下:%求零输入响应yzi=dsolve('D2y+3*Dy+2*y=0','y(0)=1,Dy(0)=1') %用符号画图函数ezplot( )画零输入响应的波形 t=0:0.01:3;ezplot(yzi,t);%零输入响应曲线 axis([0,3,-1,2]); title('零输入响应yzi'); ylabel('yzi');程序运行所得波形如图四:图四 零输入响应的波形(5)全响应:程序如下:%求全响应y=dsolve('D2y+3*Dy+2*y=2*t+2*t^2','y(0)=1,Dy(0)=1') %用符号画图函数ezplot( )画全响应响应的波形00.511.522.53-1-0.50.511.52t零输入响应yziy z it=0:0.01:3;ezplot(y,t); %全响应曲线 axis([0,3,-1,5]); title('全响应y'); ylabel('y');程序运行所得波形如图五:00.511.522.53-112345t全响应yy图五 全响应的波形(6)自由响应:程序如下:%自由响应y=dsolve('D2y+3*Dy+2*y=2*t+2*t^2','y(0)=1,Dy(0)=1'); %全响应 yht=dsolve('D2y+3*Dy+2*y=0','y(0)=1,Dy(0)=1'); % 求齐次通解yt=dsolve('D2y+3*Dy+2*y=2*t+2*t^2','y(0)=0,Dy(0)=0'); % 求非齐次通解 yp=yt-yht;yh=y-yp; % 求齐次解,即自由响应 t=0:0.01:3; ezplot(yh,t); title('自由响应yh'); ylabel('yh');程序运行所得波形如图六:0.511.522.530.511.52t自由响应yhy h图六 自由响应的波形(7)强迫响应: 程序如下:%强迫响应yht=dsolve('D2y+3*Dy+2*y=0','y(0)=1,Dy(0)=1'); % 求齐次通解yt=dsolve('D2y+3*Dy+2*y=2*t+2*t^2','y(0)=0,Dy(0)=0'); % 求非齐次通解 yp=yt-yht; % 求特解,即强迫响应 t=0:0.01:3; ezplot(yp,t); title('强迫响应yp'); ylabel('yp');程序运行所得波形如图七:0.511.522.53-112345t强迫响应ypy p图七 强迫响应的波形2.给定一个连续线性时不变系统,描述其输入输出之间关系的微分方程为:编写MATLAB 程序,绘制系统的幅频响应、相频响应、频率响应的实部和频率响应的虚部的波形,确定滤波器的类型。
随机信号分析报告实验:随机过程通过线性系统地分析报告

实验三 随机过程通过线性系统的分析实验目的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)试推导系统输出的功率谱密度、相关函数、相关时间和系统的等效噪声带宽。
第三章 随机信号通过线性系统分析

• • • • • 1、输出表达式(零状态响应,因果系统) 2、输出的均值 3、系统输入与输出之间的互相关函数 4、系统输出的自相关函数 5、系统输出的高阶距
x (t ) ► 输入为随机信号X(t)的某个实验结果的一个样本函数,则输 出为:
y (t )
h ( ) x ( t ) d
2012-6-30 3
3.1 线性系统的基本理论
系统可分为: (1)线性系统:线性放大器、线性滤波器 (2)非线性系统:限幅器、平方律检波器 对于线性系统:已知系统特性和输入信号的统计特性,可以求出系统输 出信号的统计特性
2012-6-30
4
• 下面的分析线性系统是单输入单输出(响应)的、连续或离散时不变 的、物理可实现的稳定系统。
证明:已知系统输入随机信号的自相关函数,可以求出系统 输出端的自相关函数
R Y ( t1 , t 2 ) E [ Y ( t1 ) Y ( t 2 )] h ( t1 ) h ( t 2 ) R X ( t1 , t 2 )
R Y ( t1 , t 2 ) E [Y ( t1 )Y ( t 2 )]
R Y X ( t1 , t 2 ) R X ( t1 , t 2 ) * h ( t1 )
2012-6-30 17
3.2 随机信号通过连续时间系统的分析
证明:由于系统的输出是系统输入的作用结果,因此,系统 输入输出之间是相关的,系统输入输出相关函数为
R X Y ( t1 , t 2 ) R X ( t 1 , t 2 ) * h ( t 2 )
时不变线性系统
若输入信号x(t)时移时间C, 输出y(t)也只引起一个相同 的时移,即 y(t-C) = L[x(t-C)]
第三章 随机信号通过线性系统分析讲解

第三章 随机信号通过线性系统的分析本章主要内容:● 线性系统的基本理论● 随机信号通过连续时间系统的分析 ● 随机信号通过离散时间系统的分析 ● 色噪声的产生与白化滤波器 ● 等效噪声带宽 ● 解析过程● 窄带随机过程基本概念● 窄带高斯过程包络与相位的概率密度 ● 窄带高斯过程包络平方的概率密度3.1随机信号通过连续时间系统的分析在给定系统的条件下,输出信号的某个统计特性只取决于输入信号的相应的统计特性。
分析方法:卷积积分法;频域法。
3.1.1、时域分析法1、输出表达式(零状态响应,因果系统) 输入为随机信号)(t X 某个实验结果ζ的一个样本函数),(ζt x ,则输出),(ζt y 为:对于所有的ζ,输出为一族样本函数构成随机过程Y(t):2. 输出的均值:)(*)()(t h t m t m X Y =证明:3.系统输入与输出之间的互相关函数)(*),(),(22121t h t t R t t R X XY = )(*),(),(12121t h t t R t t R X YX =证明:4、系统输出的自相关函数已知输入随机信号的自相关函数,求系统输出端的自相关函数。
显然,有:5、系统输出的高阶距输出n阶矩的一般表达式为注意:上面的分析方法是零状态响应的一般分析方法。
它既适用于输入是平稳随机信号的情况,也适用于输入是非平稳的情况。
3.1.2、系统输出的平稳性及其统计特性的计算1、双侧随机信号在这种情况下,系统输出响应在t=0时已处于稳态。
(1)若输入X(t)是宽平稳的,则系统输出Y(t)也是宽平稳的,且输入与输出联合宽平稳。
那么由于假定连续系统是稳定的,所以由于输出的均值是常数,而输出的相关函数只是 的函数,且输出均方值有界。
所以,输出随机过程为宽平稳的。
可总结如下:输出均值:输入与输出间的互相关函数为输出的自相关函数为输出的均方值即输出总平均功率为若用卷积的形式,则可分别写为(2)若输入X(t)是严平稳的,则输出Y(t)也是严平稳的。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
实验三 随机信号通过线性系统的分析
一、实验目的
1 模拟产生特定相关函数的连续随机序列或者离散的随机序列,考察其特性。
2 模拟高斯白噪声环境下信号通过系统的问题,实现低通滤波。
3 掌握系统输出信号的数字特征和功率谱密度的求解。
二、实验设备
1计算机
2 Matlab 软件
三、实验原理
随机信号通过线性系统分析的中心问题是:给定系统的输入函数(或统计特性:均值和 自相关函数)和线性系统的特性,求输出函数。
如下图所示,H 为线性变换,信号X (t )为系统输入, Y (t )为系统的输出,它也是随机信号。
图3.1 随机信号通过系统的示意图
并且满足: H [X (t )] = Y (t )
在时域:
若X(t)时域平稳,系统冲激响应为h(t),则系统输入和输出的关系为:
()()*()()()()()Y t X t h t X h t d h X t d ττττττ∞∞
-∞-∞==-=-⎰⎰ 输出期望:∑∞
===0m X
Y )m (h m )]t (Y [E m 输出的自相关函数:)(h )(h )(R )(R X Y τ*τ-*τ=τ
输出平均功率:⎰⎰∞∞-∞∞--=
τdvdu )u (h )v (h )u v (R )(R X Y 互相关:)()()()()(ττσσσττh R d h R R X X XY *=-=
⎰∞∞-
在频域:
输入与输出的关系:)(H )(X )(Y ωω=ω
输出的功率谱:2X X Y )(H )(S )(H )(H )(S )(S ωω=ωω-ω=ω
功率谱:)(H )(S )(S X XY ωω=ω
四、实验内容与步骤
1已知平稳随机过程X(n)的相关函数为:5),()(22==σδσ
m m R ; 线性系统的单位冲击响应为111,0,)(+-
=≥=实验者学号后两位r k r k h k 。
编写程序求:
1)输入信号的功率谱密度、期望、方差、平均功率;
2)利用时域分析法求输出信号的自相关函数、功率谱密度、期望、方差、平均功率;
3)利用频域分析法求输出信号的自相关函数、功率谱密度、期望、方差、平均功率;
4)利用频域分析法或时域分析法求解输入输出的互相关函数、互功率谱密度。
2 用MATLAB 的sumulink 模拟白噪声通过下图的RC 电路,用示波器观察输入和输出的波形,改变RC 的值,改变电路的时间常数,观察输入输出波形的变化。
图3.2 RC 电路
实验步骤:
步骤一 打开Simulink 仿真窗口,找到相应的模块,连接成相应仿真图;
步骤二 对各个模块进行参数设计;
步骤三 对设计好的仿真模块进行仿真,观察输入和输出的波形图。
步骤四 改变相应的RC的值重新观察结果,进行分析。
五、实验报告要求
1 写出时域分析、频域分析的必要原理,以及求上述特征的必要公式;
2 输出上述各步骤地功率谱密度和相关函数的序列波型,输出各数字特征的值;
3 附上程序和必要的注解;
4 对实验的结果做必要的分析(如时域分析法与频域分析法求解结果的对比等)。