MATLAB仿真教程分析和总结.docx

一、设计目的

通过运用MATLAB对函数进行Z域分析和单边带信号的调制与解调，使我们进一步加深对MATLAB的认识和运用，以实现以下目的：

本次试验进一步熟悉了MATLAB软件的使用方法及相关的操作。

对Z变换及其反变换函数在MATLAB中的调用有了掌握。

理论与实际的仿真相结合，更直观的看到结果。

观察了单边带信号调制与解调后的图像，加深认识。

二、设计原理

MATLAB是TheMathWorks公司在1984年推出的一种商品化软件，它提供了大量丰富的应用函数，并且具有扩充的开放性结构。目前，该软件包涵盖了控制系统应用、数字信号处理、数字图像处理、通讯、神经网络、小波理论分析、优化与统计、偏微分方程、动态系统实时仿真等多学科专业领域。

其中单边带调制信号是将双边带信号中的一个边带滤掉而形成的。根据方法的不同，产生单边带调制信号的方法有：滤波和相移法。

由于滤波法在技术上比较难实现所以在此我们将用相移法对单

边带调制与解调系统进行讨论与设计。

三、设计内容和MATLAB图像

1、数字系统的响应

H

H(Z)?

z2?2z ?1

Z3?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);[hw]=freqz(b,a,21);

subplot(424);stem(n,y1);title(零输入响应);gridon;subplot(425);stem(n,y2);title(零状态响应);gridon;subplot(426);stem(n,y3);title(系统的全响应);gridon;subplot(427);plot(w,abs(h));title(幅频特性曲线);gridon;subplot(428);plot(w,angle(h));title(相频特性曲线);gridon;

MATLAB运行结果图样：

pedutil

p

m

traP

yra

gan

ga

mI

10

5

0

系统的零极点图

21

2

0

-20

-2

0

2

A 0

5

10

15

20

RealPart

n(samples)

单位阶跃响应

ude 5

ilp

ilp

m -5

1

0

-1

单位脉冲响应

零输入响应

A 0 5 10 15 20

n(samples)

零状态响应

5

0

0 5 10 15 20

系统的全响应

5

0

-5

0 5 10 15

幅频特性曲线

10

5

-5

20 0 5 10 15 20

相频特性曲线

5

0

0

0 1 2 3

-5

0 1 2 3

本次课题分析：

在系统的零极点分布图中，符号“○”表示零点，符号“○”旁边的数字表示零点的阶数，符号“ⅹ”表示极点，图中的虚线画的是单位圆。由系统的零极点分布图可以看出，该因果系统的点全在单位圆内，故该离散时间系统是稳定的。

2、单边带调制与解调

载波频率：

载波频率：250KHZ;

? 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.20.2-1.21.2]);

title(原始信号);c=cos(2*pi*fc.*t);d=sin(2*pi*fc.*t);

u=(m.*c-imag(hilbert(m)).*d)/2;v=(m.*c+imag(hilbert(m)).