MATLAB实验报告(5-9)(3)

b6=[1]; a6=[2,0,1]; sys6=tf(b6,a6); subplot(121); pzmap(sys6); axis([-2,2,-2,2]); subplot(122); impulse(b6,a6);

第七次MATLAB实验报告

一、实验目的

1.学会运用MATLAB完成信号抽样及对抽样信号的频谱进行分析。 2.学会运用MATLAB改变抽样间隔观察抽样后的信号的频谱变化。 3.学会运用MATLAB对抽样后的信号进行重建。 4.学会运用MATLAB表示的常用离散时间信号。 5.学会运用MATLAB实现离散时间信号的基本运算。

二、实验任务

此次实验旨在让我们学习有关抽样信号的各种编程和离散时间信号的一些运算，同时明白例题的编程并举一反三完成作业习题。

三、实验内容

(1).信号抽样：从数学上讲，抽样过程就是抽样脉冲P(t)和原连续信号f(t)相乘的过程，即fs(t)=f(t)p(t);

例9-1： Ts=1; dt=0.1; t1=-4:dt:4;

ft=((1+cos(t1))/2).*(uCT(t1+pi)-uCT(t1-pi)); subplot(221) plot(t1,ft),grid on axis([-4,4,-0.5,1.2])

xlabel('Time(sec)'),ylabel('f(t)'); title('升余弦脉冲信号') N=500; k=-N:N; W=pi*k/(N*dt); Fw=dt*ft*exp(-j*t1'*W); subplot(222)

plot(W,abs(Fw)),grid on axis([-10,10,-1,1.1*pi]); title('升余弦脉冲信号的频谱'); xlabel('\\omega'),ylabel('F(w)') t2=-4:Ts:4;

fst=(1+cos(t2)/2).*(uCT(t2+pi)-uCT(t2-pi)); subplot(223) plot(t1,ft,':'),hold on stem(t2,fst),grid on axis([-4,4,-0.5,1.2]);

title('抽样后的信号'),xlabel('Time(sec)'),ylabel('fs(t)')

hold off

Fsw=Ts*fst*exp(-j*t2'*W); subplot(224)

plot(W,abs(Fsw)),grid on axis([-10,10,-1,1.1*pi]); xlabel('\\omega'),ylabel('Fs(w)'); title('抽样信号的频谱')

(2).抽样定理：若f(t)是带限信号，带宽为Wn，则信号f(t)可以用等间隔的臭氧治来唯一标示，f(t)经抽样后的频谱Fs(w)就是将f(t)的频谱F(w)在频率轴上以抽样频率ws为间隔进行周期延拓。

(3).信号重建：抽样定理表明，当抽样间隔笑语奈奎斯间隔时，可以用抽样信号fs(t)唯一的表示原信号f(t),即信号重建。为了从频谱中无失真地恢复原信号，可以采用截止频率wc>=wm的理想低通滤波器。 例9-4：