声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

12
返回列表 发新帖
楼主: jinfengshi_2005

[编程技巧] 如何给把矩阵的值赋值

[复制链接]
 楼主| 发表于 2010-8-11 09:06 | 显示全部楼层
本帖最后由 ChaChing 于 2010-8-31 14:46 编辑

我是想用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 编辑 ]

形状函数W.doc

57 KB, 下载次数: 1

形状函数W

回复 支持 反对
分享到:

使用道具 举报

 楼主| 发表于 2010-8-12 10:16 | 显示全部楼层

回复 10楼 jinfengshi_2005 的帖子

for i=1:3
      for j=1:4
     A(i,j)=i+j;
      end
   end
=====================
2        3        4        5
3        4        5        6
4        5        6        7

大家,不好意思,我把事情给搞复杂化了,其实我在10楼提的问题,想要的结果就是上面的这个程序。现在终于搞懂了,突然之间,恍然大悟,但是在搞懂这个问题的过程中还是学到了很多东西。呵呵,这就是王国维的学问四境界吧,心中有说不出的惊喜!
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-6-17 15:40 , Processed in 0.054813 second(s), 20 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表