信号与系统分析习报告 MATLAB.docVIP

信号与系统分析的MATLAB实现 学生：XXX 学号: 完成日期：2004.12.26 上机题一： 设，,对，采样形成离散序列。参数选择：采样率4ms，频率30HZ，2.5，2。编程计算： 计算 将时移,得=，计算 画出各序列图形，，，， 算法及分析 卷积函数算法与分析：由所学知识可知，连续信号的卷积是积分运算，对连续信号按采样定理进行采样，可得到离散信号或离散序列，这时卷积积分变成了卷积求和，是一种延迟加权求和运算。设有离散序列和，则有卷积公式。设，分别为序列、的长度，则。根据MATLAB语言的向量存储表示可推出，，我们用向量存储的结果，可设。 由于给出的算法涉及求和运算，于是在设计程序时，需要使用两重for循环，外重循环置的初始值，内重循环实现循环求和，而在循环之前需判断是否存在，即判断的值是否在区间内，这用一个if语句即可实现。具体程序实现如下： function f=dconv(f1,f2,n1,n2)%定义求卷积的函数 k3=length(n1)+length(n2); for n=2:k3 f(n-1)=0; for m=1:length(n1) if ((n-m)=1 (n-m)=length(n2)) f(n-1)=f(n-1)+f1(m)*f2(n-m);%计算卷积和 end end end （2）相关函数算法与分析：根据相关函数公式，其中为延迟时间变量，设，分别为序列、的长度，，可知，。用向量存储，则。 类似卷积计算，相关函数的程序设计也要使用两重for循环，外重循环置初值,内重循环实现求和。唯一不同的是，此时要判断的是的值是否在区间内。具体程序实现如下： function f=rrconv(f1,f2)%定义相关函数 k1=length(f1)-1; k2=length(f2)-1; for m=-k2:k1 f(m+k2+1)=0; for n=1:length(f1) if ((n-m)=1 (n-m)=length(f2)) f(m+k2+1)=f(m+k2+1)+f1(n)*f2(n-m); end end end （3）总体设计 2.程序列表 M=2.5; f=30;%频率 T=16;%点数 dt=0.004;%采样间隔 n1=-T*dt:dt:T*dt; f1=1-abs(n1)*f;%定义序列f1 subplot(3,2,1) plot(n1,f1,m-*)%绘制f1的图形 xlabel(t) ylabel(f1) title(f1=(1-|t|/T)[u(t+T)-u(t-T)]) grid on a=2*(f^2)*log(M); n2=0:dt:2*T*dt; f2=sin(2*pi*f*n2).*(exp(-a.*(n2.^2)));%定义序列f2 subplot(3,2,2) plot(n2,f2,r--*)%绘制f2的图形 xlabel(t) ylabel(f2) axis([0 0.1 -0.5 1]) title(f2=sin(2*pi*f*t)(e^-at^2)) grid on f3=dconv(f1,f2,n1,n2);%调用卷积函数计算f3=f1*f2 subplot(3,2,3) k3=0:(length(n1)+length(n2)-2); plot(k3,f3,m-*)%绘制卷积图形 xlabel(t) ylabel(f3) grid on f4=rrconv(f2,f2);%调用相关函数计算序列f2的自相关 subplot(3,2,5) k4=-(length(n2)-1):(length(n2)-1); plot(k4,f4,m-*)%绘制自相关函数图形 xlabel(t) ylabel(f2自相关) axis([-15 15 -1 2]) grid on n0=20; m=length(f2); for n=1:(m+n0)%求f2的时移序列f2(n-n0) if n=n0 fn2(n)=0; else fn2(n)=f2(n-n0); end end for n=(m+1):(m+n0) f2(n)=0; end subplot(3,2,4) k2=0:dt:(length(fn2)-1)*dt; plot(k2,fn2,r--*)%绘制时移图形 xlabel(t) ylabel(f2(n-n0)) axis([0.05 0.15 -0.5 1])