南邮信课程实验完整版

南邮信课程实验集团标准化办公室：[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

x 1(t )

00.51 1.52 2.53

00.5

x 2(t )

00.51 1.52 2.53

0.5

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');

00.20.40.60.8

1x 1(t )

0.511.5

2x 2(t )

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

x (t )

024********

00.5

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

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

xlabel('n');grid

-10-8-6-4-20246810

五、实验思考题

对于上机题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

-8-6-4-202468

六、实验小结

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

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

七、教师评阅

日期:

南邮广播电视工程数字视频非线性编辑制作课程设计实验报告定稿版

南邮广播电视工程数字视频非线性编辑制作课程设计实验报告精编 W O R D版 IBM system office room 【A0816H-A0912AAAHH-GX8Q8-GNTHHJ8】

通信与信息工程学院 2016 / 2017 学年第一学期课程设计实验报告模块名称数字视频非线性编辑制作专业广播电视工程学生班级 B130114 学生学号学生姓名陈超指导教师姚锡林日期： 2016 年 11 月 21 日摘要本次课程设计利用软件premiere进行数字视频非线性编辑制作。本文首先就本次实验主题归纳总结电视节目制作一般流程方法，接着对此次课程设计主要软件工具进行系统介绍，主要涉及实验相关借本操作的详细阐述；接下来两大章节部分先从取材、构思角度详细分析此次课程设计所做的主题内容与规划，并以此为指导再从具体操作上分步骤、多角度实现视频序列的制作；最后对本次课程设计的体味与收获进行思考。此次作品《再次出发》电影鉴赏是将导演约翰卡尼的一部经典音乐影视作品利用premiere软件，在制作的过程中添加了转场特效，关键帧，字幕，音频等功能，并运用多种素材，重新删减编辑，形成一部情节连贯，内容完整、主题明确的电影鉴赏短片。短片的片长时间为9分40秒，大小为720*576，AVI格式，PAL制式（48Khz）。

关键词：数字视频非线性编辑制作；premiere；视频制作；

目录第一章电视节目制作的一般流程与方法 (1) 1.1 电视节目制作一般流程 (1) 1.1.1 前期制作流程 (1) 1.1.2 后期制作工作流程 (2) 1.2 电视节目制作的一般方法 (2) 1.2.1 ENG方式 (2) 1.2.2 EFP方式 (2) 1.2.3 ESP方式 (3) 第二章 Premiere的功能介绍及操作方法 (4) 2.1 Premiere概述 (4) 2.1.1 概述 (4) 2.1.2 基本操作界面 (4) 2.2 Premiere的基本操作 (5) 2.2.1 新建项目 (5) 2.2.2 新建序列 (6)

数学实验答案-1

1.（1） [1 2 3 4;0 2 -1 1;1 -1 2 5;]+(1/2).*([2 1 4 10;0 -1 2 0;0 2 3 -2]) 2. A=[3 0 1;-1 2 1;3 4 2],B=[1 0 2;-1 1 1;2 1 1] X=(B+2*A)/2 3. A=[-4 -2 0 2 4;-3 -1 1 3 5] abs(A)>3 % 4. A=[-2 3 2 4;1 -2 3 2;3 2 3 4;0 4 -2 5] det(A),eig(A),rank(A),inv(A) 求计算机高手用matlab解决。 >> A=[-2,3,2,4;1,-2,3,2;3,2,3,4;0,4,-2,5] 求|A| >> abs(A) ans = （ 2 3 2 4 1 2 3 2 3 2 3 4 0 4 2 5 求r（A） >> rank(A) ans =

4 求A-1 《 >> A-1 ans = -3 2 1 3 0 -3 2 1 2 1 2 3 -1 3 -3 4 求特征值、特征向量 >> [V,D]=eig(A) %返回矩阵A的特征值矩阵D 与特征向量矩阵V , V = - + + - - + - + - + - + D = { + 0 0 0 0 - 0 0 0 0 + 0 0 0 0 - 将A的第2行与第3列联成一行赋给b >> b=[A(2,:),A(:,3)'] b = 《 1 - 2 3 2 2 3 3 -2

1． a=round(unifrnd(1,100)) i=7; while i>=0 i=i-1; b=input('请输入一个介于0到100的数字：'); if b==a ￥ disp('You won!'); break; else if b>a disp('High'); else if b