生产实习报告内容.doc

一． 实验内容 1.设，,将和采样形式离散序列，参数选择如下： 采样率,频率=30，=2²()， =2.5。 编程计算： ⑴； ⑵(m)； ⑶ 将时移得，计算互相关函数，并检测出时移值； ⑷ 画出各离散序列图形：，，，，。 二． 算法分析 1、在求得x(n)、h(n)的基础上计算出函数x(n)、h(n)的取点个数，在求卷积的过程中利用一个二重循环求得y(n); 2、已知h(n)的基础上，求自相关也是利用一个二重循环求得 3、因为在利用时可以发现右边函数中的自变量会出现负号，因此令n=m+ n0,这样就可以求出时移h’(n)。 相关函数计算 对序列进行相关计算的基本方法与卷积计算相似，主要差别是卷积计算中需要反褶两个序列中的一个序列，而在相关计算中则不反褶序列，但要考虑取适当的相关时长度。 设x（n）和y（n）的序列有限长，长度分别为Nx和Ny。 互相关函数公式为：rxy(n)= ∑x(m)y(m-n) 自相关函数公式为：，rxx(m)= ∑x(m)x(m-n) 式中，n为时移值。 相关函数的基本计算步骤如下： 设相关时窗长度为W，因此计算相关函数时取时移值为： n=-W/2, -W/2+1, -W/2+2,…，-2，-1，0，1，2 …，-W/2+2, -W/2+1, W/2 第一步：初始化r（n）=0； 第二步：对于时移值n的一个取值n0，计算r∑h(m)x(n0-m)，计算此和式由一个循环实现； 第三步：依次取以下时移值，计算出r（n）: n=-W/2, -W/2+1, -W/2+2,…，-2，-1，0，1，2 …，-W/2+2, -W/2+1, 三． 实现方法 .源程序如下： dt=0.004; T=30; f=30; M=2.5; n=-T*dt:dt:T*dt; x=1-abs(n)/(T*dt); subplot(3,2,1);stem(n,x,r.);grid on xlabel(x(t)); %h函数： a=2*f.^2*log(M); t=0:dt:2*T*dt; h=exp((-a)*t.^2).*sin(2*pi*f*t); subplot(3,2,2);stem(t,h,b.);grid on xlabel(h(t)); %x和h的卷积： N=length(x); M=length(h); for(n=1:N+M-1) y(n)=0; for(m=1:M) k=n-m+1; if(k=1k=N) y(n)=y(n)+h(m)*x(n-m+1); end end end subplot(3,2,3);plot(y);grid on xlabel(y(t)=x(t)*h(t)) %关于h的自相关函数： hlen = length(h); th=0:0.004:0.1; f=30; M=2.5; a=2*f.^2*log10(M); h=exp(-a*th.^2).*sin(2*pi*f*th); r=zeros(1,50); for n=1:50 for m=1:25 if(m-n=0m-n=-24) r(n)=r(n)+h(m)*h(25+m-n); end end end n=-0.096:0.004:0.1; subplot(3,2,4);stem(n,r,r.);grid on %h的延迟函数： th=0:0.004:0.1; f=30; M=2.5; a=2*f.^2*log10(M); h=exp(-a*th.^2).*sin(2*pi*f*th); th2=0.05:0.004:0.15; h2=exp(-a*(th2-0.05).^2).*sin(2*pi*f*(th2-0.05)); subplot(3,2,5);stem(th2,h2,r.);grid on %h的互相关函数： th=0:0.004:0.1; f=30; M=2.5; a=2*f.^2*log10(M); h=exp(-a*th.^2).*sin(2*pi*f*th); th2=0.05:0.004:0.15; h2=exp(-a*(th2-0.05).^2).*sin(2*pi*f*(th2-0.05)); r1=zeros(1,50); for n=1:50 for m=1:25 if(m-n=0m-n=-24) r1(n)=r1(n)+h(m)*h2(25+m-n); end end end n