wangwenting 发表于 2008-9-7 20:34

请教问题:已知振动源振动加速度,求一定距离处A计权声压级

外行请教,看了一些书和论坛上的东西,但是实际解决起问题来还是有很多没有办法下手的地方,希望大家多多指教
       需要解决的问题:根据测得的振动源的加速度,来求在一定距离处的A计权声压级
       所测振动源有不封闭的外壳,所测点位于外壳之外的一定距离处。
       为了找到振动源加速度和一定距离处A计权声压级之间的关系,给振动源固定频率100Hz的驱动力(由于驱动中含有线圈,所以驱动力实际为100Hz的整数倍),下附图为振动加速度信号和声压信号的fft图,从图中可以看出振动加速度信号和声压信号都分布在100Hz的整数倍处,但是要如何找这两者之间的关系,我还是无从下手,希望高手指点。
       另外,我想从声压信号得到A集权的声压级,下面为所编MATLAB程序
data=csvread('F0068CH2.CSV',0,3);
noi=data(:,2);
=oct3bank(noi);
plot(ff,p);
Cj=[-50.5 -44.7 -39.4 -34.6 -30.2 -26.2 -22.5 -19.1 -16.1 -13.4 -10.9 -8.6 -6.6 -4.8 -3.2 -1.9 -0.8 0 0.6 1.0 1.2 1.3 1.2 1.0 0.5 -0.1 -1.1 -2.5 -4.3 -6.6]
La=0;
Lj=zeros(1,30);
for n=1:30
Lj(n)=10^(0.1*(10*log10(p(n)/(10*exp(-12)))+Cj(n)));
end
La=10*log10(sum(Lj));
function = oct3bank(noi);
pi = 3.14159265358979;
Fs = 500000;    % Sampling Frequency
N = 3;      % Order of analysis filters.
ff=
p = zeros(1,30);
% Design filters and compute RMS powers in 1/3-oct. bands
for i = 1:30
    = oct3dsgn(ff(i),Fs,N);
   y = filter(B,A,noi);            %<-----滤波器
   p(i) = sqrt(var(y));      %得到rms
end

function = oct3dsgn(Fc,Fs,N);
pi = 3.14159265358979;
f1 = Fc/(2^(1/6));
f2 = Fc*(2^(1/6));
Qr = Fc/(f2-f1);
Qd = (pi/2/N)/(sin(pi/2/N))*Qr;
alpha = (1 + sqrt(1+4*Qd^2))/2/Qd;
W1 = Fc/(Fs/2)/alpha;
W2 = Fc/(Fs/2)*alpha;
= butter(N,);
算出来的值与用声级计测出来的值相差甚远,希望大家能帮我看看程序,问题到底处在什么地方?谢谢!

w89986581 发表于 2008-9-8 12:08

“从图中可以看出振动加速度信号和声压信号都分布在100Hz的整数倍处”这是必然的,呵呵。
根据表面振速或者加速度预报辐射声场,对于你提到的模型只能采用数值方法。如果自己不打算涉足算法的话,那就去尝试sysnoise等商业边界元软件吧。

wangwenting 发表于 2008-9-8 15:16

回复 沙发 w89986581 的帖子

数值方法不是对振动源的形状有一定要求吗?现在所测的这个振动源他不是一个简单的规则的形状,怎么来求啊?不能用实验来确定一些未知的参数吗?

w89986581 发表于 2008-9-8 18:35

数值方法是针对任意形状的,解析法是针对特殊形状的,呵呵。

wangwenting 发表于 2008-9-8 19:37

回复 地板 w89986581 的帖子

有没有推荐的资料啊?数值法是什么方法啊?
另外,用试验采取类似于黑盒的方法,没有可以推荐的方法入手点吗?
因为刚看声学振动方面的东西,知道的东西太少了
多谢指点!

w89986581 发表于 2008-9-8 19:44

呵呵,精确解只能使用数值方法,例如边界元法。这方面商业软件已经比较成熟。
对精度要求不高的话,可以尝试试验方法,获取大量试验数据的基础上从中拟合,经验很重要,呵呵。

wangwenting 发表于 2008-9-9 09:07

回复 6楼 w89986581 的帖子

如果采用试验的方法,应该从什么地方入手啊?从能量的角度还是直接从声压和振动的幅值上来进行比较啊?
另外,如果用数值方法要用哪些软件啊,新手需要多久才能上手啊

w89986581 发表于 2008-9-9 17:52

试验方法请参考
周春凯(译), 张均平(校). 潜艇自噪声监测系统. 国外舰船工程. 2002, 1: 77-81.

数值方法请参考
LMS International, SYSNOISE Users Manual, 1998

wangwenting 发表于 2008-9-9 20:01

回复 8楼 w89986581 的帖子

好的,我先找找书

wangwenting 发表于 2008-9-9 20:07

原帖由 w89986581 于 2008-9-9 17:52 发表 http://www.chinavib.com/forum/images/common/back.gif
试验方法请参考
周春凯(译), 张均平(校). 潜艇自噪声监测系统. 国外舰船工程. 2002, 1: 77-81.

数值方法请参考
LMS International, SYSNOISE Users Manual, 1998

我在我们学校的图书馆,期刊网上都没有找到,你有没有这两篇文章啊,能不能发给我啊,谢谢
wangwenting_5@hotmail.com

tj591992034 发表于 2012-5-7 20:39

{:{19}:}
页: [1]
查看完整版本: 请教问题:已知振动源振动加速度,求一定距离处A计权声压级