.. ..
一、 设计目的
通过运用MATLAB对函数进行Z域分析和单边带信号的调制与解调,使我们进一步加深对MATLAB的认识和运用,以实现以下目的: 1.本次试验进一步熟悉了MATLAB软件的使用方法及相关的操作。 2.对Z变换及其反变换函数在MATLAB中的调用有了掌握。 3.理论与实际的仿真相结合,更直观的看到结果。 4.观察了单边带信号调制与解调后的图像,加深认识。
二、 设计原理
MATLAB是The MathWorks公司在1984年推出的一种商品化软件,它提供了大量丰富的应用函数,并且具有扩充的开放性结构。目前,该软件包涵盖了控制系统应用、数字信号处理、数字图像处理、通讯、神经网络、小波理论分析、优化与统计、偏微分方程、动态系统实时仿真等多学科专业领域。
其中单边带调制信号是将双边带信号中的一个边带滤掉而形成的。根据方法的不同,产生单边带调制信号的方法有:滤波和相移法。
由于滤波法在技术上比较难实现所以在此我们将用相移法对单边带调制与解调系统进行讨论与设计。
三 、设计容和MATLAB图像
1、数字系统的响应
.. .. ..
.. ..
z2?2z?1H(Z)?3Z?0.5Z2?0.005Z?1?0.3?y(0)=-1,y(-1)=1?x(n)=e-n?试画出系统的零极点分布图、?求系统的单位脉冲与阶跃响应?求系统的零输入、零状态、全响应?求系统的频率响应?判断系统是否稳定?要求在一幅图形中给出所有结果
源代码如下:
b=[0 1 2 1 0];
a=[1 -0.5 0 0.3 -0.005]; subplot(421);zplane(b,a); title('系统的零极点图'); subplot(422);impz(b,a,21); title('单位脉冲响应'); subplot(423);stepz(b,a,21); title('单位阶跃响应'); N=21;n=0:N-1;
.. .. ..
.. ..
x=exp(-n); x0=zeros(1,N); y0=[1,-1]; xi=filtic(b,a,y0); y1=filter(b,a,x0,xi); xi0=filtic(b,a,0); y2=filter(b,a,x,xi0); y3=filter(b,a,x,xi); [h w]=freqz(b,a,21); subplot(424);stem(n,y1); title('零输入响应');grid on; subplot(425);stem(n,y2); title('零状态响应');grid on; subplot(426);stem(n,y3); title('系统的全响应');grid on; subplot(427);plot(w,abs(h)); title('幅频特性曲线');grid on; subplot(428);plot(w,angle(h)); title('相频特性曲线');grid on;
MATLAB运行结果图样:
.. .. ..
.. ..
Imaginary Part Amplitude系统的零极点图10-1-2202Real Part单位阶跃响应单位脉冲响应50-5051015n (samples)零输入响应20 Amplitude1050051015n (samples)零状态响应2010-105101520系统的全响应5050-510500123051015幅频特性曲线20-550-5051015相频特性曲线200123
本次课题分析:
在系统的零极点分布图中,符号“○”表示零点,符号“○”旁边的数字表示零点的阶数,符号“ⅹ”表示极点,图中的虚线画的是单位圆。由系统的零极点分布图可以看出,该因果系统的点全在单位圆,故该离散时间系统是稳定的。
2、单边带调制与解调
.. .. ..
.. ..
载波频率:250 KHZ;?x=sa(200t)?原始信号及频谱?调制信号及频谱?解调信号及频谱?要求在一幅图形中给出所有结果
源代码如下:
dt=0.001; fs=1/dt; t=-0.2:dt:0.2; fc=250; m=sinc(200*t); m=m.*m;
subplot(421);plot(t,m); axis([-0.2 0.2 -1.2 1.2]); title('原始信号');
.. .. ..