2024年南邮dsp实验报告.doc

南京邮電大學

实验报告

试验名称：离散時间信号与系统的時、频域表达

离散傅立叶变换和z变换

数字滤波器的频域分析和实現

数字滤波器的设计

課程名称：数字信号处理A(双語)

班级學号：B1

姓名：

開課時间：/年第二學期

试验一

试验名称：离散時间信号与系统的時、频域表达

试验目的：

熟悉Matlab基本命令和信号处理工具箱，加深理解和掌握离散時间信号与系统的時、频域表达及简朴应用。

试验任务：

在Matlab环境中，按照规定产生序列，對序列進行基本运算；對简朴离散時间系统進行仿真，计算线性時不变（LTI）系统的冲激响应和卷积输出；计算和观测序列的离散時间傅立叶变换（DTFT）幅度谱和相位谱。

试验内容：

基本序列产生和运算：Q1.1～1.3，Q1.23，Q1.30～1.33

离散時间系统仿真：Q2.1～2.3

LTI系统：Q2.19，Q2.21，Q2.28

DTFT：Q3.1，Q3.2，Q3.4

试验過程描述：

Q1.1

程序：

clf

n=-10:20;

u=[zeros(1,10)1zeros(1,20)];

stem(n,u);

xlabel(時间序列n);ylabel(振幅);

title(單位样本序列);

axis([-102001.2]);

显示的波形如下:

Q1.2

clf：清除图形；

axis：设置坐標轴范围、可讀比例等；

title：給图形加標題；

xlable：給x轴加標注；

ylable：給y轴加標注。

Q1.3

程序：

clf

n=-10:20;

u=[zeros(1,10)1zeros(1,20)];

stem(n+11,u);

xlabel(時间序列n);ylabel(振幅);

title(單位样本序列);

axis([03201.2]);

显示的波形如下:

Q1.23

程序：

n=0:50;

f=0.08;

phase=pi/2;

A=2.5;

arg=2*pi*f*n-phase;

x=A*cos(arg);

clf;

stem(n,x);

axis([050-33]);

grid;

title(正弦序列);

xlabel(時间序列n);

ylabel(振幅);

axis;

显示的波形如下:

Q1.30

加性噪声d[n]是均匀分布在-0.4和+0.4之间的随机序列

Q1.31

不能。由于d是列向量，s是行向量

Q1.32

x1是x的延時，x2和x相等，x3超前于x

Q1.33

legend用于产生图例阐明

Q1.30

未污染的信号s[n]是什么样的形式？加性噪声d[n]是什么样的形式？

答:未污染的信号s[n]:是线性增長伴伴随实指数缓慢衰減的图像

加性噪声d[n]:在-0.4和+0.4间均匀分布的自由序列

Q1.31

使用語句s=s+d能产生被噪声污染的信号吗？若不能，為何？

答:不能，由于-d是一种列向量,而s是一种行向量,需要在添加它們之前调换其中一种向量。

Q1.32

信号x1、x2、x3与x之间的关系是什么？

答:這三個信号x1,x2,和x3是x扩展的版本,左右边各一种附加的采样。x1是x延迟的版本,一种样本转移到右边并且左边补零。信号x2等于x,左右补0来填充多出的長度。最终,x3是x時间提前的版本,转移一种样本到右边，左边补0。

Q1.33

legend的作用是什么

答:thelegend命令的目的——创立图表的阐明。在P1_5,信号绘制使用不一样的颜色和线类型;阐明哪种颜色信息和行类型与每個信号有关联。

Q2.1

程序：

clf;

n=0:100;

s1=cos(2*pi*0.05*n);

s2=cos(2*pi*0.47*n);

x=s1+s2;

M=input(滤波器所需的長度=);

滤波器所需的長度=2

num=ones(1,M);

y=filter(num,1,x)/M;

subplot(2,2,1);

plot(n,s1);

axis([0,100,-2,2]);

xlabel(時间序列n);ylabel(振幅);

title(信号#1);

subplot(2,2,2);

plot(n,s2);

axis([0,100,-2,2]);

xlabel(時间序列n);ylabel(振幅);

title(信号#2);

subplot(2,2,3);

plot(n,x);

axis([0,100,-2,2]);

xlabel(時间序列n);ylabel(振幅);

title(输入信号