信号分析——导数(Java/Matlab)

信号求导，可以消除信号的背景，确定谱峰的位置，改善谱峰分辨情况，一般地，在信号分析中，常用的信号求导为一阶、二阶导数。基于S-G滤波的基础上，获取拟合参数矩阵A后，可以利用滤波窗口的2m+1个数据对窗口中心点（i=0）求该点信号处的p阶导数，对其多项式：

进行求导，当i=0时，可得：

因此，窗口中心点数据的p阶导数等于A矩阵的第p+1个元素和p的阶乘的乘积。其一阶导数和二阶导数分别如下：

由于上一篇博客中已知A：

而系数矩阵X也已经可以获取，这里不再赘述，直接上代码。

Matlab：

function [ der,D,A ] = Differential_derivation( xdata,n,k,order)
%信号求导一般方法 一阶求导，二阶求导
%Input: xdata-原始数据
%       n-窗口宽度
%       k-多项式阶数
%       order-阶数
%Output:der-导数
m=(n-1)/2;%数据步长
X=[];%S-G平滑系数矩阵
for i=0:(n-1)for j=0:(k-1)X(i+1,j+1)=power(i-m,j);end
end
D=inv(X'*X)*X';
A=zeros(k,n);%系数矩阵
temp=0;
for i=0:k-1for j=0:n-1temp = temp +D(i+1,j+1)*xdata(j+1);endA(i+1,j+1)=temp;
end
%求导
der=factorial(order)*A(order+1);
end

Java：

public static double derivation(RealMatrix xdata,int n,int k,int order) {int m = (n-1)/2;RealMatrix X_matrix = MatrixUtils.createRealMatrix(n, k );for (int i = 0; i < n; i++) {for (int j = 0; j < k; j++) {X_matrix.setEntry(i, j, Math.pow(i-m, j));}}RealMatrix D_matrix = MatrixAlgorithm.inverseMatrix(X_matrix.transpose().multiply(X_matrix)).multiply(X_matrix.transpose());System.out.println(D_matrix);RealMatrix A_matrix = D_matrix.multiply(xdata);double der = factorial(order)*A_matrix.getEntry(order+1,0);System.out.println(der);return der;
}

——参考《化学计量学方法及MATLAB实现》史永刚等编著

信号分析——导数(Java/Matlab)相关推荐

【肌电信号】脉搏信号分析(去噪+特征提取)matlab源码含GUI
一.简介 EMG,应用电子学仪器记录肌肉静止或收缩时的电活动,及应用电刺激检查神经.肌肉兴奋及传导功能的方法.英文简称EMG.通过此检查可以确定周围神经.神经元.神经肌肉接头及肌肉本身的功能状态.\ ...
matlab fftshift_matlab下实现FFT信号分析
奈奎斯特定理: 对于信号分析,我们需要掌握一个定理: 在进行模拟/数字信号的转换过程中,当采样频率 fs.max大于信号中最高频率 fmax 的 2 倍时(fs.max > 2fmax),采样之 ...
matlab语音信号分析实验报告,DSP实验二语音信号分析与处理,南京理工大学紫金学院实验报告,信号与系统...
实验二语音信号分析与处理学号姓名注:1)此次实验作为<数字信号处理>课程实验成绩的重要依据,请同学们认真.独立完成,不得抄袭. 2)请在授课教师规定的时间内完成: 3)完成作业后,请以 ...
将示波器信号完整数据导入Matlab进行分析
学习如何将麦科信示波器信号的完整数据导入matlab进行分析,我们先要学习如何完整导出示波器的信号数据.示波器可将模拟通道或数学通道波形保存到本地或者U盘,文件类型可选择WAV.CSV或BIN. WA ...
【Matlab语音分析】语音信号分析【含GUI源码 1718期】
一.代码运行视频(哔哩哔哩) [Matlab语音分析]语音信号分析[含GUI源码 1718期] 二.matlab版本及参考文献 1 matlab版本 2014a 2 参考文献 [1]韩纪庆,张磊,郑铁 ...
《MATLAB语音信号分析与合成（第二版）》：第4章语音信号的线性预测分析
<MATLAB语音信号分析与合成(第二版)>:第4章语音信号的线性预测分析前言 1. 数据与函数路径设置 2. MATLAB仿真一:LPC的频谱分析 3. MATLAB仿真二:普通预测 ...
matlab symadd,信号调理器设计—大二暑期信号分析实习报告精选多篇
第1篇第2篇第3篇第4篇第5篇更多顶部目录第一篇:信号调理器设计-大二暑期信号分析实习报告第二篇:大二暑期信号分析实习报告第三篇:传感器信号调理电路第四篇:信号调理电路第五篇:信号实习报告 ...
基于Matlab App Designer的语音信号分析与处理（二）：IIR和FIR滤波器的设计，语音信号的滤波
接上文:https://blog.csdn.net/weixin_53877178/article/details/122470759 目录一.课题的任务二.内容.步骤和要求 (1)语音信号的采集 ...
《MATLAB语音信号分析与合成（第二版）》：第8章基音周期的估算方法
<MATLAB语音信号分析与合成(第二版)>:第8章基音周期的估算方法前言 1. 数据与函数路径设置 2. MATLAB仿真一:基音周期提取的预处理 3. MATLAB仿真二:倒谱法的 ...

信号分析——导数(Java/Matlab)

信号分析——导数(Java/Matlab)相关推荐

最新文章

热门文章