数字信号处理实验报告MATLAB第二章

数字信号处理MATLAB

第二章离散时间系统的时域分析

例2.1 滑动平均系统

%程序P2_1

%一个M点滑动平均滤波器的仿真

%产生输入信号

clf;

n=0:100;

s1=cos(2*pi*0.05*n);%一个低频正弦

s2=cos(2*pi*0.47*n);%一个高频正弦

x=s1+s2;

%滑动平均滤波器的实现

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

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('振

实 验 报 告

课程名称： 数字信号处理实验 专业班级： 姓 名： 学 号：

实验名称 数字信号处理 实验地点 实验成绩 一、实验目的及任务 ? 学会运用MATLAB求解离散时间系统的零状态响应； ? 学会运用MATLAB求解离散时间系统的单位取样响应； ? 学会运用MATLAB求解离散时间系统的卷积和。 实验时间 二、实验内容与步骤 3.2.1 离散时间系统的响应 离散时间LTI系统可用线性常系数差分方程来描述，即 ?ay(n?i)??bx(n?j) （3-1） iji?0j?0NM其中，ai（i?0，1，?，N）和bj（j?0，1，?，M）为实常数。 MATLAB中函数filter可对式（13-1）的差分方程在指定时间范围内的输入序列所产生的响应进行求解。函数filter的语句格式为 y=filter(b,a,x) 其中，x为输入的离散序列；y为输出的

课程名称：数字信号处理实验

学号：实验地点：综合楼C407

专业班级：2014级生物医学工程

姓名：leifeng

指导老师：

第一次实验

第一章离散时间信号的时域分析

Q1.1运行程序P1.1，以产生单位样本序列u[n]并显示它

clf; n=-10:20;

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

xlabel('时间序号');ylabel('振幅'); title('单位样本序列'); axis([-10 20 0 1.2]);

单位样本序列10.8振幅0.60.40.20-10-505时间序号101520

Q1.2 命令clf，axis，title，xlabel和ylabel的作用是什么

clf：清除图形窗口内容； axis:规定横纵坐标的范围；

title：使图像面板上方显示相应的题目名称； xlable：定义横坐标的名字； ylable：定义纵坐标的名字。

Q1.3修改程序P1.1以产生带有延时11个样本的延迟单位样本序列ud[n]，运行修改的程序并

且显示产生的序列。

clf; n=0:30;

u=[zeros(1,11) 1 zeros(1,19)]; stem(n,u);

实 验 报 告

课程名称： 数字信号处理实验 专业班级： 姓 名： 学 号：

实验名称 数字信号处理 实验地点 实验成绩 一、实验目的及任务 ? 学会运用MATLAB求解离散时间系统的零状态响应； ? 学会运用MATLAB求解离散时间系统的单位取样响应； ? 学会运用MATLAB求解离散时间系统的卷积和。 实验时间 二、实验内容与步骤 3.2.1 离散时间系统的响应 离散时间LTI系统可用线性常系数差分方程来描述，即 ?ay(n?i)??bx(n?j) （3-1） iji?0j?0NM其中，ai（i?0，1，?，N）和bj（j?0，1，?，M）为实常数。 MATLAB中函数filter可对式（13-1）的差分方程在指定时间范围内的输入序列所产生的响应进行求解。函数filter的语句格式为 y=filter(b,a,x) 其中，x为输入的离散序列；y为输出的

数字信号处理

第二次实验报告

学 院：信息工程学院

班 级：2012级电子信息工程*班 姓 名：

学 号：20125507** 指导老师：

实验四：IIR数字滤波器设计及软件实现

一、实验目的

1、熟悉双线性变换设计IIR滤波器的原理与方法 2、掌握IIR滤波器的MATLAB实现方法 二、实验原理简述

IIR数字滤波器间接法基本设计过程：

1、将给定的数字滤波器的指标转换成过渡模拟滤波器的指标； 2、设计过渡模拟滤波器；

3、将过渡模拟滤波器系统函数转换成数字滤波器的系统函数 三、程序与图形

1、%-----------------信号产生函数mstg--------------- function st=mstg %功能函数的写法

%产生信号序列向量st,并显示st的时域波形和频谱

%st=mstg 返回三路调幅信号相加形成的混合信号，长度N=1600 N=1600 %N为信号st的长度。

Fs=10000;T=1/Fs;Tp=N*T; %采样频率Fs=10kHz，Tp为采样时间 t=0:T:(N-1)*T;k=0:N-1;f=k/Tp; fc1=Fs/10; %第1路调幅信号的载波频率fc1=1000Hz,

fm1=fc1/1

数字信号处理软件实验 ——MatLab仿真实验报告

学院:电子工程学院 班级：2013211202 姓名： 学号：

实验一：数字信号的 FFT 分析

1、实验内容及要求

(1) 离散信号的频谱分析：

? 设信号 x(n)?0.001*cos(0.45n?)?sin(0.3n?)?cos(0.302n??)4 此信号的0.3pi 和 0.302pi两根谱线相距很近，谱线 0.45pi 的幅度很小，请选择合适的序列长度 N 和窗函数，用 DFT 分析其频谱，要求得到清楚的三根谱线。 (2) DTMF 信号频谱分析

用计算机声卡采用一段通信系统中电话双音多频（DTMF）拨号数字 0～9的数据，采用快速傅立叶变换（FFT）分析这10个号码DTMF拨号时的频谱。

2、实验目的

通过本次实验，应该掌握：

(a) 用傅立叶变换进行信号分析时基本参数的选择。

(b) 经过离散时间傅立叶变换（DTFT）和有限长度离散傅立叶变换（DFT） 后信号频谱上的区别，前者 DTFT 时间域是离散信号，频率域还是连续的，而 DFT 在两个域中都是离散的。 (c) 离散傅立叶变换的基本原理、特性，以

1

第二章

2.1 判断下列序列是否是周期序列。若是，请确定它的最小周期。 （1）x(n)=Acos(6

85ππ+n ) （2）x(n)=)8(

π-n

e j

（3）x(n)=Asin(3

43π

π+n )

解 (1)对照正弦型序列的一般公式x(n)=Acos(?ω+n )，得出=ω85π。因此5

16

2=ωπ是有理数，所以是周期序列。最小周期等于N=

)5(165

16

取k k =。 （2）对照复指数序列的一般公式x(n)=exp[ωσj +]n,得出8

1

=ω。因此πωπ162=是无理数，所以不

是周期序列。

（3）对照正弦型序列的一般公式x(n)=Acos(?ω+n )，又x(n)=Asin(343ππ+n )＝Acos(-2π3

43ππ-n )＝Acos(6143-n π)，得出=ω43π。因此3

8

2=ωπ是有理数，所以是周期序列。最小周期等于

N=

)3(83

8

取k k =

2.2在图2.2中，x(n)和h(n)分别是线性非移变系统的输入和单位取样响应。计算并列的x(n)和h(n)的线性卷积以得到系统的输出y(n)，并画出y(n)的图形。

(a)

1

11

1

(b)

(c)

11

111

0 0

-1-1

-1

-1

-1

-1

-1

-1

2

2

2

222 3

3

3

3 34

44

…

…

…n

n

n

中北大学

实验报告

课 程 名： 数字信号处理I 任课教师： 陈平 专 业： 信息与计算科学 学 号： 1408024111 姓 名： 张冉

实验一 采样定理

一、实验内容

给定信号为x(t)?exp(?at)cos(100*?*at)，其中a为学号， （1）确定信号的过采样和欠采样频率

（2）在上述采样频率的条件下，观察、分析、记录频谱，说明产生上述现象的原因。 二、基本要求

验证采样定理，观察过采样和欠采样后信号的频谱变化。 三、实验结果 (1)过采样频率：

a=11; dt=0.0009; t=0:dt:0.05;

x1=exp(-a*t).*cos(100*pi*a*t); N=length(x1); k=0:(N-1); Y1=fft(x1); Y1=fftshift(Y1); subplot(2,1,1); plot(t,x1);hold on; stem(t,x1,'o'); subplot(2,1,2); plot(k,abs(Y1)); gtext('1408024111张冉');

(2)欠采样频率：

a=11; dt=

数字信号处理 实验报告

实验一 序列的傅立叶变换

一、实验目的

1.进一步加深理解DFS,DFT算法的原理; 2.研究补零问题;

3.快速傅立叶变换（FFT）的应用。 二、 实验步骤

1.复习DFS和DFT的定义，性质和应用；

2熟悉MATLAB语言的命令窗口、编程窗口和图形窗口的使用；3利用提供的程序例子编写实验用程序；4.按实验内容上机实验，并进行实验结果分析；5.写出完整的实验报告，并将程序附在后面。 三、 实验内容

1.周期方波序列的频谱

试画出下面四种情况下的的幅度频谱, 并分析补零后，对信号频谱的影响。

x(n)?cos(0.48?n)?cos(0.52?n)2.有限长序列x(n)的DFT（1）取x(n)(n=0:10)时，画出x(n)的频谱X(k) 的幅度；（2）将（1）中的x(n)以补零的方式，使x(n)加长到(n:0~100)时，画出x(n)的频谱X(k) 的幅度；

（3）取x(n)(n:0~100)时，画出x(n)的频谱X(k) 的幅度。利用FFT进行谱分析

x(t)?2sin(4?t)?5cos(8?t)3.已知：模拟信号

以t=0.01n(n=0:N-1)进行采样，求N点DFT的幅值谱。 请分别画出N=45; N=50

四川大学电气信息学院 数字信号处理实验报告

实验二 时域采样与频域采样

1. 实验结果和分析 （1）时域采样

(a)Fs=1000Hz2001000(a) FT[xa(nT)],Fs=1000Hzx1(n)0-2000204060幅度500005001000f(Hz)(b) FT[xa(nT)],Fs=300Hzn(b)Fs=300Hz200400x2(n)0-200051015幅度20000100200300n(c)Fs=200Hz200200f(Hz)(c) FT[xa(nT)],Fs=200Hzx3(n)0-20005n10幅度1000050100f(Hz)150200分析：时域采样定理：1、对模拟信号以间隔T进行时域等间隔理想采样，形成的采样信号的频谱是原模拟信号频谱以采样角频率为周期进行周期延拓。2、采样频率必须大于等于模拟信号最高频率的两倍以上，才能使采样信号的 频谱不产生频谱混叠。

由图可见，左边在时域上的采样频率逐渐降低，右边所对应的频域图样的混叠情况由微弱变得越来越大。

（2）频域采样

(a)FT[x(n)]20020(b) 三角波序列x(n)|X(ej?)|100000.5?/?(c) 16点频域采样2001x(n)100