应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » 汇编语言 » matlab 傅里叶变换 频谱分析
北京石油化工学院2026年研究生招生接收调剂公告
131/2返回列表12下一页
查看: 6574  |  回复: 12
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
本帖产生 1 个 程序强帖 ，点击这里进行查看

whfire

金虫 (正式写手)

[求助] matlab 傅里叶变换 频谱分析

我想利用matlab进行编程，进行傅里叶频谱分析，但是需要将频谱转换成功

率谱。详情见附件
回复此楼

» 收录本帖的淘帖专辑推荐

matlab 频谱

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-07-25 19:45:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

sudo

木虫 (正式写手)
★
dubo(金币+1): 欢迎常来程序语言版讨论 2011-07-31 14:25:14
唉，可能解释起来一言难尽~建议楼主阅读一下信号处理的入门书嗯

首先，你师兄的图，应该是降采样过后的（或者是后来截取了一段，不过我倾向于认为是降采样），也就是说他用来做功率谱的数据采样率只有400Hz，你给的数据的采样率是2000Hz

其次，你师兄给的图是没有取dB作为纵坐标单位的，也就是没有计算10*log10(P)，而是直接用“功率密度”来做纵坐标~

最后，其实MATLAB有现成函数可以做功率谱：
CODE:
[P,F]=periodogram(x, window, nfft, Fs)

其中输出P和F分别是功率谱和对应的频率
输入x是信号，window是数据窗口（可以写[]用默认矩形窗），nfft是进行fft变换的点数（可以写[]取默认值，一般会取到最近的2的整数次幂左右的地方），Fs是采样率，单位是Hz

然后我对你的数据裸进行功率谱处理的话，结果如下图：
赞一下(8人) 回复此楼
2楼2011-07-25 20:28:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sudo

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
xzhdty(金币+2): 欢迎常来 2011-07-25 20:57:22
dubo(金币+5, 程序强帖+1): 感谢应助 2011-07-31 14:26:41
另外说明：MATLAB的功率谱实现函数是像这样子的：
CODE:
function [P,f]=fftpsd(x,Fs)
    m = length(x);          % 整个数据长度
    n = pow2(nextpow2(m));  % 扩展到合适长度，2的整数次幂
    y = fft(x,n);           
    f = (0:n-1)*(Fs/n);     %频率范围
    P = y.*conj(y)/n;   %功率谱密度

    f=f(1:n/2);
    P=P(1:n/2);
   
    %plot(f,10*log10(P)); %把数据转换为dB
    plot(f,P);
    grid;
    xlabel('Frequency (Hz)');
    ylabel('Power(dB)');
    title('{\bf Periodogram}');

为了使用FFT算法（速度快），MATLAB把数据扩展了，本来降采样你的数据之后，一共有4000个点，但是实际上使用了4096个点。但是由于数据是实数序列，所以，仅有采样率一半的数据是有效的（后半部分其实是镜像），所以在图上，F和P都是2048个点（periodogram取了2049个，这个无所谓）
赞一下(7人) 回复此楼
4楼2011-07-25 20:53:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sudo

木虫 (正式写手)

【答案】应助回帖

whfire(金币+20): 哥们，谢谢。太谢谢了 2011-07-27 09:13:49
对了，还有一些想补充的~

首先就是直流分量的问题，如果不先对数据都减去一个均值，那么接近0Hz部分的功率会非常非常大~这是有害的，一般做数据处理的时候如果只关注交流量的话，就要先使得整个序列的均值为0，方法就是给每个数减去该序列的一个均值

然后就是降采样的方法，其实很简单啦，在MATLAB里面的话，如果y是你原来的数据的话（注意是测量数据，你给的数据的第2列）.....还是给出命令的全部吧：
CODE:
ty = y(1:5:end); %将原来2000Hz的数据降采样到400Hz
ty = ty - mean(ty); %减去直流分量
[P,F]=periodogram(ty,[],[],400); %求功率谱
plot(F, P) %画图

结果就是
赞一下(6人) 回复此楼

» 本帖已获得的红花（最新10朵）

独孤_逍遥
3楼2011-07-25 20:44:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mogul

金虫 (正式写手)
引用回帖:
696000楼: Originally posted by sudo at 2011-07-25 20:28:50
唉，可能解释起来一言难尽~建议楼主阅读一下信号处理的入门书嗯

首先，你师兄的图，应该是降采样过后的（或者是后来截取了一段，不过我倾向于认为是降采样），也就是说他用来做功率谱的数据采样率只有400Hz，你 ...

[P,F]=periodogram(x, window, nfft, Fs)
[Pxx,w] = periodogram(x,window,nfft) uses the modified periodogram to estimate the PSD while specifying the length of the FFT with the integer nfft.
"If you set nfft to the empty vector [], it takes the default value
for N listed in the previous syntax. "(ref, matlab help)
"If you set nfft to the empty vector [], nfft defaults to min(256,length(x)" (ref, 薛年喜. MATLAB在数字信号处理中的应用(第2版). p286)
赞一下(1人) 回复此楼
6楼2013-03-13 22:04:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

sudo

木虫 (正式写手)
★
dubo(金币+1): 欢迎常来程序语言版讨论 2011-07-31 14:26:48
还有一个需要补充~

可能你会疑惑，第一列数据（时间）为什么计算中没用到？

实际情况是，我只需要知道，你的第二列序列是等间隔采样出来的，以及他们之间的间隔（采样周期），这就足够了。虽然我也是从第一列算出来的，不过其实第一列并不需要全部给出来
赞一下(5人) 回复此楼
5楼2011-07-25 20:56:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

夏天summer

铜虫 (小有名气)
看起来很厉害。。
回复此楼
时间如拉车的牛。
7楼2013-07-10 10:33:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤_逍遥

铁虫 (初入文坛)
送红花一朵
引用回帖:
3楼: Originally posted by sudo at 2011-07-25 20:44:28
对了，还有一些想补充的~

首先就是直流分量的问题，如果不先对数据都减去一个均值，那么接近0Hz部分的功率会非常非常大~这是有害的，一般做数据处理的时候如果只关注交流量的话，就要先使得整个序列的均值为0，方 ...

sudo是大神啊，直流分量和交流量的问题真是让我茅塞顿开啊！！
我不是学信号处理的，所以国外文献给出的结果跟我算出来的死活不对，原来他们都做的是交流量的处理，都减了一个平均值。
赞一下 回复此楼
8楼2014-02-28 16:18:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mallory19

新虫 (初入文坛)
引用回帖:
3楼: Originally posted by sudo at 2011-07-25 20:44:28
对了，还有一些想补充的~

首先就是直流分量的问题，如果不先对数据都减去一个均值，那么接近0Hz部分的功率会非常非常大~这是有害的，一般做数据处理的时候如果只关注交流量的话，就要先使得整个序列的均值为0，方 ...

sudo，你好。
我在用Mathematica中的Periodogram和PeriodogramArray对同一个数据列表计算时，得到的结果却不同。我在网上查了一下，也有别人在问这个问题。但没有回答。
不知道你对Mathematica是否了解，是否可以交流一下。
赞一下 回复此楼
9楼2014-04-03 23:53:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

whfire

金虫 (正式写手)
引用回帖:
9楼: Originally posted by mallory19 at 2014-04-03 23:53:13
sudo，你好。
我在用Mathematica中的Periodogram和PeriodogramArray对同一个数据列表计算时，得到的结果却不同。我在网上查了一下，也有别人在问这个问题。但没有回答。
不知道你对Mathematica是否了解，是否可以 ...

不好意思，我也不懂
赞一下 回复此楼
10楼2014-04-04 08:12:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 whfire 的主题更新
131/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 302分求调剂 一志愿安徽大学085601 +8 zyx上岸！ 2026-04-04 8/400 2026-04-05 23:08 by lbsjt
[考研] 320分人工智能调剂 +8 振—TZ 2026-04-03 8/400 2026-04-05 22:33 by 范式思维
[考研] 材料专硕283求调剂 +14 试试看呗 2026-04-04 15/750 2026-04-05 19:21 by jkddd
[考研] 288求调剂 +7 没有答案_ 2026-04-05 7/350 2026-04-05 17:01 by yulian1987
[考研] 一志愿北交大材料工程总分358求调剂 +6 cs0106 2026-04-05 6/300 2026-04-05 16:34 by imissbao
[考研] 346分的生物与医药08600求调剂 +5 常雨阳上岸 2026-04-05 6/300 2026-04-05 13:42 by imissbao
[考研] 271分求调剂学校 +12 zph158488！ 2026-04-02 13/650 2026-04-05 10:13 by lqwchd
[考研] 求调剂 +3 小沢 2026-04-03 3/150 2026-04-05 09:10 by sihailian3
[考研] 求生物学调剂 +14 15172915737 2026-04-01 14/700 2026-04-04 20:13 by babysonlkd
[考研] 本科211，专业085404，293分请求调剂 +5 莲菜就是藕吧 2026-04-04 5/250 2026-04-04 14:08 by 这是一个无聊的�
[考研] 求调剂 +3 wos666 2026-04-03 3/150 2026-04-04 05:16 by gswylq
[考研] 化工调剂303分，过四级 +28 栖梧待风 2026-04-02 28/1400 2026-04-03 21:40 by qlm5820
[考研] 本科985，专业0812分336求调剂 +4 莫莫很行 2026-04-03 4/200 2026-04-03 21:31 by zhq0425
[考研] 303求调剂 +9 DLkz1314. 2026-03-30 9/450 2026-04-03 18:34 by ls刘帅
[考研] 081200-11408-276学硕求调剂 +5 崔wj 2026-04-03 5/250 2026-04-03 15:06 by arrow8852
[考研] 一志愿山东大学化学与化工学院材料与化工专硕，360分求调剂 +4 不愿透露姓名的� 2026-04-02 4/200 2026-04-03 09:29 by 遗忘消失的灆
[考研] 学硕化学工程与技术，一志愿中国海洋大学320+求调剂 +8 披星河 2026-04-02 8/400 2026-04-02 14:12 by oooqiao
[考研] 372求调剂 +3 jj涌77 2026-04-02 3/150 2026-04-02 09:57 by olim
[考研] 材料专业调剂 +5 啦啦啦哭 2026-03-31 6/300 2026-04-01 16:48 by JourneyLucky
[考研] 0710生物学求调剂！ +6 叙述文 2026-03-31 6/300 2026-04-01 09:39 by JourneyLucky
信息提示
关闭
请填处理意见
关闭