声强测量问题
有个朋友在论坛上写了个关于声强计算的例子。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行吗? 对原始信号进行FFT后再按上面的程序计算肯定错了,因为频域的乘积对应时域的卷积了,这样结果肯定不对的。你如果想得到频域的结果,可以对上面程序得到的结果进行频域的分析啊
页:
[1]