求ＭＡＴＬＡＢ做傅立叶频谱图的实例

sssssxxxxx921

在ＭＡＴＬＡＢ的帮助中看了半天也没看懂如何画频谱图，哪位能发一个带注释的实例啊

无水1324

在matlab help--〉search--〉FFT
在那里面出来一个实例：
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);% FFT变换

Pyy = Y.* conj(Y) / 512;

f = 1000*(0:256)/512;
plot(f,Pyy(1:257))
title('Frequency content of y')
xlabel('frequency (Hz)')



再看不懂，你就把头伸过来！:lol :lol :lol :lol

wenzihui

2楼给的信号是定周期采集的FFT，敢问，如果是变周期采集的数据如何进行FFT

wy558558558

无水１３２４给的程序幅值的物理意义是什么？

wy558558558

原帖由 无水1324 于 2007-6-21 18:14 发表
在matlab help--〉search--〉FFT
在那里面出来一个实例：
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  ...

本人觉得稍微改下要好点点
y=abs(fft(y,n)*2/n);
要不然幅值太大．

sssssxxxxx921

改动和不改动有什么区别呀，它们分别代表什么意思啊

sssssxxxxx921

在用数值算法（如龙格库塔法）计算得出：时间tc=[tc,t_span(8000:n)];
位移 yc=[yc,q(8000:n,2)];
我用以下ＭＡＴＬＡＢ语句怎么显示不出来啊　　帮忙看看
Y = fft(yc,1024);
Pyy = Y.* conj(Y) / 1024;
f = 1000*(0:1023)/1024;
plot(f(1:512),Pyy(1:512))

liyaohua522

wy558558558的意思可能是改成：abs(fft(y,n)).^2/n; 这样算出来和 Pyy = Y.* conj(Y) /n;其实是一样的，都是周期图法求功率谱密度，没有谁大谁小的问题。
至于sssssxxxxx921在若干个帖子里问的显示不出来的几条语句，从fft那条开始的四句没有问题，如果显示不出来肯定有出错信息，自己多分析一下吧。


个人看法。

sssssxxxxx921

我已经知道了　但是出来的总感觉到不太对　　

请问：这样出来的是不是功率谱　　而不是所谓的幅值谱啊

wy558558558

这个出来的必然是功率谱密度，决不是幅值．看数值也能看出．

sssssxxxxx921

那你说幅值谱该怎么得出呢　　请大虾指点

swk_sun

不懂,等待高手指点,学习! 期待中!!!!!!!!!!!!!!!!

中原

受不了这类似的问题了，我干脆推荐一本书，里面例子多的是
电子工业出版社 吴新余 周井泉 沈元隆<信号与系统——时域、频域分析及matlab软件的应用>的第三章 信号的频谱分析与傅立叶变换，绝对管用:handshake

zhangnan3509

这个在信号版和MATLAB版的例子有很多，可以搜索一下就到手了