总声压幅值为:p(x,y,z)=ω*ρ*|v(x,y)*∑_(nx=0)^∞▒∑_(ny=0)^∞▒〖RP〗_(nx,ny) |
式中〖RP〗_(nx,ny)={█((2*c0*sin〖(π*(c-z/lamda)*√(4-[(nx/a)^2+(ny/b)^2 ] )〗)/(√(4-[(nx/a)^2+(ny/b)^2])*cos〖(π*c0*√(4-[(nx/a)^2+(ny/b)^2])〗 )@(2*c0*〖(c-z/lamda)*sinh〗〖[π*√([(nx/a)^2+(ny/b)^2 ]-4)〗)/(√([(nx/a)^2+(ny/b)^2 ]-4)*cosh〖(π*√([(nx/a)^2+(ny/b)^2 ]-4)〗 ))┤
(nx,ny)对应一组波,所以不同的(nx,ny)对应不同的波。矩形槽的尺寸Lx=a*lamda=280mm;Ly=b*lamda=200mm;Lz=c*lamda=200mm;
由总声压幅值可以看出:声压幅值p(x,y,z)沿(x,y)方向的分布取决于v(x,y),沿z方向的相对变化可表示为RP=|∑_(nx=0)^∞▒∑_(ny=0)^∞▒〖RP〗_(nx,ny) |
(在程序中,我已经根据计算取nx=0:100,ny=0:100)
老师,程序是独立的v(x,y)和RP,这里我想绘出的是总声压的图(用三维表四维的图,x,y,z轴依旧是空间位置位置坐标,用颜色来代表对应声压的大小)。还有在编程中用RP乘以v(x,y)是否能绘出空间中总声压p的变化呢?您帮我看看
已退回1积分