数字波导网络算法研究.pdfVIP

数字波导网络算法研究 田玲爱 朴胜春 (哈尔滨工程大学水声工程学院哈尔滨150001) The Network researchofthe DigitalWaveguide Tian Piao Lingai Shengchun ofUnderwaterAcoustic (HarbinEngineeringUniversity,College Engineering，150001) 1引言 Smith教授于上世纪 80年代中期首次提出的一种波动方程求解方法…，它由时域有限差分法演变发展得到，具有算法简单、各 参数物理量意义直观、清晰的优点。此方法己成功应用于音乐的数字合成等领域。 2数字波导网络法的基本原理 数字波导网络法的基础就是将振动系统或声学系统在时间域和空间域进行离散化，如图1所示，数字 波导中离散的点称为散射结点，结点上的物理量称为波变量”1。数字波导网络法就是根据类比法将振动系 统或声学系统通过电网络求解，计算波变量在多 维空间的散射和传播。该网络的基本组成是双向 时延线路和散射结点，声波在传播过程中被分解 为各向行波的迭加”。，当入射到散射结点时发生 散射。 在阻抗类比中，电压U类比于振动系统中的 力，电流，类比于振动系统中的速度；而对于声 学系统，电压V类比于声压，电流，类比于体积 速度。对于导纳类比则反之。另外还需引入物理 量电阻抗Z，对应于力阻抗和声阻抗。在振动系 图l数字波导网络和结点结构图 统中，质点的各向振动速度相等，当达到力平衡 时，各向的应力之和为零，所以有 v1=V2=…=vⅣ 石+五+…+^=0 (1) 导纳类比所得结点称为并联结点，如图I中右图所示，并联结点中各端口的电压与结点电压相等，电流之 和为零(根据并联电路中的基尔霍夫电流定律)，那么对于一个无损耗的有N个端口的并联结点有： U“=U“∑lq。l=0 (2) 其中q表示各端口，如图1中二维网络的结点有上下左右4个端口，因此对于一维网络而言有端口r，x一， 并有： U“=u主i+uif19．；=，主i十，i， (3) 式(3)中上标“+”表示输入结点的量，一”表示输出结点的量，说明端口电压和端口电流等于该端口输 入输出电压和电流之和。而且电压电流之间满足 ，0=■，，u0 』0=一■，；u0 (4) 其中匕．1是电导纳，对应电阻抗有：Z“=l／‘，，·联立式(2)、(3)、(4)得结点电压： 。 75 u’j-≠l∑‘，，％√ (5) 1J』＼目 ／ 数字波导中连接相邻结点是双向时延线路，散射结点的输入是前一步相邻结点的输出，即： ％．如)=晖j+z(n一1)啦．。(，z)=％,I--i(”一1) (6) 数字波导网络中，所有的延迟时间必须