我是想用Matlab编程计算响度，而在计算响度的过程中，必须要计算一个形状函数W（也叫加权函数），W(g(i))=(1+p(i)g(i))*e^(-p(i)*g(i));
因为可听声的频率是20Hz-20000Hz（也就是附件公式中的i),对应于每一个频率，都有个g(i)和p(i)，而g(i)和p(i)就是通过附件中的公式来确定。
1. 附件公式中的F表示的是每个带宽的中心频率，总共有372个中心频率。而我已经把这372个中心频率已经求好放在数组里了。也就是对于1个中心频率，可以求出19981个P(i)和g(i).
2. ERB为对应中心频率F的等效矩形带宽，总共也有372个，也已经放在另外一个数组里了。
3.X为对应频率i处的ERB级，X=（20,30,40，50， 60， 70， 80，90，100），X总共有9个数值，也就是X=20的时候，求出372个中心频率所对应的19981个p(i)和g(i).然后X=30时在求出372个中心频率所对应的19981个p(i)和g(i)，依次类推！
而我在13楼写的程序只是计算了第一个中心频率F1，和X=20时所对应的19981个p(i)和g(i)。
===================================================================
不知道我上面讲的清楚吗？还有就是我这个帖子的题目是不是起的有点文不对题，版主看看呢，如果不适合，帮忙改一下吧！呵呵。
[ 本帖最后由 jinfengshi_2005 于 2010-8-11 09:16 编辑 ]
回复 10楼 jinfengshi_2005 的帖子
for i=1:3for j=1:4
A(i,j)=i+j;
end
end
=====================
2 3 4 5
3 4 5 6
4 5 6 7
大家，不好意思，我把事情给搞复杂化了，其实我在10楼提的问题，想要的结果就是上面的这个程序。现在终于搞懂了，突然之间，恍然大悟，但是在搞懂这个问题的过程中还是学到了很多东西。呵呵，这就是王国维的学问四境界吧，心中有说不出的惊喜！
页:
1
[2]