测试信号实验三作业.docVIP

WORD完美整理版 范文范例 参考指导 实验三作1. 已知信号f1(t)和f2(t)的时域波形如下，试用Matlab绘制卷积积分的波形。 输入采样时间间隔分别取0.01、0.05和0.1，分别获得的三个波形，画出三个波形并比较它们，说明卷积为什么产生相异的效果。（需调用Heaviside与gggfconv函数） 代码： t1=-1:0.01:3; f1=Heaviside(t1)-Heaviside(t1-2); % 定义信号 t2=t1; f2=0.5*t2.*(Heaviside(t2)-Heaviside(t2-2)); % 定义信号 [t,f]=gggfconv(f1,f2,t1,t2); % 计算卷积积分并绘出时域波形 %%%%%%%%% function f=Heaviside(t) f=(t0); %%%%%%%% function [f,t]=gggfconv(f1,f2,t1,t2) d=input(请输入采样时间间隔：); f=conv(f1,f2); f=f*d; ts=t1(1)+t2(1); l=length(t1)+length(t2)-2; t=ts:d:(ts+l*d); subplot(2,2,1);plot(t1,f1) subplot(2,2,2);plot(t2,f2) subplot(2,2,3);plot(t,f) 当间隔t=0.01时 当时间间隔t=0.05时 当时间间隔t=0.1时 采样时间间隔Ts越小，影响t,从而影响f(t1)和f(t2),而绘制的f(t)波形就越接近实际波形。 2. 某离散系统差分方程为，激励信号为，序列数取为40，用MATLAB求该系统的单位脉冲响应和零状态响应，按三行一列的格式在一张图中分别画出、以及波形，并给出波形名。（可用impz及filter函数） 解：matlab程序如下 n = -20:20; a = [2 -3 1]; b = [2]; figure(1); subplot(311); impz(b,a,n); title(单位脉冲响应h(n)); xlabel(n); k=0:40; xn = cos(k*pi/3); y1 = filter(b,a,xn); subplot(312); stem(k,xn,filled); title(激励信号x(n)=cos(n*pi/3); xlabel(n); subplot(313); stem(k,y1,filled); title(零状态响应); xlabel(n); 图形如下 3. 已知信号，利用MATLAB实现采样间隔时的时域采样，分别按照一行两列的格式画出与其频谱，以及3个抽样信号及其频谱，并给出波形名。比较四幅图像的异同，说明采样间隔不同对频谱生成的影响。 解： matlab程序如下： display(奈奎斯特抽样间隔1.5秒，Ts1.5，过采样，Ts1.5，欠采样); display(Please input the value of sample period); Ts = input(Ts = ); % 绘制有限长余弦信号 y = cos(2/3*pi*t) t = 0:0.01:40; y = cos(2/3*pi*t); subplot(221); plot(t,y); axis([0 6 -1.1 1.1]); xlabel(t 单位 :s,Fontsize,8); title(f(t)); line([0 6],[0 0],color,[0 0 0]); % 数值求解余弦信号的频谱 N = 300; W = 2*pi*5; k = -N:N; w = k*W/N; Y = 0.01*y*exp(-j*t*w); % 求 的傅里叶变换 Y = abs(Y); subplot(222); plot(w/pi,Y) axis([-2,2,0,pi*7+0.2]); title(F(j\omega)); xlabel(\omega 单位:pi); % 抽样后的余弦信号 subplot(223); plot(t,y,b:); % 绘制包络 hold on t2=0:Ts:40; y2=cos(2/3*pi*t2); stem(t2,y2); axis([0 6 -1.1 1.1]); xlabel(t 单位:s,Fontsize,8); title(fs(t