声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1996|回复: 0

[FFT] 请教如何将功率谱转换为幅值谱

[复制链接]
发表于 2010-8-6 15:44 | 显示全部楼层 |阅读模式

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

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

x
请教各位采用功率谱模型函数生成的功率谱如何正确的转换为傅里叶幅值谱,我转换后的好像不正确,与实际幅值谱幅值之间相差很大。程序如下:
clc
clear all
close all
direc='Data';
records=fname(direc);
finame=char(records);
[xg,dt,NbPt]=peer2acc(direc,finame);%读取信号,dt为时间步长,xg为信号幅值,n为采样点数
EFz=1/dt;%采样频率
nfft=2^nextpow2(NbPt) ;
nbf=nfft-NbPt;
nzer=zeros(1,nbf);
xg=[xg,nzer];
n=nfft;
Fgen=EFz/nfft:EFz/nfft:EFz/2;%产生离散频率
SaCgen=Respec(Fgen,xg,dt);%计算反应谱
Wgen=2*pi*Fgen;
Z=0.05;%阻尼比
S=15;%强震持时
% 通过功率谱密度函数模型将反应谱转换为功率谱
GSUM=0.0;
for I=1:length(Fgen),
    POW=2*Z*Wgen(I)*S;
    TRANS=1-exp(-POW);
    ZS=Z/TRANS;
    QYS=sqrt(4*ZS/pi);
    XSP=Wgen(I)*S/2/pi/(-log(exp(-1)));
    RSTAR=sqrt(2*log(2*XSP));
    ET=-RSTAR*QYS*sqrt(pi/2);
    ARG=2*XSP*(1-exp(ET));
    RSP=sqrt(2*log(ARG));
    GW(I)=(4*ZS/(Wgen(I)*pi))*((SaCgen(I)/RSP)^2-GSUM);
    if GW(I)<=0.0, GW(I)=0.0;end
    if I==1,
        GSUM=0.5*GW(1)*Wgen(1);
    else GSUM=GSUM+GW(I)*(Wgen(I)-Wgen(I-1));
    end
end
%
% 将功率谱转换为傅里叶幅值谱
AA(1)=sqrt(GW(1)*Wgen(1)*1);
for I=2:length(Fgen)-1
    AA(I)=sqrt(GW(I)*(Wgen(I+1)-Wgen(I-1))*1);
end
I=length(Fgen);
AA(I)=sqrt(GW(I)*(Wgen(I)-Wgen(I-1))*1);

[ 本帖最后由 wdecai 于 2010-8-6 15:57 编辑 ]
回复
分享到:

使用道具 举报

您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

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

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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