matlab做点声源辐射模型

难渡秋心

建立Matlab模型，下图中移动声源为某一具体车辆，汽车按行驶中心线行驶的距离为20m，传感器距离车辆行驶中心线7.5m，传感器距离地面高度为1.2m。

      模型建立思想相当简单，为简单计算考虑，只取一个点:发动机附近，取其等效几何中心位置视为该等效点源位置测量车辆处于静止位置，即驻车工况车辆前端处于起始线位置处，发动机附近声压频谱图，测量其近场测量距离，以及车辆匀加速运行过程中，任一时刻时该点声源与传声器之间的距离。有关参数会在下面给出

按照半自由场空间点声源辐射模型理论，设某声源声功率级为Lw，则在距离声源某位置处的声压级满足Lw=Lp1+20logr1+8=Lp2+20logr2+8，从而求出车辆运行状态下，发动机辐射噪声随时间变化即随车辆运行位置变化时辐射到远场传声器位置的声压级。

图形如下图所示，移动声源为车辆，下面给出相应参数：

S=20;                               %总加速距离/m

L=13.5;                             %车辆长度

a=1;                                %加速度大小/m.s-2

p0=2e-5;                            %基准声压/Pa

h=1.2;                              %接收器高度/m

HF=1.6;                             %发动机高度/m

xf=2.76;                            %发动机近场测点及等效几何中心离车前端的距离/m

yf=1.4;                             %发动机近场测点离车纵向中心线的距离/m

disp('请输入速度v/km.h-1')

v=input('')

g=sqrt((v/3.6/a)^2+2*(S+L)/a);

tm=-v/3.6/a+g;                      %总匀加速运行时间/t

s=v/3.6*t+0.5*a*t.^2;               %t时刻车辆运行距离/m

rf=sqrt((10-s+xf).^2+7.5^2+(1.2-HF)^2);     %t时刻发动机距离传声器距离/m  

r1=yf;                             %近场声压测试距离

程序一直写不出来，还请懂的朋友帮忙看看

我们以发动机为例，在发动机等效点声源加速通过测试路段时，所用的时间差不多是2.236秒，而测试数据采样频率为51200，频率步长为1Hz,但是最后取的数值点只有8193个，为保证频率步长为1，频率点数为8193，则在fft的时候，时间步长将为1/(8193*2),而总时间为2.236，那么时间步长将远大于频率步长，无法对整个测试路段进行fft分析，在满足频率步长与频率点的情况下，只能求得车辆差不多开到一半位置时的情况，这个。。。
可能是我太笨了，理解上有错误，反正就是不会写，如果可以，请帮忙