山东大学信号与系统MATLAB实验7-8.docxVIP

实验七1.编程实现对的抽样信号（参考例 7-2） 1）画出抽样间隔分别为 1, 1.5, 2 时抽样后的信号时域图和频谱图；源程序（我还加上了f（t）的时域图和频谱图）：dt=0.1;t1=-4:dt:4; %设置自变量范围ft=(1+cos(t1)).*(heaviside(t1+pi)-heaviside(t1-pi)); %输入连续函数f（t）=1+cost1，|t1|πsubplot(421) %将多个图画到一个平面，排成4行2列，现在第1个图上绘制plot(t1,ft) %画图grid on %在画图时添加网格线axis([-4,4,-0.2,2.2]) %图像上横坐标区间定义为【-4，4】，纵坐标区间定义为【-0.2，2.2】xlabel(Time(sec)),ylabel(f(t)); %分别显示横纵坐标标题title(升余弦脉冲信号) %图像名称为“升余弦脉冲信号”N=500;k=-N:N; %设置自变量范围，k为-500到500，间隔为1的数组W=pi*k/(N*dt); %W为范围为-10π到10π，间隔为π/50的数组Fw=dt*ft*exp(-j*t1*W); %输入函数F(w)=∫f(t)e^(-jt1w)dt，这个是对f（t）的傅里叶变换subplot(422) %将多个图画到一个平面，排成4行2列，现在第2个图上绘制plot(W,abs(Fw)) %绘制信号的频谱图grid on %在画图时添加网格线axis([-10,10,-0.2,2.2*pi]) %图像上横坐标区间定义为【-10，10】，纵坐标区间定义为【-0.2，2.2π】xlabel(\omega),ylabel(F(w)); %“\”为转义字符，\omega用来显示希文字母wtitle(升余弦脉冲信号的频谱) %图像名称为“升余弦脉冲信号的频谱”Ts=[1,1.5,2]; %Ts为一行三列矩阵for m=1:3 %for循环t2=-4:Ts(m):4; %设置自变量范围fst=(1+cos(t2)).*(heaviside(t2+pi)-heaviside(t2-pi)); %输入抽样函数fs（t）=1+cost2，|t2|πsubplot(4,2,4*Ts(m)-1) %将多个图画到一个平面，排成4行2列，现在第4*Ts(m)-1个图上绘制plot(t1,ft,:) %绘制t1与ft的关系图像，plot(t1,ft,:)中‘：’表示对应线用的线条（也就是线条是....这样的）hold on %当前轴（坐标系）及图形保持而不被刷新，准备接受此后的绘制（即多图共存）stem(t2,fst)%绘制t2与fst关系的离散图形grid on %在画图时添加网格线axis([-4,4,-0.2,2.2]) %图像上横坐标区间定义为【-4，4】，纵坐标区间定义为【-0.2，2.2】xlabel(Time(sec)),ylabel(fs(t)); %分别显示横纵坐标标题title([抽样后的信号,Ts=,num2str(Ts(m))]) %图像名称为“抽样后的信号，Ts=Ts（m）”hold off %hold off使当前轴及图形不再具备被刷新的性质，与hold on功能相反，在当前图的轴（坐标系）中画了一幅图，再画另一幅图时，原来的图就看不到了，在轴上绘制的是新图，原图被替换了Fsw=Ts(m)*fst*exp(-j*t2*W); %输入函数Fs(w)=∫fs(t)e^(-jt2w)dt，这个是对fs（t）的傅里叶变换subplot(4,2,4*Ts(m)) %将多个图画到一个平面，排成4行2列，现在第4*Ts(m)个图上绘制plot(W,abs(Fsw)) %绘制信号的频谱图grid on %在画图时添加网格线axis([-10,10,-0.2,2.2*pi])%图像上横坐标区间定义为【-10，10】，纵坐标区间定义为【-0.2，2.2π】xlabel(\omega),ylabel(Fs(w));%分别显示横纵坐标标题title([抽样后的信号的频谱,Ts=,num2str(Ts(m))])%图像名称为“抽样后的信号的频谱，Ts=Ts（m）”end %for循环结束程序运行结果：结论分析：要想抽样样本值包含原连续信号全部信息，则Ws=2π／Ts,Ws2Wm,已知Wm=2,故Ws4,Ts＜π/2,故当Ts=1或1.5时抽样信号包含原连续信号的全部信息，而当Ts=2时，抽样信号丢失了原连续信号的部分信息。 2）对三个抽样后的信号用合适的低通滤波器处理，比较输出后的信号，说明抽样定理的准确性。源程序：wc=[2.4,2,1.6]; %Ts=1时得到wm=2，取