请教问题：已知振动源振动加速度，求一定距离处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,);
算出来的值与用声级计测出来的值相差甚远，希望大家能帮我看看程序，问题到底处在什么地方？谢谢！ “从图中可以看出振动加速度信号和声压信号都分布在100Hz的整数倍处”这是必然的，呵呵。
根据表面振速或者加速度预报辐射声场，对于你提到的模型只能采用数值方法。如果自己不打算涉足算法的话，那就去尝试sysnoise等商业边界元软件吧。
回复 沙发 w89986581 的帖子
数值方法不是对振动源的形状有一定要求吗？现在所测的这个振动源他不是一个简单的规则的形状，怎么来求啊？不能用实验来确定一些未知的参数吗？ 数值方法是针对任意形状的，解析法是针对特殊形状的，呵呵。回复 地板 w89986581 的帖子
有没有推荐的资料啊？数值法是什么方法啊？另外，用试验采取类似于黑盒的方法，没有可以推荐的方法入手点吗？
因为刚看声学振动方面的东西，知道的东西太少了
多谢指点！ 呵呵，精确解只能使用数值方法，例如边界元法。这方面商业软件已经比较成熟。
对精度要求不高的话，可以尝试试验方法，获取大量试验数据的基础上从中拟合，经验很重要，呵呵。
回复 6楼 w89986581 的帖子
如果采用试验的方法，应该从什么地方入手啊？从能量的角度还是直接从声压和振动的幅值上来进行比较啊？另外，如果用数值方法要用哪些软件啊，新手需要多久才能上手啊 试验方法请参考
周春凯(译), 张均平(校). 潜艇自噪声监测系统. 国外舰船工程. 2002, 1: 77-81.
数值方法请参考
LMS International, SYSNOISE Users Manual, 1998
回复 8楼 w89986581 的帖子
好的，我先找找书 原帖由 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 {:{19}:}
页:
[1]