课程设计：信号分析与处理C语言编程.docxVIP

???????????????????????最新资料推荐??????????????????? 勘查技术课程设计：信号分析与处理基础 （西南石油大学 --- 资源与环境学院） 对于勘查技术与工程专业的学生来说， 《信号分析与处理基础》是一门专业基础课，我是 2010 级的， 我们是在大三第一学期上的，这门课数学与物理知识要求比较高，不过一开认真仔细学的话，也会学的好 的，起码要比那空洞、生奥、蛋疼《弹性波动力学》好学些。 随着课程的结束， 《信号分析与处理基础》的课程设计也随之而来，我们是老师布置了 4 个题目，分 单号与双号各自做 2 道，我是单号，做的是滤波与相关。 这次课程设计， 注意考验大家的编程能力，目前我们学过得就只有 C 语言，可以用 Fortran,Matlab 等 等， Matlab 可以现学现用，上手快。但是大家也可以挑战下自己的 C 语言，提高下自己的编程能力，这是 一次很好地机会，真正实用的时刻。我就是用 C 语言编的。 其余 2 题，我也把程序与结果图收集到了这里，以供学弟、学妹们参考之用！我的 QQ：593066480， 有什么不懂的，或者好的见教，欢迎来信息交流！ 题目如下： 1、滤波 已知原始地震记录 x(t), 要求：设计滤波器，消除 x(t)中 10Hz 以下， 80Hz 以上的干扰信号。 x(t) s(t) n(t ) 2 2 s(t ) A1e t sin(2 f1t ) A2 e t sin(2 f2t ) n(t ) A3[sin( 2 f3t ) cos(2 f 4t)] 建议参数： A 1=1，A 2=0.8， A 3=0.5 f 1=25Hz ， f2=45Hz ， f3=5Hz ， f 4=80 取样点数： N=200 抽样间隔： =0.004，?=100 ○1时域滤波 ： 由 y(t) x(t ) * h(t ) (h(t) 为滤波因子 ) h(t ) 2 cos( 0 t) sin( t) t 建议参数：第一参数： ff1=20Hz ， ff2=50Hz 第二参数： ff1=25Hz ， ff2=45Hz 抽样点数： M=60 抽样间隔： =0.004 12 ff1 ， 2 2 ff 2 2 1 ， 0 21 2 2 2 单独计算： h(0) 要求：画出 x(t),h(t),y(t) 图形，为了分析方便，也可以画出有效波 s(t),干扰波 n(t)及其频 谱进行分析，如下图 : 1 ???????????????????????最新资料推荐??????????????????? 2 ???????????????????????最新资料推荐??????????????????? 3 ???????????????????????最新资料推荐??????????????????? 最后就是答辩，老师问问题，学生回答。主要注意几点就行了 1， 熟悉课本滤波部分知识； 2， 第二参数要比第一参数滤波效果好，因为门第一参数开大了，进来的干扰波也多了，从第一参数： Fy、第二参数 Fy 图形上可以看出来，干扰波频谱被压小了。 第二参数 压制了干扰波，突显了有效波， 所有好。 ○2频域滤波 由公式： Y (k ) X (k ) H (k ) 1， 对 x(t) 进行补 0， 28 或者 29 总之必须是 2 的次方，因为要用到 FFT 公式与 IFFT 公式， 进行 FFT 变换得到 X （ k）； 2， H（ k）的求法： H（ k）也必须与 X （ k）点数相同 1 ， N 单门： m1 f 双门： m1 f  20, m2 f 50 20, m2 f 50, m3 =N- m2 , m4 =N- m1 , 3, Y(k) X (k ) H (k) 在用 IFFT 反变换得 y(t), 存在实部与虚部，需要分析与处理 要求：画出 H(k) 、 X(k) 、Y(k) 图形，并且分析 X(k) 、Y(k) 的区别，还有开单双门的区别与差异，时域与频域滤波谁好、为什么？ 4 ???????????????????????最新资料推荐??????????????????? 5 ???????????????????????最新资料推荐??????????????????? 3， 分析： 从单双门实部图形看出，与有效波是完全一样的，但是幅值变大，这体现了 滤波 突显有效波 的特性；从图形看出虚部对结果没有影响； 单双门虚部完全不同，这是由于开单双门效果不同引起的。单门虚部变化大，而 且幅值与起伏变化也大，而双门幅值很小，小到可以忽略不计。按理说反变化 IFFT 后 应该只有实部，没有虚部，至于为什么会产生虚部，希望读者自己下去研究下，希望大家相