weipeng 发表于 2007-10-14 16:28

关于信号频率问题

(1)
N=2048;
=wavread('d:\11.wav'); %读取语音信号的数据,赋给变量x1
y1=fft(x1,N);         %对信号做2048点FFT变换
y1 = fftshift(y1);
f=fs*(0:1024)/2048;
pyy=y1.*conj(y1)/2048;
plot(f,pyy(1:1025));
(2)t = 0:0.001:0.6;
x = sin(2*pi*50*t)+sin(2*pi*120*t);
y = x + 2*randn(size(t));
plot(1000*t(1:50),y(1:50))

title('Signal Corrupted with Zero-Mean Random Noise')
xlabel('time (milliseconds)')
Y = fft(y,512);
Pyy = Y.* conj(Y) / 512;
f = 1000*(0:256)/512;
plot(f,Pyy(1:257))

程序(1)(2)中的红色部分,一个用fs乘,一个用1000乘。请问f=fs*(0:1024)/2048;中的fs位置的数值是怎么取值的?画出的图横坐标是频率还是点数?我想知道一个声音信号的频率范围应该怎么做?画功率谱图还是画其它的?谢谢!

[ 本帖最后由 eight 于 2007-10-14 17:07 编辑 ]

eight 发表于 2007-10-14 17:07

原帖由 weipeng 于 2007-10-14 16:28 发表 http://www.chinavib.com/forum/images/common/back.gif
(1)
N=2048;
=wavread('d:\11.wav'); %读取语音信号的数据,赋给变量x1
y1=fft(x1,N);         %对信号做2048点FFT变换
y1 = fftshift(y1);
f=fs*(0:1024)/2048;
pyy=y1.*conj(y1)/2048; ...

请移步信号处理方法板块,然后自己搜索一下,估计就能找到答案了
页: [1]
查看完整版本: 关于信号频率问题