信号分析处理课程设计报告.docVIP

信 号 分 析 与 处 理 课 程 设 计 报 告 学院：油气资源学院 班级：物探0901 学号：5 姓名：杨继东 指导老师：宁忠华老师 本次课程设计用了一周时间，从十八周周二上午到十八周周五上午，实际上机操作为两个上午，在西安石油大学校本部第一实验楼四楼机房，由宁忠华老师带领、指导完成。这次课程设计的主要目的是巩固我们所学的信号分析与处理基本理论知识，掌握用计算机对信号进行采集、处理基本方法。所应用的主要软件为数学编程软件MATLAB7.0.，我们这次重要应用其对基本信号的画图、滤波、快速傅里叶变换和快速傅里叶反变换及形成“.dat”格式文件的功能。 这次课程设计的主要内容有熟悉MATLAB编程环境，掌握基本的变量定义、赋值、基本简谐波函数调用、画图、傅里叶函数fft( )和ifft( )的调用、求绝对值或模或振幅函数abs( )、取整函数int16( )的编程语句。掌握简单数字滤波器（低通滤波器、高通滤波器、带通滤波器）设计方法，学会输出基本简谐波的时域、频域波形图的方法和步骤及简谐波叠加后复杂波滤波后时域、频域波形图的方法和步骤。学会利用所编写的滤波器程序对实际地震信号进行处理并且学会分析滤波前后地震波波形的变化，更深一步了解滤波器的实际物理意义以便以后更好地进行地震数据资料的处理和解释。并且通过这次课程设计了解地震数据基本的记录方式，即要了解单炮和多炮地震数据不同特点以及距震源不同距离每道地震记录的不同特征。这次课程设计的具体内容主要有以下四点，主要图形有十四幅。 了解Matlab软件特点，熟悉Matlab编程环境。………………… 2 对单频简谐波和叠加波的时域频域波形的输出（4幅图）。……3 对叠加波进行低通高通滤波处理以及时域频域波形的输出（4幅图）。 ……………………………………………………………………… 5 利用带通滤波器对实际的地震数据进行处理及对其滤波后特点的分析（8幅图）。…………………………………………………………7 五、本次课程设计所遇到的问题。………………………………………11 详细内容在下面一一详细介绍 一、了解Matlab软件特点，熟悉Matlab编程环境 ㈠ 本部分的目的 本部分是这次课程设计的基础，要求我们掌握MATLAB7.0中基本的编程语句，了解编程的环境以及运行的方法，学会对所编写程序及运行后成果图的保存和调用。 ㈡ 本部分的任务及内容 这部分主要涉及编程的基本语句，具体如下： 1. 定义，赋值： 给变量赋值： f1=20; 给数组赋值：i=1:N; 或tt=0:deltT:(N-1)*deltT; 2. 函数调用： ①sin( )函数或者cos( )函数的调用；x1(i)=sin(tmp*f1*i);x2(i)= cos(tmp*f2*i); ②画图函数plot( )的调用：figure(1); figure(2); plot(x1);plot(tt,x2); ③title(原始信号的时间序列);xlabel(时间(ms));ylabel(振幅); ④傅里叶函数fft( )和ifft( )的调用：XX=fft(xx,N);yy=ifft(XX,N); ⑤求绝对值或模或振幅函数abs( ) ：XXabs=abs(XX); ⑥取复数中的实部X=real(Z ); 取复数中的虚部Y = imag(Z); 3. 画频谱图：横坐标为频率；纵坐标为振幅。 4. 滤波：XX(1:k)=0.0; 5. 取整函数int16( )：lp=int16(lf/f0); 6. 其它： ①注释：%============人工合成纪录 ②每个语句后面用分号结束“；” ③π应输入pi ④抓图：Save或Alt+PrScrn 二 、对单频简谐波和叠加波的时域频域波形的输出 ㈠ 本部分的目的 这部分主要要求我们会利用数学编程软件MATLAB7.0编写相应的程序对单频简谐波和单频简谐波叠加后复杂波输出时域和频域的波形图，了解基本简谐波的波形特点以及叠加后波形的变形特点。这部分主要完成四幅图，即单频谐波时域图、单频谐波振幅谱图、双频信号时域图、双频信号振幅谱图。 ㈡ 任务及方法实现 对于单频简谐波我所选的是：x1(i)=cos(2*pi*70*i*0.001)(间隔为0.001s 单频波频率为70hz)。 对于叠加波我所选的是：x1(i)=cos(2*pi*70*i*0.001)+ cos(2*pi*120*i*0.001)(间隔全为0.001s 低频波频率f1为70hz 高频波频率f2为120hz)。 对于时域信号的两张图主要是利用plot（）输出的，而频域的谱线图先经过fft（）实现傅里叶变换再经过取绝对值abs（）和画图plot（）得出。 ㈢ 结果及分析 这部分主要由四张成果图，详细如下： 单频简