利用matlab对音频做FFT进行频谱分析
利用matlab对音频的左右声道做频谱分析。
syms y fs;
syms in_put;
in_put= 'C:\Users\getao\Desktop\立体声.wav'
%in_put 为音频文件所在的位置,.wav文件为文件的名称。
[y,fs] = audioread(in_put);
info=audioinfo(in_put) ;
sound(y,fs);
T=1/fs;
t=(0:length(y)-1)*T;
f=(0:length(y)-1)*fs/length(y);
figure(1);
yz=y(:,1);
%提取音频信号的左声道数据
subplot(2,2,1);
plot(t,yz);
title('原始信号时域(左声道)');
xlabel('时间');
ylabel('振幅');
yz2=y(:,2);
%提取音频信号的右声道数据
subplot(2,2,2);
plot(t,yz2);
title('原始信号时域(右声道)');
xlabel('时间');
ylabel('振幅');
subplot(2,2,3);
n=length(yz);
y1=fft(yz,n);
F=fs/length(yz);
plot(f,abs(y1));
title('原始信号频谱(左声道)');
xlabel('F(Hz)');
ylabel('H(jw)');
subplot(2,2,4);
n2=length(yz2);
y1=fft(yz2,n2);
F=fs/length(yz2);
plot(f,abs(y1));
title('原始信号频谱(右声道)');
xlabel('F(Hz)');
ylabel('H(jw)');
syms y fs;syms in_put;
in_put= 'C:\Users\getao\Desktop\右声道.wav';[y,fs] = audioread(in_put);info=audioinfo(in_put) ;
sound(y,fs);
T=1/fs;
t=(0:length(y)-1)*T;
f=(0:length(y)-1)*fs/length(y);
figure(1);
yz=y(:,1);
subplot(2,2,1);
plot(t,yz);
title('原始信号时域(左声道)');
xlabel('时间');
ylabel('振幅');yz2=y(:,2);
subplot(2,2,2);
plot(t,yz2);
title('原始信号时域(右声道)');
xlabel('时间');
ylabel('振幅');subplot(2,2,3);
n=length(yz);
y1=fft(yz,n);
F=fs/length(yz);plot(f,abs(y1));
title('原始信号频谱(左声道)');
xlabel('F(Hz)');
ylabel('H(jw)');subplot(2,2,4);
n2=length(yz2);
y1=fft(yz2,n2);
F=fs/length(yz2);plot(f,abs(y1));
title('原始信号频谱(右声道)');
xlabel('F(Hz)');
ylabel('H(jw)');
下图分别为三段音频的分析,三段音频文件对同一首音乐的某一段进行录音得到。立体声左右声道都录制了声音,左右声道的文件则通过录音设备的两端进行分别录制。
立体声 音频分析
左声道 音频分析
右声道 音频分析
从上图可以看出,该段声音的频率在0-10KHz范围内,在未录制声音的音道有杂音的干扰使其频谱有尖峰的存在。
利用matlab对音频做FFT进行频谱分析相关推荐
- matlab录音函数怎么保存在哪里,利用matlab 录音
利用matlab的音频信号处理工具箱,可以实现声音的录制和播放. 录音函数wavrecord语法为: y=wavrecord(n,fs,channel,dataType); 其中 n为采样点数,fs为 ...
- 利用matlab做dsp实验,华工DSP实验一利用Matlab进行频谱分析
实验一利用Matlab进行频谱分析 一.对无限长时域离散信号进行频谱分析 1.代码: function test1(N) N = input('Type in the length = '); k = ...
- Matlab如何进行利用离散傅里叶变换DFT (快速傅里叶变换FFT)进行频谱分析
文章目录 1. 定义 2. 变换和处理 3. 函数 4. 实例演示 例1:单频正弦信号(整数周期采样) 例2:单频正弦信号(非整数周期采样) 例3:含有直流分量的单频正弦信号 例4:正弦复合信号 例5 ...
- matlab相关性分析频谱_利用Matlab绘制正弦信号的频谱图并做相关分析[共6页]
<利用Matlab绘制正弦信号的频谱图并做相关分析[共6页]>由会员分享,可在线阅读,更多相关<利用Matlab绘制正弦信号的频谱图并做相关分析[共6页](6页珍藏版)>请在人 ...
- 快速傅里叶变换FFT进行频谱分析(matlab)
快速傅里叶变换FFT进行频谱分析(matlab) 本章摘要:FFT是离散傅立叶变换的快速算法,可以将一个信号变换到频域.有些信号在时域上是很难看出什么特征的,但是如果变换到频域之后,就很容易看出特征了 ...
- 利用MATLAB实现WAV音频的LSB信息隐藏
0. 背景 音频信息隐藏技术就是将秘密数据直接嵌入音频载体中,但不影响原载体的使用价值,也不容易被人的知觉系统觉察或注意到. 音频信息隐藏技术作为信息隐藏技术的一个重要分支,受到越来越多的关注. 近年 ...
- 利用matlab怎样进行频谱分析
图像的频率是表征图像中灰度变化剧烈程度的指标,是灰度在平面空间上的梯度.如:大面积的沙漠在图像中是一片灰度变化缓慢的区域,对应的频率值很低:而对于地表属性变换剧烈的边缘区域在图像中是一片灰度变化剧烈的 ...
- matlab 频谱图例子_利用matlab怎样进行频谱分析
图像的频率是表征图像中灰度变化剧烈程度的指标,是灰度在平面空间上的梯度.如:大面积的沙漠在图像中是一片灰度变化缓慢的区域,对应的频率值很低:而对于地表属性变换剧烈的边缘区域在图像中是一片灰度变化剧烈的 ...
- matlab adc仿真,[转载]关于ADC仿真做FFT的设置和结果分析
这些都是比较基础和实用的东西,故做分享~ 以前在实验室跑ADC动态性能仿真,会跑不少点,浪费了不少时间.主要是由于当时没有用相干采样的频率,而是直接选一个和采样频率比较互质的输入频率,然后再进行加窗处 ...
最新文章
- kafka 异步发送阻塞_Kafka学习一
- Boosting、Adaboost、AdaBoost模型的优缺点、提升树、梯度提升树GBDT
- 评审关上了你CVPR的门?这还有Rebuttal的窗,7个小技巧送上
- PCL点云库用Poisson网格化实现点云的表面重建
- Java+MyEclipse+Tomcat (三)配置MySQL及查询数据显示在JSP网页中
- es中的高效文件读取方式
- 360推出国内首个工业互联网安全态势感知系统
- 干货| 掌握这种逻辑思维,大厂面试成功率可提升90%
- springmvc前后端传值总结
- UNIX环境高级编程习题——第七章
- 修改PDF文件软件中文问题,福昕软件Foxit PDF Editor
- 系统架构图 云架构案例
- P,NP,NPC,NP-HARD 图片基于P!=NP
- 恶意代码分析实战——反汇编
- 万年历程序例题(农历阴历转换)
- mac访达中显示隐藏的文件夹和文件
- tensorflow和cuda以及cudnn版本对齐
- 阿里实人认证java_GitHub - lap888/react-native-face-aliverify: 基于阿里金融级实人认证sdk封装分享...
- ICPC-TOOLs resolver滚榜单工具使用指南
- VQA系列论文(二)
热门文章
- Latex编译报错 “Undefined control sequence. \chapter“
- 咕咕机显示服务器请求异常,咕咕机 云服务器异常
- 【研究】周耀旗写好英语科技论文的诀窍
- 应用软件中如何实现多币种自动换算
- 电商项目:拼团活动设计及表结构设计
- oracle导出dnf文件,ORACLE 常见故障恢复-dnf文件损坏
- BugkuCTF-web-秋名山车神 writeup
- Nginx命令stop或者quit停止服务器无效失效
- 关于博客园官方博客技术改造升级的启发
- 关于计算机的网络作文,关于网络的作文:网络_450字