南邮信课程实验完整版
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
南邮信课程实验集团标准化办公室:[VV986T-J682P28-JP266L8-68PNN]
课程实验报告
( 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');
00.51 1.52 2.53
00.5
1
x 1(t )
00.51 1.52 2.53
00.5
1
x 2(t )
00.51 1.52 2.53
0.5
1
y (t )=x 1*x 2
修改后
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');
5
00.20.40.60.8
1x 1(t )
05
0.511.5
2x 2(t )
05
0.511.5
2y (t )=x 1*x 2
上机题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');
024********
00.5
1
x (t )
024********
00.5
1
h (t )
024********
0.51
---->t/s
y =x *h
(二)信号的频域分析
上机题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);