南邮信号课程实验

南邮信号课程实验 TYYGROUP system office room 【TYYUA16H-TYY-TYYYUA8Q8-

课程实验报告

（ 2014 / 2015 学年第 1 学期）

实验名称：

用Matlab进行信号与系统的时频域分析

专业电子信息工程

学生姓名

班级学号

指导教师

指导单位通信与信息工程学院

日期 2014 年 11月 6 日

一、实验目的

1.了解连续信号、离散信号的波形特点；

2.掌握连续信号、离散信号的Matlab实现；

3.熟悉Matlab中plot、stem等函数的应用；

4.掌握利用matlab函数表示常见信号波形。

二、实验任务

1.用matlab软件画出连续时间信号的卷积的波形并修改其时间区间改成不同脉宽观察其卷积波形；

2.用matlab软件计算卷积结果并画出卷积波形；

3.用matlab软件画出矩形脉冲频谱图并修改其时间常数和脉宽和幅度，观察比较实验结果。

三、主要仪器设备

硬件：微型计算机

软件：Matlab?

四、实验内容

（一）连续时间信号的卷积

上机题1.已知两个信号)2()1()(1---=t t t x εε和)1()()(2--=t t t x εε，试分别画出)(),(21t x t x 和卷积)()()(21t x t x t y *=的波形。

（上机原程序及所画出的波形图）

T=;

t1=1;t2=2;

t3=0;t4=1;

t=0:T:t2+t4;

x1=ones(size(t)).*((t>t1)-(t>t2));

x2=ones(size(t)).*((t>t3)-(t>t4));

y=conv(x1,x2)*T;

subplot(3,1,1),plot(t,x1);

ylabel('x1(t)');

subplot(3,1,2),plot(t,x2);

ylabel('x2(t)');

subplot(3,1,3),plot(t,y(1:(t2+t4)/T+1));

ylabel('y(t)=x1*x2');

xlabel(----→t/s');

修改后

T=;

t1=2;t2=4;

t3=0;t4=1;

t=0:T:t2+t4;

x1=ones(size(t)).*((t>t1)-(t>t2));

x2=2*ones(size(t)).*((t>t3)-(t>t4));

y=conv(x1,x2)*T;

subplot(2,5,1),plot(t,x1);

ylabel('x1(t)');

subplot(2,5,2),plot(t,x2);

ylabel('x2(t)');

subplot(2,5,3),plot(t,y(1:(t2+t4)/T+1));

ylabel('y(t)=x1*x2');

xlabel(----→t/s');

上机题2.已知两个信号)()(t e t x t ε-=和)()(2/t te t h t ε-=，试用数值计算法求卷积，并分别画出)(),(t h t x 和卷积)()()(t h t x t y *=的波形。

（上机原程序及所画出的波形图）

syms tao ;

t=sym( 't','positive');

xt=exp(-t);

ht=t*exp(-t/2);

xh_tao=subs(xt,t,tao)*subs(ht,t,t-tao);

yt=int(xh_tao,tao,0,t)

源程序及波形图

t2=3;t4=11;

T=;

t=0:T:t2+t4;

x=exp(-t).*((t>0)-(t>t2));

h=t.*exp(-t/2).*((t>0)-(t>t4));

y=conv(x,h)*T;

yt=4*exp(-t)+2*t.*exp(-1/2*t)-4*exp(-1/2*t);

subplot(3,1,1),plot(t,x);

ylabel('x(t)');

subplot(3,1,2),plot(t,h);

ylabel('h(t)');

subplot(3,1,3),plot(t,y(1:(t2+t4)/T+1),t,yt,'--r');

legend('by numerical','Theoretical');

ylabel('y=x*h');

xlabel('---->t/s');

（二）信号的频域分析

上机题3.求周期矩形脉冲信号的频谱图，已知s T s A 5.0,1.0,1===τ

a=1;tao=;t=;

n0=t/tao;

n=0:2*n0;

fn_p=a*tao/t*(sin(n*pi*tao/t+eps*(n==0)))./(n*pi*tao/t+eps*(n==0));

fn_pabs=abs(fn_p);

fn_pang=angle(fn_p);

fn_mabs=fliplr(fn_pabs(2:11));

fn_mang=-fliplr(fn_pang(2:11));

fnabs=[fn_mabs fn_pabs];

fnang=[fn_mang fn_pang];

subplot(2,1,1),stem((-2*n0:2*n0),fnabs);

text(4,,'amplitude spectrum');

subplot(2,1,2),stem((-2*n0:2*n0),fnang);

text(-2,2,'phase spectrum');

xlabel('n');grid

五、实验思考题

对于上机题3，改变信号的周期，比较他们的计算结果。

a=8;tao=;t=;

n0=t/tao;

n=0:2*n0;

fn_p=a*tao/t*(sin(n*pi*tao/t+eps*(n==0)))./(n*pi*tao/t+eps*(n==0)); fn_pabs=abs(fn_p);

fn_pang=angle(fn_p);

fn_mabs=fliplr(fn_pabs(2:9));

fn_mang=-fliplr(fn_pang(2:9));

fnabs=[fn_mabs fn_pabs];

fnang=[fn_mang fn_pang];

subplot(2,1,1),stem((-2*n0:2*n0),fnabs);

text(4,,'amplitude spectrum');

subplot(2,1,2),stem((-2*n0:2*n0),fnang);

text(-2,2,'phase spectrum');

xlabel('n');grid

六、实验小结

这次信号与系统实验用Matlab进行信号与系统的时频域分析，我再次熟悉了matlab软件的使用方法，并且通过此了解了连续信号、离散信号的波形特点，

熟悉Matlab中plot、stem等函数的应用，掌握利用matlab函数表示常见信号波形。在这次实验中，我也遇到了一些困难，在解决这些困难的同时也加深了我对信号与系统课程中许多知识的理解与掌握，同时我也在向同学请教的过程中明白了合作的重要性，最终终于在老师的悉心不厌其烦的讲解与同学帮忙查找错误后圆满的完成了这次实验，我不但再次熟悉了matlab软件的使用方法，还加深理解了信号与系统中的许多知识，通过这次实验使我受益匪浅。

七、教师评阅

日期: