答案-信号与系统实验报告 下载本文

大连理工大学 信息与通信工程学院

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页)