interp,interp

插值是采样率增加整数倍。

库::`工程师`

语法

函数调用

[参数:y]=interp(<参数:x>>,<参数:r>>) -增加输入信号的采样率 [参数:x] 在 [参数:r] 有一次。

[参数:y]=interp(<参数:x>>,<参数:r>>,<参数:n>>,<参数:截止>>) -定义两个附加值:
- [参数:n] -用于对扩展信号进行插值的初始采样值的数量的一半。
- [参数：截止] -输入信号的归一化截止频率，设置为奈奎斯特频率的一小部分。

[参数:y],[参数:b]=interp(<参数:x>>,<参数:r>>,<参数:n>>,<参数:截止>>) -还返回用于插值的滤波器系数。

争论

输入参数

xref/# x — 输入信号

+ 向量资料

输入信号设置为矢量。

数据类型: 漂浮64, 漂浮物32</无翻译>

xref/# *r*是的插值系数

+ 一个正整数

的插值系数，设为正整数。

数据类型: 漂浮64, 漂浮物32</无翻译>

xref/# *n*是用于插值的输入样本值数量的一半

+ 4 （默认情况下）| 一个正整数

用于插值的输入样值个数的一半，设为正整数。为获得最佳效果，请使用 n 不超过 10. 低通插值滤波器的长度为 2×n×<参数:r>> + 1.

数据类型: 漂浮64, 漂浮物32</无翻译>

xref/# 截止 — 归一化截止频率

+ 0.5 （默认情况下）| 正标量

输入信号的归一化截止频率，设为不超过 1，这是奈奎斯特频率的一部分。意义 1 这意味着信号占据整个奈奎斯特区间。

数据类型: 漂浮64, 漂浮物32</无翻译>

输出参数

xref/# y — 内插信号

+ 向量资料

作为矢量返回的内插信号。论点 y 在 [参数:r] 比原始输入向量长倍 [参数:x].

数据类型: 漂浮64, 漂浮物32</无翻译>

xref/# b — 低通内插滤波器的系数

+ 列向量

返回的低通内插滤波器的系数作为列向量。

数据类型: 漂浮64, 漂浮物32</无翻译>

例子：

信号插值

让我们创建一个采样频率的正弦信号 1 千赫。我们用系数进行插值 4.

import EngeeDSP.Functions: interp

t = 0:1/1e3:1
x = sin.(2π*30*t) .+ sin.(2π*60*t)
y = interp(x, 4).y

让我们绘制原始信号和插值信号。

plot(
    plot(0:30, x[1:31], seriestype = :stem, markershape = :circle, markerstrokewidth = 0, markersize = 3,
         xlabel = "Sample Number", ylabel = "Original", grid = true),
    plot(0:120, y[1:121], seriestype = :stem, markershape = :circle, markerstrokewidth = 0, markersize = 3,
         xlabel = "Sample Number", ylabel = "Interpolated", grid = true),
    layout = (2, 1)
)

interp

算法

插值增加了序列的初始采样率。这与变薄相反。功能 interp 将零插入到原始信号中，然后将低通内插滤波器应用于扩展序列。该函数使用[1]中描述的低通插值算法8.1:

通过在原始数据值之间插入零，将输入向量扩展到所需的长度。
设计一个特殊的对称FIR滤波器，它将原始数据不变地传递并对其进行插值，以最小化插值点与其理想值之间的均方根误差。函数使用的过滤器 interp，匹配返回的过滤器 [医]内.
对展开的输入向量应用滤波器以获得输出值。

文学作品

Ieee Acoustics,Speech,And Signal Processing Society,eds的数字信号处理委员会。 _用于数字信号处理的程序。_纽约：IEEE出版社，1979年。
Oetken，G.，Thomas W.Parks和H.W.Schüssler。 "数字内插器设计的新成果。"_ieee^®关于声学、语音和信号处理的事务。_卷。 ASSP-23，第3号，1975年6月，第301-309页。