如何从PSD频谱求总功率

jacky19831102

想从声场PSD频谱图中计算总声功率AOSWL。可能还有些基本概念不是很清楚，下面有描述错误的恳请指出。

假设声场的PSD频谱图(1/3倍频图)如下：

我的计算过程是把各频带的声功率叠加得到总声功率
公式是：AOSWL=10*log10(∑10^(SWL(i)/10))
其中SWL(i)是各频带的纵坐标值，如上图，分别为55.2、61.4、58.5、56.1、46.5......30.5
用FORTRAN简单编程：
PARAMETER (N=30)
DIMENSION SWL(N)
DATA SWL/55.2,61.4,58.5,56.1,46.5,51.0,57.6,56.7,57.6,59.3,57.9,57.2,55.8,&
                   53.7,54.4,53.1,52.1,51.5,48.7,46.7,44.0,41.0,37.3,34.8,32.9,31.5,&
                   30.6,30.4,30.4,30.5/
sum=0.0
DO i=1,N
sum=sum+10**(0.1*SWL(i))
END DO
AOSWL=10*ALOG10(sum)
WRITE(*,*)AOSWL
END
程序计算结果为：AOSWL=69.12808(dB)
与上图中显示的70.5(dB)略有不同，所以不知道这个计算方法对不对？

[ 本帖最后由 jacky19831102 于 2009-12-4 12:57 编辑 ]

fendou22

不错啊，先学习一下
股票知识 秦皇岛人才网 奇迹私服  RO私服  天堂2私服

网页游戏  秦皇岛电脑维修  混凝土泵车配件   

Jeremiah

算法是对的，将公式调整如下：
AOSWL=10*log10(∑10^(SWL(i)/10)/1.5)
你再试试

cologate

Jeremiah 发表于 2010-12-14 15:47

                               
登录/注册后可看大图

算法是对的，将公式调整如下：
AOSWL=10*log10(∑10^(SWL(i)/10)/1.5)
你再试试

你好，请问一下这个公式的来源是什么？关于频谱和声功率处理转换这方面有没有什么资料可以推荐参考的？