Question

我有一段声音信号，想用matlab来分析，求它的最高频率等。请高手帮忙下。

我是自己用录音笔采集的声音信号，大概40几秒，请问如何用matlab来分析，求它最高频率，并画出图形，希望高手帮忙下，一定重谢，高分追加。如果你需要我的声音信号，求留下邮箱，我发给你，谢谢了。

Answer 1

直接FFT分析就可以了，各个频率成分全部出来

Answer 2

网上盛传的《如何从按键音中听出360总裁周鸿祎 的手 机 号 码》你百度后参考下，呵呵

Answer 3

你得知道录音笔的采样频率

追问

谢谢，我知道采样频率后接下来要怎么处理？谢谢了

追答

知道采样频率fs后，由这个录音笔取得的信号的最大频率将会是 fmax = fs/2 (采样定理)。
假如 你的声音信号保存在变数V里。
他的频谱就是 y=abs(fft(V));
假设y的长度是n = length(y);

这个y会是一个左右对称的图形所以只取左半部分 y=y(1:n/2,:);
接下来就是要做出频率轴，也就是频谱的x轴。
x轴 x=linspace(0,fmax,n/2);
意思是在0Hz到fmaxHz之中插入n/2个空当，可以和y一一对应上。
最后plot(x,y)将会看到你的声音信号的频谱。

追问

大哥，你是对的。我按你的方法弄了，得到了你说的图，麻烦再问下，进行傅里叶变化时候取的点数N有什么规定，我有时候取的挺大的，他就比较密。我在取fs的时候，取大点，X轴也跟着变化，比如我FS取12000，信号集中在x轴的6000hz左右，当我去5000的时候，信号集中x轴的2000左右，我要怎么区分，大哥教教我，这个问题困扰了我挺久的。

追答

fft 后面的点数n越多 图应该越清晰 越圆滑，详细查一下　傅立叶变换的 zero padding就应该会知道。但是你声音信号自身的频谱是不会变的，按照我这个方法的话。

我不明白的是，你的fs为什么会变？系统的fs是固定的，不是你随便定的。就算变了，信号的位置是不会变的。

还有一点就是 你说 你fs取12000，信号在6000。。。明显根据采样定理12000的fs是足够看到6000里存在的信号的 因为 fs>2fmax，各种奇怪啊。。。

追问

大哥很感谢你的回答，以后我在我的论文上一定写上你的网名。大哥，我的程序是这样弄的，就是把fs变化了下你帮我看看。

clear

fs=12000;

v=wavread('0072.wav');

y=abs(fft(v));

n = length(y);

y=y(1:n/2,:);

x=linspace(0,fs,n/2);

plot(x,y)

得到的图形是这样的

然后我把上面的fs改成6000，得到下面图形

感觉图形压缩了，这样要怎么看，谢谢了。帮我分析下，感谢。

追答

x=linspace(0,fs,n/2);这句不对，应该是x=linspace(0,fs/2,n/2);你仔细看看我之前说的。。。

还有那个录音笔 是可以变换采样的么？也就是这两个data 一个是在12000的fs下取得的信号，另一个是在6000fs下取得的信号么？
如果你给我看的这两个信号是同一个信号的话，是不可以变换fs的，以多少fs取得的，就得用多少fs来看

假设你的fs是12000，那个用这个信号能观测到的最大频率就是6000，这时如果你想只显示到3000的频率的话
就在plot 下面加一句 xlim([0 3000])这个命令就好了。。
但是变化fs绝对是错的。。

看了看之前你的回答与提问，我就直接了当的说吧，不知道你是什么目的想要观测什么分析什么，如果你这是通过某种器材，机器，在实验中取得的信号的话，针对同一个信号手动变化fs就是错，fs是定数，不是变量。如果想要以不同的fs来观测信号的话，只能是调整机器的fs，重新取得一遍信号，然后进行分析，而不是在用同一个信号进行分析时变化fs。

追问

感谢你这么详细的回答，大哥可以留个邮箱吗，希望以后可以再咨询你下，谢谢你了。

追答

t2mca004@eng.kitakyu-u.ac.jp

Answer 4

频谱分析

Answer 5

把波段立体化

追问

大侠，详细点，谢谢了。能不能帮我处理下，我愿奉献全部的200多财富。

追答

悟

追问

留个邮箱，我把信号发给你吧。谢谢了。

追答

否