xiaolin7200476 发表于 2009-4-5 18:53

声强测量问题

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

fp=fopen('dai12.txt');
T=fscanf(fp,'%f',);
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)=*p1*3;   %质点速度
    P(i)=u(i)*p(i)';   %求声强
end
这个程序是关于时域的计算,如果我想把他变成频域的,我该怎么办?把原始信号先FFT行吗?

mao 发表于 2009-5-6 15:21

对原始信号进行FFT后再按上面的程序计算肯定错了,因为频域的乘积对应时域的卷积了,这样结果肯定不对的。你如果想得到频域的结果,可以对上面程序得到的结果进行频域的分析啊
页: [1]
查看完整版本: 声强测量问题