大连理工大学 信息与通信工程学院
tao=2;T=10;w3=2*pi/T; c=n.^2;
x=n*pi*tao/(2*T); d=sin(x); e=d.^2;
运行结果如下：
tao=1,T=10.5fn=8*e./(tao*c*4*pi*pi/T); subplot(414)
stem(n*w3,fn),grid on; title('tao=2,T=10'); hold on
stem(0, 0.1);
0-200.05-15-10-50tao=1,T=1051015200-200.10.050-200.10.050-20-15-10-50tao=1,T=55101520-15-10-50tao=2,T=105101520-15-10-505101520
从图中可以看出，脉冲宽度τ 越大，信号的频谱带宽越小；而周期越小，谱线之间间隔越大.
3. 试用MATLAB命令求下列信号的傅里叶变换，并绘出其幅度谱和相位谱。
解：
调试程序如下：
ft1=sym('sin(2*pi*(t-1))/(pi*(t-1))'); ft2=sym('(sin(pi*t)/(pi*t))^2'); Fw1=fourier(ft1); Fw2=fourier(ft2); subplot(411);
ezplot(abs(Fw1));grid on; title('f1幅度谱');
phase=atan(imag(Fw1)/real(Fw1)); subplot(412);
ezplot(phase);grid on; title('f1相位谱'); subplot(413);
ezplot(abs(Fw2));grid on; title('f2幅度谱');
phase=atan(imag(Fw2)/real(Fw2)); subplot(414); ezplot(phase);grid on; title('f2相位谱');
运行结果如下：
第 4页 /（共 28页）
1-6-4-20wf1相位谱24610-1-6-4大连理工大学 信息与通信工程学院
-2f1幅度谱111-6-4-20wf1相位谱24610-1-6-4-210.50-6-4-20wf2幅度谱2460wf2相位谱24610-1-6-4-20wf2幅度谱2460w246
4. 试用MATLAB命令求下列信号的傅里叶反变换，并绘出其时域信号图。
100解： w调试程序如下： f2相位谱1-6-4-22460.5
clear; syms t1; 0syms omega; -1Fw1 = -6-4-20246fourier((10/(3+j*omega))-(4/(5+j*omega))); wft1 = ifourier(Fw1,t1); syms t2;
Fw2 = fourier(exp(-4*omega^2));
ft2 = ifourier(Fw2,t2);
subplot(211);
ezplot(t1,ft1);grid on; title('f1时域信号'); subplot(212);
ezplot(t2,ft2);grid on; title('f2时域信号');
f1时域信号2.52y1.510.50123xf2时域信号4561.510.50-0.50123x456y
第 5页 /（共 28页）
大连理工大学 信息与通信工程学院
5. 试用MATLAB数值计算方法求门信号的傅里叶变换，并画出其频谱图。门信号即
dt = 0.005; t = -5:dt:5; y1 = [t>=-0.5]; y2 = [t>=0.5]; ft = y1 - y2; N = 2000; k = -N:N;
W = 2*pi*k/((2*N+1)*dt); F = dt * ft*exp(-j*t'*W); plot(W,F), grid on; xlabel('W'), ylabel('F(W)'); axis([-20*pi 20*pi -0.3 1.2]); title('频谱图');
频谱图10.80.6F(W)0.40.20-0.2-60-40-200W204060
第 6页 /（共 28页）
大连理工大学 信息与通信工程学院
实验2 连续时间系统分析
一、实验目的
1 建立系统的概念；
2 掌握连续时间系统的单位冲激响应的求解； 3 掌握连续时间系统单位阶跃响应的求解； 4 掌握连续时间系统零极点的求解； 5 分析系统零极点对系统幅频特性的影响； 6 分析零极点对系统稳定性的影响； 7 介绍常用信号处理的MATLAB工具箱；
二、实战演练
1. 已知系统的微分方程为
1.41.210.8h(t)0.60.40.20012345time678910y''(t)+3yⅱ(t)+2y(t)=x(t)+4x(t)，计算该系统的单位冲激响应和单位阶跃响应。
单位冲激响应：
21.81.61.41.2a=[1 3 2];
b=[1 4]; sys=tf(b,a); t=0:0.1:10;
y=impulse(sys,t); plot(t,y);
xlabel('time'); ylabel('h(t)');
单位阶跃响应：
h(t)*u(t)10.80.60.40.20012345time678910第 7页 /（共 28页）