图Q3-2
你编写的程序Q3_2抄写如下:
clc, close all T =0.01; dw =0.1; t=0:T:40;
w=-4*pi:dw:4*pi; b=input('b='); a = input('a=');
x = input('Type in the expression of the input signal x(t)'); X=x*exp(-j*t'*w)*T; X1=abs(X); y=lsim(b,a,x,t);
Y=y'*exp(-j*t'*w)*T; Y1=abs(Y);
h=impulse(b,a,40); [H,w]= freqs(b,a); Hm= abs(H); subplot(324) plot(w,Hm),
axis([0 4*pi -0 1.2]); grid on,
title('Magnitude response of the system') subplot(321) plot(t,x),
axis([0 40 -3 3]);
grid on,
title('Input sihnal x(t)') subplot(323) impulse(b,a,40), axis([0 40 -0.2 1]); grid on,
title('Impulse response h(t) of the system') subplot(325) lsim (b,a,x,t), axis([0 40 -1 1]);
grid on,title('Output signal y(t)') w=-4*pi:dw:4*pi; subplot(322) plot(w,X1),
axis([-4*pi 4*pi 0 20]); grid on,
title('Magnitude response of input signal x(t)') subplot(326) plot(w,Y1),
axis([-4*pi 4*pi 0 20]); grid on,
title('Magnitude response of output signal y(t)'); xlabel('Frequency in rad/sec')
执行程序Q3_2,输入信号x(t) = sin(t) + sin(8t),输入由Eq.3.3描述的系统。得到的图形如下:
在键盘上输入
请手工绘制出信号x(t) = sin(t) + sin(8t) 的幅度频谱图如下:
你手工绘制的信号x(t) = sin(t) + sin(8t) 的幅度频谱图与执行程序Q3_2得到的x(t) = sin(t) + sin(8t) 的幅度频谱图是否相同?如不同,是何原因造成的?
答:不相同。手工绘制的幅度频谱图是无法做到那么密集的取值,所以看起来比
较离散,而用程序绘制的图比较密集。
执行程序Q3_2得到的x(t) = sin(t) + sin(8t) 的幅度频谱图实际上是另外一个信号x1(t)的幅度频谱,这个信号的时域数学表达式为
X1(t)=
请利用傅里叶变换的相关性质计算并绘制信号x1(t)的幅度频谱图。 计算过程:x1(t)=sin(t)+sin(8t)
由sin(w0t)----->πj[ょ(w+w0)-ょ(w-w0)]
所以F(x1(t))=πj [ょ(w+1)-ょ(w-1) +ょ(w+8)-ょ(w-8)] % ょ(t)为冲激函数
结合所学的有关滤波的知识,根据上面所得到的信号的时域和频域图形,请从时域和频域两个方面解释滤波的概念。 答:
四、实验结论与体会
本次试验中,在做第一个试验时遇到了困难,在绘制系统的频率图时,遇到了麻烦,在查书和查看课本以前试验没有解决问题后,通过查资料和问同学,终于解决国问题。 。
实验四 通信系统仿真
一、实验目的
1、理解信号的抽样及抽样定理以及抽样信号的频谱分析; 2、掌握和理解信号抽样以及信号重建的原理; 3、掌握傅里叶变换在信号调制与解调中的应用。
基本要求:掌握并理解“抽样”定理及其重要意义,理解抽样信号的频谱特征。一般理解信号重建的物理过程以及内插公式所描述的信号重建原理。理解调制与解调的基本概念,理解信号调制过程中的频谱搬移。掌握利用MATLAB仿真正弦幅度调制与解调的方法。
二、实验要求
学会利用MAYLAB求解信号与系统的相关问题;
三、实验内容
Q4-1 在1/2—1/10之间选择若干个不同Ts值,反复执行执行范例程序
Program4_1,保存执行程序所得到的图形。
Ts = 1/2时的信号时域波形和频谱
Ts = 1/4时的信号时域波形和频谱图
Ts = 1/8时的信号时域波形和频谱图