地震动功率谱与反应谱的转换关系_赵凤新

基于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 模型和谐波模型（噪声中含有复指数）。一旦模型选择好后，下一步就是计算模型的参数。最后将计算得到的参数带

通信系统建模与仿真课程设计

通信系统建模与仿真课程设计2011 级通信工程专业1113071 班级 题目基于SIMULINK的基带传输系统的仿真姓名学号 指导教师胡娟 2014年6月27日

1任务书 试建立一个基带传输模型，采用曼彻斯特码作为基带信号，发送滤波器为平方根升余弦滤波器，滚降系数为0.5，信道为加性高斯信道，接收滤波器与发送滤波器相匹配。发送数据率为1000bps，要求观察接收信号眼图，并设计接收机采样判决部分，对比发送数据与恢复数据波形，并统计误码率。另外，对发送信号和接收信号的功率谱进行估计。假设接收定时恢复是理想的。 2基带系统的理论分析 1.基带系统传输模型和工作原理 数字基带传输系统的基本组成框图如图1 所示，它通常由脉冲形成器、发送滤波器、信道、接收滤波器、抽样判决器与码元再生器组成。系统工作过程及各部分作用如下。 g T(t) n 定时信号 图 1 ：数字基带传输系统方框图 发送滤波器进一步将输入的矩形脉冲序列变换成适合信道传输的波形g T(t)。这是因为矩形波含有丰富的高频成分，若直接送入信道传输，容易产生失真。 基带传输系统的信道通常采用电缆、架空明线等。信道既传送信号，同时又因存在噪声n(t)和频率特性不理想而对数字信号造成损害，使得接收端得到的波形g R(t)与发送的波形g T(t)具有较大差异。 接收滤波器是收端为了减小信道特性不理想和噪声对信号传输的影响而设置的。其主要作用是滤除带外噪声并对已接收的波形均衡，以便抽样判决器正确判决。 抽样判决器首先对接收滤波器输出的信号y(t)在规定的时刻（由定时脉冲cp控制）进行抽样，获得抽样信号{r n}，然后对抽样值进行判决，以确定各码元是“1”码还是“0”码。 2.基带系统设计中的码间干扰和噪声干扰以及解决方案

功率谱估计

功率谱估计及其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;

功率谱估计方法的比较

功率谱估计方法的比较 摘要： 本文归纳了信号处理中关键的一种分析方法, 即谱估计方法。概述了频谱估计中的周期图法、修正的协方差法和伯格递推法的原理，并且对此三种方法通过仿真做出了对比。 关键词：功率谱估计；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

振型分解反应谱法

结构设计系列之振型分解反应谱法 苏义

前言 我国规范对于常规结构设计有两个方法：底部剪力法和振型分解反应谱法。其中，底部剪力法视多质点体系为等效单质点体系，且其地震作用沿高度呈倒三角形分布，当结构层数较高或体系较复杂时，其计算假再用，因部剪时，其计算假定不再适用，因此规范规定底部剪力法仅适用于高度不超过40m、以剪切变形为主且质量和刚度沿高度分布比较均匀的结构。因此，一般结构均采用振型分解反应谱法。

振型分解反应谱法的基本步骤： 通过体系的模态分析，求出多自由度体系的振型通过体系的模态分析求出多自由度体系的振型向量、参与系数等等；然后把每个振型看作单自由度体系，求出其在规定反应谱的地震加速度作用下产生的地震效应；最后把所有振型的地震效应式进行叠，得到体系震应应按一定方式进行叠加，就会得到体系地震效应的解。 注意 注意： 振型分解反应谱法只适用于弹性分析，对于弹塑性体系，由于力与位移不再具有对应关系，性体系，由于力与位移不再具有一一对应关系， 该法不再适用。

目录 一模态分析二 反应谱分析 三 振型组合方法 四 方向组合方法

一、模态分析 模态分析也被称作振型叠加法动力分析，是线性体系地震分析中最常用且最有效的方法。它最主要的 优势在于其计算一组正交向量之后，可以将大型 整体平衡方程组缩减为相对数量较少的解耦二阶平解阶微分方程，这样就明显减少了用于数值求解这些 方程的计算时间。模态分析为结构相关静力分析 提供相关结构性能，包括结构静力地震作用分析 和静力风荷载分析。 模态分析是其它动力分析的基础，包括反应谱分析和时程分析。

《随机信号处理》课程设计

《随机信号处理》课程设计

————————————————————————————————作者：————————————————————————————————日期：

华北水利水电大学 随机信号处理上机实验报告 学院：数学与信息科学 专业：信息与计算科学 姓名：孙志攀 学号：201216511 指导老师：蒋礼 日期：2015年10月20日

实验一 1、熟悉并练习使用下列Matlab 的函数，给出各个函数的功能说明和内部参数的意义，并给出至少一个使用例子和运行结果 1.rand() (1)Y = rand(n) 生成n×n 随机矩阵，其元素在（0，1）内 (2)Y = rand(m,n) 生成m×n 随机矩阵 (3)Y = rand([m n]) 生成m×n 随机矩阵 (4)Y = rand(m,n,p,…) 生成m×n×p×…随机矩阵或数组 (5)Y = rand([m n p…]) 生成m×n×p×…随机矩阵或数组 (6)Y = rand(size(A)) 生成与矩阵A 相同大小的随机矩阵 选择（3）作为例子，运行结果如下： 2.randn() 产生随机数数组或矩阵，其元素服从均值为0，方差为1的正态分布 （1）Y = randn 产生一个伪随机数 （2）Y = randn(n) 产生n×n的矩阵，其元素服从均值为0，方差为1的正态分布（3）Y = randn(m,n) 产生m×n的矩阵，其元素服从均值为0，方差为1的正态分布（4）Y= randn([m n]) 产生m×n的矩阵，其元素服从均值为0，方差为1的正态分布选择（3）作为例子，运行结果如下： 3．normrnd() 产生服从正态分布的随机数 （1）R = normrnd(mu,sigma) 产生服从均值为mu，标准差为sigma的随机数，mu和sigma 可以为向量、矩阵、或多维数组。 （2）R = normrnd(mu,sigma,v) 产生服从均值为mu 标准差为sigma的随机数，v是一个行向量。如果v是一个1×2的向量，则R为一个1行2列的矩阵。如果v是1×n的，那么R 是一个n维数组 （3）R = normrnd(mu,sigma,m,n) 产生服从均值为mu 标准差为sigma的随机数，标量m和n是R的行数和列数。

回归分析课程设计(最终版)

回归分析课程设计 （题目） （副标题） 指导教师 学院名称专业名称 设计提交日期年月

目录 1.课程设计简述-------------------------------------------------------2 2.多元线性回归-------------------------------------------------------3 3.违背基本假设的情况------------------------------------------------5 3.1 异方差性-------------------------------------------------------5 3.2 自相关性-------------------------------------------------------6 3.3 异常值检验-----------------------------------------------------6 4.自变量的选择与逐步回归--------------------------------------------7 4.1 所有子集回归---------------------------------------------------7 4.2 逐步回归--------------------------------------------------------8 5.多重共线性的情形及其处理-----------------------------------------10 5.1 多重共线性诊断------------------------------------------------10 5.2 消除多重共线性------------------------------------------------11 6.岭回归--------------------------------------------------------------12 7.主成分回归----------------------------------------------------------14 8.含定性变量的回归模型------------------------------------------------ 9.附录（程序代码）-----------------------------------------------------

参数法功率谱估计

参数法功率谱估计 一、信号的产生 （一）信号组成 在本实验中，需要事先产生待估计的信号，为了使实验结果较为明显，我产生了由两个不同频率的正弦信号（频率差相对较大）和加性高斯白噪声组成的信号。 （二）程序 N=1024;n=0:N-1; xn=2*cos(2*pi*0.2*n)+ cos(2*pi*0.213*n)+randn(1,1024); 这样就产生了加有白噪声的两个正弦信号 其波形如下

0100200300400500600 -8-6 -4 -2 2 4 6 8 10 (a) 两个正弦信号与白噪声叠加的时域波形 二、参数模型法功率谱估计 （一）算法原理简介 1．参数模型法是现代谱估计的主要内容，思路如下： ① 假定所研究的过程)(n x 是由一个白噪声序列)(n 激励一个因果稳定的可逆线性系统)(z H 的输出； ② 由已知的)(n x ，或其自相关函数)(m r x 估计)(z H 的参数； ③ 由)(z H 的参数来估计)(n x 的功率谱。 2．自回归模型，简称AR 模型，它是一个全极点的模型。“自回归”的含义是：该模型现在的输出是现在的输入和过去p 个输出的加权和。此模型可以表现

为以下三式：

① ∑=+--=p k k n u k n x a n x 1 )()()(； ② ∑=-+==p k k k z a z A z H 111)(1)(； ③ 212 1)(∑=-+=p k jwk k jw x e a e P σ。 3．AR 模型的正则方程建立了参数k a 和)(n x 的自相关函数的关系，公式如下： =)(m r x ∑=--p k x k k m r a 1)( 1≥m 时，=)(m r x 21)(σ+-∑=k r a p k x k 0=m 时。 （二）两种AR 模型阶次的算法 1．Yule-Walker 算法(自相关法) （1）算法主要思想 Yule-Walker 算法通过解Yule-Walker 方程获得AR 模型参数。从低阶开始递推，直到阶次p ，给出了在每一个阶次时的所有参数。公式如下： ① 11 11/])()()([--=-∑+--=m m k x x m m m r k m r k a k ρ； ② )()()(11k m a k k a k a m m m m -+=--；

振型分解反应谱法知识讲解

振型分解反应谱法

振型分解反应谱法 振型分解反应谱法是用来计算多自由度体系地震作用的一种方法。该法是利用单自由度体系的加速度设计反应谱和振型分解的原理，求解各阶振型对应的等效地震作用，然后按照一定的组合原则对各阶振型的地震作用效应进行组合，从而得到多自由度体系的地震作用效应。振型分解反应谱法一般可考虑为计算两种类型的地震作用：不考虑扭转影响的水平地震作用和考虑平扭藕联效应的地震作用。 适用条件 （1）高度不超过40米，以剪切变形为主且质量和刚度沿高度分布比较均匀的结构，以及近似于单质点体系的结构，可采用底部剪力法计算。（此为底部剪力法的适用范围） （2）除上述结构以外的建筑结构，宜采用“振型分解反应谱法”。 （3）特别不规则的建筑、甲类建筑和规范规定的高层建筑，应采用时程分析法进行补充计算。 刚重比 刚重比是指结构的侧向刚度和重力荷载设计值之比，是影响重力二阶效应的主要参数 刚重比=Di*Hi/Gi Di-第i楼层的弹性等效刚度，可取该层剪力与层间位移的比值Hi-第i楼层层高

Gi-第i楼层重力荷载设计值 刚重比与结构的侧移刚度成正比关系；周期比的调整将导致结构侧移刚度的变化，从而影响到刚重比。因此调整周期比时应注意，当某主轴方向的刚重比小于或接近规范限值时，应采用加强刚度的方法；当某主轴方向刚重比大于规范限值较多时，可采用削弱刚度的方法。同样，对刚重比的调整也可能影响周期比。特别是当结构的周期比接近规范限值时，应采用加强结构外围刚度的方法规范上限主要用于确定重力荷载在水平作用位移效应引起的二阶效应是否可以忽略不计。见高规5.4.1和5.4.2及相应的条文说明。刚重比不满足规范上限要求，说明重力二阶效应的影响较大，应该予以考虑。规范下限主要是控制重力荷载在水平作用位移效应引起的二阶效应不致过大，避免结构的失稳倒塌。见高规5.4.4及相应的条文说明。刚重比不满足规范下限要求，说明结构的刚度相对于重力荷载过小。但刚重比过分大，则说明结构的经济技术指标较差，宜适当减少墙、柱等竖向构件的截面面积。 长细比 长细比=计算长度/回转半径。 所以很显然，减小计算长度或者加大回转半径即可。 这里需要注意的是，计算长度并非实际长度，而是实际长度乘以长度系数，长度系数则与柱子两端的约束刚度有关。说白了就是

数字信号处理课程设计正余弦信号的谱分析

指导教师： 日期： 《数字信号处理》课程设计 题目：正余弦信号的谱分析 姓名： 院系：电子信息工程系 专业：通信工程 班级：通信091 学号： 指导教师： 2012年6 月

正余弦信号的谱分析 （电子信息工程学系 指导教师：留黎钦） 中文摘要：使用MATLAB 软件，通过编写程序，对正余弦信号进行傅里叶变换。用DFT 和FFT 实现对正余弦信号的 谱分析，并且分析DFT 长度对频谱的影响。 关键词：matlab ； 傅里叶变换； DFT; FFT; 一、概述 数字信号处理方法的一个重要用途是在离散时间域中确定一个连续时间信号的频谱，通常称为频谱分析，更具体的说，它也包括确定能量谱或功率谱。数字频谱分析可以应用在很广阔的领域。 二、设计目的 1.用DFT 实现对正余弦信号的谱分析； 2.观察DFT 长度和窗函数长度对频谱的影响； 3.对DFT 进行谱分析中的误差现象获得感性认识。 三、设计原理 1、谱分析原理 频谱分析方法是基于以下的观测：如果连续时间信号)(t g a 是频带有限的，那么对其离散时间等效信号)(n g 的DFT 进行谱分析。它的离散时间等效物g(n)应当能给出a g (t)频谱的一个很近似的估计两者之间只差一个带数因子T 。然而，在大多数情况下，)(t g a 是在∞<<∞-t 范围内定义的，因此)(n g 也就定义在∞<<∞-n 的无线范围内，要估计一个无限长信号的频谱是不可能的。实用的方法是：先让模拟连续信号)(t g a 通过一个抗混叠的模拟滤波器，然后把它采样成一个离散序列)(n g 。假定反混叠滤波器的设计是正确的，则混叠效应可以忽略，又假设A/D 变换器的字长足够长，则A/D 变换中的量化噪声也可忽略。 假定表征正余弦信号的基本参数，如振幅、频率和相位不随时间改变，则此信号的傅立叶变换 )(ωj e G 可以用计算它的DTFT 得到 ∑∞ -∞ =-= n n j j e n g e G ωω )()( （1） 实际上无限长序列)(n g 首先乘以一个长度为M 的窗函数)(n w ，使它变成一个长为M 的有限长序列， )()()(1n w n g n g =，对)(1n g 求出的DTFT )(1ωj e G 应该可以作为原连续模拟信号)(t g a 的频谱估计， 然后求出)(1ω j e G 在πω20≤≤区间等分为N 点的离散傅立叶变换DFT 。为保证足够的分辨率，DFT 的 长度N 选的比窗长度M 大，其方法是在截断了的序列后面补上N －M 个零。计算采用FFT 算法。 更详细地考察一下上面的方法。这样才能了解它的限制，并正确利用它所得出的结果。特别要分析加窗的效果，以及和由DFT 样本来估计DTFT 频率采样值的问题。 在讨论由)(1k G 来估计频谱)(1jw e G 和)(jw e G 时，需要重新探讨一下这些变换和它们所对应的频

参数法功率谱估计

参数法功率谱估计 一、 信号的产生 （一）信号组成 在本实验中，需要事先产生待估计的信号，为了使实验结果较为明显，我产生了由两个不同频率的正弦信号（频率差相对较大）和加性高斯白噪声组成的信号。 （二）程序 N=1024;n=0:N-1; xn=2*cos(2*pi*0.2*n)+ cos(2*pi*0.213*n)+randn(1,1024); 这样就产生了加有白噪声的两个正弦信号 其波形如下 0100200300400500600 -8 -6-4-202468 10(a) 两个正弦信号与白噪声叠加的时域波形

二、参数模型法功率谱估计 （一）算法原理简介 1．参数模型法是现代谱估计的主要内容，思路如下： ① 假定所研究的过程)(n x 是由一个白噪声序列)(n ω激励一个因果稳定的可逆线性系统)(z H 的输出； ② 由已知的)(n x ，或其自相关函数)(m r x 估计)(z H 的参数； ③ 由)(z H 的参数来估计)(n x 的功率谱。 2．自回归模型，简称AR 模型，它是一个全极点的模型。“自回归”的含义是：该模型现在的输出是现在的输入和过去p 个输出的加权和。此模型可以表现为以下三式： ① ∑=+--=p k k n u k n x a n x 1)()()(； ② ∑=-+== p k k k z a z A z H 111) (1 )(； ③ 2 12 1)(∑=-+= p k jwk k jw x e a e P σ。 3．AR 模型的正则方程建立了参数k a 和)(n x 的自相关函数的关系，公式如下： =)(m r x ∑=--p k x k k m r a 1 )( 1≥m 时，=)(m r x 21 )(σ+-∑=k r a p k x k 0=m 时。

反应谱和傅里叶谱+地震波选取

【拓展知识1-2】功率谱，反应谱和傅里叶谱，地震波选取，地震持续时间确定功率谱 功率谱是功率谱密度函数的简称。对于一般情况的随机振动，其时间历程具有明显的非周期性，具有连续的多种频率成分，每种频率有对应的功率或能量，用图像来表示这种关系，称为功率在频率域内的函数，简称功率谱密度。加速度功率谱是对地震动加速度时程进行快速傅里叶变换（FFT）得到的[1]。 对于非平稳随机过程，功率谱密度的单位是G的平方/频率。G指的是随机过程。 对于加速度功率谱，加速度的单位是m/s2，则功率谱密度的单位是（m/s2）2/Hz，Hz的单位是1/s，故加速度功率谱密度的单位为m2/s3。加速度功率谱密度函数曲线下方的面积代表随机加速度的总方差，即加速度功率谱可以理解为“随机加速度方差的密度分度”。 参考文献 [1] 庄表中. 随机振动入门.科学出版社，1981. 反应谱和傅里叶谱 反应谱（earthquake response spectrum），是单自由度弹性系统对于某个实际地震加速度的最大反应（可以是加速度、速度和位移）和体系的自振特征（自振周期或频率和阻尼比）之间的函数关系。 反应谱是地震工程中分析结构和设备在地震中的性能的非常有用的工具，因为许多主要表现为简单的振荡器(也称为单自由度系统)。因此，如果能找出结构的固有频率，那么建筑的峰值响应可以通过从地面响应谱中读取相应频率的值来估计。在地震区域的大多数建筑规范中，这个值构成了计算结构必须抵抗的力的基础(地震分析)。如前所述，地面响应谱是在地球自由表面所做的响应图。如果建筑物的响应与地面运动(共振)的组成部分“协调”，可能会发生重大的地震破坏，这些成分可以从响应谱中识别出来。 傅里叶谱，全称为傅里叶振幅谱。地震波是在时间上连续的随机过程，地震

通信原理课程设计报告2

￥ 课程设计报告? < 课程名称通信原理 设计题目 DSB与2ASK调制与解调 专业通信工程 班级 学号 姓名 完成日期 …

课程设计任务书 设计题目：DSB与2ASK调制与解调 设计内容与要求： 设计内容： 1.根据DSB的调制原理设计线路，进行仿真模拟调制DSB的调制和解调过程，并通过仿真软件观察信号以及的调制过程中信号波形和频谱的变化。 2. 根据ASK的调制原理设计线路，进行仿真模拟调制DSB的调制和解调过程，并通过仿真软件观察信号以及的调制过程中信号波形和频谱的变化。 3.在设计过程中分析信号变化的过程和思考仿真过程的设计原理。 ; 设计要求： 1.独立完成DSB与ASK的调制与解调； 2.运用仿真软件设计出DSB与ASK的调制线路 3.分析信号波形和频谱 指导教师：范文 2012年12月16日 课程设计评语 （ 成绩： 指导教师：_______________

年月日

一．调制原理： 调制: 将各种数字基带信号转换成适于信道传输的数字调制信号(已调信号或频带信号)； 时域定义：调制就是用基带信号去控制载波信号的某个或几个参量的变化，将信息荷载在其上形成已调信号传输，而解调是调制的反过程，通过具体的方法从已调信号的参量变化中将恢复原始的基带信号。 频域定义：调制就是将基带信号的频谱搬移到信道通带中或者其中的某个频段上的过程，而解调是将信道中来的频带信号恢复为基带信号的反过程. 根据所控制的信号参量的不同，调制可分为： 调幅，使载波的幅度随着调制信号的大小变化而变化的调制方式。 调频，使载波的瞬时频率随着调制信号的大小而变，而幅度保持不变的调制方式。 调相，利用原始信号控制载波信号的相位。 调制的目的是把要传输的模拟信号或数字信号变换成适合信道传输的信号，这就意味着把基带信号（信源）转变为一个相对基带频率而言频率非常高的代通信号。该信号称为已调信号，而基带信号称为调制信号。调制可以通过使高频载波随信号幅度的变化而改变载波的幅度、相位或者频率来实现。调制过程用于通信系统的发端。在接收端需将已调信号还原成要传输的原始信号，也就是将基带信号从载波中提取出来以便预定的接受者（信宿）处理和理解的过程。该过程称为解调。

通信原理课程设计

通信原理课程设计

通信原理课程设计 专业 学号 学生姓名 指导教师

完成日期 1 月 5 日 一、课程设计目的 在本课程设计中使用的软件工具是MATLAB。目的是希望在以下几方面有所收获： 1.会MATLAB软件的最基本运用。 MATLAB是一种很实用的数学软件，它易学易用。MATLAB对于许多的通信仿真类问题来说是比较合适的。 2.了解计算机仿真的基本原理及方法，知道怎样经过仿真的方 法去研究通信问题。 3.加深对信号与系统和通信原理及其相关课程内容的理解。 二软件实现特点 与硬件实验相比，软件实验具如下一些特点： 1.软件实验具有广泛的实用性和极好的灵活性。在硬件实验中改变系统参数可能意味着要重做硬件，而在软件实验中这只是

该一两个数据，或者只是在屏幕上按几下鼠标。 2.软件实验更有助于我们较为全面地研究通信系统。有许多问题，经过硬件试验来研究可能非常困难，但在软件实验中却易于解决。 3.硬件实验的精确度取决于元器件及工艺水平，软件实现的精确度取决于CPU的运算速度或者说是程序的运算量。 4.软件实验开发周期短，成本低。 三基本要求 1掌握matlab的基本操作及了解基本的仿真方法。 2按以下要求编制仿真程序并调试运行 （1）基本信号的仿真 （2）数字基带传输码型的仿真 （3）调制解调系统设计及仿真 （4）数字基带系统设计及仿真 四课程设计内容 1、编程实现基本信号的仿真

（1）产生并绘出以下信号 单位阶跃序列 k= -30:30; uk=[zeros(1,30),ones(1,31)]; stem(k,uk) 图1 单位阶跃序列图 周期方波square() t=-2*pi/100:pi/1024:2*pi/100; y=square(2*pi*30*t,50); plot(t,y); ylim([-1.5 1.5])

功率谱密度估计方法的MATLAB实现

功率谱密度估计方法的MATLAB实现 在应用数学和物理学中，谱密度、功率谱密度和能量谱密度是一个用于信号的通用概念，它表示每赫兹的功率、每赫兹的能量这样的物理量纲。在物理学中，信号通常是波的形式，例如电磁波、随机振动或者声波。当波的频谱密度乘以一个适当的系数后将得到每单位频率波携带的功率，这被称为信号的功率谱密度（power spectral density, PSD）或者谱功率分布（spectral power distribution, SPD）。功率谱密度的单位通常用每赫兹的瓦特数（W/Hz）表示，或者使用波长而不是频率，即每纳米的瓦特数（W/nm）来表示。信号的功率谱密度当且仅当信号是广义的平稳过程的时候才存在。如果信号不是平稳过程，那么自相关函数一定是两个变量的函数，这样就不存在功率谱密度，但是可以使用类似的技术估计时变谱密度。信号功率谱的概念和应用是电子工程的基础，尤其是在电子通信系统中，例如无线电和微波通信、雷达以及相关系统。因此学习如何进行功率谱密度估计十分重要，借助于Matlab工具可以实现各种谱估计方法的模拟仿真并输出结果。下面对周期图法、修正周期图法、最大熵法、Levinson递推法和Burg法的功率谱密度估计方法进行程序设计及仿真并给出仿真结果。 以下程序运行平台：Matlab R2015a（8.5.0.197613） 一、周期图法谱估计程序 1、源程序 Fs=100000; %采样频率100kHz N=1024; %数据长度N=1024 n=0:N-1; t=n/Fs; xn=sin(2000*2*pi*t); %正弦波,f=2000Hz Y=awgn(xn,10); %加入信噪比为10db的高斯白噪声 subplot(2,1,1); plot(n,Y) title('信号') xlabel('时间');ylabel('幅度');

高层建筑地震作用计算的时域显式随机模拟法

文章编号 5777=4>4?!6758"75=7753=57!"##57<5@774A B587 华南理工大学亚热带建筑科学国家重点实验室项目 6753e I75 作者简介 苏成 5?4> ! 男 广东潮阳人 工学博士 教授 X=O'+( 0K0J2$^20$/:1,$:0& 收稿日期 675@年3月 35

Matlab课程设计正文部分

1 Matlab概述 1.1 简介 MATLAB是矩阵实验室（Matrix Laboratory）的简称，是美国MathWorks公司出品的商业数学软件，用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境，主要包括MATLAB和Simulink两大部分。 1.2 主要功能 MATLAB是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵运算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如C、Fortran）的编辑模式，代表了当今国际科学计算软件的先进水平。 MATLAB的基本数据单位是矩阵，它的指令表达式与数学、工程中常用的形式十分相似，故用MATLAB来解算问题要比用C，FORTRAN等语言完成相同的事情简捷得多，并且MATLAB 也吸收了像Maple等软件的优点，使MATLAB成为一个强大的数学软件。在新的版本中也加入了对C，TORTRAN，C++，JAVA的支持。可以直接调用,用户也可以将自己编写的实用程序导入到MATLAB函数库中方便自己以后调用。 1.3 应用 MATLAB 产品族可以用来进行以下各种工作： ●数值分析 ●数值和符号计算 ●工程与科学绘图 ●控制系统的设计与仿真 ●数字图像处理技术 ●数字信号处理技术 ●通讯系统设计与仿真

●财务与金融工程 MATLAB 的应用范围非常广，包括信号和图像处理、通讯、控制系统设计、测试和测量、财务建模和分析以及计算生物学等众多应用领域。附加的工具箱（单独提供的专用MATLAB 函数集）扩展了MATLAB 环境，以解决这些应用领域内特定类型的问题。 2 试题（第一套） 1．计算y1= 5 1)3.0sin(2+π和y2= 5 1)3.0cos(2+π； 2．画出衰减震荡曲线t e y t 3sin )3/(-=及其他的包络线)3/(0t e y -=,t 的取值范围是[0,4π]. 3．画出2 2 2 2sin y x y x z ++= 所表示的三维曲面。x,y 的取值范围是[-8，8]。 4．分析下面每条指令的功能并运行，观察执行结果。 (1) X=0:0.1:1; Y=X.*exp(-X); plot(X,Y),xlabel(‘x’), ylabel(‘y’),title(‘y=x*exp(-x)’); (2) A=zeros(2,5) A(:)=-4:5 L=abs(A)>3 islogical(L) X=A(L) (3) A=[1:4;5:8] pow2(A) (4) A=zeros(2,3) A(:)=1:6 A=A*(1+i) A1=A.’; B1=A’; (5) A=ones(1,2)

反应谱理论与人工模拟地震波技术简介

第３３卷第２６期?１０６?２００７年９月山西建筑 ＳＨＡＮＸＩＡＲＣＨＩＴＥＣＴＵＲＥ Ｖｄ３３Ｎｏ．２６ Ｓｅｐ．２００７ 文章编号：１００９—６８２５｛２００７）２６—０１０６—０３ 反应谱理论与人工模拟地震波技术简介 邱玉国王玉富 摘要：介绍了反应谱理论的发展历程和国内外研究现状，分析了研究问题的思路，指出了利用反应谱理论来解决实际工程时遇到的问题，并简单介绍了国外对人工模拟地震波技术的应用和研究，为抗震理论提供了参考依据。 关键词：反应谱理论，地震波，随机振动，非弹性地震波 中图分类号：ＴＵ３５２文献标识码：Ａ １概述 反应谱理论是建筑结构抗震设计的重要理论基础之一。从２０世纪５０年代开始，反应谱理论逐渐成为结构抗震设计的重要方法，经过５０多年的发展，目前这种方法已经为世界上大多数国家的设计规范所采用。但是，由于地震产生机理和作用效果的复杂性，采用反应谱理论进行分析和设计与工程实践还存在很多与实际不相符合之处。此外，对于反应地震重要特性的时间问题，反应谱法也无能为力。 人工模拟地震波技术是近年来才发展起来的一项新的结构抗震设计的技术手段，目前主要用于计算机模拟和特别重要结构模型的振动台试验。它能够通过模拟地震波的特性来用于对结构进行时程分析，是～种新兴的、具有革命性意义的试验手段。 图２数值模拟结果２．３计算结果分析 通过数值模拟和试验得到瓦斯管承载力等数值如表２所示。 表２数值模拟和试验结果 Ｉ研究方法承载力仆但ａ最大应变／％最大剪应力／ＳＰａＩ数值模拟７．１４Ｏ．０８４２１６０室内试验６．６２０．０９６４ ３结语 通过对丁集煤矿瓦斯管材质和整体抗外压的试验研究以及数值模拟分析，可以获得如下重要结论： １）通过对管材材质的试验研究表明：工作管材质采用Ｑ３４５，尺寸为柘３０ｒｆｌｌＴｌ×１４ｉｎｌｎ，能够满足强度和稳定性要求。 ２）瓦斯管整体抗外压试验结果表明：工作管抗外压承载力为６，６２ＭＰａ；通过大变形有限元数值计算，采用变形稳定性控制其承载力，结果为７．１４ＭＰａ，两者数值十分接近，说明用文中方法模拟大直径瓦斯管的承载力是可行的。 参考文献： ［１］李正来．瓦斯抽排钻孔定向技术的改进［Ｊ］．安徽科技，２００６（３）：４９—５０． ［２］汪东生．瓦斯抽排技术治理本煤层采空区瓦斯涌出的实践［Ｊ］．煤矿安全，２００６（１）：１３—１５． ［３］张敦伍，任胜杰．瓦斯抽排钻孔防偏斜实践［Ｊ］．矿业安全与环保，２００５（８）：６７—６８． ［４］刘克功，范再良，赵新华．采空区瓦斯抽排法治理综放面瓦斯超限［Ｊ］．煤，１９９８（２）：４８—５０． Ｓｔｕｄｙｉｎｇｏｎｒａｄｉａｌｓｔａｂｉｌｉｔｙｎｕｍｅｒｉｃａｌｓｉｍｕｌａｔｉｏｎｏｆｌａｒｇｅｐｉｐｅｉｎｍｉｎｅ ＴＯＮＧＷｅｎ－ｌｉｎ Ａｂｓｔｒａｃｔ：ＴｈｅｅｘｐｅｒｉｍｅｎｔａｌａｎｄｖａｌｕｅｓｉｍｕｌａｔｉｏｎｍｅｔｈｏｄｓｈａｖｅｓｔｕｄｉｅｄｔｈｅＤｉｎｇｉｉｃｏａｌｍｉｎｅｌａｒｇｅｄｉａｍｅｔｅｒｇａｓｔｕｂｅｕｎｄｅｒｍｅｃｈａｎｉｃｓｃｈａｒａｃｔｅｒ—ｉｓｔｉｅ．Ｒｅｓｕｌｔｉｎｄｉｃａｔｅｄ：ｔｈｅｌａｒｇｅｄｉａｍｅｔｅｒｇａｓｔｕｂｅｉｓｐｒｅｓｅｎｔｅｄｓｔａｂｉｌｉｔｙｆａｉｌｕｒｅｍｏｄｅｌｉｎｅｎｃｉｒｃｌｅｓｐｒｅｓｓｅｓｓｈａｐｅ，ｉｔｓｓａｆｅｔｙｆａｃｔｏｒｒｅａｃｈｅｓ３．０，ｉｔｉｓｄｅｓｉｇｎｔｈｅｌａｒｇｅｄｉａｍｅｔｅｒｇａｓｔｕｂｅａｎｄｔｈｅｃｏｎｓｔｒｕｃｔｐｍｖｉｄｅｓｔｈｅｒｅｆｅｒｅｎｃｅ． Ｋｅｙｗｏｒｄｓ：ｌａｒｇｅｄｉａｍｅｔｅｒｇａｓｔｕｂｅ，ｅｘｐｅｒｉｍｅｎｔａｌｉｎｌａｂ，ｎｕｍｅｒｉｃａｌｓｉｍｕｌａｔｉｏｎ，ｓｔａｂｉｌｉｔｙｆａｉｌｕｒｅｍｏｄｅｌ 收稿日期：２００７．０４．０６ 作者简介：邱玉国（１９７３。），男，工程师，辽宁工程技术大学软件学院，辽宁阜新１２３０００ 王玉富（１９７０．），男，工程师，中铁十九局集团第三工程有限公司，辽宁辽阳１１１０００