【数字信号处理】傅里叶变换性质 ( 傅里叶变换频移性质示例 | PCM 音频信号处理 | 使用 matlab 进行频移操作 )

举报
韩曙亮 发表于 2022/03/13 22:36:23 2022/03/13
【摘要】 文章目录 一、PCM 音频信号处理二、要点说明 一、PCM 音频信号处理 给定一段 PCM 音频数据 , 封装到 .WAV 文件 中 , 此时的信号 是 实信号 ...





一、PCM 音频信号处理



给定一段 PCM 音频数据 , 封装到 .WAV 文件 中 ,

此时的信号 是 实信号 , 每个 PCM 音频采样都是一个 实数值 ,

绘制该 音频数据 实信号 的 " 幅频特性 " ,


将该 PCM 音频的 频率 移动 6 k H z 6kHz 6kHz ,

频移后的 音频信号 变成了 复信号 ,

绘制该 复信号的 " 幅频特性 " ;



准备一个 wav 格式的音频文件 , 该 music.wav 文件是一个

16 16 16

单声道

44100 H z 44100Hz 44100Hz 采样率

的音频文件 , 其 采样波形 与 频谱 如下 :

在这里插入图片描述

代码示例 :

% wave 文件采样
pcmFileName='music.wav';

% 读取  music.wav 文件
% s1 是读取的信号向量
% sampleRate 是采样率 44100Hz 表示每秒有 44100 个音频样本
[signal,sampleRate] = audioread(pcmFileName);

% 上面的 sampleRate 计算出来是 44100 单位是 Hz
% 这里需要将 sampleRate 的单位变为 kHz , 除以 1000 
sampleRate = sampleRate/1000;

% 矩阵转置 
%1x132300 矩阵转为 132300x1 矩阵
signalT=signal';

% 取 信号 的前 2048 个点进行分析
x=signalT(1:2048);

% 求 signalT 实信号的 频谱
xf=20 * log10(abs(fft(x,4096)));
xf1(1:2048)=xf(2049:4096);
xf1(2049:4096)=xf(1:2048);

% 绘制实信号的频谱图
figure;
plot((-2048:2047) * sampleRate/4096,xf1);
grid on;

% 根据 傅里叶变换 频移性质, 进行频谱搬移
% 频率搬移 6kHz
y=x.*exp(2*j*pi*6/32*(0:2047));

% 求 signalT 频移后的 复信号的 频谱
yf=20*log10(abs(fft(y,4096)));
yf1(1:2048)=yf(2049:4096);
yf1(2049:4096)=yf(1:2048);

% 绘制频率搬移 6kHz 后的复信号的频谱图
figure;
plot((-2048:2047)*fs/4096,yf1);
grid on;

  
 
  • 1
  • 2
  • 3
  • 4
  • 5
  • 6
  • 7
  • 8
  • 9
  • 10
  • 11
  • 12
  • 13
  • 14
  • 15
  • 16
  • 17
  • 18
  • 19
  • 20
  • 21
  • 22
  • 23
  • 24
  • 25
  • 26
  • 27
  • 28
  • 29
  • 30
  • 31
  • 32
  • 33
  • 34
  • 35
  • 36
  • 37
  • 38
  • 39
  • 40
  • 41
  • 42

在这里插入图片描述





二、要点说明



根据 图1 所示 , 实信号的 " 幅频特性 " 是 偶对称的 ;
在这里插入图片描述
图2 复信号 的 " 幅频特性 " , 向右平移了 6 kHz , 此时就不是 偶对称 了 ;

在这里插入图片描述

傅里叶变换时移 只影响 " 相频特性 " , 不影响 " 幅频特性 " ;

傅里叶变换频移 同时影响 " 相频特性 " 和 " 幅频特性 " ;

文章来源: hanshuliang.blog.csdn.net,作者:韩曙亮,版权归原作者所有,如需转载,请联系作者。

原文链接:hanshuliang.blog.csdn.net/article/details/123396457

【版权声明】本文为华为云社区用户转载文章,如果您发现本社区中有涉嫌抄袭的内容,欢迎发送邮件进行举报,并提供相关证据,一经查实,本社区将立刻删除涉嫌侵权内容,举报邮箱: cloudbbs@huaweicloud.com
  • 点赞
  • 收藏
  • 关注作者

评论(0

0/1000
抱歉,系统识别当前为高风险访问,暂不支持该操作

全部回复

上滑加载中

设置昵称

在此一键设置昵称,即可参与社区互动!

*长度不超过10个汉字或20个英文字符,设置后3个月内不可修改。

*长度不超过10个汉字或20个英文字符,设置后3个月内不可修改。