经典谱估计算法研究与实现

毕业设计（论文）

论文题目：经典谱估计算法研究与实现

教学中心：电子科技大学网络教育学院苏州学习中心指导老师：职称：

学生姓名：学号：

专业：通信工程

毕业设计（论文）任务书

题目：经典谱估计算法研究与实现

任务与要求：

探讨我国经典谱估计算法运用和影响下所面临的机遇与挑战以及经典谱估计算法研究与实现,结合所知识,理论联系实际,

写出毕业论文。

时间：2014年 1 月25日至 2014年 4 月 14日共 12 周教学中心：电子科技大学网络教育学院苏州学习中心

学生姓名：学号：

专业：通信工程

指导单位或教研室：电子科技大学网络教育学院苏州学习中心

指导教师：职称：

毕业设计(论文)进度计划表

电子科技大学毕业设计(论文)中期检查记录表

注:此表同学生毕业设计(论文)一起存档

1 随机信号的经典谱估计方法

估计功率谱密度的平滑周期图是一种计算简单的经典方法。它的主要特点是

与任何模型参数无关，是一类非参数化方法[4]。它的主要问题是：由于假定信号的自相关函数在数据观测区以外等于零，因此估计出来的功率谱很难与信号的真实功率谱相匹配。在一般情况下，周期图的渐进性能无法给出实际功率谱的一个满意的近似，因而是一种低分辨率的谱估计方法。本章主要介绍了周期图法、相关法谱估计（BT ）、巴特利特(Bartlett)平均周期图的方法和Welch 法这四种方法。

2.1 周期图法

周期图法又称直接法。它是从随机信号x(n)中截取N 长的一段，把它视为能量有限x(n)真实功率谱)(jw x e S 的估计)(jw x e S 的抽样.

周期图这一概念早在1899年就提出了，但由于点数Ｎ一般比较大，该方法的计算量过大而在当时无法使用。只是1965年FFT 出现后，此法才变成谱估计的一个常用方法。周期图法[5]包含了下列两条假设：

1．认为随机序列是广义平稳且各态遍历的，可以用其一个样本x(n)中的一段)(n x N 来估计该随机序列的功率谱。这当然必然带来误差。

2．由于对)(n x N 采用DFT ，就默认)(n x N 在时域是周期的，以及)(k x N 在频域是周期的。这种方法把随机序列样本x(n)看成是截得一段)(n x N 的周期延拓，这也就是周期图法这个名字的来历。与相关法相比，相关法在求相关函数)(m R x 时将)(n x N 以外是数据全都看成零，因此相关法认为除)(n x N 外x(n)是全零序列，这种处理方法显然与周期图法不一样。

但是，当相关法被引入基于FFT 的快速相关后，相关法和周期图法开始融合。通过比较我们发现：如果相关法中M=N ，不加延迟窗，那么就和补充（N-1）个零的周期图法一样了。简单地可以这样说：周期图法是M=N 时相关法的特例。因此相关法和周期图法可结合使用。

2.2 相关法谱估计（BT ）法

这种方法以相关函数为媒介来计算功率谱，所以又叫间接法。它是1958年由Blackman 和Tukey 提出。这种方法的具体步骤是：

第一步：从无限长随机序列x(n)中截取长度N 的有限长序列列)(n x N 第二步：由N 长序列)(n x N 求（2M-1）点的自相关函数)(m R x ∧

序列。即

)()(1)(1

m n x n x N m R N n N N x +=∑-=∧

(2-1)

这里，m=-(M-1)…,-1,0,1…,M-1，

M N ，)(m R x 是双边序列，但是由自相关函数的偶对称性式，只要求出m=0,。。。,M-1的傅里叶变换，另一半也就知道了。

第三步：由相关函数的傅式变换求功率谱。即

jwm M M m X

jw

x e m R

e S ----=∧∧

∑=

)()(1)

1(

(2-2)

以上过程中经历了两次截断，一次是将x(n)截成N 长，称为加数据窗，一次是将x(n)截成（2M-1）长，称为加延迟窗。因此所得的功率谱仅是近似值，也叫谱估计，式中的)(jw x e S 代表估值。一般取M<

当FFT 问世后，情况有所变化。因为截断后的)(n x N 可视作能量信号，由相关卷积定理可得

)(m R x ∧

)]()([1

m x m x N

N N -*=

(2-3)

这就将相关化为线性卷积，而线性卷积又可以用快速卷积来实现。我们可对上式两边取（2N-1）点DFT ，则有

2

121212)(1)()([1)(K X N

K x k x N m R N N N x ---∧

=*=

(2-4)

于是将时域卷积变为频域乘积，用快速相关求)(m R x ∧

的完整方案如下： 1.

对N 长)(n x N 的补充（N-1）个零，成为（2N-1）长的。

2.

求（2N-1）点的FFT ，得∑-=----=2

20

121

212)()(N N mk N N N W n x

K X 。

3.

求

2

12)(1K X N

N -。由DFT 性质，)(12n x N -是纯实的，)(12k x N -满足共轭偶对称，而2

12)(1K X N

N -一定是实偶的，且以（2N-1）为周期。 4. 求（2N-1）点的IFFT ：

mk N N N k N x W K X N N m R -----=-∧

∑-=121

)1(2

12)(1121)(

(2-5)

这里

2

12)(1K X N

N -是实偶的，m=-(N-1)...0...N-1。本来IFFT 求和范围是0至2N-2，由于2

12)(1K X N

N -的实偶性与周期性，求和范围改为-（N-1）至（N-1）

不影响计算结果。同理可将m 的范围改为-（N-1）至（N-1）。

上述的快速相关中，补充零的目的是为了能用圆周卷积代替线性卷积，以便进一步采用快速卷积算法。快速相关输出是-（N-1）至（N-1）的2N-1点，加)(m W M 窗后截取的是-（M-1）至（M-1）的频段，最后作（2M-1）点FFT ，得)(k S x ∧

。

我们注意到：如果数据点数与自相关序列点数相同即M=N ，则（2N-1）点的IFFT 后紧跟一个（2N-1）点的FFT ，利用)(m R x ∧

的对称性，FT 运算框的计算式变为

=

)(K S X ∑---=∧

1)

1()(N N m X

m R

1212---N mk N X W

(2-6)

由于N=M 并假设窗形状是矩形的，第二次()m W M 的截断就不需要了。比较

式（2-5）和式（2-6）， ](m)R FFT[k S x x ∧

=）

（,])(1[)( 2

12K X N

IFFT m R N x -∧

=正反傅氏变换可以抵消，直接得

）

（k S x =2

12)(1K X N

N - (2-7)

为了实行基2FFT ，也可将（2N-1）点换成2N 点，这样做不影响结果的正确性。

2.3 巴特利特(Bartlett)平均周期图法

首先让我们来看一下为什么周期图经过某种平均(或平滑)后会使它的方差当∞→N 时趋于零，达到一致估计的目的。如果L x x x , , ,21 是不相关的随机变量，每一个具有期望值μ，方差2σ，则可以证明它们的数学平均

L x x x x l /)...(21+++=的期望值等于μ，数学平均的方差等于L /2σ，即：

[][]μμ=?=+++=

L L

x x x E L x E L 1

121 [][]

[]222)(][))((x E x E x E x E x Var -=-=

[]

22212)(1

μ-+++=

L x x x E L

()[]

211

2

222121μ-??

???

??

???++++=∑∑=≠=L j L j i i j

i L x x E x x x E L [][][]∑∑∑∑=≠==≠=?=L j L

j

i i i

j

L j j

i i j

i

x E x E x x E 1

111

222)1(μμμμL L L L -=-=

所以

[][][]

??

????-=-??????-+=∑∑==21222

22122211μμμμL x E L L L x E L x Var L i i L i i

[][][][][][]{}

22

222221212])[(])[(])[(1L L x E x E x E x E x E x E L

-++-+-=

L L L

2221σσ=

=

(2-8)

由(2-8)可见，L 个平均的方差比每个随机变量的单独方差小L 倍。当

[]0 →∞→x Var L ，可达到一致谱估计的目的。因而降低估计量的方差的一种有

效方法是将若干个独立估计值进行平均。把这种方法应用于谱估计通常归功于Bartlett 。

Bartlett 平均周期图的方法是将序列)10( )(-≤≤N n n x 分段求周期图再平均。设将)(n x 分成L 段，每段有M 个样本，因而LM N =，第i 段样本序列可写成

L i M n M iM n x n x i ≤≤-≤≤-+=1 ,10 )()(

第i 段的周期图为

2

10

)(1)(∑-=-=

M n n j j

i M e n x

M

I ωω

如果)( ,m M m xx φ>很小，则可假定各段的周期图)(ωi

M I 是互相独立的。对

功率谱密度的概念的讨论，谱估计可定义为L 段周期图的平均，即

∑==L i i

M xx I L P 1

)(1)(?ωω

(2-9)

于是它的期望值为

[]

[][]

∑===L i i M i M xx I E I E L P E 1

)()(1)(?ωωω []

?-

-=π

πθωθθπ

ωd e W P

P E j B xx

xx

)()(21)(?)(

()[()[?-??

????--=

π

πθθωθωθπd M P M

xx 2

2/sin 2sin )(21 (2-10) 这里L N M /=，因此Bartlett 估计的期望值是真实谱)(ωxx P 与三角窗函数的卷积。由于三角窗函数不等于δ函数，所以Bartlett 估计也是有偏估计即0≠Bias ，但当∞→N 时，0→Bias 。

由于我们假定各段周期图是相互独立的，所以可按式(2-8)得到下式：

[]

[])(1)(?ωωM

xx I Var L

P Var = ???

?

??????? ??+≈2

2s i n )s i n (1)(1ωωωM M P L xx

(2-11)

由此可见，随着L 的增加[]

)(?ωxx P Var 是下降的，当∞→L 时，[]

0)(?→ωxx

P Var 。因此Bartlett 估计是一致估计。

比较式(2-10)的[]

)(?ωxx

P E 与式(2-1)的[])(ωN I E 可见在二种情况的估计量的期望值都是真值)(ωxx P 与窗口函数)(ωj B e W 的卷积形式，但后者将前者WB 中N 改为M ，

N L N M <=/。因而使)(ωj B e W 主瓣的宽度增大。由于主瓣的宽度愈窄愈接近δ函

经典功率谱和Burg法的功率谱估计

现代信号处理作业 实验题目： 设信号)()8.0cos(25.0)47.0cos()35.0cos()(321n v n n n n x ++++++=θπθπθπ，其中321,,θθθ是[]ππ,-内的独立随机变量，v(n)是单位高斯白噪声。 1.利用周期图法对序列进行功率谱估计。数据窗采用汉明窗。 2.利用BT 法对序列进行功率谱估计，自相关函数的最大相关长度为M=64,128,256,512采用BARTLETT 窗。 3.利用Welch 法对序列进行功率谱估计，50%重叠，采用汉明窗，L=256,128,64。 4.利用Burg 法对序列进行AR 模型功率谱估计，阶数分别为10，13. 要求每个实验都取1024个点，fft 作为谱估计，取50个样本序列的算术平均，画出平均的功率谱图。 实验原理： 1）。周期图法： 又称间接法，它把随机信号的N 个观察值x N (n)直接进行傅里叶变换，得到X N (e jw ),然后取其幅值的平方，再除以N ，作为对x （n ）真实功率谱的估计。 2^ )(1)(jw e X N w P N per = ， 其中∑-=-=1 )()(N n jwn N jw N e n x e X 2）。BT 法： 对于N 个观察值x(0),x(1),。。。,x(N-1),令x N (n)=a(n)x(n)。计算r x （m ）为

∑--=-≤+= m N n N N x N m m n x n x N m r 10 1),()(1 )(,计算其傅里叶变换 ∑-=--≤= M M m jwm x BT N M e m r m v w P 1 ,)()()(^ ^ ，作为观察值的功率谱的估计。 其中v(m)是平滑窗。 3)。Welch 法： 假定观察数据是x(n),n=0,1,2...,N-1，现将其分段，每段长度为M,段与段之间的重叠为M-K,第i 个数据段经加窗后可表示为 1,...,1,0 )()()(-=+=M i iK n x n a n x i M 其中K 为一整数，L 为分段数，该数据段的周期图为 2)(1)(^w X MU w P i M i per =，其中∑-=-=1 0)()(M n j w n i M i M e n x w X 。由此得到平均周期图为 ∑-==10 ^_ )(1)(L i i per w P L w P 。其中归一化U 取∑-== 10 2 )(1M n n a M U 。 4）。Burg 法： 在约束条件下，使得)(2 1^^^ b f ρρρ+=极小化，其中，约束条件是它所得到的 各阶模型解要求满足Levison 递归关系。 仿真结果： 1.周期图法

功率谱估计方法的比较

功率谱估计方法的比较 摘要： 本文归纳了信号处理中关键的一种分析方法, 即谱估计方法。概述了频谱估计中的周期图法、修正的协方差法和伯格递推法的原理，并且对此三种方法通过仿真做出了对比。 关键词：功率谱估计；AR 模型；参数 引言： 谱估计是指用已观测到的一定数量的样本数据估计一个平稳随机信号的谱。由于谱中包含了信号的很多频率信息，所以分析谱、对谱进行估计是信号处理的重要容。谱估计技术发展 渊源很长，它的应用领域十分广泛，遍及雷达、声纳、通信、地质勘探、天文、生物医学工程等众多领域，其容、方法都在不断更新，是一个具有强大生命力的研究领域。谱估计的理论和方法是伴随着随机信号统计量及其谱的发展而发展起来的，最早的谱估计方法是建 立在基于二阶统计量, 即自相关函数的功率谱估计的方法上。功率谱估计的方法经历了经典谱估计法和现代谱估计法两个研究历程，在过去及现在相当长一段时间里，功率谱估计一直占据着谱估计理论里的核心位置。经典谱估计也成为线性谱估计，包括BT 法、周期图法。现代谱估计法也称为非线性普估计，包括自相关法、修正的协方差法、伯格（Burg ）递推法、特征分解法等等。 原理： 经典谱估计方法计算简单，其主要特点是谱估计与任何模型参数无关，是一类非参数化的方法。它的主要问题是：由于假定信号的自相关函数在数据的观测区间以外等于零，因此估计出来的功率谱很难与信号的真实功率谱相匹配。在一般情况下，经典法的渐进性能无法给出实际功率谱的一个满意的近似，因而是一种低分辨率的谱估计方法。现代谱估计方法使用参数化的模型，他们统称为参数化功率谱估计，由于这类方法能够给出比经典法高得多的频率分辨率，故又称为高分辨率方法。下面分别介绍周期图法、修正的协方差法和伯格递推法。修正的协方差法和伯格递推法采用的模型均为AR 模型。 （1）周期图法 周期图法是先估计自相关函数, 然后进行傅里叶变换得到功率谱。假设随机信号x(n)只观测到一段样本数据，n=0, 1, 2, …, N-1。根据这一段样本数据估计自相关函数，如公式（1） 对（1）式进行傅里叶变换得到（2）式。 ∑--=+=1||0 *) ()(1 )(?m N n xx m n x n x N m r

材料现代分析方法试题及答案1

一、单项选择题（每题 2 分，共10 分） 3．表面形貌分析的手段包括【 d 】 （a）X 射线衍射（XRD）和扫描电镜（SEM）(b) SEM 和透射电镜（TEM） (c) 波谱仪（WDS）和X 射线光电子谱仪（XPS）(d) 扫描隧道显微镜（STM）和 SEM 4．透射电镜的两种主要功能：【b 】 （a）表面形貌和晶体结构（b）内部组织和晶体结构 （c）表面形貌和成分价键（d）内部组织和成分价键 二、判断题（正确的打√，错误的打×，每题2 分，共10 分） 1．透射电镜图像的衬度与样品成分无关。（×）2．扫描电镜的二次电子像的分辨率比背散射电子像更高。（√）3．透镜的数值孔径与折射率有关。（√）4．放大倍数是判断显微镜性能的根本指标。（×）5．在样品台转动的工作模式下，X射线衍射仪探头转动的角速度是样品转动角 速度的二倍。（√） 三、简答题（每题5 分，共25 分） 1. 扫描电镜的分辨率和哪些因素有关？为什么? 和所用的信号种类和束斑尺寸有关，因为不同信号的扩展效应不同，例如二次电子产生的区域比背散射电子小。束斑尺寸越小，产生信号的区域也小，分辨率就高。 1．透射电镜中如何获得明场像、暗场像和中心暗场像？ 答：如果让透射束进入物镜光阑，而将衍射束挡掉，在成像模式下，就得到明场象。如果把物镜光阑孔套住一个衍射斑，而把透射束挡掉，就得到暗场像，将入射束倾斜，让某一衍射束与透射电镜的中心轴平行，且通过物镜光阑就得到中心暗场像。 2．简述能谱仪和波谱仪的工作原理。 答：能量色散谱仪主要由Si(Li)半导体探测器、在电子束照射下，样品发射所含元素的荧光标识X 射线，这些X 射线被Si(Li)半导体探测器吸收，进入探测器中被吸收的每一个X 射线光子都使硅电离成许多电子—空穴对，构成一个电流脉冲，经放大器转换成电压脉冲，脉冲高度与被吸收的光子能量成正比。最后得到以能量为横坐标、强度为纵坐标的X 射线能量色散谱。 在波谱仪中，在电子束照射下，样品发出所含元素的特征x 射线。若在样品上方水平放置一块具有适当晶面间距 d 的晶体，入射X 射线的波长、入射角和晶面间距三者符合布拉格方程时，这个特征波长的X 射线就会发生强烈衍射。波谱仪利用晶体衍射把不同波长的X 射线分开，即不同波长的X 射线将在各自满足布拉格方程的2θ方向上被检测器接收，最后得到以波长为横坐标、强度为纵坐标的X射线能量色散谱。 3．电子束与试样物质作用产生那些信号？说明其用途。 （1）二次电子。当入射电子和样品中原子的价电子发生非弹性散射作用时会损失其部分能量(约30～50 电子伏特)，这部分能量激发核外电子脱离原子，能量大于材料逸出功的价电子可从样品表面逸出，变成真空中的自由电子，即二次电子。二次电子对试样表面状态非常敏感，能有效地显示试样表面的微观形貌。 （2）背散射电子。背散射电子是指被固体样品原子反射回来的一部分入射电子。既包括与样品中原子核作用而形成的弹性背散射电子，又包括与样品中核外电子作用而形成的非弹性散射电子。利用背反射电子作为成像信号不仅能分析形貌特征，也可以用来显示原子序数衬度，进行定性成分分析。 （3）X 射线。当入射电子和原子中内层电子发生非弹性散射作用时也会损失其部分能量(约

经典功率谱估计方法实现问题的研究

1 随机信号的经典谱估计方法 估计功率谱密度的平滑周期图是一种计算简单的经典方法。它的主要特点是与任 何模型参数无关，是一类非参数化方法[4]。它的主要问题是：由于假定信号的自相关函数在数据观测区以外等于零，因此估计出来的功率谱很难与信号的真实功率谱相匹配。在一般情况下，周期图的渐进性能无法给出实际功率谱的一个满意的近似，因而是一种低分辨率的谱估计方法。本章主要介绍了周期图法、相关法谱估计（BT ）、巴特利特(Bartlett)平均周期图的方法和Welch 法这四种方法。 2.1 周期图法 周期图法又称直接法。它是从随机信号x(n)中截取N 长的一段，把它视为能量有限x(n)真实功率谱)(jw x e S 的估计)(jw x e S 的抽样. 周期图这一概念早在1899年就提出了，但由于点数Ｎ一般比较大，该方法的计算量过大而在当时无法使用。只是1965年FFT 出现后，此法才变成谱估计的一个常用方法。周期图法[5]包含了下列两条假设： 1．认为随机序列是广义平稳且各态遍历的，可以用其一个样本x(n)中的一段 )(n x N 来估计该随机序列的功率谱。这当然必然带来误差。 2．由于对)(n x N 采用DFT ，就默认)(n x N 在时域是周期的，以及)(k x N 在频域是周期的。这种方法把随机序列样本x(n)看成是截得一段)(n x N 的周期延拓，这也就是周期图法这个名字的来历。与相关法相比，相关法在求相关函数)(m R x 时将 )(n x N 以外是数据全都看成零，因此相关法认为除)(n x N 外 x(n)是全零序列，这种处 理方法显然与周期图法不一样。 但是，当相关法被引入基于FFT 的快速相关后，相关法和周期图法开始融合。通过比较我们发现：如果相关法中M=N ，不加延迟窗，那么就和补充（N-1）个零的周期图法一样了。简单地可以这样说：周期图法是M=N 时相关法的特例。因此相关法和周期图法可结合使用。 2.2 相关法谱估计（BT ）法

现代谱估计

现代谱估计实验报告 1实验目的 功率谱估计在实际工程中有重要应用价值。如在语音信号识别、雷达杂波分析、波达方向估计、地震勘探信号处理、水声信号处理、系统辨识中非线性系统识别、物理光学中透镜干涉、流体力学的内波分析、太阳黑子活动周期研究等许多领域发挥了重要作用。 本次实验的目的主要是深入理解现代谱估计的基本理论，包括ARMA模型、ARMA谱估计。掌握现代谱估计的基本方法，包括SVD-TLS算法等。利用ARMA 功率谱估计中Cadzow谱估计子和Kaveh谱估计子来进行谱估计。 2实验原理 冃景 若离散随机过程｛x(n)｝服从线性差分方程 p q x(n) a i X(n i) e(n) b j e(n j) (1) i 1 j 1 式中e (n)是一离散白噪声，则称｛x(n)｝为ARMA过程，而式(1)所示的差分方程称为ARMA模型。系数a1,a2??p,和b1,b2?…b q,分别称为自回归参数和滑动平均参数，而p和q分别叫做AR阶数和MA阶数。式(1)所示的ARMA过程，其功率谱密度为 jw I / 、2〔B(z)| 2 B(e)l w両z e jw |B(e jw)| (2) Px ARMA谱估计的目的是使用N个已知的观测数据x(0)，x(1)：.x(N-1)计算出ARMA过程｛x(n)｝的功率谱密度估计。 在实际中，可以运用cadzow谱估计子和kaveh谱估计子来估计，cadzow谱估计子秩序确定AR阶数p和估计AR参数，而kaveh谱估计子也只需要确定AR 阶数p和估计AR参数以及MA阶数。

相关算法 AR阶数p的确定用奇异值分解(SVD,AR参数的估计用总体最小二乘法(TLS), 即应用(SVD-TLS算法来完成ARMA谱估计。 SVD-TLS 算法： 步骤1计算增广矩阵B的SVD,并储存奇异值和矩阵V; 步骤 2 确定增广矩阵 B 的有效秩p； 步骤 3 计算矩阵S; 步骤4求S的逆矩阵S--并计算出未知参数的总体最小二乘估计。 3 实验内容 仿真的观测数据由下式给出： xn = square(W*n)+*randn(1,N) ( 3) 其中, fs = 20000, n = 0:1/fs: , N = length(n), W = 2000*pi 。 1、采样周期图法进行谱估计 2、假设AR阶数未知，用SVD-TLS^法确定AR阶数和参数，然后使用Cadzow 谱估计子进行谱估计。 4 Matlab 仿真 仿真的观测数据时域信号如图 1 所示

基于Burg算法的AR模型功率谱估计简介

基于Burg 算法的AR 模型功率谱估计简介 摘要：在对随机信号的分析中，功率谱估计是一类重要的参数研究，功率谱估计的方法分为经典谱法和参数模型方法。参数模型方法是利用型号的先验知识，确定信号的模型，然后估计出模型的参数，以实现对信号的功率谱估计。根据wold 定理,AR 模型是比较常用的模型，根据Burg 算法等多种方法可以确定其参数。 关键词：功率谱估计；AR 模型；Burg 算法 随机信号的功率谱反映它的频率成分以及各成分的相对强弱, 能从频域上揭示信号的节律, 是随机信号的重要特征。因此, 用数字信号处理手段来估计随机信号的功率谱也是统计信号处理的基本手段之一。在信号处理的许多应用中, 常常需要进行谱估计的测量。例如, 在雷达系统中, 为了得到目标速度的信息需要进行谱测量; 在声纳系统中, 为了寻找水面舰艇或潜艇也要对混有噪声的信号进行分析。总之, 在许多应用领域中, 例如, 雷达、声纳、通讯声学、语言等领域, 都需要对信号的基本参数进行分析和估计, 以得到有用的信息, 其中, 谱分析就是一类最重要的参数研究。 1 功率谱估计简介 一个宽平稳随机过程的功率谱是其自相关序列的傅里叶变换，因此功率谱估计就等效于自相关估计。对于自相关各态遍历的过程，应有： )()()(121lim *k r n x k n x N N x N N n =? ?????++∞→∑-= 如果所有的)(n x 都是已知的，理论上功率谱估计就很简单了，只需要对其自相关序列取傅里叶变换就可以了。但是，这种方法有两个个很大的问题：一是不是所有的信号都是平稳信号，而且有用的数据量可能只有很少的一部分；二是数据中通常都会有噪声或群其它干扰信号。因此，谱估计就是用有限个含有噪声的观测值来估计)(jw x e P 。 谱估计的方法一般分为两类。第一类称为经典方法或参数方法，它首先由给定的数据估 计自相关序列)(k r x ，然后对估计出的)(?k r x 进行傅里叶变换获得功率谱估计。第二类称为非经典法，或参数模型法，是基于信号的一个随机模型来估计功率谱。非参数谱估计的缺陷是其频率分辨率低，估计的方差特性不好, 而且估计值沿频率轴的起伏甚烈，数据越长, 这一现象越严重。 为了改善谱分辨率，研究学者对基于模型的参数方法进行了大量研究。参数方法的第一步是对信号选择一个合适的模型，这种选择可能是基于有关信号如何产生的先验知识，也可能是多次试验后获得的结果。通常采用的模型包括AR 、MA 、ARMA 模型和谐波模型（噪声中含有复指数）。一旦模型选择好后，下一步就是计算模型的参数。最后将计算得到的参数带

材料现代分析方法练习题及答案

8. 什么是弱束暗场像？与中心暗场像有何不同？试用Ewald图解说明。 答：弱束暗场像是通过入射束倾斜，使偏离布拉格条件较远的一个衍射束通过物镜光阑，透射束和其他衍射束都被挡掉，利用透过物镜光阑的强度较弱的衍射束成像。 与中心暗场像不同的是，中心暗场像是在双光束的条件下用的成像条件成像，即除直射束外只有一个强的衍射束，而弱束暗场像是在双光阑条件下的g/3g的成像条件成像，采用很大的偏离参量s。中心暗场像的成像衍射束严格满足布拉格条件，衍射强度较强，而弱束暗场像利用偏离布拉格条件较远的衍射束成像，衍射束强度很弱。采用弱束暗场像，完整区域的衍射束强度极弱，而在缺陷附近的极小区域内发生较强的反射，形成高分辨率的缺陷图像。图：PPT透射电子显微技术1页 10. 透射电子显微成像中，层错、反相畴界、畴界、孪晶界、晶界等衍衬像有何异同？用什么办法及根据什么特征才能将它们区分开来？ 答：由于层错区域衍射波振幅一般与无层错区域衍射波振幅不同，则层错区和与相邻区域形成了不同的衬度，相应地出现均匀的亮线和暗线，由于层错两侧的区域晶体结构和位相相同，故所有亮线和暗线的衬度分别相同。层错衍衬像表现为平行于层错面迹线的明暗相间的等间距条纹。 孪晶界和晶界两侧的晶体由于位向不同，或者还由于点阵类型不同，一边的晶体处于双光束条件时，另一边的衍射条件不可能是完全相同的，也可能是处于无强衍射的情况，就相当于出现等厚条纹，所以他们的衍衬像都是间距不等的明暗相间的条纹，不同的是孪晶界是一条直线，而晶界不是直线。 反相畴界的衍衬像是曲折的带状条纹将晶粒分隔成许多形状不规则的小区域。 层错条纹平行线直线间距相等 反相畴界非平行线非直线间距不等 孪晶界条纹平行线直线间距不等 晶界条纹平行线非直线间距不等 11.什么是透射电子显微像中的质厚衬度、衍射衬度和相位衬度。形成衍射衬度像和相位衬度像时，物镜在聚焦方面有何不同？为什么？ 答：质厚衬度：入射电子透过非晶样品时，由于样品不同微区间存在原子序数或厚度的差异，导致透过不同区域落在像平面上的电子数不同，对应各个区域的图像的明暗不同，形成的衬度。 衍射衬度：由于样品中的不同晶体或同一晶体中不同部位的位向差异导致产生衍射程度不同而形成各区域图像亮度的差异，形成的衬度。 相位衬度：电子束透过样品，试样中原子核和核外电子产生的库伦场导致电子波的相位发生变化，样品中不同微区对相位变化作用不同，把相应的相位的变化情况转变为相衬度，称为相位衬度。 物镜聚焦方面的不同：透射电子束和至少一个衍射束同时通过物镜光阑成像时，透射束和衍射束相互干涉形成反应晶体点阵周期的条纹成像或点阵像或结构物象，这种相位衬度图像的形成是透射束和衍射束相干的结果，而衍射衬度成像只用透射束或者衍射束成像。

功率谱估计

功率谱估计及其MATLAB仿真 詹红艳 （201121070630控制理论与控制工程） 摘要:从介绍功率谱的估计原理入手分析了经典谱估计和现代谱估计两类估计方法的原理、各自特点及在Matlab中的实现方法。 关键词:功率谱估计;周期图法;AR参数法;Matlab Power Spectrum Density Estimation and the simulation in Matlab Zhan Hongyan （201121070630Control theory and control engineering） Abstract:Mainly introduces the principles of classical PSD estimation and modern PSD estimation,discusses the characteristics of the methods of realization in Matlab.Moreover,It gives an example of each part in realization using Matlab functions. Keywords:PSDPstimation,Periodogram method,AR Parameter method,Matlab 1引言 现代信号分析中,对于常见的具有各态历经的平稳随机信号,不可能用清楚的数学关系式来描述,但可以利用给定的N个样本数据估计一个平稳随机信号的功率谱密度叫做功率谱估计(PSD)。它是数字信号处理的重要研究内容之一。功率谱估计可以分为经典功率谱估计(非参数估计)和现代功率谱估计(参数估计)。 功率谱估计在实际工程中有重要应用价值,如在语音信号识别、雷达杂波分析、波达方向估计、地震勘探信号处理、水声信号处理、系统辨识中非线性系统识别、物理光学中透镜干涉、流体力学的内波分析、太阳黑子活动周期研究等许多领域,发挥了重要作用。 Matlab是MathWorks公司于1982年推出的一套高性能的数值计算和可视化软件,人称矩 阵实验室,它集数值分析、矩阵运算、信号处理和图形显示于一体,构成了一个方便的、界面友好的用户环境,成为目前极为流行的工程数学分析软件。也为数字信号处理进行理论学习、工程设计分析提供了相当便捷的途径。本文的仿真实验中,全部在Matlab6.5环境下调试通过;随机序列由频率不同的正弦信号加高斯白噪声组成。 2经典功率谱估计 经典功率谱估计是将数据工作区外的未知数据假设为零,相当于数据加窗。经典功率谱估计方法分为:相关函数法(BT法)、周期图法以及两种改进的周期图估计法即平均周期图法和平滑平均周期图法,其中周期图法应用较多,具有代表性。 1.1相关函数法(BT法) 该方法先由序列x(n)估计出自相关函数R(n),然后对R(n)进行傅立叶变换,便得到x(n)的功率谱估计。当延迟与数据长度相比很小时,可以有良好的估计精度。 Matlab代码示例1: Fs=500;%采样频率 n=0:1/Fs:1;

材料现代分析方法北京工业大学

材料现代分析方法北京工业大学 篇一：13103105-材料现代分析方法 《材料现代分析方法》课程教学大纲 一、课程基本信息 课程编号：13103105 课程类别：专业核心课程 适应专业：材料物理 总学时：54学时 总学分:3 课程简介： 本课程介绍材料微观形貌、结构及成分的分析与表面分析技术主要方法及基本技术，简单介绍光谱分析方法。包括晶体X射线衍射、电子显微分析、X射线光电子谱仪、原子光谱、分子光谱等分析方法及基本技术。 授课教材：《材料分析测试方法》，黄新民解挺编，国防工业出版社，20XX年。 参考书目： [1]《现代物理测试技术》，梁志德、王福编，冶金工业出版社，20XX 年。 [2]《X射线衍射分析原理与应用》，刘粤惠、刘平安编，化学工业出

版社，20XX年。 [3]《X射线衍射技术及设备》，丘利、胡玉和编，冶金工业出版社，20XX年。 [4]《材料现代分析方法》，左演声、陈文哲、梁伟编，北京工业大学出版社，20XX年。 [5]《材料分析测试技术》，周玉、武高辉编，哈尔滨工业大学出版社，2000年。 [6]《材料结构表征及应用》，吴刚编，化学工业出版社，20XX年。 [7]《材料结构分析基础》，余鲲编，科学出版社，20XX年。 二、课程教育目标 通过学习，了解X射线衍射仪及电子显微镜的结构，掌握X-射线衍射及电子显微镜的基本原理和操作方法，了解试样制备的基本要求及方法，了解材料成分的分析与表面分析技术的主要方法及基本技术，了解光谱分析方法，能够利用上述相关仪器进行材料的物相组成、显微结构、表面分析研究。学会运用以上技术的基本方法，对材料进行测试、计算和分析，得到有关微观组织结构、形貌及成分等方面的信息。 三、教学内容与要求 第一章X射线的物理基础 教学重点：X射线的产生及其与物质作用原理 教学难点：X射线的吸收和衰减、激发限 教学时数：2学时

利用经典谱估计法估计信号的功率谱(随机信号)

随机信号 利用经典谱估计法估计信号的功率谱

作业综述： 给出一段信号“asd.wav”,利用经典谱估计法的原理，通过不同的谱估计方法，求出信号的功率谱密度函数。采用MATLAB语言，利用MATLAB语言强大的数据处理和数据可视化能力，通过GUI的对话框模板，使操作更为简便！在一个GUI界面中，同时呈现出不同方法产生出的功率谱。 这里给出了几种不同的方法：BT法，周期图法，平均法以及Welch法。把几种不同方法所得到的功率谱都呈现在一个界面中，便于对几种不同方法得到的功率谱作对比。 一．题目要求 给出一段信号及采样率，利用经典谱估计法估计出信号的功率谱。 二．基本原理及方法 经典谱估计的方法，实质上依赖于传统的傅里叶变换法。它是将数据工作区外的未知数据假设为零，相当于数据加窗，主要方法有BT法，周期图法，平均法以及Welch法。 1. BT法(Blackman-Tukey) ●理论基础： (1)随机序列的维纳-辛钦定理 由于随机序列{X(n)}的自相关函数Rx(m)=E[X(n)X(n+m)]定义在离散点m上，设取样间隔为，则可将随机序列的自相关函数用连续时间函数表示为 等式两边取傅里叶变换，则随机序列的功率谱密度 (2)谱估计 BT法是先估计自相关函数Rx(m)(|m|=0,1,2…,N-1),然后再经过离散傅里叶变换求的功率谱密度的估值。即 其中可有式得到。 2. 周期图法 ●理论基础： 周期图法是根据各态历经随机过程功率谱的定义来进行谱估计的。在前面我们已知，各态历经的连续随机过程的功率谱密度满足

式中 是连续随机过程第i 个样本的截取函数 的频谱。对应在随机序列中则有 由于随机序列中观测数据 仅在 的点上存在，则 的N 点离散傅里叶变换为： 因此有随机信号的观测数据 的功率谱估计值（称“周期图”）如下： 由于上式中的离散傅里叶变换可以用快速傅里叶变换计算，因此就可以估计出功率 谱。 3.平均法： 理论基础： 平均法可视为周期图法的改进。周期图经过平均后会使它的方差减少，达到一致估计的目的，有一个定理：如果 ， ， ， 是不相关的随机变量，且都有个均值 及其方差 ，则可以证明它们的算术平均的均值为 ，方差为 。 由定理可见：具有 个独立同分布随机变量平均的方差，是单个随机变量方差的 ， 当 时，方差 ，可以达到一致估计的目的。因此，将 个独立的估计量经过算术 平均后得到的估计量的方差也是原估计量方差的 。 平均图法即是将数据 ， ， 分段求周期图法后再平均。例如，给定N=1000个数据样本（平均法适用于数据量大的场合），则可以将它分成10个长度为100的小段，分别计算每一段的周期图 ()()2 1001100,100(1) 1 ,1,2,```,10100 l j l n l G w X e l ω-=-= =∑ 然后将这10个周期图加以平均得谱估计值： ()() 10 100100,1 110l l G w G w ==∑ 由于这10小段的周期图取决于同一个过程，因而其均值相同。若这10个小段的周期图是统计独立的，则这10个小段平均之后的方差却是单段方差的 。

基于经典谱估计的多普勒频移算法

仿真程序说明文档 1．平坦衰落信道仿真 1.1仿真设计 1.产生频率为c ω的01N +路正弦余弦信号； 2.产生在[0,2)π均匀分布的随机相位，产生在[,]m m f f -均匀分布的频率偏移，将随机相位和频率偏移加入到第一步产生的每路正余弦信号中； 3.产生每路信号的衰减cos n β，将第二步产生的每路信号乘以衰减cos n β； 4.将加入随机相位频率偏移以及乘以衰减以后的正弦信号叠加，得到同相分量，将加入随机相位频率偏移以及乘以衰减以后的余弦信号叠加，得到正交分量； 5.由同相和正交信号分别作为实部()c T t 和虚部()s T t 得到平坦衰落信道的输出； 1.2程序流程图

1.3仿真结果 下面给出了仿真的平坦衰落信道的统计特性 图1通过信道后的接收信号包络 在图1中用信道的最大增益对衰落进行了归一化。可以看出，仿真的数据流能够较好的符合典型的Rayleigh衰落信号。信道在某些点会引起深度衰落。 图2通过仿真信道信号的包络概率密度函数

在图2中，绘出了当N =34时的包络分布。可以看出当N =34时，包络分布与标准的瑞利分布基本吻合。随着N 的增加，包络更加趋向瑞利分布，且分布函数与时间t 无关，这一点满足广义平稳过程的要求。 图3通过仿真信道信号的自相关函数 可以看出自相关函数趋近贝赛尔(Bessel)函数。 以上我们讨论了仿真信道产生的随机过程是广义平稳的，并且其信号包络、包络概率分布、自相关性等统计特性，都能与Clarke 模型较好的吻合，因此能够较真实的反映信道。 2．LCR 仿真 2.1仿真设计 1.信号采样：首先对接收信号进行采样，得到输入信号的离散序列()g n ； 2.计算包络：根据输入的()g n 序列，计算相应的包络()()n g n α=； 3.确定电平R ：计算()n α的均方根包络电平rms R ，使rms R R =； 4.估计LCR L ：根据第二、第三步计算()n α和rms R 估计()n α每秒通过rms R 的次数LCR L ； 5.求解v ：进行数值计算，求LCR 法的速度估计值v ；

材料现代分析方法实验报告

力学与材料学院 材料现代分析方法实验报告二 XRD图谱分析 专业年级：1 姓名：1 指导老师：1 学号：1 2016年12月 中国南京 目录 实验名称：XRD图谱分析…………………………………………… 一、实验目的……………………………………………………

二、实验要求…………………………………………………… 三、操作过程…………………………………………………… 四、结果分析与讨论……………………………………………… 实验名称：XRD图谱分析 一、实验目的 了解XRD基本原理及其应用，不同物相晶体结构XRD图谱的区别，熟练掌握如何来分析利用X射线测试得到的XRD图谱。 二、实验要求

1、熟练掌握如何来利用软件打开、分析XRD图谱，以及输出分析结果。 2、明确不同物质的XRD图谱，掌握XRD图谱包含的晶体结构的关系，通过自己分析、数据查找和鉴别的全过程，了解如何利用软件正确分析和确定不同物相的XRD图谱，并输出分析结果。 3、实验报告的编写，要求报告能准确的反映实验目的、方法、过程及结论。 三、操作过程 1、启动Jade 6.0，并打开实验数据。 2、点击图标使图谱平滑后，再连续两次点击图标扣除背景影响。 3、右击工具栏中的图标，全选左侧的项目，取消选择右侧中的Use Chemistry Filter，最后在下方选择S/M Focus on Major Phases（如图一），并点击OK。 图一

4、得到物相分析，根据FOM值（越小，匹配性越高）可推断出该物相为以ZnO为主，可能含有CaF2、Al2O3、Mg(OH)2混合组成的物质（如图二），双击第一种物质可以得到主晶相的PDF卡片（如图三），点击图三版面中的Lines可以观察到不同角度处的衍射强度（如图四）。 图二

(完整版)功率谱估计性能分析及Matlab仿真

功率谱估计性能分析及Matlab 仿真 1 引言 随机信号在时域上是无限长的，在测量样本上也是无穷多的，因此随机信号的能量是无限的，应该用功率信号来描述。然而，功率信号不满足傅里叶变换的狄里克雷绝对可积的条件，因此严格意义上随机信号的傅里叶变换是不存在的。因此，要实现随机信号的频域分析，不能简单从频谱的概念出发进行研究，而是功率谱[1]。 信号的功率谱密度描述随机信号的功率在频域随频率的分布。利用给定的 N 个样本数据估计一个平稳随机信号的功率谱密度叫做谱估计。谱估计方法分为两大类：经典谱估计和现代谱估计。经典功率谱估计如周期图法、自相关法等，其主要缺陷是描述功率谱波动的数字特征方差性能较差，频率分辨率低。方差性能差的原因是无法获得按功率谱密度定义中求均值和求极限的运算[2]。分辨率低的原因是在周期图法中，假定延迟窗以外的自相关函数全为0。这是不符合实际情况的，因而产生了较差的频率分辨率。而现代谱估计的目标都是旨在改善谱估计的分辨率，如自相关法和Burg 法等。 2 经典功率谱估计 经典功率谱估计是截取较长的数据链中的一段作为工作区，而工作区之外的数据假设为0，这样就相当将数据加一窗函数，根据截取的N 个样本数据估计出其功率谱[1]。 周期图法( Periodogram ) Schuster 首先提出周期图法。周期图法是根据各态历经的随机过程功率谱的定义进行的谱估计。 取平稳随机信号()x n 的有限个观察值(0),(1),...,(1)x x x n -，求出其傅里叶变换 1 ()()N j j n N n X e x n e ω ω---==∑ 然后进行谱估计

(完整版)材料现代分析方法考试试卷

班级学号姓名考试科目现代材料测试技术A 卷开卷一、填空题（每空1 分，共计20 分；答案写在下面对应的空格处，否则不得分） 1. 原子中电子受激向高能级跃迁或由高能级向低能级跃迁均称为_辐射跃迁__ 跃迁或_无辐射跃迁__跃迁。 2. 多原子分子振动可分为__伸缩振动_振动与_变形振动__振动两类。 3. 晶体中的电子散射包括_弹性、__与非弹性___两种。 4. 电磁辐射与物质（材料）相互作用，产生辐射的_吸收_、_发射__、_散射/光电离__等，是光谱分析方法的主要技术基础。 5. 常见的三种电子显微分析是_透射电子显微分析、扫描电子显微分析___和_电子探针__。 6. 透射电子显微镜（TEM）由_照明__系统、_成像__系统、_记录__系统、_真空__系统和__电器系统_系统组成。 7. 电子探针分析主要有三种工作方式，分别是_定点_分析、_线扫描_分析和__ 面扫描_分析。 二、名词解释（每小题3 分，共计15 分；答案写在下面对应的空格处，否则不得分） 1. 二次电子二次电子：在单电子激发过程中被入射电子轰击出来的核外电子. 2. 电磁辐射：在空间传播的交变电磁场。在空间的传播遵循波动方程，其波动性表现为反射、折射、干涉、衍射、偏振等。 3. 干涉指数：对晶面空间方位与晶面间距的标识。 4. 主共振线：电子在基态与最低激发态之间跃迁所产生的谱线则称为主共振线 5. 特征X 射线：迭加于连续谱上，具有特定波长的X 射线谱，又称单色X 射线谱。 三、判断题（每小题2 分，共计20 分；对的用“√”标识，错的用“×”标识） 1．当有外磁场时，只用量子数n、l 与m 表征的原子能级失去意义。（√） 2．干涉指数表示的晶面并不一定是晶体中的真实原子面，即干涉指数表示的晶面上不一定有原子分布。（√） 3．晶面间距为d101／2 的晶面，其干涉指数为（202）。（×） 4．X 射线衍射是光谱法。（×） 5．根据特征X 射线的产生机理，λKβ<λK α。 （√ ） 6．物质的原子序数越高，对电子产生弹性散射的比例就越大。（√ ） 7．透射电镜分辨率的高低主要取决于物镜。（√ ）8．通常所谓的扫描电子显微镜的分辨率是指二次电子像的分辨率。（√）9．背散射电子像与二次电子像比较，其分辨率高，景深大。（× ）10．二次电子像的衬度来源于形貌衬度。（× ） 四、简答题（共计30 分；答案写在下面对应的空格处，否则不得分） 1. 简述电磁波谱的种类及其形成原因？（6 分）答：按照波长的顺序，可分为：（1）长波部分，包括射频波与微波。长波辐射光子能量低，与物质间隔很小的能级跃迁能量相适应，主要通过分子转动能级跃迁或电子自旋或核自旋形成；（2）中间部分，包括紫外线、可见光核红外线，统称为光学光谱，此部分辐射光子能量与原子或分子的外层电子的能级跃迁相适应；（3）短波部分，包括X 射线和γ射线，此部分可称射线谱。X 射线产生于原子内层电子能级跃迁，而γ射线产生于核反应。

材料现代分析方法试题 6

材料现代分析方法试题 一、基本概念题（共10题，每题5分） 1．什么是光电效应？光电效应在材料分析中有哪些用途？ 2．当波长为λ的X射线在晶体上发生衍射时，相邻两个（hkl）晶面衍射线的 波程差是多少？相邻两个HKL干涉面的波程差又是多少？ 3．测角仪在采集衍射图时，如果试样表面转到与入射线成30 0角，则计数管 与入射线所成角度为多少？能产生衍射的晶面，与试样的自由表面是何种几何关 系？ 4．宏观应力对X射线衍射花样的影响是什么？衍射仪法测定宏观应力的方法 有哪些？ 5．薄膜样品的基本要求是什么? 具体工艺过程如何? 双喷减薄与离子减薄 各适用于制备什么样品? 6．图说明衍衬成像原理，并说明什么是明场像、暗场像和中心暗场像。 7．说明透射电子显微镜成像系统的主要构成、安装位置、特点及其作用。 8．何为晶带定理和零层倒易截面? 说明同一晶带中各晶面及其倒易矢量与 晶带轴之间的关系。 9．含苯环的红外谱图中，吸收峰可能出现在哪4个波数范围? 10．陶瓷纳米/微米颗粒的红外光谱的分析样品该如何制，为什么？ 二、综合及分析题（共5题，每题10分） 1．请说明多相混合物物相定性分析的原理与方法？ 2．对于晶粒直径分别为100，75，50，25nm的粉末衍射图形，请计算由于晶粒细化引起的衍射线条宽化幅度B（设θ=450，λ=0.15nm）。对于晶粒直径为25nm的粉末，试计算θ=100、450、800时的B 值。 3．二次电子像和背散射电子像在显示表面形貌衬度时有何相同与不同之处? 4．何为波谱仪和能谱仪？说明其工作的三种基本方式及其典型应用，并比较波谱仪和能谱仪的优缺点。要分析钢中碳化物成分和基体中碳含量，应选用哪种电子探针仪? 为什么? 5．分别指出谱图中标记的各吸收峰所对应的基团? 材料现代分析方法试题（参考答案） 一、基本概念题（共10题，每题5分） 1．什么是光电效应？光电效应在材料分析中有哪些用途？ 答：光电效应是指：当用X射线轰击物质时，若X射线的能量大于物质原子 对其内层电子的束缚力时，入射X射线光子的能量就会被吸收，从而导致其内层 电子被激发，产生光电子。材料分析中应用光电效应原理研制了光电子能谱仪和 荧光光谱仪，对材料物质的元素组成等进行分析。 2．什么叫干涉面？当波长为λ的X射线在晶体上发生衍射时，相邻两个（hkl） 晶面衍射线的波程差是多少？相邻两个HKL干涉面的波程差又是多少？ 答：晶面间距为d’/n、干涉指数为nh、nk、nl的假想晶面称为干涉面。当波 长为λ的X射线照射到晶体上发生衍射，相邻两个（hkl）晶面的波程差是nλ， 相邻两个（HKL）晶面的波程差是λ。

作业——现代谱估计法

现代谱估计法 （殷恒刚 107010254） 1. 现代谱估计简介 经典谱估计法可以利用FFT 计算，因而有计算效率高的优点，在谱分辨力要求不是太高的地方常用这种方法。但频率分辨率地是经典谱估计的一个无法回避的缺点。如周期图法在计算中把观测到的有限长的N 个数据以外的数据认为是零，而BT 法仅利用N 个有限的观测数据作自相关函数估计，实质上也就是假设除已知数据外的自相关函数全为零，这些显然都是与事实不符的。为了克服以上缺点，人们提出了平均，加窗平滑等方法，在一定程度上改善了经典谱估计的性能。但是，经典谱估计，始终无法解决，频率分辨率与谱估计稳定性之间的矛盾，特别是在数据记录长度比较短时，这一矛盾尤其突出。 现代谱估计理论也就是在这种背景下产生的，以1967年Burg 提出的最大熵谱分析法为代表的现代谱估计法，不认为在观察到的N 个数据以外的数据全为零。因此克服了经典法的这个缺点，提高了谱估计的分辨率。后来发现线性预测自回归模型法(简称AR 模型法)与Burg 的最大熵谱分析法是等价的，它们都可归结为通过Yule-Walker 方程求解自回归模型的系数问题。目前常用的求自回归模型系数的算法有三种：①为Levinson 递推算法；②为Burg 递推算法；③为正反向线性预测最小二乘算法。 2.现代谱估计的三种模型 由信号与系统相关知识可知，任何具有有理功率谱密度的随机信号都可以看成是由一白噪声激励一物理网络所形成。如图一所示。我们可以先假设一个模型，然后根据已记录数据估计参数值，这样就不用假设N 以外的所有数据全为零，这就克服了经典谱估计的缺点。 图1 一个系统的Z 域传递函数的一般形式如下： 00 ()() b a n j j j n i i i b z Y z X z a z -=-== ∑∑ (1.1) 参数建模的任务也就是如何确定阶数a n 和b n 以及系统数组(1,,)i a a i n = 和

材料现代分析方法知识点

材料现代分析方法知识点 1.什么是特征X射线？ 当管压增至与阳极靶材对应的特定值U k时，在连续谱的某些特定波长位置上出现一系列陡峭的尖峰。该尖峰对应的波长λ与靶材的原子序数Z存在着严格的对应关系，尖峰可作为靶材的标志或特征，故称尖峰为特征峰或特征谱。 2.什么是电子探针的点分析、线分析、面分析？ ①点分析：将电子束作用于样品上的某一点，波谱仪分析时改变分光晶体和探测器的位置，收集分析点的特征X射线，由特征X射线的波长判定分析点所含的元素；采用能谱仪工作时，几分钟内可获得分析点的全部元素所对应的特征X射线的谱线，从而确定该点所含有的元素及其相对含量。②线分析：将探针中的谱仪固定于某一位置，该位置对应于某一元素特征X射线的波长或能量，然后移动电子束，在样品表面沿着设定的直线扫描，便可获得该种元素在设定直线上的浓度分布曲线。改变谱仪位置则可获得另一种元素的浓度分布曲线。③面分析：将谱仪固定于某一元素特征X射线信号（波长或能量）位置上，通过扫描线圈使电子束在样品表面进行光栅扫描（面扫描），用检测到的特征X射线信号调制成荧光屏上的亮度，就可获得该元素在扫描面内的浓度分布图像。 3. XRD对样品有何要求？ 粉末样品应干燥，粒度一般要求约10~80μm，应过200目筛子（约0.08mm），且避免颗粒不均匀。块状样品应将其处理成与窗孔大小一致，可扫描宽度宜大于5mm，小于30mm，至少保证一面平整。 4.电子探针分析原理？ 电子探针是一中利用电子束作用样品后产生的特征X射线进行微区成分分析的仪器。其结构与扫描电竞基本相同，所不同的只是电子探针检测的是特征X射线，而不是二次电子或背散射电子。 5.结构因子的计算？P68 （1）简单点阵：简单点阵的晶胞仅有一个原子，坐标为（0，0，0），即X=Y=Z=0，设原子的散射因子为f，则(公式3-69) （2）底心点阵：底心点阵的晶胞有两个原子，坐标分别为（0，0，0），（1/2,1/2,0）各原子的散射因子为f，则（公式3-70） （3）体心点阵：体心点阵的晶胞有两个原子，坐标分别为（0，0，0），（1/2,1/2,1/2）各原子的散射因子为f，则（公式3-71） （4）面心点阵：面心点阵的晶胞有4个原子，坐标分别为（0，0，0），（1/2,1/2,0），（1/2,0,1/2），（0,1/2,1/2）各原子的散射因子为f，则（公式3-72） 6.X射线衍射与电子衍射的关系（比较）？P150 （1）电子波的波长短，远远小于X射线，同等衍射条件下，它的衍射半角很小，衍射束集中在前方额，而x射线的衍射半角可接近90度。 （2）电子衍射反射球半径大 （3）电子衍射散射强度高，物质对电子的散射比对x射线散射强约1000000倍 (4) 电子衍射不仅可以进行微区结构分析，还可以进行形貌观察，而x射线衍射却无法进行形貌分析 （5）薄晶样品的倒易点阵为沿厚度方向的倒易杆，大大增加了反射球与倒易杆相截的机会，即使偏离布拉格方

经典谱估计算法研究与实现

毕业设计（论文） 论文题目：经典谱估计算法研究与实现 教学中心：电子科技大学网络教育学院苏州学习中心指导老师：职称： 学生姓名：学号： 专业：通信工程

毕业设计（论文）任务书 题目：经典谱估计算法研究与实现 任务与要求： 探讨我国经典谱估计算法运用和影响下所面临的机遇与挑战以及经典谱估计算法研究与实现,结合所知识,理论联系实际, 写出毕业论文。 时间：2014年 1 月25日至 2014年 4 月 14日共 12 周教学中心：电子科技大学网络教育学院苏州学习中心 学生姓名：学号： 专业：通信工程 指导单位或教研室：电子科技大学网络教育学院苏州学习中心 指导教师：职称：

毕业设计(论文)进度计划表

电子科技大学毕业设计(论文)中期检查记录表 注:此表同学生毕业设计(论文)一起存档

1 随机信号的经典谱估计方法 估计功率谱密度的平滑周期图是一种计算简单的经典方法。它的主要特点是 与任何模型参数无关，是一类非参数化方法[4]。它的主要问题是：由于假定信号的自相关函数在数据观测区以外等于零，因此估计出来的功率谱很难与信号的真实功率谱相匹配。在一般情况下，周期图的渐进性能无法给出实际功率谱的一个满意的近似，因而是一种低分辨率的谱估计方法。本章主要介绍了周期图法、相关法谱估计（BT ）、巴特利特(Bartlett)平均周期图的方法和Welch 法这四种方法。 2.1 周期图法 周期图法又称直接法。它是从随机信号x(n)中截取N 长的一段，把它视为能量有限x(n)真实功率谱)(jw x e S 的估计)(jw x e S 的抽样. 周期图这一概念早在1899年就提出了，但由于点数Ｎ一般比较大，该方法的计算量过大而在当时无法使用。只是1965年FFT 出现后，此法才变成谱估计的一个常用方法。周期图法[5]包含了下列两条假设： 1．认为随机序列是广义平稳且各态遍历的，可以用其一个样本x(n)中的一段)(n x N 来估计该随机序列的功率谱。这当然必然带来误差。 2．由于对)(n x N 采用DFT ，就默认)(n x N 在时域是周期的，以及)(k x N 在频域是周期的。这种方法把随机序列样本x(n)看成是截得一段)(n x N 的周期延拓，这也就是周期图法这个名字的来历。与相关法相比，相关法在求相关函数)(m R x 时将)(n x N 以外是数据全都看成零，因此相关法认为除)(n x N 外x(n)是全零序列，这种处理方法显然与周期图法不一样。 但是，当相关法被引入基于FFT 的快速相关后，相关法和周期图法开始融合。通过比较我们发现：如果相关法中M=N ，不加延迟窗，那么就和补充（N-1）个零的周期图法一样了。简单地可以这样说：周期图法是M=N 时相关法的特例。因此相关法和周期图法可结合使用。