编程实现预测系数与反射系数间的转换和反射系数与声管面积间的转换。函数定义如下:1 由预测系数求反射系数名称:lpcar2rf功能:已知预测系数计算出反射系数调用格式:rf=lpcar2rfar说明:ar 是预测系数;rf 是反射系数。2 由反射系数求出预测系数名称:lpcrf2ar功能:已知反射系数计算出预测系数调用格式:ararparug=lpcrf2arrf说明:rf 是反射系数;ar 是预测
以下是Python实现的代码:
import numpy as np
def lpcar2rf(ar): p = len(ar) rf = np.zeros(p) for i in range(p): k = ar[i] for j in range(i): k += rf[j] * ar[i-j-1] rf[i] = -k / (1 + k) return rf
def lpcrf2ar(rf): p = len(rf) ar = np.zeros(p) arp = np.zeros(p) aru = np.zeros(p) g = np.zeros(1) ar[0] = rf[0] arp[0] = 1 aru[0] = 1 g[0] = np.sqrt(1 - rf[0] * rf[0]) for i in range(1, p): k = rf[i] for j in range(i): k += rf[j] * ar[i-j-1] ar[i] = -k for j in range(i): arp[j+1] = arp[j] + ar[i] * arp[i-j-1] aru[j+1] = aru[j] + ar[i] * aru[i-j-1] arp[0] = -ar[i] aru[0] = 1 return ar, arp, aru, g
def lpcrf2ao(rf): return np.abs((1 + rf) / (1 - rf))
def lpcao2rf(ao): return (ao - 1) / (ao + 1)
原文地址: https://www.cveoy.top/t/topic/b99G 著作权归作者所有。请勿转载和采集!