检测原理第二次作业.docxVIP

过程参数检测作业 一、 超声波顺流传播信号采样 ? 实验程序 f=2*10^6;a=10^12;tao=5*10^-6;fs=10^9;T=1/fs;t=0:T:2*tao; %初始化参数 s=sin(2*pi*f*t).*exp(-a*(t-tao).^2/2); %函数关系式 plot(t,s) %按两个数组画图 xlabel(t/s); title(超声波顺流信号采样图) %图表信息 ? 实验结果 图1. 超声波顺流信号采样图 二、 超声波逆流传播信号 ? 实验程序 f=2*10^6;a=10^12;tao=5*10^-6;fs=10^9;T=1/fs;t=0:T:2*tao; %初始化参数 t1=2.05*10^-7;t2=2.10*10^-7;t3=2.15*10^-7; s=sin(2*pi*f*t).*exp(-a*(t-tao).^2/2); %函数关系式 s1=sin(2*pi*f*(t-t1)).*exp(-a*(t-t1-tao).^2/2); %函数关系式 s2=sin(2*pi*f*(t-t2)).*exp(-a*(t-t2-tao).^2/2); %函数关系式 s3=sin(2*pi*f*(t-t3)).*exp(-a*(t-t3-tao).^2/2); %函数关系式 plot(t,s); hold on; plot(t,s1,y); plot(t,s2,c); plot(t,s3,r); xlabel(t/s); title(超声波逆流信号采样图) %图表信息 legend(顺流信号,dt=205ns,dt=210ns,dt=215ns); f=2*10^6;a=10^12;tao=5*10^-6;fs=10^9;T=1/fs;t=0:T:2*tao; %初始化参数 ? 实验结果 图2. 超声波逆流信号采样整体图 图3. 超声波逆流信号采样局部图 三、 相关法估计时差 ? 实验程序 由于超声波逆流波形与顺利波形相差的只是一个时延，所以采用相关法分析两者的波形，从而得到波形的时延，当相关性分析所得的最大值对应的t值即为两者的时延。下面分析了时延为205ns的情况，其他的同理。 f=2*10^6;a=10^12;tao=5*10^-6;fs=10^9;T=1/fs;t=0:T:2*tao; %初始化参数 t1=2.05*10^-7;t2=2.10*10^-7;t3=2.15*10^-7; s=sin(2*pi*f*t).*exp(-a*(t-tao).^2/2); %函数关系式 s1=sin(2*pi*f*(t-t1)).*exp(-a*(t-t1-tao).^2/2); %函数关系式 s2=sin(2*pi*f*(t-t2)).*exp(-a*(t-t2-tao).^2/2); %函数关系式 s3=sin(2*pi*f*(t-t3)).*exp(-a*(t-t3-tao).^2/2); %函数关系式 [a,b]=xcorr(s,s1); plot(b*T,a) xlabel(‘t/s’); title(‘信号相关性分析’) %图表信息 [maxy,maxi]=max(a); N=(maxi-length(s))*T N = -2.0500e-007 ? 实验结果 图4. 信号相关性分析效果图 四、 噪声影响 ? 实验程序 f=2*10^6;a=10^12;tao=5*10^-6;fs=10^9;T=1/fs;t=0:T:2*tao; %初始化参数 t1=2.05*10^-7;t2=2.10*10^-7;t3=2.15*10^-7; s=sin(2*pi*f*t).*exp(-a*(t-tao).^2/2); %函数关系式 s1=sin(2*pi*f*(t-t1)).*exp(-a*(t-t1-tao).^2/2); %函数关系式 s2=sin(2*pi*f*(t-t2)).*exp(-a*(t-t2-tao).^2/2); %函数关系式 s3=sin(2*pi*f*(t-t3)).*exp(-a*(t-t3-tao).^2/2); %函数关系式 g(1)=0; for i=1:10000 noise1=g(i)*randn(1,length(s)); noise2=g(i)*randn(1,length(s)); saddnoise1=s+noise1; s1addnoise2=s1+noise2; [a,b]=xcorr(saddnoise1,s1addnoise2); [maxy,maxi]=max(a); n(i)=(maxi-length(s))*T; if i10000 g(i+1)=g(i)+0.0001; else g(i+1)=g(i)+0.0001; end en