实验案例七.docVIP

实验案例七 实验序号： 08-7 实验组数及人数：30组/1人 实验名称：Matlab信号处理分析 实验内容： 用Matlab信号分析工具箱，实现对地球物理观测数据的滤波（如高通、低通或带通滤波等） 一、实习要求 根据给定的某种地球物理观测信号，要求学生用Matlab实现观测数据的导入，并绘制相应的地球物理观测信号时程曲线。 利用Matlab中的信号处理库函数，实现对输入观测信号的滤波处理，包括高通、低通、带通等。并将滤波后的数据重新绘制。 通过完成上述任务，掌握基本的Matlab编程技能，并加深对地球物理观测数据信号处理的认识。 图7.1为滤波前的地球物理观测信号数据（*.dat），该文件可课程网址上下载。 图7.1 数据（采样间隔：0.1ns） -206 -160 -195 -147 -174 -180 -129 -207 -165 -142 -121 -213 -190 -191 -198 -189 -199 -174 -197 -186 -207 -182 -173 -150 -188 -183 -167 -193 -185 -196 -156 -182 -184 -187 -180 -156 -159 -177 -168 -186 -210 -238 -404 -533 -689 -932 -1729 -3526 -5520 -4813 526 9297 15632 14999 6681 -3872 -15083 -20228 -22323 -23378 -21578 -14496 -1629 10311 17179 15953 8767 1990 793 4597 9312 11208 10496 7958 4641 1961 1473 3297 6202 7452 5859 2040 -1831 -3388 -2642 -911 419 1120 1736 1980 1640 776 -254 -1174 -2195 -2878 -2939 -2214 -1189 -480 -305 -629 -1169 -1806 -2585 -3093 -3110 -2620 -2170 -2241 -2813 -3207 -2748 -1311 565 2054 2607 1942 384 -1505 -3154 -3301 -2287 -780 37 -502 -1802 -2694 -2512 -1384 465 2108 3223 3497 2912 1691 92 -1045 -1576 -1735 -1744 -1752 -1773 -1659 -1239 -601 230 904 1243 1032 537 -117 -632 -752 -772 -539 -378 -215 -114 -20 139 202 290 275 126 -62 -228 -304 -182 -38 100 119 113 86 42 16 3 -70 -59 -89 -202 -259 -279 -303 -250 -217 -170 -166 -82 -155 -141 -171 -86 23 88 94 118 60 -9 -60 -64 -93 -21 24 24 25 -40 -73 -124 -154 -79 -45 -28 -4 2 -86 -116 -232 -238 -253 -270 -282 -304 -309 -313 -289 -285 -311 -311 -304 -326 -331 -299 -287 -282 -241 -239 -216 -236 -221 二、编程提示 1. 数据文件的读取，用load(*.dat)函数载入原始观测数据，并将各采样点的观测值用一矩阵保存。 2. 用plot（data）函数绘制观测信号的时程曲线图。 3. 用FFT(data)函数对原始数据进行谱分析； 4. 设置相应的高通、低通或带通截止频率值，进行频率的切除； 5. 用IFFT（）函数进行逆变换，获得滤波后的数据； 6. 绘制滤波后的时程曲线； 7. 输入数据文件及滤波后的时程曲线图像文件。 三、提交图件 1. Matlab源程序（*.m）； 2. 不同类型滤波后的数据文件（*.dat）； 3. 滤波前后的时程曲线的对比图，并附相关的文件说明。