数据处理的基础知识.ppt

强震动观测数据处理分析的 基础知识 周雍年 2013.11 1 傅里叶变换 2 连续函数的数字化 3 强震动观测记录常规处理 2 连续函数的数字化 2.1 数字化方法与误差 2.2 采样定理 2.3 窗函数 2.4 有限记录长度对频谱的影响 2.5 混淆问题 3.5 反应谱计算 3.5.1 反应谱的基本概念 2.2 采样定理 数字化后得到的一系列离散数据能在多大程度上代表原来的连续函数？ 采样定理 设函数 的傅里叶变换为 ，且 满足 则函数 可唯一地由一系列等间距采样数值 确定，即 离散采样序列 与连续函数 有相同的傅里叶谱，即 的最高频率 由下式确定 通常称为奈奎斯特(Nyquist)频率 采样间隔 (采样步长) 想从记录中获取的最高频率是 , 理论上采样步长应取 实际采 用 采样步长决定了可以从记录中获取的最高频率成分 2.3 窗函数 用有限长度矩形窗截断后的函数的频谱与原无限长函数的频谱是否完全一样？ 矩形窗函数可以定义为 傅里叶谱 采用矩形窗后的傅里叶频谱已不等于原来的真实傅里叶频谱，而是真实频谱在每一点以 为权重的平均 的主瓣宽1/T，同时，因为除了主瓣外还有一系列边瓣，每个边瓣宽度为1/2T，这种边瓣带进了虚假的高频。边瓣中还有负的边瓣，这就造成了能量的损失。这种现象称为能量的漏泄 三角形窗 汉宁窗 时间窗为 谱窗为 2.4 有限记录长度对频谱的影响 设有一个无限长的记录 ，它的真实频谱为 。这个函数也可以看作 相乘的结果 是和无限长矩形窗 当矩形窗无限长时,截得的记录的谱是它的真实谱 当矩形窗为有限长 时，截得的记录的谱已不再是它的真实谱 矩形窗为有限长 时，它的谱窗 是一个较宽的频带， 并带有漏泄 。 越长，这个谱窗越窄，计算得到的频谱 越接近真实谱 一个强震动记录的长度 T 确定了频谱的分辨率 T越长，频谱分辨得越细，越接近真实频谱。 对长度为T的记录，其频谱中的频率间隔为1/T 最小的可分辨频率也为1/T 矩形窗的主瓣宽度为1/T。某一点的频率分量是在这一频率左右 之间光滑化的结果，这决定了可以分辨的频率间隔为 由此也可以推断，频谱中最小的可分辨频率为 2.5 混淆问题 根据采样定理，如果信号的傅里叶变换满足 当 只要采样步长取到 ， 就可唯一地由一系列采样值 确定 而采样值 的谱 与原连续记录的谱在 范围内是相同的， 即 时 在 范围内谱 与连续记录的谱 是否还是相同呢？ 中不仅有 时的 信息, 还包含了 在 时的高频信息。这种高频信息叠加到 的离散低频谱上的现象 就称为混淆或高频混淆。 如果 fc 2fc 3fc 4fc 5fc 6fc fc 2fc 4fc 6fc 3fc 5fc 0 0 (a) (b) 减小混淆现象的途径: 尽量小,以提高奈奎斯特频率 使余下的 (2) 采用不等步长采样，减少混淆现象。 对数字强震动仪来说，是不现实的 (3) 在谱分析前先进行数字滤波，滤去不需要的高频成分。 高频分量尽可能小 (1) 加大采样率,使采样步长 3 强震动观测记录常规处理 3.1 强震动观测记录的一般特征 3.2 强震动观测记录的常规处理分析 3.3 强震动加速度记录误差分析 3.4 强震动记录的校正处理 3.5 反应谱计算 3.6 傅里叶谱计算 3.1 强震动观测记录的一般特征 最大加速度值、主要频率、 持续时间 3.2 强震动观测记录的常规处理分析 常规处理内容 零基线调整 误差校正 计算校正加速度时程、速度和位移时程 反应谱计算 傅里叶谱分析 （1）零基线调整 模拟强震动仪记录 数字强震动仪记录 （2）仪器误差校正