声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2418|回复: 5

[FFT] 二维信号如何求功率谱

[复制链接]
发表于 2007-12-28 13:47 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
:@) 大家好,我现在有一个二维的离散信号,我要求它的径向对数功率谱曲线,请问我该怎么求啊,是直接幅度平方再取对数么?
回复
分享到:

使用道具 举报

发表于 2007-12-28 16:00 | 显示全部楼层
理论上进行自相关得到的就是功率谱了
发表于 2008-3-5 11:08 | 显示全部楼层

借贵地一用

我现在同样在做这方面的工作,也很想知道怎么去求二维径向功率谱怎么去求?如果你知道了方法希望告诉哈!谢谢
 楼主| 发表于 2008-3-18 09:59 | 显示全部楼层
不知道你做的是什么样的数据处理,我做得是重磁异常的数据处理,径向平均对数功率谱其实是一个统计平均值,假如你有256*256个数据,那你要把频率分成128份,然后统计每一圈内的数据点数,算功率谱的平均值再取对数,就行了
发表于 2009-2-17 16:31 | 显示全部楼层

re::::::::::::::::

自相关的fft2就是功率谱函数!!(矩阵)
%% 由功率谱求自相关
Z=importfile('112.mat');   %%原始数据
Z=Z-mean(Z(:));
Pzz=fft2(Z,256,256);
Pzz=Pzz.*conj(Pzz)/(256*256);
Pzz1=fftshift(Pzz);
R=real(ifft2(Pzz));
R=R/(max(R(:)));
R=fftshift(R);

%% 由自相关求功率谱
Z=importfile('112.mat');   %%原始数据
Z=Z-mean(Z(:));
Rr=xcorr2(Z);
Rr=Rr/max(Rr(:));
Pzzr=fft2(Rr);
[m,n]=size(Pzzr);
Pzzr=Pzzr.*conj(Pzzr)/(m*n);
Pzzr=fftshift(Pzzr);
发表于 2013-3-16 12:00 | 显示全部楼层

求分享源代码
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-25 11:46 , Processed in 0.071539 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表