SG算法解析

Savitzky-Golay 滤波器的核心代码主要集中在计算投影矩阵B并使用这个矩阵对输入信号进行滤波。这部分核心代码包括计算B矩阵、处理边界效应和进行实际滤波操作。以下是对核心代码的一点解释:

① 计算 Savitzky-Golay 投影矩阵B

复制代码
B = sgolay(order, framelen, weights); % 计算 Savitzky-Golay 投影矩阵

Savitzky-Golay 投影矩阵B是通过多项式拟合计算出来的。这个矩阵用于对输入信号进行平滑处理。函数 sgolay 负责生成这个矩阵。

② 处理数据的维度

复制代码
if isempty(dim)
    [x, nshifts] = shiftdim(x); % 沿第一个非单维度工作
else
    perm = [dim, 1:dim-1, dim+1:ndims(x)];
    x = permute(x, perm); % 将 DIM 置于第一维度
end

if size(x, 1) < framelen, error(message('signal:sgolayfilt:InvalidDimensionsTooSmall')), end

这段代码确保信号x的长度至少为 framelen,并将数据沿指定维度进行处理。

③ 计算滤波结果

前端过渡

复制代码
ybegin = B(end:-1:(framelen-1)/2+2,:) * x(framelen:-1:1,:);

这部分代码计算信号的前端过渡部分。通过倒序取出 framelen 个数据点并与矩阵B相乘,计算出平滑后的前端部分。

稳态输出

复制代码
ycenter = filter(B((framelen-1)/2+1,:), 1, x);

这部分代码计算信号的稳态输出部分。通过应用 filter 函数进行卷积操作,实现对信号中心部分的平滑处理。

后端过渡

复制代码
yend = B((framelen-1)/2:-1:1,:) * x(end:-1:end-(framelen-1),:);

这部分代码计算信号的后端过渡部分。通过倒序取出信号末尾的 framelen 个数据点并与矩阵B相乘,计算出平滑后的后端部分。

合并结果

复制代码
y = [ybegin; ycenter(framelen:end,:); yend];

将前端过渡部分、稳态输出部分和后端过渡部分合并,得到完整的平滑信号。

④ 恢复原始数据的形状

复制代码
if isempty(dim)
    y = shiftdim(y, -nshifts); % 恢复原始维度
else
    y = ipermute(y, perm); % 恢复原始维度顺序
end

将平滑后的信号恢复到与输入信号相同的维度顺序。

相关推荐
Zheng照邻、4 分钟前
VLM Prompt优化之 DynaPrompt(ICLR 2025)论文总结
人工智能·算法·语言模型·prompt·aigc
CoovallyAIHub5 分钟前
机器人“大脑”遭遇认知冻结攻击!复旦等提出FreezeVLA,一张图片即可瘫痪多模态大模型
深度学习·算法·计算机视觉
无敌最俊朗@31 分钟前
死锁 (Deadlock) 深度解析
算法
西阳未落1 小时前
欧拉路径与欧拉回路
算法·深度优先
Swift社区1 小时前
LeetCode 390 消除游戏
算法·leetcode·游戏
橘颂TA2 小时前
【剑斩OFFER】优雅的解法——三数之和
算法
我爱工作&工作love我2 小时前
2024-CSP-J T3 小木棍
算法·动态规划
DatGuy2 小时前
Week 18: 深度学习补遗:Stacking和量子运算Deutsch算法
人工智能·深度学习·算法
Nie_Xun4 小时前
ROS1 go2 vlp16 局部避障--3 篇
算法
Da Da 泓7 小时前
LinkedList模拟实现
java·开发语言·数据结构·学习·算法