东华大学信号与线性系统课程设计答案

东华大学 信号与线性系统 课程设计

班级：

学号： 姓名：

东华大学信息学院

目 录

实验一 连续信号的时域分析 …………………………… 3 实验二 连续时间系统的时域分析 ……………………… 8 实验三 连续信号的频域分析 ……………………………14 实验四 连续系统的频域分析 ……………………………22 实验五 信号采样与重建 …………………………………33 实验六 离散时间信号和系统分析 ………………………39 实验总结 ………………………………………………… 46

- 2 -

实验一 连续信号的时域分析 一、实验目的

1、熟悉MATLAB 软件。

2、掌握常用连续信号与离散信号的MATLAB 表示方法。

二、实验设备

安装有matlab6.5 以上版本的PC 机一台。 三、实验原理

四、实验内容

1、用MATLAB表示连续信号：Aet，Acos(ω0 t +ϕ )，Asin(ω0 t +ϕ )。 源程序：

clc clear close syms t;

f1=2*exp(t);

f2=2*cos(3*t+4); f3=2*sin(3*t+4); subplot(2,2,1);

ezplot(f1,[-10,2]); xlabel('t');

- 3 -

title('f(t)=2e^t'); grid on;

subplot(2,2,2); ezplot(f2,[-5,5]); xlabel('t');

title('f(t)=2cos(3t+4)'); grid on;

subplot(2,2,3); ezplot(f3,[-5,5]); xlabel('t');

title('f(t)=2sin(3t+4)'); grid on

f(t)=2et20-10-5tf(t)=2sin(3t+4)210-1-2-50t5 f(t)=2cos(3t+4)210-1-20-50t52、用MATLAB表示抽样信号(sinc(t))、矩形脉冲信号(rectpuls(t, width)) 及三角脉冲信号(tripuls(t, width, skew))。 源程序：

clc clear close

t=-5:0.01:5; f1=sinc(t);

f2=3*rectpuls(t,4);

- 4 -

f3=3*tripuls(t,4,0); subplot(2,2,1); plot(t,f1); xlabel('t');

title('f(t)=sinc(t)'); grid on;

subplot(2,2,2) plot(t,f2); xlabel('t');

title('f(t)=3rectpuls(t,4)'); grid on;

axis([-5,5,-1,4]); subplot(2,2,3); plot(t,f3); xlabel('t');

title('f(t)=3rectpuls(t,4,0)'); grid on;

axis([-5,5,-1,4]);

f(t)=sinc(t)1430.5210-0.5-50tf(t)=3rectpuls(t,4,0)5-1-50t5f(t)=3rectpuls(t,4)043210-1-50t5 3、编写如图3 的函数并用MATLAB 绘出满足下面要求的图形。 (1) f (−t); (2) f (t − 2); (3) f (1− 2t); (4) f (0.5t +1).

- 5 -

源程序：

clc clear close

t1=-14:0.05:2; t2=0:0.05:16; t=-6.5:0.05:1.5; t4=-6:0.05:26;

f1=4*rectpuls(-t1-6,12)+3*tripuls(-t1-6,4,0); f2=4*rectpuls(t2-8,12)+3*tripuls(t2-8,4,0);

f3=4*rectpuls(1-2*t-6,12)+3*tripuls(1-2*t-6,4,0);

f4=4*rectpuls(0.5*t4+1-6,12)+3*tripuls(0.5*t4+1-6,4,0); subplot(2,2,1); plot(t1,f1); xlabel('t'); title('f(-t)'); subplot(2,2,2); plot(t2,f2); xlabel('t');

title('f(t-2)'); subplot(2,2,3); plot(t,f3); xlabel('t');

title('f(1-2t)'); subplot(2,2,4); plot(t4,f4); xlabel('t');

title('f(0.5t+1)');

- 6 -

f(-t)820-15820f(t-2)-10-5tf(1-2t)050510tf(0.5t+1)1520820-10-5t05820-10010t2030 - 7 -

实验二 连续时间系统的时域分析

一、实验目的

1、掌握卷积计算方法。

2、掌握函数lsim，impulse，step 的用法，lsim 为求取零状态响应， impulse 为求取单位脉冲响应，step 为求取单位阶跃响应。 3、运用课堂上学到的理论知识，从RC、RL 一阶电路的响应中正确 区分零输入响应、零状态响应、自由响应与受迫响应。

二、实验设备

安装有matlab6.5 以上版本的PC 机一台。 三、实验原理

- 8 -

四、实验内容

1. 分别用函数lsim 和卷积积分两种方法求如图7 所示系统的零状态 响应。其中L=1，R=2，e(t) = etε(t)，i(0− ) = 2。

 e(t)  L i(t) R

- 9 -

源程序： 方法一：

clear close clc

t=0:0.01:10; f=exp(-t); a=[1 2]; b=[1];

y=lsim(b,a,f,t); plot(t,y);

xlabel('Time(sec)'); ylabel('i(t)');

axis([0,10,-0.025,0.275]); grid on;

0.250.20.15i(t)0.10.0500123456Time(sec)710 方法二：

易求得系统的冲激响应为e2tε(t)

clear; clc; close;

- 10 -

syms t x; e=exp(-x);

h=exp(-2.*(t-x)); i=int(e.*h,x,0,t); ezplot(i,[0,10]); xlabel('Time(sec)'); ylabel('i(t)');

title('exp(-t)*exp(-2t)'); grid on;

2. 求上述系统的冲激响应与阶跃响应。 冲激响应源程序：

clear; close; clc;

a=[1 2]; b=[1];

impulse(b,a,10);

xlabel('Time(sec)'); ylabel('i(t)');

- 11 -

axis([0,10,-0.1,1]); grid on;

Impulse Response0.90.80.70.60.5i(t)0.40.30.20.10-0.10123456710Time(sec) (sec) 阶跃响应源程序：

clear; close; clc;

a=[1 2]; b=[1];

step(b,a,10);

xlabel('Time(sec)'); ylabel('i(t)');

axis([0,10,0,0.55]); grid on;

- 12 -

Step Response0.50.40.3i(t)0.20.100123456710Time(sec) (sec) 五、思考题

1. 为什么连续时间系统的零状态响应为激励与冲击响应的卷积？

答：根据卷积的定义，函数e(t)与函数h(t）相卷积后，就是在变量由负无穷到正无穷范围内，对于某一t值时乘积e(τ)h(t-τ)曲线下的面积，也就是：r(t)=e(t)*h(t），又零状态响应与系统的特性和外加激励有关，所以如问题。

2. 利用卷积积分法计算系统响应应从几个方面进行?

答：利用卷积积分法先要将系统的冲击响应求出，之后再将其与激励卷积即可

- 13 -

实验三 连续信号的频域分析 一、实验目的

1. 掌握周期信号的频谱—— Fourier 级数的分析方法。

2. 深入理解信号频谱的概念，掌握典型信号的频谱以及Fourier 变换 的主要性质。

3. 掌握调制与解调的基本原理及滤波器的使用。

二、实验设备

安装有matlab6.5 以上版本的PC 机一台。

三、实验原理

- 14 -

- 15 -

- 16 -

四、实验内容

1. 求如图所示周期矩形脉冲信号的Fourier级数表达式，画出频谱图，并用前N次谐波合成的信号近似。

f(t)1o-1TT/2t 源程序：

close clear clc

syms t n;

T=2*pi; %设T为2*pi N=9; %设N为9

f=heaviside(t)-2*heaviside(t-T/2)+heaviside(t-T); subplot(2,2,1); ezplot(f,[0,2.1*pi]);

title('原函数');

h=exp(-j*n*2*pi/T*t); %用指数傅里叶级数表示，角频率2*pi/T A1=int(f.*h,t,0,T); A=2/T*A1; for n1=-N:-1

C(n1+10)=subs(A,n,n1); end

for n1=1:N

C(n1+10)=subs(A,n,n1); end

C(10)=0;%N=0时系数 subplot(2,2,3); k=-N:N;

stem(k,abs(C)); %作出幅度谱 ylabel('Cn 的幅度'); xlabel('\\Omega'); subplot(2,2,4);

stem(k,angle(C)); %作出相位谱 ylabel('Cn的相位'); xlabel('\\Omega');

- 17 -

f1=0;

for m=-N:N

f1=f1+1/2*C(m+10)*exp(j*m*t); %前N次谐波合成的信号近似 end

subplot(2,2,2);

ezplot(f1,[0,2.01*pi]);

title('前N次谐波合成的信号近似');

2、试用fourier()函数求下列信号的傅里叶变换F( jω) ,并画出F( jω) (1) f (t) = te−3tε (t) (2) f (t) = sgn(t)

(1)源程序：

clear close clc

syms t x ;

x=fourier(t*exp(-3*t)*heaviside(t)); x

- 18 -

z=abs(x);

ezplot(z);%符号函数作图函数 xlabel('Time(sec)'); ylabel('|F(jω)|'); grid on

得F(jω)= 1/(3+ j*ω)^2

(2) 源程序：

clear close clc

syms t x ;

x=fourier(2*heaviside(t)-1);%2*heaviside(t)-1即为sgn(t) x

z=abs(x);

ezplot(z);%符号函数作图函数 xlabel('Time(sec)'); ylabel('|F(jω)|'); grid on

得F(jω)= -2*j/ω

- 19 -

3、调制信号为一取样信号，利用MATLAB 分析幅度调制(AM)产生 的信号频谱，比较信号调制前后的频谱并解调已调信号。设载波信号 的频率为100Hz。 源程序：

clc; clear; close; Fm=10;

t1=0:0.00002:0.2; syms t v;

x=sin(2.0*pi*Fm*t)/(2.0*pi*Fm*t); subplot(3,2,1); ezplot(x,[0,0.2]);

title('原函数'); Fx=fourier(x,v); subplot(3,2,2);

ezplot(Fx,[-50*pi,50*pi]);

axis([-50*pi,50*pi,-0.05,0.1]);

- 20 -

title('频谱');

y=x*cos(200*pi*t); subplot(3,2,3); b=subs(y,t,t1); plot(t1,b);

title('调制后');

axis([0,0.2,-1,1]); Fy=fourier(y,v); subplot(3,2,4);

ezplot(Fy,[-250*pi,250*pi]);

axis([-250*pi,250*pi,-0.05,0.1]); title('频谱');

z=y*cos(200*pi*t); Fz=fourier(z,v);

G=-heaviside(v-20*pi)+heaviside(v+20*pi);%门函数 Fx1=Fz*G;

x1=2*ifourier(Fx1,v);%滤波过程中幅度减半且反向。 subplot(3,2,5);

ezplot(x1,[0,0.2]);

title('解调后'); subplot(3,2,6);

ezplot(2*Fx1,[-50*pi,50*pi]); axis([-50*pi,50*pi,-0.05,0.1]); title('频谱');

原函数10.5000.050.1t调制后0.150.20.10.050-0.05-1000v频谱100频谱10-10.10.050-0.0500.050.1解调后0.150.20.10.05000.050.1v0.150.2-0.05-1000v100-5000v频谱50010.50 - 21 -

五、思考题

1、根据试验1 的结果，解释Gibbs 现象。

答：因为对于具有不连续点的函数，即使级数的项无限增大，在不连续处，级数之和不收敛于函数f(t)；在跃变点附近的波形，总是不可避免的存在有起伏震荡，从而使跃变点的值超过一形成过冲，造成吉布斯现象。

2、比较周期信号与非周期信号的频谱。

区别：

1．周期信号的频谱为离散频谱，非周期信号的频谱为连续频谱。

2.周期信号的频谱为Fn的分布，表示每个谐波分量的复振幅；而周期信号的频谱为F(jω)的分布，(F(jω)/2)ω表示合成谐波分量的复振幅，所以也将称为频谱密度函数。 联系：

1.都是反映将时域信号表示为正弦类信号时各谐波分量的分布特性。 2.若周期信号是连续非周期信号的周期延拓，则两者的关系为 F(jω)= limT0Fn； Fn=

T0F(j)| T03、调制与解调的基本原理是什么？为什么要进行调制？

调制：调制就是使信号f(t)控制载波的某一个或某些参数（如振幅、频率、相位等），是这些参数按照信号f(t)的规律变化的过程。载波可以是正弦波或脉冲序列。以正弦型信号作载波的调制叫做连续波调制。调制后的载波就载有调制信号所包含的信息，称为已调波。对于连续波调制，已调信号可以表示为：

(t)A(t)cosot(t)

它有振幅频率和相位三个参数构成。改变三个参数中的任何一个都可以携带

- 22 -

同样的信息。因此连续波的调制可分为调幅、调相、和调频。

解调：解调是调制的逆过程，它的作用是从已调波信号中取出原来的调制信号。调制过程是一个频谱搬移的过程，它将低频信号的频谱搬移到载频位置。如果要接收端回复信号，就要从已调信号的频谱中，将位于载频的信号频谱再搬回来。

之所以进行解调，是因为无线电通信是通过空间辐射方式传输信号的，调制过程可以将信号的频谱搬移到容易以电磁波形势辐射的较高频率范围；此外，调制过程可以将不同的信号通过频谱搬移托付至不同频率的载波上，实现多路复用，不至于互相干扰。

- 23 -

实验四 连续系统的频域分析 一、实验目的：

掌握连续时间系统变换域分析的基本方法。

二、实验设备：

安装有matlab6.5以上版本的PC机一台。

三、实验原理

- 24 -

四、实验内容：

如图所示系统：

+-Rf(t)Cy(t)

(a) 对不同的RC值，用freqs函数画出系统的幅频曲线。

易求得H(jω)=1/（1+jωRC），RC的取值依次为100、10、……、0.00001时的幅频曲线。 源程序：

close clear clc

b=[0 1]; for c=-5:2 RC=10^c; a=[RC 1]; freqs(b,a);

axis([10^(-2),10^(5),0.1,1]); hold on end

得到一系列幅频曲线，从左到右依次为RC的取值依次为100、10、……、0.00001时的幅频曲线。图中褐色虚线表示纵坐标取值为0.707，红线表示横坐标为100，绿线横坐标为2000。

- 25 -

(b) 信号f(t)cos(100t)cos(2000t)包含了一个低频分量和一个高频分量。确定适当的RC值，滤除信号中的高频分量并画出信号f(t)和y(t)在

t0~0.2s范围内的波形。提示： |H( jω)|为最大值的2/2处对应

的频率为通带截止频率ωc,首先求取|H( jω)|并找到ωc和RC关系，然后根据题意选定ωc即可确定RC值。 由（a）中的图可知，当RC=-2时符合题意。 源程序：

close clear clc

t=0:0.001:0.2;

f=cos(100*t)+cos(2000*t); subplot(2,1,1); plot(t,f);

y1=cos(100*t)/(1+j*100*10^(-2))+cos(2000*t)/(1+j*2000*10^(-2)); subplot(2,1,2); plot(t,y1)

- 26 -

2、信号任选，分析以下几种情况下信号的频谱和波形变化： （1） 系统满足线性不失真条件时； （2） 系统只满足恒定幅值条件时； （3） 系统只满足相位条件时； （4）系统两个条件均不满足时。

提示：利用fourier求取信号的傅立叶变换E(jω)，然后设计

H( jω) = H( jω) eφ (ω )使之满足不同条件，计算R( jω)= E(jω) H(jω)并画频

谱图。 （1）源程序：

clc clear close

syms t v;

e=exp(-2*abs(t));

- 27 -

subplot(2,3,1); ezplot(e,[-3,3]); axis([-3,3,-0.2,2]); Fe=fourier(e,v); subplot(2,3,2); ezplot(Fe,[-3,3]); title('幅度谱'); axis([-3,3,0,2]); i=1;

for a=-3:0.02:3

R11=subs(Fe,v,a); C(i)=angle(R11); i=i+1; end

b=-3:0.02:3; subplot(2,3,3); plot(b,C);

title('相位谱'); axis([-3,3,-1,1]); H1=2*exp(-j*v*1); R1=Fe*H1;

r1=ifourier(R1,t); subplot(2,3,4); ezplot(r1,[-3,3]);

title('满足线性不失真条件'); axis([-3,3,-0.2,2]) subplot(2,3,5);

ezplot('abs(8/(4+v^2)*exp(-i*v))'); title('幅度谱');

axis([-3,3,0,2.2]); i=1;

for a=-3:0.02:3

R11=subs(R1,v,a); C(i)=angle(R11); i=i+1; end

b=-3:0.02:3; subplot(2,3,6); plot(b,C);

title('相位谱'); axis([-3,3,-3,3]);

- 28 -

exp(-2 abs(t))21.510.5002t满足线性不失真条件21.510.50-20t2-221.510.50幅度谱10.50-0.5-20v幅度谱2-1相位谱-20相位谱221.510.50-20v220-2-202 （2）源程序：

clc clear close

syms t v;

e=exp(-2*abs(t)); subplot(2,3,1); ezplot(e,[-3,3]); axis([-3,3,-0.2,2]); Fe=fourier(e,v); subplot(2,3,2); ezplot(Fe,[-3,3]); title('幅度谱'); axis([-3,3,0,2]); i=1;

for a=-3:0.02:3

R11=subs(Fe,v,a); C(i)=angle(R11); i=i+1; end

b=-3:0.02:3; subplot(2,3,3); plot(b,C);

- 29 -

title('相位谱'); axis([-3,3,-1,1]); H1=(1-j*v)/(1+j*v); R1=Fe*H1; D=abs(R1);

r1=ifourier(R1,t); subplot(2,3,4); ezplot(r1,[-3,3]);

title('只满足恒定幅值'); axis([-3,3,-1,2]); subplot(2,3,5);

ezplot('4*abs(1/(4+v^2)*(1-i*v)/(1+i*v))'); title('幅度谱'); axis([-3,3,0,2]); subplot(2,3,6) i=1;

for a=-3:0.02:3

R11=subs(R1,v,a); C(i)=angle(R11); i=i+1; end

a=-3:0.02:3; plot(a,C);

title('相位谱'); axis([-3,3,-3,3]);

- 30 -

exp(-2 abs(t))21.510.5002t只满足恒定幅值21-221.510.50幅度谱10.50-0.5-20v幅度谱2-1相位谱-20相位谱221.5120-2-20v2-20200.5-20t20-1 （3）源程序：

clc clear close

syms t v;

e=exp(-2*abs(t)); subplot(2,3,1); ezplot(e,[-3,3]); axis([-3,3,-0.2,2]); Fe=fourier(e,v); subplot(2,3,2); ezplot(Fe,[-3,3]); title('幅度谱'); axis([-3,3,0,2]); i=1;

for a=-3:0.02:3

R11=subs(Fe,v,a); C(i)=angle(R11); i=i+1; end

b=-3:0.02:3; subplot(2,3,3); plot(b,C);

- 31 -

title('相位谱'); axis([-3,3,-1,1]); H1=v^2*exp(-j*v*1); R1=Fe*H1;

r1=ifourier(R1,t); subplot(2,3,4); ezplot(r1,[-3,3]); title('满足相位条件'); axis([-3,3,-4,0.2]) subplot(2,3,5); ezplot('R1');

title('幅度谱'); axis([-3,3,-3,3]); i=1;

for a=-3:0.02:3

R11=subs(R1,v,a); C(i)=angle(R11); i=i+1; end

b=-3:0.02:3; subplot(2,3,6); plot(b,C);

axis([-3,3,-3,3]); title('相位谱');

exp(-2 abs(t))21.510.5002t满足相位条件02-1-2-3-4-20t20-2-20R120-2-2022-221.510.50幅度谱10.50-0.5-20v幅度谱2-1相位谱-20相位谱2 - 32 -

（4）源程序：

clc clear close

syms t v;

e=exp(-2*abs(t)); subplot(2,3,1); ezplot(e,[-3,3]); axis([-3,3,-0.2,2]); Fe=fourier(e,v); subplot(2,3,2); ezplot(Fe,[-3,3]); title('幅度谱'); axis([-3,3,0,2]); i=1;

for a=-3:0.02:3

R11=subs(Fe,v,a); C(i)=angle(R11); i=i+1; end

b=-3:0.02:3; subplot(2,3,3); plot(b,C);

title('相位谱'); axis([-3,3,-1,1]);

H1=v^2*(1-j*v)/(1+j*v); R1=Fe*H1; D=abs(R1);

r1=ifourier(R1,t); subplot(2,3,4); ezplot(r1,[-3,3]);

title('两个条件均不满足'); axis([-3,3,-0.5,7]); subplot(2,3,5);

ezplot('4*abs(1/(4+v^2)*v^2*(1-i*v)/(1+i*v))'); title('幅度谱'); axis([-3,3,0,2]); subplot(2,3,6) i=1;

for a=-3:0.02:3

R11=subs(R1,v,a); C(i)=angle(R11); i=i+1;

- 33 -

end

a=-3:0.02:3; plot(a,C);

title('相位谱'); axis([-3,3,-3,3]);

exp(-2 abs(t))21.510.5002t两个条件均不满足61.54120-20t20.50-20v20-2-202-221.510.50幅度谱10.50-0.5-20v幅度谱2-1相位谱-20相位谱222 五、思考题

1、连续系统频域与复频域分析的基本方法是什么？

频域分析方法：将激励信号分解为正弦分量，找出联系响应与激励的系统函数H(jw)，求每一频率分量的响应，从响应的频谱函数R(jw)求傅里叶反变换从而求得响应r(t)。

复频域分析方法：通过拉普拉斯变换将时域中的积分微分方程变成复频域中的代数方程，在复频域中进行代数运算后则可得到系统响应的复频域解，将此解在经反变换则得到最终的时域解。

2、若信号经过系统不发生失真，则对系统频响有何要求？

系统频响的幅频特性在整个频率范围中为一常数，即系统具有无限频宽的响应均

- 34 -

匀的通频带，系统的相频特性是经过原点的直线。

- 35 -

实验五 信号采样与重建 一、实验目的

1、 深刻理解采样定理的内容。 2、掌握信号采样后的频谱。

二、实验设备

安装有matlab6.5 以上版本的PC 机一台。 三、实验原理

- 36 -

- 37 -

四、实验内容：

1、已知f (t) = Sa(2t)，以fs为采样频率，对f (t)进行采样得到fa(t) ，观 察随着fs 由小变大， fa(t)频谱的变化，最后得出fa(t)与f (t)两者频谱 之间的关系。 源程序：

clc close clear

syms t w;

f=sin(2*t)/(2*t); subplot(3,2,1); ezplot(f,[-5,5]); title('时域信号'); ff=fourier(f); subplot(3,2,2); ezplot(ff,[-5,5]);

title('频谱信号'); fs=1;fc=0; for n=-6:6

fc=fc+(pi/2)*fs*(-heaviside(w-2-n*2*pi*fs)+heaviside(w+2-n*2*pi*fs)); end

subplot(3,2,4);

ezplot(fc,[-40,40]);

title('采样信号频谱，采样频率1'); fs=3;fc=0; for n=-3:3

fc=fc+(pi/2)*fs*(-heaviside(w-2-n*2*pi*fs)+heaviside(w+2-n*2*pi*fs)); end

subplot(3,2,5);

ezplot(fc,[-40,40]);

title('采样信号频谱，采样频率3'); fs=5;fc=0; for n=-3:3

fc=fc+(pi/2)*fs*(-heaviside(w-2-n*2*pi*fs)+heaviside(w+2-n*2*pi*fs)); end

- 38 -

subplot(3,2,6);

ezplot(fc,[-40,40]);

title('采样信号频谱，采样频率5'); fs=0.5;fc=0; for n=-3:3

fc=fc+(pi/2)*fs*(-heaviside(w-2-n*2*pi*fs)+heaviside(w+2-n*2*pi*fs)); end

subplot(3,2,3);

ezplot(fc,[-40,40]);

title('采样信号频谱，采样频率0.5');

时域信号10.5005t采样信号频谱，采样频率0.51.510.50-4002040w采样信号频谱，采样频率3-20-51.510.50-5频谱信号05w采样信号频谱，采样频率11.510.50-4002040w采样信号频谱，采样频率5-20420-40-200w2040820-40-200w2040 2、 由实验1 中采样得到的离散信号重建对应的连续时间信号：情况 一、f s ≥ 2f m；情况二、f s < 2f m。提示根据内插公式重建连续时间信 号，或者根据采样后频谱通过低通滤波器的形式重建连续时间信号。 源程序：

clc close clear

syms t w;

f=sin(2*t)/(2*t);

- 39 -

subplot(4,2,1); ezplot(f,[-5,5]); title('时域信号'); ff=fourier(f); subplot(4,2,2); ezplot(ff,[-5,5]);

title('频谱信号'); fs=3;fc=0; for n=-3:3

fc=fc+(pi/2)*fs*(-heaviside(w-2-n*2*pi*fs)+heaviside(w+2-n*2*pi*fs)); end

subplot(4,2,3);

ezplot(fc,[-40,40]);

title('采样信号频谱，fs≥2fm(fs=3)'); G=-heaviside(w-4)+heaviside(w+4); Fc=fc*G;

f1=ifourier(Fc)/fs; subplot(4,2,5);

ezplot(Fc,[-40,40]);

title('重建信号频谱,fs≥2fm(fs=3)'); subplot(4,2,7); ezplot(f1,[-5,5]);

title('重建时域信号,fs≥2fm(fs=3)'); fs=0.5;fc=0; for n=-3:3

fc=fc+(pi/2)*fs*(-heaviside(w-2-n*2*pi*fs)+heaviside(w+2-n*2*pi*fs)); end

subplot(4,2,4);

ezplot(fc,[-40,40]);

title('采样信号频谱，fs<2fm(fs=0.5)'); G=-heaviside(w-4)+heaviside(w+4); Fc=fc*G;

f1=ifourier(Fc)/fs; subplot(4,2,6);

ezplot(Fc,[-40,40]);

title('重建信号频谱,fs<2fm(fs=0.5)'); subplot(4,2,8); ezplot(f1,[-5,5]);

title('重建时域信号,fs<2fm(fs=0.5)');

- 40 -

时域信号10.50-505t采样信号频谱，fs≥2fm(fs=3)1.510.50-5频谱信号420-4002040w重建信号频谱,fs≥2fm(fs=3)-2050-4002040w重建时域信号,fs≥2fm(fs=3)-2010.50-50x505w采样信号频谱，fs<2fm(fs=0.5)1.510.50-40-2002040w重建信号频谱,fs<2fm(fs=0.5)210-40-2002040w重建时域信号,fs<2fm(fs=0.5)1.510.50-0.5-505x 五、思考题

1、随着采样频率 fs 从小到大变化，fa(t) 的频谱发生怎样的变化，与

f (t)频谱之间的关系如何？与理论计算结果之间是否完全一致？如

果不一致，请分析可能是什么原因导致的？

随着采样频率 fs 从小到大变化，fa(t) 的频谱刚开始时出现混叠现象，之后会出现不相互叠加的频谱，并且相邻频谱之间的间隔越来越大，幅度会增大一些。

2、 采样频率f s 分别满足情况一与情况二时，由fa (t)重建的信号f (t)是否完全相同？如果不相同，试分析原因。

情况一：重建的图像衰减一定幅度得到的图像和原图像完全相同。

f(t)与原

情况二：不一样。因为fa (t) 出现混叠现象，使其解调后的图像与原图像不同。

- 41 -

实验六 离散时间信号和系统分析 一、 实验目的

1. 掌握常用离散信号的MATLAB 表示方法。 2、掌握用MATLAB 计算卷积和及零状态响应的方法。 二、 实验设备

安装有matlab6.5 以上版本的PC 机一台。 三、实验原理

- 42 -

四、实验内容

1、用MATLAB表示离散信号：ak，Asin(k)。 源程序：

clc close clear k=-3:3; fk=2.^(k);

subplot(1,2,1);

stem(k,fk,'filled'); axis([-4,4,0,9]); h=-7:7;

fh=2*sin(h); subplot(1,2,2);

stem(h,fh,'filled');

axis([-7.5,7.5,-2.5,2.5]);

9876504-0.53210-4-1-1.5-2-2024-2.5-5052.521.510.5 2、有两离散序列f1(k)={-2,-1,0,1,2}，f2(k)={1,1,1}用MATLAB 绘

k0k0出它们的波形及f1(k)+f2(k)及f1(k)*f2(k)的波形。

- 43 -

源程序：

clc close clear k=-2:2;

f1=[-2,-1,0,1,2]; h=-1:1;

f2=[1,1,1];

subplot(2,2,1);

stem(k,f1,'filled'); title('f1(k)'); subplot(2,2,2);

stem(h,f2,'filled'); title('f2(k)'); f3=[-2,0,1,2,2]; subplot(2,2,3);

stem(k,f3,'filled'); title('f1(k)+f2(k)'); f4=conv(f1,f2); subplot(2,2,4); g=-3:3;

stem(g,f4,'filled'); title('f1(k)*f2(k)');

f1(k)2100.4-1-2-20.2-10f1(k)+f2(k)210-1-2-2-1012420-2-4-4-2024120-1-0.50f1(k)*f2(k)0.5110.80.6f2(k)

- 44 -

3、已知离散序列波形如图11 所示，试用MATLAB 绘出满足下列要 求的序列波形。

(1)f(k-2)ε(k) (2)f(-k) (3)f(-k+2) (4)f(k-2)ε(k-2) 源程序：

clc close clear k=-6:5;

f=[0,1,3,6,10,15,14,12,9,5,0,0]; subplot(3,2,1);

stem(k,f,'filled'); title('f(k)'); subplot(3,2,3); h=0:7;

f1=[10,15,14,12,9,5,0,0]; stem(h,f1,'filled');

title('f(k-2)\\epsilon(k)'); subplot(3,2,4)

stem(-k,f,'filled'); title('f(-k)'); subplot(3,2,5)

stem(-k+2,f,'filled'); title('f(-k+2)'); subplot(3,2,6); h=2:7;

f1=[14,12,9,5,0,0]; stem(h,f1,'filled');

title('f(k-2)\\epsilon(k-2)');

- 45 -

f(k)151050-1020100-5f(k-2)(k)20100-52010005f(-k)024f(-k+2)680510f(k-2)(k-2)20100-505102468 4、试用MATLAB 的conv()函数计算实验2 中第1 题的结果。 源程序：

clc; close; clear;

p=0.0005;%时间间隔 t1=0:p:10;

h=exp(-2.*t1); t2=0:p:10; e=exp(-t2); y=conv(h,e); y=p*y;

c=length(t1)+length(t1)-2; t=0:p:c*p; plot(t,y);

xlabel('Time(sec)'); ylabel('i(t)');

axis([0,10,-0.025,0.275]); grid on;

- 46 -

0.250.20.15i(t)0.10.0500123456Time(sec)710 5、假设某系统的单位函数响应h(k)=(0.8)k(k)，系统激励信号 e(k) ={1,1,1,1}，求系统的零状态响应。 源程序：

clc close clear k=0:3;

h=(0.8).^k; e=[1 1 1 1]; y=conv(h,e); t=0:6;

stem(t,y,'filled'); title('零状态响应');

- 47 -

零状态响应32.521.510.500123456 - 48 -

实验总结：

- 49 -

- 50 -