信息分析与处理

1、信息分析与处理第1页，共60页，2022年，5月20日，16点32分，星期一Filtering奥姆斯比 滤波器第2页，共60页，2022年，5月20日，16点32分，星期一Line 127 With Spiking Decon and 10-80 Bandpass Filter第3页，共60页，2022年，5月20日，16点32分，星期一第4页，共60页，2022年，5月20日，16点32分，星期一第5页，共60页，2022年，5月20日，16点32分，星期一第6页，共60页，2022年，5月20日，16点32分，星期一第7页，共60页，2022年，5月20日，16点32分，星期一第8页，共6

2、0页，2022年，5月20日，16点32分，星期一第9页，共60页，2022年，5月20日，16点32分，星期一第10页，共60页，2022年，5月20日，16点32分，星期一第11页，共60页，2022年，5月20日，16点32分，星期一第12页，共60页，2022年，5月20日，16点32分，星期一第13页，共60页，2022年，5月20日，16点32分，星期一第14页，共60页，2022年，5月20日，16点32分，星期一第15页，共60页，2022年，5月20日，16点32分，星期一第16页，共60页，2022年，5月20日，16点32分，星期一第17页，共60页，2022年，5月20

3、日，16点32分，星期一第18页，共60页，2022年，5月20日，16点32分，星期一第19页，共60页，2022年，5月20日，16点32分，星期一第20页，共60页，2022年，5月20日，16点32分，星期一第21页，共60页，2022年，5月20日，16点32分，星期一第22页，共60页，2022年，5月20日，16点32分，星期一第23页，共60页，2022年，5月20日，16点32分，星期一第24页，共60页，2022年，5月20日，16点32分，星期一第25页，共60页，2022年，5月20日，16点32分，星期一第26页，共60页，2022年，5月20日，16点32分，星期一

4、第27页，共60页，2022年，5月20日，16点32分，星期一第28页，共60页，2022年，5月20日，16点32分，星期一第29页，共60页，2022年，5月20日，16点32分，星期一第30页，共60页，2022年，5月20日，16点32分，星期一第31页，共60页，2022年，5月20日，16点32分，星期一第32页，共60页，2022年，5月20日，16点32分，星期一第33页，共60页，2022年，5月20日，16点32分，星期一第34页，共60页，2022年，5月20日，16点32分，星期一第35页，共60页，2022年，5月20日，16点32分，星期一T1Ap (dB)As

5、(dB)0.50.3854300.40.6467400.30.9059360.380.698744加一个样本值T1T1T2Ap(dB)As(dB)0.590.110.3176650.600.110.2918590.580.110.343359加两个样本值T1和T2H(k)=DFTh(n)第36页，共60页，2022年，5月20日，16点32分，星期一第37页，共60页，2022年，5月20日，16点32分，星期一第38页，共60页，2022年，5月20日，16点32分，星期一第39页，共60页，2022年，5月20日，16点32分，星期一第40页，共60页，2022年，5月20日，16点32分

6、，星期一第41页，共60页，2022年，5月20日，16点32分，星期一Hrs=zeros(1,7),T1,T2,ones(1,14),T2,T1,zeros(1,6)k1=0:floor(N-1)/2)k2=floor(N+1)/2):N-1angH=-LL*(2*pi)/N*k1, LL*(2*pi)/N*(N-k2);Hdk=Hrs.*exp(j*angH);h=real(ifft(Hdk,N);K=067 8915162223 242530H= 7*0T1 T27*1 7*1T2 T1 6*1第42页，共60页，2022年，5月20日，16点32分，星期一第43页，共60页，2022年

7、，5月20日，16点32分，星期一K=056713141519202125262733343540H=6*0T17*1T1 5*0 0 5*0T1 7*1T1 6*0Hrs=zeros(1,6),T1,ones(1,7),T1,zeros(1,11),T1,ones(1,7),T1,zeros(1,6)k1=0:floor(LL)k2=floor(LL)+1:N-1angH=pi/2-LL*(2*pi)/N*(k1+0.5),-pi/2+LL*(2*pi)/N*(N-k2-0.5);Hdk=Hrs.*exp(j*angH);h1=ifft(Hdk,N);n=0:1:N-1;h=real(h1.

8、*exp(j*pi*n/N);第44页，共60页，2022年，5月20日，16点32分，星期一第45页，共60页，2022年，5月20日，16点32分，星期一第46页，共60页，2022年，5月20日，16点32分，星期一第47页，共60页，2022年，5月20日，16点32分，星期一第48页，共60页，2022年，5月20日，16点32分，星期一白噪声是均值为零，自相关函数含有冲激函数，功率谱为常数的随机信号。第49页，共60页，2022年，5月20日，16点32分，星期一第50页，共60页，2022年，5月20日，16点32分，星期一b=10,3,0,0; 3,10,3,0; 0,3,10

9、,3; 0,0,3,10;c=3.25;0.85;0.31;0.03;b1=inv(b);a3=b1*cx=3,1;y=conv(a3,x)第51页，共60页，2022年，5月20日，16点32分，星期一第52页，共60页，2022年，5月20日，16点32分，星期一第53页，共60页，2022年，5月20日，16点32分，星期一第54页，共60页，2022年，5月20日，16点32分，星期一第55页，共60页，2022年，5月20日，16点32分，星期一实例:井震反射系数分析比较井深: 5536.255642.29m深度厚度: 106.04m双程时: 3534.73579.9ms单程时:1767.361789.98时间厚度: 45.2ms (双程时) 22.6ms (单程时)rtv第56页，共60页，2022年，5月20日，16点32分，星期一图1 期望输出为1,0,.0第57页，共60页，2022年，5月20日，16点32分，星期一图