地震波观测系统的MATLAB仿真报告.docVIP

地震波观测系统的MATLAB仿真 课程名称 数字信号处理 实验项目 题目6 地震波观测系统的MATLAB仿真 指导教师 赵双琦 学 院 光电信息与通信工程 _ 专 业 电子信息工程 班级/学号 学生姓名 课设时间 2011-12- 28至2012-1-5 09级“数字信号处理课程设计”任务书 题目主要 内容 设计 要求 主要仪器设备 1、计算机1台，安装 主要参考文献 [M].西安：西安交通大学出版社，2002. 课程设计进度计划（起止时间、工作内容） 本课程设计共安排，这是其中题目之一。整个课程设计共学时，分周安排，具体进度如下： ； 用语言实现 4学时 进一步完善功能，现场检查、答辩； 完成课程设计报告。 课程设计开始日期 20.12.26 课程设计完成日期 20.1.6 课程设计实验室名称 实验室 地 点 实验楼03、605 资料下载地址 数字信号处理课程设计 目录 摘要 - 4 - 正文 - 4 - 一、目的 - 4 - 二、原理 - 4 - 三、要求 - 5 - 四、步骤 - 5 - 五、程序实现 - 6 - 实验结果 - 12 - 六、体会 - 15 - 参考文献 - 15 - 摘要 本文的目的是实现地震波观测系统的MATLAB仿真。一个线性系统y(t)=h(t)*x(t)，x(t)为地面运动，h(t)为系统的冲击响应，y(t)为系统输出。根据卷积定理，有Y（ω）=H(ω)X(ω)。由地震波观测文件数据y（t），再设计一个宽频带滤波器h（t），就可以恢复地面运动x（t）。对于短周期地震仪，其系统函数为H1（w），对于输入地面运动x（t），有Y1(ω)=H1(ω)X（ω），我们可以推导出Y1（w）=H1(w)Y(w)/H(w)，再对Y1（w）作ifft就可以实现宽频带仪器到短周期窄带仪器的仿真。同样，对长周期地震仪，其系统函数为H2（w），我们也可以得到Y2（w）=H1(w)Y(w)/H(w)，然后对Y2（w）作ifft实现仿真。椭圆滤波器、巴特沃斯滤波器和切比雪夫滤波器的设计都很简单，只要滤波器的指标没问题，调用相应的函数就能实现。仿真的结果请参考本文的正文部分。 正文 一、目的 运用所学数字信号处理的基本知识，掌握地震波观测系统的数字信号处理方法。实现宽频带系统的输出仿真到窄频带输出及地面运动恢复。 二、原理 对于一个线性系统，可以用它的系统函数或脉冲响应来表示 y(t)=h(t)*x(t) ① 式中，x(t)为输入信号，相当于地震观测系统的地面运动；y(t)为系统的输出，相当于地震观测系统的地震记录；h(t)为系统的冲击响应。在频率域内，根据卷积定理，该式可以表示为 Y（ω）=H(ω)X(ω) ② 式中，H(ω)为系统的传递函数， X(ω)、 Y（ω）为x(t)、y(t)的傅里叶变换。 设想一个频带范围很宽的线性系统，如宽频带地震仪，其系统函数为H(ω);另一个频带较窄的系统，如短周期地震仪，其系统函数为H1(ω)，对于同样的输入X(ω)有 Y(ω)=H(ω)X(ω), Y1(ω)=H1(ω)X（ω） ③ 式中，Y1(ω)为频带较窄的系统记录的频谱；H1(ω为频带较窄系统的传递函数。由式③可得 H1(ω) Y(ω) Y1(ω)= H(ω) ④ 将上式变换到时间域就得到频带较窄系统的输出y1(t)。也就是说，如果知道宽频带和窄频带系统的传递函数H(ω)和 H1(ω)，原则上可以从宽频带系统的输出推测出窄频带系统的输出。但如果我们知道窄频带系统输出及其两种系统的传递函数，却无法得到宽频带系统的输出。这样就使得我们在记录某种信号时采用宽频带记录，然后仿真到各种窄频带的记录仪器上对信号进行分析。 如果已知地震仪的输出和地震仪的传递函数，我们可以求出地面运动为 X（ω）= Y(ω)/ H(ω) ⑤ 三、要求 以某地震台站记录的地震观测文件为例，选择合适滤波器揭示地面运动恢复和仿真的概念