最后得:
)1()3()(2)(---+=n u n u n x n n
3、如图)(),(n x n h 的波形,求他们的线性卷积和L=6的循环卷积,并画出他们的图形。 解:
()()()()()
k y n x n h n x k h n k ∞
=-∞
=*=
-∑
y(0)=
)
0(
)
(n
x
n
h
k
-
∑∞
-∞
=
=1 y(1)=
)
1(
)
(n
x
n
h
k
-
∑∞
-∞
=
=2
y(2)=
)
2(
)
(n
x
n
h
k
-
∑∞
-∞
=
=3 y(3)=
)
3(
)
(n
x
n
h
k
-
∑∞
-∞
=
=4
y(4)=
)
4(
)
(n
x
n
h
k
-
∑∞
-∞
=
=4 y(5)=
)
5(
)
(n
x
n
h
k
-
∑∞
-∞
=
=3
y(6)=
)
6(
)
(n
x
n
h
k
-
∑∞
-∞
= =2 y(7)=
)
7(
)
(n
x
n
h
k
-
∑∞
-∞
= =1
当n>7或n<0时y(n)都等于0图形如图
为求h(n)和x(n)在L=6时的循环卷积应将h(n)补两个零点,将x(n)补一个零点,再进行循环卷积
yc(0)=3,yc(1)=3,y(2)=3,yc(3)=4,yc(4)=4,yc(5)=3
数字信号处理期末试卷!
数字信号处理模拟试题一 一、单项选择题(本大题共10小题,每小题2分,共20分) 在每小题列出的四个备选项中只有一个是符合题目要求的,请将其代码填写在题后的括号内。错选、多选或未选均无分。 1.在对连续信号均匀采样时,要从离散采样值不失真恢复原信号,则采样角频率Ωs与信号最高截止频率Ωc应满足关系(A ) A.Ωs>2Ωc B.Ωs>Ωc C.Ωs<Ωc D.Ωs<2Ωc 2.下列系统(其中y(n)为输出序列,x(n)为输入序列)中哪个属于线性系统?(D) A.y(n)=y(n-1)x(n) B.y(n)=x(n)/x(n+1) C.y(n)=x(n)+1 D.y(n)=x(n)-x(n-1) 3.已知某序列Z变换的收敛域为5>|z|>3,则该序列为(D ) A.有限长序列 B.右边序列 C.左边序列 D.双边序列 4.实偶序列傅里叶变换是(A ) A.实偶序列 B.实奇序列 C.虚偶序列 D.虚奇序列 5.已知x(n)=δ(n),其N点的DFT[x(n)]=X(k),则X(N-1)=(B) A.N-1 B.1 C.0 D.-N+1 6.设两有限长序列的长度分别是M与N,欲通过计算两者的圆周卷积来得到两者的线性卷积,则圆周卷积的点数至少应取(B ) A.M+N B.M+N-1 C.M+N+1 D.2(M+N) 7.下面说法中正确的是(C) A.连续非周期信号的频谱为周期连续函数 B.连续周期信号的频谱为周期连续函数 C.离散非周期信号的频谱为周期连续函数 D.离散周期信号的频谱为周期连续函数 8.下列各种滤波器的结构中哪种不是IIR滤波器的基本结构?(C ) A.直接型 B.级联型 C.频率抽样型 D.并联型 9.下列关于FIR滤波器的说法中正确的是(C) A.FIR滤波器容易设计成线性相位特性
随机信号分析实验报告
一、实验名称 微弱信号的检测提取及分析方法 二、实验目的 1.了解随机信号分析理论如何在实践中应用 2.了解随机信号自身的特性,包括均值、方差、相关函数、频谱及功率谱密度等 3.掌握随机信号的检测及分析方法 三、实验原理 1.随机信号的分析方法 在信号与系统中,我们把信号分为确知信号和随机信号。其中随机信号无确定的变化规律,需要用统计特新进行分析。这里我们引入随机过程的概念,所谓随机过程就是随机变量的集合,每个随机变量都是随机过程的一个取样序列。 随机过程的统计特性一般采用随机过程的分布函数和概率密度来描述,他们能够对随机过程作完整的描述。但由于在实践中难以求得,在工程技术中,一般采用描述随机过程的主要平均统计特性的几个函数,包括均值、方差、相关函数、频谱及功率谱密度等来描述它们。本实验中算法都是一种估算法,条件是N要足够大。 2.微弱随机信号的检测及提取方法 因为噪声总会影响信号检测的结果,所以信号检测是信号处理的重要内容之一,低信噪比下的信号检测是目前检测领域的热点,而强噪声背景下的微弱信号提取又是信号检测的难点。 噪声主要来自于检测系统本身的电子电路和系统外空间高频电磁场干扰等,通常从以下两种不同途径来解决 ①降低系统的噪声,使被测信号功率大于噪声功率。 ②采用相关接受技术,可以保证在信号功率小于噪声功率的情况下,人能检测出信号。 对微弱信号的检测与提取有很多方法,常用的方法有:自相关检测法、多重自相法、双谱估计理论及算法、时域方法、小波算法等。 对微弱信号检测与提取有很多方法,本实验采用多重自相关法。 多重自相关法是在传统自相关检测法的基础上,对信号的自相关函数再多次做自相关。即令: 式中,是和的叠加;是和的叠加。对比两式,尽管两者信号的幅度和相位不同,但频率却没有变化。信号经过相关运算后增加了信噪比,但其改变程度是有限的,因而限制了检测微弱信号的能力。多重相关法将 当作x(t),重复自相关函数检测方法步骤,自相关的次数越多,信噪比提高的越多,因此可检测出强噪声中的微弱信号。
西电随机信号大课后复习
随机信号大作业 班级:02xxxx 姓名:xx
学号:02xxxxx 第一章 1.23上机题:设有随机初相信号X(t)=5cos(t+φ),其中相位φ是在区间(0,2π)上均匀分布的随机变量。试用Matlab编程产生其三个样本函数。 解:程序: clc clear m=unifrnd(0,2*pi,1,10); for k=1:3 t=1:0.1:10; X=5*cos(t+m(k)); plot(t,X); hold on
end title('其三个样本函数'); xlabel('t');ylabel('X(t)'); grid on ;axis tight ; 由 Matlab 产生的三个样本函数如下图所示: 第二章 2.22 上机题:利用Matlab 程序设计一正弦型信号加高斯白噪声的复合信号。 (3)分析复合信号通过理想低通系统后的功率谱密度和相应的幅度分布特性。 1 2 3 4 5 6 7 8 9 10 -4-3-2-101 23 4其三个样本函数 t X (t )
解:取数据如下: 正弦信号的频率为:fc=10HZ,抽样频率为:fs=100HZ; 信号:x=sin(2*pi*fc*t); 高斯白噪声产生复合信号y: y=awgn(x,10); 复合信号y通过理想滤波器电路后得到信号y3 ,通过卷积计算可以得到y3 即:y3=conv2(y,sin(10*t)/(pi*t)); y3的幅度分布特性可以通过傅里叶变换得到Y3(jw)=fft(y3),y3的功率谱密度:G3(w)=Y3(jw).*conj(Y3(jw)/length(Y3(jw)))。 程序: clear all; fs=100; fc=10; n=201; t=0:1/fs:2; x=sin(2*pi*fc*t); y=awgn(x,10); m=50; i=-0.49:1/fs:0.49; for j=1:m R(j)=sum(y(1:n-j-1).*y(j:199),2)/(n-j); Ry(49+j)=R(j);
测试信号处理实验
实验一 离散时间系统的时域分析 一、实验目的 1. 运用MATLAB 仿真一些简单的离散时间系统,并研究它们的时域特性。 2. 运用MATLAB 中的卷积运算计算系统的输出序列,加深对离散系统的差分方程、冲激响应和卷积分析方法的理解。 二、实验原理 离散时间系统其输入、输出关系可用以下差分方程描述: ∑=∑=-=-M k k N k k k n x p k n y d 00] [][ 当输入信号为冲激信号时,系统的输出记为系统单位冲激响应 ][][n h n →δ,则系统响应为如下的卷积计算式: ∑∞ -∞=-= *=m m n h m x n h n x n y ][][][][][ 当h[n]是有限长度的(n :[0,M])时,称系统为FIR 系统;反之,称系统为IIR 系统。在MA TLAB 中,可以用函数y=Filter(p,d,x) 求解差分方程,也可以用函数 y=Conv(x,h)计算卷积。 例1 clf; n=0:40; a=1;b=2; x1= 0.1*n; x2=sin(2*pi*n); x=a*x1+b*x2; num=[1, 0.5,3]; den=[2 -3 0.1]; ic=[0 0]; %设置零初始条件 y1=filter(num,den,x1,ic); %计算输入为x1(n)时的输出y1(n) y2=filter(num,den,x2,ic); %计算输入为x2(n)时的输出y2(n) y=filter(num,den,x,ic); %计算输入为x (n)时的输出y(n) yt= a*y1+b*y2; %画出输出信号 subplot(2,1,1) stem(n,y); ylabel(‘振幅’); title(‘加权输入a*x1+b*x2的输出’);
数字信号处理期末试卷(含答案)全..
数字信号处理期末试卷(含答案) 一、单项选择题(在每小题的四个备选答案中,选出一个正确答案,并将正确答案的序号填在括号内。 1.若一模拟信号为带限,且对其抽样满足奈奎斯特采样定理,则只要将抽样信号通过( )即可完全不失真恢复原信号。 A.理想低通滤波器 B.理想高通滤波器 C.理想带通滤波器 D.理想带阻滤波器 2.下列系统(其中y(n)为输出序列,x(n)为输入序列)中哪个属于线性系统?( ) A.y(n)=x 3(n) B.y(n)=x(n)x(n+2) C.y(n)=x(n)+2 D.y(n)=x(n 2) 3..设两有限长序列的长度分别是M 与N ,欲用圆周卷积计算两者的线性卷积,则圆周卷积的长度至少应取( )。 A .M+N B.M+N-1 C.M+N+1 D.2(M+N) 4.若序列的长度为M ,要能够由频域抽样信号X(k)恢复原序列,而不发生时域混 叠现象,则频域抽样点数N 需满足的条件是( )。 A.N ≥M B.N ≤M C.N ≤2M D.N ≥2M 5.直接计算N 点DFT 所需的复数乘法次数与( )成正比。 A.N B.N 2 C.N 3 D.Nlog 2N 6.下列各种滤波器的结构中哪种不是FIR 滤波器的基本结构( )。 A.直接型 B.级联型 C.并联型 D.频率抽样型 7.第二种类型线性FIR 滤波器的幅度响应H(w)特点( ): A 关于0=w 、π、π2偶对称 B 关于0=w 、π、π2奇对称 C 关于0=w 、π2偶对称 关于=w π奇对称 D 关于0=w 、π2奇对称 关于=w π偶对称 8.适合带阻滤波器设计的是: ( ) A )n N (h )n (h ---=1 N 为偶数 B )n N (h )n (h ---=1 N 为奇数 C )n N (h )n (h --=1 N 为偶数
随机信号分析课程设计报告
随机信号分析课程设计 报告 标准化工作室编码[XX968T-XX89628-XJ668-XT689N]
随机信号分析课程设计报告 题目 学院信息电子技术 专业电子信息工程 班级 15级1班 学籍号 1 姓名朱李伟 指导教师刘文科 信息电子技术学院 2018年6月18日
实验二随机过程的模拟与数字特征 一、实验目的 1. 学习利用MATLAB模拟产生随机过程的方法。 2. 熟悉和掌握特征估计的基本方法及其MATLAB 实现。 二、实验原理 1.正态分布白噪声序列的产生 MATLAB提供了许多产生各种分布白噪声序列的函数,其中产生正态分布白 噪声序列的函数为randn。 函数:randn 用法:x = randn(m,n) 功能:产生 m×n 的标准正态分布随机数矩阵。 如果要产生服从N (,) 分布的随机序列,则可以由标准正态随机序列 产生。如果X ~ N(0,1),则N (,)。 2.相关函数估计 MATLAB提供了函数 xcorr用于自相关函数的估计。 函数:xcorr 用法:c= xcorr (x,y) c= xcorr (x) c= xcorr (x,y ,'opition') c= xcorr (x, ,'opition') 功能:xcorr(x,y) 计算X (n ) 与Y (n)的互相关,xcorr(x)计算X (n )的自相关。 option 选项可以设定为: 'biased' 有偏估计。 'unbiased' 无偏估计。 'coeff' m = 0 时的相关函数值归一化为1。 'none' 不做归一化处理。 3.功率谱估计 对于平稳随机序列X(n),如果它的相关函数满足那么它的功率谱定义为自相关函数R X(m)的傅里叶变换: 功率谱表示随机信号频域的统计特性,有着重要的物理意义。我们实际所 能得到的随机信号的长度总是有限的,用有限长度的信号所得的功率谱只是真实功率谱的估计,称为谱估计或谱分析。功率谱估计的方法有很多种,这里我们介绍基于傅里叶分析的两种通用谱估计方法。 (1)自相关法
工程教育专业认证标准
工程教育专业认证标准(讨论稿) (2011年11月) 1.总则 (1)本标准适用于普通高等学校工程教育本科专业认证。 (2)本标准提供工程教育本科培养层次的基本质量要求。 (3)本标准由通用标准和专业补充标准组成。
2.通用标准 2.1 专业目标 2.1.1 专业设置 专业设置适应国家和地区、行业经济建设的需要,适应科技进步和社会发展的需要,符合学校自身条件和发展规划,有明确的服务面向和人才需求。申请认证或重新认证的专业必须具有: 1.明确充分的专业设置依据和论证,有相应学科作依托,专业口径、布局符合学校的定位。 2.明确的、可衡量、公开的人才培养目标。根据经济建设和社会发展的需要、自身条件和发展潜力,确定在一定时期内培养人才的层次、类型和人才的主要服务面向。 3.至少已有3届毕业生。 2.1.2 毕业生能力 专业必须证明所培养的毕业生达到如下知识、能力与素质的基本要求: 1.具有较好的人文社会科学素养、较强的社会责任感和良好的工程职业道德; 2.具有从事工程工作所需的相关数学、自然科学知识以及一定的经济管理知识; 3.掌握扎实的工程基础知识和本专业的基本理论知识,了解本专业的前沿发展现状和趋势; 4.具有综合运用所学科学理论和技术手段分析并解决工程问题的基本能力; 5.掌握文献检索、资料查询及运用现代信息技术获取相关信息的基本方法; 6.具有创新意识和对新产品、新工艺、新技术和新设备进行研究、开发和设计的初步能力; 7.了解与本专业相关的职业和行业的生产、设计、研究与开发的法律、法规,熟悉环境保护和可持续发展等方面的方针、政策和法津、法规,能正确认识工程对于客观世界和社会的影响; 8.具有一定的组织管理能力、较强的表达能力和人际交往能力以及在团队
数字信号处理实验五
实验五:FIR数字滤波器设计与软件实现 信息学院 10电本2班王楚炘 2010304224 10.5.1 实验指导 1.实验目的 (1)掌握用窗函数法设计FIR数字滤波器的原理和方法。 (2)掌握用等波纹最佳逼近法设计FIR数字滤波器的原理和方法。 (3)掌握FIR滤波器的快速卷积实现原理。 (4)学会调用MATLAB函数设计与实现FIR滤波器。 2.实验内容及步骤 (1)认真复习第七章中用窗函数法和等波纹最佳逼近法设计FIR数字滤波器的原理; (2)调用信号产生函数xtg产生具有加性噪声的信号xt,并自动显示xt及其频谱,如图10.5.1所示; 图10.5.1 具有加性噪声的信号x(t)及其频谱如图(3)请设计低通滤波器,从高频噪声中提取xt中的单频调幅信号,要求信号幅频失真小于0.1dB,将噪声频谱衰减60dB。先观察xt的频谱,确定滤波器指标参数。 (4)根据滤波器指标选择合适的窗函数,计算窗函数的长度N,
调用MATLAB函数fir1设计一个FIR低通滤波器。并编写程序,调用MATLAB快速卷积函数fftfilt实现对xt的滤波。绘图显示滤波器的频响特性曲线、滤波器输出信号的幅频特性图和时域波形图。 (4)重复(3),滤波器指标不变,但改用等波纹最佳逼近法,调用MATLAB函数remezord和remez设计FIR数字滤波器。并比较两种设计方法设计的滤波器阶数。 提示:MATLAB函数fir1和fftfilt的功能及其调用格式请查阅本书 第7章和第?章; 采样频率Fs=1000Hz,采样周期T=1/Fs; 根据图10.6.1(b)和实验要求,可选择滤波器指标参数:通带截止频率fp=120Hz,阻带截至频率fs=150Hz,换算成数字频率,通带截止频率,通带最大衰为0.1dB,阻带截至频率,阻带最小衰为60dB。]实验程序框图如图10.5.2所示,供读者参考。 Fs=1000,T=1/Fs xt=xtg 产生信号xt, 并显示xt及其频谱 用窗函数法或等波纹最佳逼近法 设计FIR滤波器hn 对信号xt滤波:yt=fftfilt(hn,xt) 1、计算并绘图显示滤波器损耗函数 2、绘图显示滤波器输出信号yt End 图10.5.2 实验程序框图 4.思考题 (1)如果给定通带截止频率和阻带截止频率以及阻带最小衰减,如何用窗函数法设计线性相位低通滤波器?请写出设计步骤. 答:用窗函数法设计线性相位低通滤波器的设计步骤: a.根据对阻带衰减及过渡带的指标要求,选择窗函数的类型,并估计窗口的长度N; b.构造希望逼近的频率响应函数; c.计算h d(n); d.加窗得到设计结果h(n)=h d(n)w(n)。 (2)如果要求用窗函数法设计带通滤波器,且给定通带上、下截止频率为和,阻带上、下截止频率为和,试求理想带通滤波器的截止频率。 答:希望逼近的理想带通滤波器的截止频率分别为:
数字信号处理试题及答案
数字信号处理试题及答案 一、填空题:(每空1分,共18分) 1、 数字频率ω是模拟频率Ω对采样频率s f 的归一化,其值是 连续 (连续还是离散?)。 2、 双边序列z 变换的收敛域形状为 圆环或空集 。 3、 某序列的 DFT 表达式为∑-==1 0)()(N n kn M W n x k X ,由此可以看出,该序列时域的长度为 N ,变换后数字频域上相邻两个频率样点之间的间隔是 M π 2 。 4、 线性时不变系统离散时间因果系统的系统函数为2 52) 1(8)(2 2++--=z z z z z H ,则系统的极点为 2,2 1 21-=-=z z ;系统的稳定性为 不稳定 。系统单位冲激响应)(n h 的初值 4)0(=h ;终值)(∞h 不存在 。 5、 如果序列)(n x 是一长度为64点的有限长序列)630(≤≤n ,序列)(n h 是一长度为128点 的有限长序列)1270(≤≤n ,记)()()(n h n x n y *=(线性卷积),则)(n y 为 64+128-1=191点 点的序列,如果采用基FFT 2算法以快速卷积的方式实现线性卷积,则FFT 的点数至少为 256 点。 6、 用冲激响应不变法将一模拟滤波器映射为数字滤波器时,模拟频率Ω与数字频率ω之间的 映射变换关系为T ω = Ω。用双线性变换法将一模拟滤波器映射为数字滤波器时,模拟频率Ω 与数字频率ω之间的映射变换关系为)2 tan(2ω T =Ω或)2arctan(2T Ω=ω。 7、当线性相位 FIR 数字滤波器满足偶对称条件时,其单位冲激响应)(n h 满足的条件为 )1()(n N h n h --= ,此时对应系统的频率响应)()()(ω?ω ωj j e H e H =,则其对应的相位函数 为ωω?2 1 )(-- =N 。 8、请写出三种常用低通原型模拟滤波器 巴特沃什滤波器 、 切比雪夫滤波器 、 椭圆滤波器 。 二、判断题(每题2分,共10分) 1、 模拟信号也可以与数字信号一样在计算机上进行数字信号处理,只要加一道采样的工序就可 以了。 (╳) 2、 已知某离散时间系统为)35()]([)(+==n x n x T n y ,则该系统为线性时不变系统。(╳)
北理工随机信号分析实验报告
本科实验报告实验名称:随机信号分析实验
实验一 随机序列的产生及数字特征估计 一、实验目的 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)均匀分布随机数按上式进行变
西电随机信号分析大作业
随机信号分析大作业 学院:电子工程学院 班级:021151 学号:02115037 姓名:隋伟哲
第一题:设有随机信号X(t)=5cos(t+a),其中相位a是在区间(0,2π)上均匀分布的随机变量,使用Matlab编程产生其三个样本函数。 解: 源程序如下: clc;clear; C=2*pi*rand(1,3);%在[0,2π]产生均匀分布的相位角 t=1:.1:80; y1=5*cos(t+C(1)); %将产生的随机相位角逐一代入随机过程中 y2=5*cos(t+C(2)); %将产生的随机相位角逐一代入随机过程中 y3=5*cos(t+C(3)); %将产生的随机相位角逐一代入随机过程中 plot(t,y1,'r-'); hold on; plot(t,y2,'g--'); hold on; plot(t,y3,'k-'); xlabel('t');ylabel('X(t)'); grid on;axis([0 30 -8 8]); title('随机相位的三条样本曲线'); 产生的三条样本曲线:
第二题:利用Matlab程序设计一正弦型信号加高斯白噪声的复合信号。(1)分析复合信号的功率谱密度、幅度分布特性; (2)分析复合信号通过RC积分电路后的功率谱密度和相应的幅度分布特性; (3)分析复合信号通过理想低通系统后的功率谱密度和相应的幅度分布特性。 解:设定正选信号的频率为10HZ,抽样频率为100HZ x=sin(2*pi*fc*t)
(1)正弦函数加上高斯白噪声: y=awgn(x,10) y 的幅度分布特性可以通过傅里叶变换得到: Y(jw)=fft(y) y 的功率谱密度: G(w)=Y(jw).*conj(Y(jw)/length(Y(jw))) 随机序列自相关函数的无偏估计公式为: 1 01()()()N m xx n R m x n x n m N m --==+-∑ 01m N ≤≤- (2)复合信号 y 通过RC 积分电路后得到信号y2 通过卷积计算可以得到y2 即:y2= conv2(y,b*pi^-b*t) y2的幅度分布特性可以通过傅里叶变换得到: Y2(jw)=fft(y2) y2的功率谱密度: G2(w)=Y2(jw).*conj(Y2(jw)/length(Y2(jw))) (3)复合信号 y 通过理想滤波器电路后得到信号y3 通过卷积计算可以得到y3 即:y3=conv2(y,sin(10*t)/(pi*t)) y3的幅度分布特性可以通过傅里叶变换得到: Y3(jw)=fft(y3) y3的功率谱密度: G3(w)=Y3(jw).*conj(Y3(jw)/length(Y3(jw)))
测试信号分析与处理作业实验一二
王锋 实验一:利用FFT 作快速相关估计 一、实验目的 a.掌握信号处理的一般方法,了解相关估计在信号分析与处理中的作用。 b.熟悉FFT算法程序;熟练掌握用FFT作快速相关估计的算法。 c.了解快速相关估计的谱分布的情况。 二、实验内容 a.读入实验数据[1]。 b.编写一利用FFT作相关估计的程序[2]。 c.将计算结果表示成图形的形式,给出相关谱的分布情况图。 注[1]:实验数据文件名为“Qjt.dat”。 实验数据来源:三峡前期工程 “覃家沱大桥” 实测桥梁振动数据。 实验数据采样频率:50Hz。 可从数据文件中任意截取几段数据进行分析,数据长度N 自定。 注[2]:采用Matlab 编程。 三、算法讨论及分析 算法为有偏估计,利用FFT计算相关函数 Step 1: 对原序列补N个零,得新序列x2N(n) Step2: 作FFT[x2N(n)]得到X2N(k) Step 3: 取X2N(k)的共轭,得 Step 4: 作 Step 5: 调整与的错位。 四、实验结果分析 1. 该信号可以近似为平稳信号么? 可以近似为平稳信号,随机过程的统计特性不随样本的采样时刻而发生变化。取N=8192,分别取间隔m=500,m=700,m=1000,所得到的均值均为0.5366,方差为47369,与时间无关。
图1-1 自相关函数图 (上图表示的R0,下图为调整后的R0) 2. 该信号是否具有周期性,信噪比如何? >> load Qjt.dat; %加载数据 N=32768; %数据长度 i=1:1:N; %提取数据 plot(i,Qjt(i)); 抛去几个极值点,从图1-2可以看出,数据具有一定的周期性,杂音比较少,说明信噪比较高。 图1-2 数据图
随机信号分析大作业
随机信号分析实验报告 信息25班 2120502123 赵梦然
作业题三: 利用Matlab 产生一个具有零均值、单位方差的的高斯白噪声随机序列X(n),并通过一脉冲响应为 (0.8)(0)0 n n h n else =≥??? 的线性滤波器。 (1) 产生一个具有零均值、单位方差的的高斯白噪声随机序列X(n),检验其一维概率密度函 数是否与理论相符。 (2) 绘出输入输出信号的均值、方差、自相关函数及功率谱密度的图形,讨论输出信号服从 何种分布。 (3) 试产生在[-1,+1]区间均匀分布的白噪声序列,并将其替换高斯白噪声通过上述系统。 画出此时的输出图形,并观察讨论输出信号服从何种分布。 作业要求 (1) 用MATLAB 编写程序。最终报告中附代码及实验结果截图。 (2) 实验报告中必须有对实验结果的分析讨论。 提示: (1) 可直接使用matlab 中已有函数产生高斯白噪声随机序列。可使用hist 函数画出序列的 直方图,并与标准高斯分布的概率密度函数做对比。 (2) 为便于卷积操作,当N 很大时,可近似认为h(N)=0。卷积使用matlab 自带的conv 函 数。 (3) 分析均值、方差等时,均可使用matlab 现有函数。功率谱密度和自相关函数可通过傅 里叶变换相互获得。傅里叶变换使用matlab 自带的fft 函数。 (4) 作图使用plot 函数。
一、作业分析: 本题主要考察的是加性高斯白噪声相关问题,因此构造一个高斯白噪声十分重要,故在本题中使用randn函数随机生成一个个符合高斯分布的数据,并由此构成高斯白噪声;而且由于白噪声是无法完全表示的,故此根据噪声长度远大于信号长度时可视为高斯白噪声,构造了一个长度为2000的高斯白噪声来进行试验。 二、作业解答: (1)matlab程序为: x-1000:1:1000; k=1*randn(1,length(x));% 生成零均值单位方差的高斯白噪声。 [f,xi]=ksdensity(x);%利用ksdensity函数估计样本的概率密度。 subplot(1,2,1); plot(x,k); subplot(1,2,2); plot(xi,f); 实验结果为:
测试信号分析与处理作业实验五
王锋 实验五:多种功率谱估计的比较 一、实验目的 a.了解功率谱估计在信号分析中的作用; b.掌握随机信号分析的基础理论,掌握参数模型描述形式下的随机信 号的功率谱的计算方法; c.掌握在计算机上产生随机信号的方法; d.了解不同的功率谱估计方法的优缺点。 二、实验准备 有三个信号源,分别代表三种随机信号(序列)。 信号源1: 123()2cos(2)2cos(2)2cos(2)()x n f n f n f n z n πππ=+++ 其中,1230.08,=0.38,0.40f f f == z(n)是一个一阶 AR 过程,满足方程: ()(1)(1)()z n a z n e n =--+ (1)0.823321a =- e(n)是一高斯分布的实白噪声序列,方差20.1σ= 信号源2和信号源3: 都是4阶的AR 过程,它们分别是一个宽带和一个窄带过程,满足方程: ()(1)(1)(2)(2)(3)(3)(4)(4)()x n a x n a x n a x n a x n e n =--------+ e(n)是一高斯分布的实白噪声序列,方差2σ,参数如下: 三、实验内容 a. 描绘出这三个实验信号的真实功率谱波形。 b. 在计算机上分别产生这个三个信号,令所得到的数据长度 N= 256 。 注意:产生信号的时候注意避开起始瞬态点。例如,可以产生长度为512 的信号序列,然后取后面256 个点作为实验数据。 c. 分别用如下的谱估计方法,对三个信号序列进行谱估计。 1、经典谱估计 周期图法 自相关法 平均周期图法(Bartlett 法)
Welch法(可选每段64 点,重叠32 点,用Hamming 窗)2、现代谱估计 Yule - Walker方程(自相关法) 最小二乘法 注:阶次p可在3-20之间,由自己给定。 四、实验结果分析 生成的信号源
信号分析与处理答案第二版完整版
信号分析与处理答案第 二版 HEN system office room 【HEN16H-HENS2AHENS8Q8-HENH1688】
第二章习题参考解答 求下列系统的阶跃响应和冲激响应。 (1) 解当激励为时,响应为,即: 由于方程简单,可利用迭代法求解: ,, …, 由此可归纳出的表达式: 利用阶跃响应和冲激响应的关系,可以求得阶跃响应: (2) 解 (a)求冲激响应 ,当时,。 特征方程,解得特征根为。所以: …(2.1.2.1) 通过原方程迭代知,,,代入式(2.1.2.1)中得:解得,代入式(2.1.2.1): …(2.1.2.2) 可验证满足式(2.1.2.2),所以: (b)求阶跃响应 通解为 特解形式为,,代入原方程有,即 完全解为 通过原方程迭代之,,由此可得 解得,。所以阶跃响应为: (3)
解 (4) 解 当t>0时,原方程变为:。 …(2.1.3.1) …(2.1.3.2) 将(2.1.3.1)、式代入原方程,比较两边的系数得: 阶跃响应: 求下列离散序列的卷积和。 (1) 解用表 格法求 解 (2) 解用表 格法求 解 (3) 和 如题图2.2.3所示 解用表 格法求 解
(4) 解 (5) 解 (6) 解参见右图。 当时: 当时: 当时: 当时: 当时: (7) , 解参见右图: 当时: 当时: 当时: 当时: 当时: (8) ,解参见右图
当时: 当时: 当时: 当时: (9) , 解 (10) , 解 或写作:
求下列连续信号的卷积。 (1) , 解参见右图: 当时: 当时: 当时: 当时: 当时: 当时: (2) 和如图2.3.2所示 解当时: 当时: 当时: 当时: 当时: (3) , 解 (4) , 解 (5) , 解参见右图。当时:当时: 当时:
数字信号处理试卷及答案
A 一、 选择题(每题3分,共5题) 1、)6 3()(π-=n j e n x ,该序列是 。 A.非周期序列 B.周期6 π = N C.周期π6=N D. 周期π2=N 2、序列)1()(---=n u a n x n ,则)(Z X 的收敛域为 。 A.a Z < B.a Z ≤ C.a Z > D.a Z ≥ 3、对)70()(≤≤n n x 和)190()(≤≤n n y 分别作 20 点 DFT ,得)(k X 和)(k Y , 19,1,0),()()( =?=k k Y k X k F ,19,1,0)],([)( ==n k F IDFT n f , n 在 围时,)(n f 是)(n x 和)(n y 的线性卷积。 A.70≤≤n B.197≤≤n C.1912≤≤n D.190≤≤n 4、)()(101n R n x =,)()(72n R n x =,用DFT 计算二者的线性卷积,为使计算量尽可能的少,应使DFT 的长度N 满足 。 A.16>N B.16=N C.16实验一-LabVIEW中的信号分析与处理
实验一 LabVIEW中的信号分析与处理 一、实验目的: 1、熟悉各类频谱分析VI的操作方法; 2、熟悉数字滤波器的使用方法; 3、熟悉谐波失真分析VI的使用方法。 二、实验原理: 1、信号的频谱分析是指用独立的频率分量来表示信号;将时域信号变换到频域,以显示在时域无法观察到的信号特征,主要是信号的频率成分以及各频率成分幅值和相位的大小,LabVIEW中的信号都是数字信号,对其进行频谱分析主要使用快速傅立叶变换(FFT)算法:·“FFT Spectrum(Mag-Phase).vi”主要用于分析波形信号的幅频特性和相频特性,其输出为单边幅频图和相频图。 ·“FFT.vi”以一维数组的形式返回时间信号的快速傅里叶运算结果,其输出为双边频谱图,在使用时注意设置FFT Size为2的幂。 ·“Amplitude and Phase Spectrum .vi”也输出单边频谱,主要用于对一维数组进行频谱分析,需要注意的是,需要设置其dt(输入信号的采样周期)端口的数据。 2、数字滤波器的作用是对信号进行滤波,只允许特定频率成份的信号通过。滤波器的主要类型分为低通、高通、带通、带阻等,在使用LabVIEW中的数字滤波器时,需要正确设置滤波器的截止频率(注意区分模拟频率和数字频率)和阶数。 3、“Harmonic Distortion Analyzer .vi”用于分析输入的波形数据的谐波失真度(THD),该vi还可分析出被测波形的基波频率和各阶次谐波的电平值。 三、实验容: (1) 时域信号的频谱分析 设计一个VI,使用4个Sine Waveform.vi(正弦波形)生成频率分别为10Hz、30Hz、50Hz、100Hz,幅值分别为1V、2V、3V、4V的4个正弦信号(采样频率都设置为1kHz,采样点数都设置为1000点),将这4个正弦信号相加并观察其时域波形,然后使用FFT Spectrum(Mag-Phase).vi对这4个正弦信号相加得出的信号进行FFT频谱分析,观察其幅频和相频图,并截图保存。
(完整版)数字信号处理试卷及答案
江 苏 大 学 试 题 课程名称 数字信号处理 开课学院 使用班级 考试日期
江苏大学试题第2A页
江苏大学试题第3A 页
江苏大学试题第页
一、填空题:(每空1分,共18分) 8、 数字频率ω是模拟频率Ω对采样频率s f 的归一化,其值是 连续 (连续还是离散?)。 9、 双边序列z 变换的收敛域形状为 圆环或空集 。 10、 某序列的DFT 表达式为∑-== 10 )()(N n kn M W n x k X ,由此可以看出,该序列时域的长度为 N , 变换后数字频域上相邻两个频率样点之间的间隔是 M π 2 。 11、 线性时不变系统离散时间因果系统的系统函数为2 52) 1(8)(22++--=z z z z z H ,则系统的极点为 2,2 1 21-=-=z z ;系统的稳定性为 不稳定 。系统单位冲激响应)(n h 的初值4)0(=h ; 终值)(∞h 不存在 。 12、 如果序列)(n x 是一长度为64点的有限长序列)630(≤≤n ,序列)(n h 是一长度为128点的有限长 序列)1270(≤≤n ,记)()()(n h n x n y *=(线性卷积),则)(n y 为 64+128-1=191点 点的序列,如果采用基FFT 2算法以快速卷积的方式实现线性卷积,则FFT 的点数至少为 256 点。 13、 用冲激响应不变法将一模拟滤波器映射为数字滤波器时,模拟频率Ω与数字频率ω之间的映射变换 关系为T ω = Ω。用双线性变换法将一模拟滤波器映射为数字滤波器时,模拟频率Ω与数字频率ω之 间的映射变换关系为)2tan(2ωT = Ω或)2 arctan(2T Ω=ω。 当线性相位FIR 数字滤波器满足偶对称条件时,其单位冲激响应)(n h 满足的条件为)1()(n N h n h --= ,
随机信号实验报告
随机信号分析 实验报告 目录 随机信号分析 (1) 实验报告 (1) 理想白噪声和带限白噪声的产生与测试 (2) 一、摘要 (2) 二、实验的背景与目的 (2) 背景: (2) 实验目的: (2) 三、实验原理 (3) 四、实验的设计与结果 (4) 实验设计: (4) 实验结果: (5) 五、实验结论 (12) 六、参考文献 (13) 七、附件 (13) 1
理想白噪声和带限白噪声的产生与测试一、摘要 本文通过利用MATLAB软件仿真来对理想白噪声和带限白噪声进行研究。理想白噪声通过低通滤波器和带通滤波器分别得到低通带限白噪声和帯通带限白噪声。在仿真的过程中我们利用MATLAB工具箱中自带的一些函数来对理想白噪声和带限白噪声的均值、均方值、方差、功率谱密度、自相关函数、频谱以及概率密度进行研究,对对它们进行比较分析并讨论其物理意义。 关键词:理想白噪声带限白噪声均值均方值方差功率谱密度自相关函数、频谱以及概率密度 二、实验的背景与目的 背景: 在词典中噪声有两种定义:定义1:干扰人们休息、学习和工作的声音,引起人的心理和生理变化。定义2:不同频率、不同强度无规则地组合在一起的声音。如电噪声、机械噪声,可引伸为任何不希望有的干扰。第一种定义是人们在日常生活中可以感知的,从感性上很容易理解。而第二种定义则相对抽象一些,大部分应用于机械工程当中。在这一学期的好几门课程中我们都从不同的方面接触到噪声,如何的利用噪声,把噪声的危害减到最小是一个很热门的话题。为了加深对噪声的认识与了解,为后面的学习与工作做准备,我们对噪声进行了一些研究与测试。 实验目的: 了解理想白噪声和带限白噪声的基本概念并能够区分它们,掌握用MATLAB 或c/c++软件仿真和分析理想白噪声和带限白噪声的方法,掌握理想白噪声和带限白噪声的性质。
信号分析与处理实验
信号分析与处理实验实验一时间信号的产生 班级:自动化1101班 姓名:陈宝平 学号: 成绩:
1. 实验目的 数字信号处理系统中的信号都是以离散时间形态存在的。研究离散时间信号,首先 需要产生出各种离散时间信号。使用MATLAB 软件可以方便的产生各种常见的离散时间信号,还具有强大的绘图功能,便于用户直观地处理输出结果。 通过本实验,学习用MATLAB 产生一些常见的离散时间信号,并通过MATLAB 绘图工具对产生的信号进行观察加深对常见离散信号和信号卷积和运算的理解。 2. 实验原理 离散时间信号用x(n)来表示,自变量n 必须是整数;连续时间信号用x(t)来表示。常见的时间信号如下: (1) 单位冲激序列δ(n)=???≠=0,00,1n n ; 如图(a); 单位冲激信号δ(t)=? ??≠=0,00 ,1t t ;如图(b): (a) (b) 如果δ(n)在时间轴上延迟了k 个单位,得到δ(n-k)=? ??≠=k n k n ,0,1。 (2) 单位阶跃序列u(n)=???<≥0,00,1n n ; 如图(c):单位阶跃信号u(t)=???<≥0 ,00 ,1t t ; 如图(d): (c) (d) 如果u(n)在时间轴上延迟了k 个单位,得到u(n-k)=???<≥k n k n ,0,1。
(3) 矩形序列R N (n)=???≥<-≤≤),0(,0) 10(,1N n n N n ,矩形序列有一个重要的参数,就是序列宽度 N 。R N (n)与u(n)之间的关系为R N (n)= u(n)- u(n-N)。如图(e): 单位矩形信号R (t)=? ??≥<≤≤),0(,0) 0(,1T t t T t , R(t)=u(t)-u(t-T),如图(f): (e) (f) (4)正弦序列x(n)=Acos(ω0n+?)。只有当 2ωπ 为有理数时,正弦序列具有周期性, 如图(g): 正弦信号x(t)=Acos(?ω+t 0),如图(h); (g) (h) (5)单边实指数序列x(n)=a n u(n),当a>0时,该序列均取正值,当a<0时,序列在正负摆动。如图分别为:x(n)=1.2n , x(n)=(-1.2)n , x(n)=0.8n , x(n)=(-0.8)n 的图。
数字信号处理试卷大全..
北京信息科技大学 2010 ~2011 学年第一学期 《数字信号处理》课程期末考试试卷(A) 一、填空题(本题满分30分,共含4道小题,每空2分) 1.两个有限长序列x1(n),0≤n≤33和x2(n),0≤n≤36,做线性卷积 后结果的长度是,若对这两个序列做64点圆周卷积,则圆周卷积结果中n= 至为线性卷积结果。 W的、和三个固有特性来实现2.DFT是利用nk N FFT快速运算的。 3.IIR数字滤波器设计指标一般由、、和等 四项组成。 4.FIR数字滤波器有和两种设计方法,其结构 有、和等多种结构。 二、判断题(本题满分16分,共含8道小题,每小题2分,正 确打√,错误打×) 1.相同的Z变换表达式一定对应相同的时间序列。() 2.Chirp-Z变换的频率采样点数M可以不等于时域采样点数N。() 3.按频率抽取基2 FFT首先将序列x(n)分成奇数序列和偶数序列。() 4.冲激响应不变法不适于设计数字带阻滤波器。() 5.双线性变换法的模拟角频率Ω与数字角频率ω成线性关系。() 6.巴特沃思滤波器的幅度特性必在一个频带中(通带或阻带)具有等
波纹特性。( ) 7. 只有FIR 滤波器才能做到线性相位,对于IIR 滤波器做不到线性相 位。( ) 8. 在只要求相同的幅频特性时,用IIR 滤波器实现其阶数一定低于 FIR 阶数。( ) 三、 综合题(本题满分18分,每小问6分) 若x (n)= {3,2,1,2,1,2 },0≤n≤5, 1) 求序列x(n)的6点DFT ,X (k)=? 2) 若)()]([)(26k X W n g DFT k G k ==,试确定6点序列g(n)=? 3) 若y(n) =x(n)⑨x(n),求y(n)=? 四、 IIR 滤波器设计(本题满分20分,每小问5分) 设计一个数字低通滤波器,要求3dB 的截止频率f c =1/π Hz ,抽样频率f s =2 Hz 。 1. 导出归一化的二阶巴特沃思低通滤波器的系统函数H an (s)。 2. 试用上述指标设计一个二阶巴特沃思模拟低通滤波器,求其系 统函数H a (s),并画出其零极点图。 3. 用双线性变换法将H a (s)转换为数字系统的系统函数H(z)。 4. 画出此数字滤波器的典范型结构流图。 五、 FIR 滤波器设计(本题满分16分,每小问4分)