cvxEDA 用法

三十八画生

已于 2022-03-06 14:37:03 修改

阅读量409

点赞数

分类专栏： Python自学遇到的问题文章标签： python

于 2021-08-29 15:39:09 首次发布

本文链接：https://blog.csdn.net/weixin_42764972/article/details/119981213

版权

该博客介绍了一种基于凸优化的EDA信号处理方法cvxEDA，用于从EDA数据中提取时域的交感神经活动（phasic）和静息（tonic）成分。通过读取CSV文件中的EDA数据，应用cvxEDA算法，最终得到的输出包括交感神经和静息成分的数值。这个算法结合了ARMA模型、样条插值和凸优化技术，为生物信号处理提供了一个有效的工具。

摘要由CSDN通过智能技术生成

import numpy as np
import cvxopt as cv
import cvxopt.solvers
import matplotlib.pyplot as plt


#从csv文件里读取EDA信号，添加到列表中
import matplotlib.pyplot as plt
import csv
import pandas as pd

filename='F:\\E4手环\data\EDA.csv'
un=pd.read_csv('F:\\E4手环\data\EDA.csv')

with open(filename) as f:
    reader=csv.reader(f)
    header_row=next(reader)
    index=[]
    highs=[]
    for i in reader:
        high=float(i[0])
        highs.append(high)



def cvxEDA(eda, sampling_rate=1000, tau0=2., tau1=0.7, delta_knot=10., alpha=8e-4, gamma=1e-2, solver=True,
           verbose=False, options={'reltol': 1e-9}):
    """
    A convex optimization approach to electrodermal activity processing (CVXEDA).
    This function implements the cvxEDA algorithm described in "cvxEDA: a
    Convex Optimization Approach to Electrodermal Activity Processing" (Greco et al., 2015).
    Parameters
    ----------
       eda : list or array
           raw EDA signal array.
       sampling_rate : int
           Sampling rate (samples/second).
       tau0 : float
           Slow time constant of the Bateman function.
       tau1 : float
           Fast time constant of the Bateman function.
       delta_knot