当前位置：文档之家 › 用MATLAB分析离散信号的频谱与信号的采样

用MATLAB分析离散信号的频谱与信号的采样

  • 格式：doc
  • 大小：476.50 KB
  • 文档页数：8

.

精品 实验六 用MATLAB分析离散信号的频谱与信号的采样

一、 实验目的

1、了解离散时间信号频谱的分析方法；

2、了解相关函数的调用格式及作用；

3、掌握用MATLAB分析信号的采样过程与原理。

二、涉及的MATLAB函数

1、fft函数：可用来计算离散周期信号频谱

X[m] = fft(x)

x：是离散周期信号0~N-1 一个周期的序列值

X[m] 是离散周期信号的频谱

函数fft还可用来计算离散非周期信号频谱、连续周期信号和连续非周期信号的频谱。

2、rectpuls函数：表示矩形脉冲信号

y=rectpuls(t,width)

产生宽度为0.4,幅度为1,以零点对称的矩形波1P(t)

三、实验内容

1、用MATLAB实现下图所示周期矩形序列的频谱

x[k]的频谱函数为：X[m]=)(sin)]12([sinNmMNm x[k]k0MM.

精品 %Program 6_1计算离散周期矩形序列的频谱

N=32; M=4; %定义周期矩形序列的参数

x=[ones(1,M+1),zeros(1,N-2*M-1),ones(1,M)]; %产生周期矩形序列

X=fft(x); %计算DFS系数

m=0:N-1;

stem(m,real(X)); %画出频谱X的实部

title('X[m]的实部');xlabel('m')

figure;

stem(m,imag(X)); %画出频谱X的虚部

title('X[m]的虚部');xlabel('m');

xr=ifft(X);

figure;

stem(m,real(xr));

xlabel('k');

title('重建的x[k]');

仿真的结果如下： .

精品

.

精品

2、用MATLAB实现脉冲采样抽样函数的过程。

)()()(tptStfa，要求：求出抽样函数Sa(t)、矩形脉冲波P(t)与抽样结果f(t)的波形。其中)(tp的波形如下：

p(t)-0.20.20.40.811.4-0.4-0.8-1-1.4t0…………1

程序如下：

t=-3*pi:0.01:3*pi; % 定义时间范围向量

s=sinc(t/pi); % 计算Sa(t)函数

subplot(3,1,1),plot(t,s); % 绘制Sa(t)的波形

title('抽样函数Sa(t)');

p=zeros(1,length(t)); % 预定义p(t)的初始值为0

for i=16:-1:-16

p=p+rectpuls(t+0.6*i,0.4); % 利用矩形脉冲函数rectpuls的平移来产生宽度为0.4,幅度为1的矩形波p(t)

end

subplot(3,1,2),stairs(t,p); % 用阶梯图形表示矩形脉冲

title('矩形脉冲');

axis([-10 10 0 1.2]);

f=s.*p;

subplot(3,1,3),plot(t,f); % 绘制f(t)=Sa(t)*p(t)的波形

xlabel('t');

ylabel('Sa(t)*p(t)');

仿真结果如下：.

精品

3、用MATLAB分析如图所示三角波信号的采样过程。

低通4c)(tf)(1ty)(tT)(ty

f(t)t0-0.50.51…………-TsTs0tTs=0.2s)(tT||21)(ttf.

精品 要求：

（1）、 求)(tf的频谱图)(F

（2）、求)(1ty的频谱图)(1Y

（3）、求)(ty的频谱图)(Y

)4(21))21cos(1(4)(22aSF

Matlab程序如下：

syms t w f; % 定义符号变量

f=(1-2*abs(t))*exp(-j*w*t); % 计算被积函数

F=int(f,t,-1/2,1/2); % 计算傅立叶系数F(w)

F=simple(F);F % 化简

subplot(3,1,1), % 绘制三角波的幅频特性曲线F(w)

low=-26*pi;high=-low; % 设置w的上界和下界

ezplot(abs(F),[low:0.01:high]);

axis([low high -0.1 0.5]); xlabel('');

title('三角波的频谱');

subplot(3,1,2), % 绘制经过截止频率为4*pi低通滤波器后的频谱Y1(w)

ezplot(abs(F),[-4*pi:0.01:4*pi]);

axis([low high -0.1 0.5]); xlabel('');

title('低通滤波后的频谱');% 采样信号的频谱是原信号频谱的周期延拓，延拓周期为(2*pi)/Ts

% 利用频移特性F[f(t)*exp(-j*w0*t)]=F(w+w0)来实现

subplot(3,1,3); % 绘制采样后的频谱Y(w) .

精品 Ts=0.2; % 采样信号的周期.

精品 w0=(2*pi)/Ts; % 延拓周期10*pi

for k=-2:2

ft=f*exp(-j*w0*k*t);

FT=int(ft,t,-1/2,1/2);

ezplot((1/Ts)*abs(FT),[(-4*pi-k*w0):0.01:(4*pi-k*w0)]);

hold on

end

axis([low high -0.1 2.5]); xlabel('');

title('采样后的频谱');

如有侵权请联系告知删除，感谢你们的配合！

合集下载

数字信号处理实验报告

数字信号处理实验报告

《数字信号处理》实验报告

学 院：信息科学与工程学院

专业班级 ：通信1303

姓 名

学 号：

实验一 常见离散时间信号的产生和频谱分析

一、 实验目的

（1） 熟悉MATLAB应用环境，常用窗口的功能和使用方法；

（2） 加深对常用离散时间信号的理解；

（3） 掌握简单的绘图命令；

（4） 掌握序列傅里叶变换的计算机实现方法，利用序列的傅里叶变换对离散信号进行频域分析。

二、 实验原理

（1） 常用离散时间信号

a）单位抽样序列

01)(n00nn

如果)(n在时间轴上延迟了k个单位，得到)(kn即：

01)(kn0nkn

b）单位阶跃序列

01)(nu00nn

c）矩形序列

01)(nRN其他10Nn

d）正弦序列

)sin()(wnAnx

e）实指数序列

f）复指数序列

()()jwnxne

（2）离散傅里叶变换：

设连续正弦信号()xt为

0()sin()xtAt

这一信号的频率为0f，角频率为002f，信号的周期为00012Tf。如果对此连续周期信号()xt进行抽样，其抽样时间间隔为T，抽样后信号以()xn表示，则有0()()sin()tnTxnxtAnT，如果令w为数字频率，满足000012ssfwTff，其中sf是抽样重复频率，简称抽样频率。

为了在数字计算机上观察分析各种序列的频域特性，通常对)(jweX在2,0上进行M点采样来观察分析。 对长度为N的有限长序列x(n), 有

10)()(NnnjwjwkkenxeX

其中 1,,1,02MkkMwk，

通常M应取得大一些，以便观察谱的细节变化。取模|)(|kjweX可绘出幅频特性曲线。

如何利用Matlab进行音频特征提取

如何利用Matlab进行音频特征提取

如何利用Matlab进行音频特征提取

引言：

音频特征提取是音频信号处理领域的重要环节，它可以提取出音频信号的特征参数，进而用于音频分类、识别、检索等应用。在实际应用中，利用Matlab进行音频特征提取是一种常见且有效的方法。本文将介绍如何利用Matlab进行音频特征提取，并逐步详解其步骤和原理。

一、背景知识

在进行音频特征提取之前，我们需要了解一些背景知识。音频信号通常是以时域波形的方式表示的，它是一段连续的时间信号。然而，时域波形无法直接提取到音频信号的有用信息。因此，我们需要将音频信号转换为另一种表示形式，即频域表示。

二、音频信号的频域分析

对音频信号进行频域分析是提取音频特征的第一步。主要包括以下几个步骤：

1. 预处理：首先，将音频信号进行预处理，包括去除噪声、进行均衡化等。这样可以提高信号的质量，减少对后续处理的干扰。

2. 采样和量化：接下来，对预处理后的音频信号进行采样和量化。采样是将连续时间信号转换为离散时间信号，量化是将连续值转换为离散值。

3. 快速傅里叶变换：然后，利用快速傅里叶变换（FFT）将时域信号转换为频域信号。FFT可以将时域波形转换为频谱图，显示出频率和振幅之间的关系。

三、常用的音频特征参数

在频域表示下，我们可以提取出各种音频特征参数。常用的音频特征参数包括： 1. 频谱特征：频谱特征描述了音频信号在频域上的分布情况。常见的频谱特征包括功率谱密度、谱平均值、频谱质心等。

2. 感知特征：感知特征是描述音频信号的听觉感知属性的参数。如音量、音调、音色等。

3. 时域特征：时域特征描述了音频信号在时间域上的变化情况。如时长、均方根能量、过零率等。

四、利用Matlab进行音频特征提取的步骤

现在，我们来详细了解如何利用Matlab进行音频特征提取。以下是具体的步骤：

1. 导入音频文件：首先，使用Matlab的音频文件导入函数，如audioread()或wavread()，将音频文件读取到Matlab的工作空间中。

基于MATLAB的信号采样与重构的实现

基于MATLAB的信号采样与重构的实现

！　＝　ＣＮｌｌ一２０３４／Ｔ　实验技术与管理　Ｅｘｐｅｒｉｍｅｎｔａｌ　Ｔｅｃｈｎｏｌｏｇｙ　ａｎｄ　Ｍａｎａｇｅｍｅｎｔ　第２４卷第２期２００７年２月　Ｖ０１．２４　Ｎｏ．２　Ｆｅｂ．２ｏ０７　

基于ＭＡＴＬＡＢ的信号采样与重构的实现　

周祥才，杨铮　

（常州工学院质量工程学院，江苏常州　２１３００２）　

摘要：该文基于ＭＡＴＬＡＢ再现了连续信号的采样频谱和重构的时域信号，给出了信号重构的误差，并分　析了采样周期对采样频谱叠加和信号重构精度的影响。　关键词：信号采样；ＭＡＴＬＡＢ；采样周期；频谱；信号重构　中图分类号：Ｎ９１１６．６　文献标识码：Ａ　文章编号：１００２－．４９５６（２００７）０２—００７９—０２　

Ｒｅａｌｉｚａｔｉｏｎ　ｏｆ　ｓｉｇｎａｌ　ｓａｍｐｌｉｎｇ　ａｎｄ　ｒｅｃｏｎｓｔｒｕｃｔｉｏｎ　ｂａｓｅｄ　ｏｎ　ＭＡＴＬＡＢ　

ＺＨＯＵ　Ｘｉａｎｇ—ｃａｉ，ＹＡＮＧ　Ｚｈｅｎｇ　

（Ｃｈａｎｇｚｈｏｕ　Ｉｎｓｔｉｔｕｔｅ　ｏｆ　Ｔｅｃｈｎｏｌｏｇｙ，Ｃｈａｎｇｚｈｏｕ　２１３００２，Ｃｈｉｎａ）　

Ａｂｓｔｒａｃｔ：Ｔｈｉｓ　ｔｅｘｔ　ｒｅｃｏｎｓｔｒｕｃｔｅｄ　ｆｒｅｑｕｅｎｃｙ　ｓｐｅｃｔｒｕｍ　ａｎｄ　ｔｈｅ　ｓｉｇｎａｌ　ｏｆ　ａ　ｃｏｎｔｉｎｕｏｕｓ　ｓｉｇｎａｌ　ｂａｓｅｄ　ｏｎ　ＭＡＴＬＡＢ，ｇａｖｅ　ｔｈｅ　ｅｎｌ０ｒ　ｏｆ　ｒｅｃｏｎｓｔｒｕｃｔｉｎｇ　ｓｉｇｎａ１．ａｎｄ　ａｎａｌｙｚｅｄ　ｔｈｅ　ｉｎｆｌｕｅｎｃｅ　ｏｆ　ｓａｍｐｌｅ　ｐｅｒｉｏｄ　ｏｎ　ｓｕｐｅｒｐｏｓｉｔｉｏｎ　ｏｆ　ｔｈｅ　ｆｒｅｑｕｅｎｃｙ　ｓｐｅｃ—　ｔｒｕｍ　ａｎｄ　ｐｒｅｃｉｓｉｏｎ　ｏｆ　ｓａｍｐｌｉｎｇ　ｓｉｇｎａｌ　ｒｅｃｏｎｓｔｒｕｃｔｉｏｎ．　Ｋｅｙ　ｗｏｒｄｓ：ｓｉｇｎａｌ　ｓａｍｐｌｉｎｇ；ＭＡＴＬＡＢ；ｓａｍｐｌｉｎｇ　ｐｅｒｉｏｄ　ｆｒｅｑｕｅｎｃｙ；ｓｐｅｃｔｒｕｍ；ｓｉｇｎａｌ　ｒｅｃｏｎｓｔｒｕｃｔｉｏｎ　

[整理]matlab离散信号的DFT及其快速算法fft以及ifft

[整理]matlab离散信号的DFT及其快速算法fft以及ifft

实验四 离散信号的DFT及其快速算法

一、 实验目的

1．在学习DFT理论的基础上，通过本实验，加深对FFT的理解，体会二者之间的关系。

2．熟悉应用FFT实现两个序列的线性卷积的方法。

二、 实验原理

N点序列x[n] 的DFT和IDFT定义：

可以用函数U=fft（u，N）和u=ifft（U，N）计算N点序列的DFT正、反变换。

三、 实验内容

1. x(n)=R5(n),求N分别取8，32时的X(k)，最后绘出图形。

离散傅立叶变换函数的MATLAB实现如下：

N=8;

x=[ones(1,5),zeros(1,N-5)];

n=0:N-1;

X=dft(x,N);

magX=abs(X); phaX=angle(X)*180/pi;

k=(0:length(magX)'-1)*N/length(magX);

subplot(2,2,1);stem(n,x);title('x(n)—8点');

subplot(2,2,2);stem(k,magX);axis([0,8,0,6]);title('|X(k)|--8点');

N=32; x=[ones(1,5),zeros(1,N-5)];

n=0:N-1;

X=dft(x,N);

magX=abs(X);

phaX=angle(X)*180/pi;

k=(0:length(magX)'-1)*N/length(magX);

subplot(2,2,3);stem(n,x);title('x(n)—32点');

subplot(2,2,4);stem(k,magX);axis([0,32,0,5]);title('|x(k)|--32点');

% dft函数 function[Xk]=dft(xn,N)

n=[0:1:N-1];

k=[0:1:N-1]; WN=exp(-j*2*pi/N);

nk=n'*k;

WNnk=WN.^nk;

Xk=xn*WNnk;