实验三-诱发响应的提取

实验三诱发响应的提取学生姓名刘冰学院名称精密仪器与光电子工程专业生物医学工程学号3010202286实验目的了解并掌握诱发响应的提取方法。实验原理由于对刺激的滞后响应具有随机性，诱发响应具有不同的潜伏期，所以不能够简单地对诱发响应累加求平均，而应该设法求出各次响应的潜伏期后，在原始诱发响应信号中除去潜伏期的影响后，再累加求平均。具体的步骤为求出多次实验的平均响应作为初步估计求平均响应与诱发响应的时间互相关求出互相关函数的最大值对应的时间t求出响应的最终估计三．Matlab程序clc;clearall;closeall;%*****产生一个幅度衰减的模拟信号fs=1000;t=0:1/fs:1;f=10;s=sin(2*pi*f*t).*exp(-15*t);%---画出理想的诱发响应figure,plot(t,s);title('模拟诱发响应');xlabel('时间/s')ylabel('幅值')%**多次实验N取1000，%***模拟诱发潜伏期的时间大概在50ms左右N=1000;time=50;%****用自带的随机数产生N个50左右的随机数，这些数字符合均匀分布tx=time+round(10*rand(1,N));L=length(s);%叠加一个符合正态分布的白噪声fori=1:N;l=zeros(1,tx(i));%用tx产生一个长度为50左右的0数组s1=[ls(1:L-tx(i))];%将0数组与模拟信号连起来，模拟含有潜伏期的信号e=0.1*randn(1,L);%产生随机数模拟噪声x(i,:)=s1+e;%叠加模拟信号与白噪声，并将多有信号存入矩阵end;figure,plot(t,x');title('N次实验的模拟原始信号');p=mean(x);%叠加后求平均，这里求得的平均结果会丢失高频信息figure,plot(t,p');title('N次实验结果平均');x2=zeros(N,L);fori=1:N;R(i,:)=xcorr(p,x(i,:));%求平均响应与每个诱发响应的互相关[MI(i)]=max(R(i,:));%找出互相关函数的最大值对应的时间点min1=min(I);%以最小的时间点为标准作为起点%每个响应需要移动的距离lx(i)=I(i)-min1;%移动诱发响应并将结果存储x2(i,1:L-lx(i))=x(i,1+lx(i):1001);end;%求出响应的最终估计p2p2=mean(x2);figure,plot(t,p2');title('最终估计结果');xlabel('时间/s')ylabel('幅值')%以下程序验证课堂PPT上的一句话“潜伏期的影响使迭加平均结果丢失高频信息”fft1=abs(fft(s));figure,plot(t.*fs,fft1);fft3=abs(fft(p2));figure,plot(t.*fs,fft3);做fft的结果：实验结论：1.对于原始信号的fft的高频极值：对于最终估计的FFT的高频极值：对比可知其相对结果高频