声强测量问题

xiaolin7200476

有个朋友在论坛上写了个关于声强计算的例子。
clear
fp=fopen('dai11.txt');
Y=fscanf(fp,'%f',[1,75000]);
fclose(fp);
sig1=Y;   %采集信号

fp=fopen('dai12.txt');
T=fscanf(fp,'%f',[1,75000]);
fclose(fp);
sig2=T;   %采集信号

q=1.29;     %空气密度
d=2.5e-2;   %传声器中心距
for i=1:75000;
    p(i)=(sig1(i)+sig2(i))/2;
    p1=sig2(i)-sig1(i);
    u(i)=[1/(q*d)]*p1*3;   %质点速度
    P(i)=u(i)*p(i)';   %求声强
end
这个程序是关于时域的计算，如果我想把他变成频域的，我该怎么办？把原始信号先FFT行吗？

mao

对原始信号进行FFT后再按上面的程序计算肯定错了，因为频域的乘积对应时域的卷积了，这样结果肯定不对的。你如果想得到频域的结果，可以对上面程序得到的结果进行频域的分析啊