将“DelayFactor”参数用于`理性`

打开生活的脚本

这个例子展示了如何使用'DelayFactor'参数来改善输出的质量理性.

这个理性函数选择一个与频域数据匹配的有理函数。如果该数据包含一个显著的“时间延迟”，这将在频域表现为相移，那么使用合理数量的极点进行拟合可能会非常困难。

在这些情况下，当输入数据包含较大的负斜率（即具有足够大时间延迟的数据）时，我们可以询问理性首先从数据中删除一些延迟，然后找到一个最适合其余“未延迟”数据的有理函数。这个理性函数通过将其存储在输出的' delay '参数中来解释被删除的延迟。默认情况下,理性不会从数据中删除任何延迟。

首先，从4端口背板s参数创建差分传递函数数据。的默认设置来拟合数据理性最后，使用“DelayFactor”参数来提高理性.

创建传递函数

从'default.s4p'读取4端口背板s参数数据。

S = sparameters ('default.s4p')；fourportdata=S.参数；freq=S.频率；fourportZ0=S.阻抗；

将4端口单端s参数转换为2端口差分s参数

diffdata=s2sdd（四端口数据）；diffZ0=2*fourportZ0；

从差异的2端口数据创建一个传输函数

tfdata = s2tf (diffdata diffZ0、diffZ0 diffZ0);

使用“DelayFactor”的默认值时，分析rationalfit的输出

使用频率响应用于计算输出响应的函数理性.

defaultfit=rationalfit（频率、tfdata）

警告：只有48个极点，而不是40分贝只达到10.2分贝精度。考虑使用“NoLead”参数指定更大数量的极点。

defaultfit = rfmodel。A: [48x1 double] C: [48x1 double] D: 0 Delay: 0 Name: '有理函数'

respfreq=0:4e6:20e9；defaultresp=freqresp（defaultfit，respfreq）；

请注意，“延迟”参数为零（没有从数据中删除延迟）。

绘制原始数据和默认输出理性.

图subplot(2,1,1) tfdataDB = 20*log10(abs(tfdata));情节(频率、tfdataDB'.-')举行在情节(respfreq, 20 * log10 (abs (defaultresp)))从包含(‘频率（Hz）’) ylabel (‘震级（dB）’) defaultnpoles = numel(defaultfit.A);defstr = ['默认延迟因子（使用'，num2str（defaultnpoles），“两极)”];标题(defstr)传说(“原始数据”,“默认理性IT”,“位置”,“最好的”)子批次（2,1,2）tfdataphase=180*展开（角度（tfdata））/pi；绘图（频率、tfdataphase、，'.-')举行在打印（respfreq，180*展开（角度（默认Resp））/pi）保持从包含(‘频率（Hz）’) ylabel (‘角度（度）’)传奇(“原始数据”,“默认理性IT”,“位置”,“最好的”)

图中包含2个轴对象。标题为“默认延迟因子”（使用48个极点）的轴对象1包含2个line类型的对象。这些对象表示原始数据，默认为rationalfit。轴对象2包含2个line类型的对象。这些对象表示原始数据，默认为rationalfit。

注意，使用默认设置时的结果理性很穷。由于原始数据的相位有一个非常大的负斜率，因此可以使用“DelayFactor”参数来提高有理函数的精度。

使用“DelayFactor”的自定义值时，分析rationalfit的输出

“DelayFactor”必须设置为介于0和1之间的值。选择哪个值是一个反复试验的过程。对于某些数据集（其相位整体向上倾斜的数据集），更改“DelayFactor”的值不会对结果产生影响。

保存所有其他可能的参数理性常数，0.98被发现可以创建良好的拟合。

customfit=rationalfit（频率、tfdata、，“延迟因素”,0.98)

customfit = rfmodel。C: [31x1 double] D: 0 Delay: 6.5521e-09 Name: '有理函数'

customresp=freqresp（customfit，respreq）；

请注意，“延迟”参数不是零(理性从数据中删除了一些延迟）。

绘制原始数据与的自定义输出理性.

子地块（2,1,1）图（频率、tfdataDB、，'.-')举行在绘图（respfreq，20*log10（abs（customresp）））保持从包含(‘频率（Hz）’) ylabel (‘震级（dB）’)customnpoles=numel（customfit.A）；customstr=['延迟系数=0.98（使用'，num2str（客户文件），“两极)”];标题(customstr)传说(“原始数据”,“自定义rationalfit”,“位置”,“最好的”)子地块（2,1,2）图（频率、tfdataphase、，'.-')举行在打印（respfreq，180*展开（角度（customresp））/pi）保持从包含(‘频率（Hz）’) ylabel (‘角度（度）’)传奇(“原始数据”,“自定义rationalfit”,“位置”,“最好的”)