先有下面的程序  
clear all  
y=wavread('E:\GLASS.WAV');%读取声音信号  
size(y) %获取信号的长度  
y=y(1:1:11488);%读取从 1到11488的信号,间隔为1  
y在这里表示瞬时声压,求用MATLAB计算这个信号的有效声压和声压级并画出声压级谱图  
有效声压Y即在一段时间T内的瞬时声压y对时间取的均方根值,即瞬时声压y对时间的积分除以时间T再开平方.  
声压级:L=20lg(Y/y0) y0--0.00002pa  
麻烦懂的帮下忙,先谢谢了!!!
主要思路是利用输入的信号计算出有效声压,然后利用声压级公式计算声压级,然后有fft对声压级进行变换,最后得出声压级的频谱图,小弟在计算上面2个公式的编程中遇到了困难望达人帮忙,谢谢!