matlab ellip传递函数,地震波观测系统的MATLAB仿真报告.doc

您所在位置：网站首页 > 海量文档

&nbsp>&nbsp建筑/施工&nbsp>&nbsp岩土与地质勘察

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

本文档一共被下载：次,您可全文免费在线阅读后下载本文档。

下载提示

1.本站不保证该用户上传的文档完整性，不预览、不比对内容而直接下载产生的反悔问题本站不予受理。

2.该文档所得收入(下载+内容+预览三)归上传者、原创者。

3.登录后可充值，立即自动返金币，充值渠道很便利

地震波观测系统的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(ω) ⑤

三、要求

以某地震台站记录的地震观测文件为例，选择合适滤波器揭示地面运动恢复和仿真的概念

发表评论

请自觉遵守互联网相关的政策法规，严禁发布色情、暴力、反动的言论。

用户名:

验证码:

匿名?

发表评论