双管声带模型：计算在声门处激发并在嘴唇处终止的双管模型的频率响应。-matlab开发

在此MATLAB练习中，用户指定了无损，人声声道的两管模型的面积和长度，并且MATLAB练习计算了从声门到嘴唇的体积速度的传递函数，其中反射系数为假定声门为 r_G（在范围-1 ≤ r_G ≤ 1 内），假定嘴唇处的反射系数为 r_L，同样在范围-1 ≤ r_L ≤1 内。 用户指南的图 2.1 显示了这种双管模型（管长度和面积值为 l1 = 8 cm、A1 = 1 cm^2、l2 = 9 cm、A2 = 7 cm2），以及流程图从声门到嘴唇的复合系统的结构如图 2.2 所示。 这个 MATLAB 练习计算从声门到嘴唇的体积速度的传递函数，绘制体积速度传递函数（对数幅度谱），尽可能确定共振峰频率（中心频率和带宽）的位置，并打印它们到指定的ascii文本文件。 对 2 管模型频谱进行采样并转换为 2 管模型脉冲响应，该脉冲响应与周期性声门序列卷积并播放出来，以便用户可以听到使用双管模型产生的