diff --git a/proc/Untitled.ipynb b/proc/Untitled.ipynb index 97facfc..b01e90b 100644 --- a/proc/Untitled.ipynb +++ b/proc/Untitled.ipynb @@ -1,692 +1,757 @@ { "cells": [ { "cell_type": "code", "execution_count": null, "metadata": { "ExecuteTime": { - "end_time": "2020-02-13T17:16:55.434600Z", - "start_time": "2020-02-13T17:16:54.920370Z" + "end_time": "2020-02-14T16:59:52.224327Z", + "start_time": "2020-02-14T16:59:51.683530Z" } }, "outputs": [], "source": [ "import numpy as np\n", "import os\n", "import pandas as pd \n", "import sys\n", "import matplotlib.pyplot as plt\n", "import eventBased as eb\n", "\n", "dataFolder = '../data'\n", "experiment = 4\n", "sub = 72\n", "dataRwLb = os.path.join(dataFolder, \"dataLabeld\")\n", "dataEBLb = os.path.join(dataFolder, \"dataEB/sub\"+str(sub))\n", "baseFreq = 50\n", "dT = 1/baseFreq\n", "\n", "start = 10000\n", "stop = int(start+baseFreq*2.56)" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "ExecuteTime": { - "end_time": "2020-02-13T17:16:55.437416Z", - "start_time": "2020-02-13T17:16:55.435669Z" + "end_time": "2020-02-14T16:59:52.231439Z", + "start_time": "2020-02-14T16:59:52.225375Z" + } + }, + "outputs": [], + "source": [ + "fileHigh = './highPass.fir'\n", + "tapsHigh = list(pd.read_csv(fileHigh,header = None)[0])\n", + "fileLow = './lowPass.fir'\n", + "tapsLow = list(pd.read_csv(fileLow,header = None)[0])" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "metadata": { + "ExecuteTime": { + "end_time": "2020-02-14T16:59:52.236961Z", + "start_time": "2020-02-14T16:59:52.233237Z" } }, "outputs": [], "source": [ "plt.close('all')" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "ExecuteTime": { - "end_time": "2020-02-13T17:16:55.474704Z", - "start_time": "2020-02-13T17:16:55.438876Z" + "start_time": "2020-02-14T16:59:51.679Z" } }, "outputs": [], "source": [ "%matplotlib auto\n", "np.set_printoptions(precision=3)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## load data " ] }, { "cell_type": "code", "execution_count": null, "metadata": { "ExecuteTime": { - "start_time": "2020-02-13T17:16:54.916Z" + "start_time": "2020-02-14T16:59:51.680Z" } }, "outputs": [], "source": [ "def getExpNo(s):\n", " idx = 0\n", " for i in range(len(s)):\n", " idx = i\n", " if s[i].isalpha():\n", " break\n", " return int(s[0:idx])\n", "def pik(ls, s):\n", " for e in ls:\n", " if s in e:\n", " return e\n", "def generateFreqScale(dt,N, correctCoeff = None):\n", " if correctCoeff:\n", " dt *= correctCoeff\n", " fs = np.arange(-1/(2*dt),1/(2*dt)-dT/10,1/(dt*N))\n", " return fs \n", "\n", "pathToRwFiles = sorted([os.path.join(dataRwLb,x) for x in os.listdir(dataRwLb) if getExpNo(x) == experiment])\n", "pathToEBFiles = sorted([os.path.join(dataEBLb,x) for x in os.listdir(dataEBLb) if getExpNo(x) == experiment])\n", "\n", "#Original data ---------------------------------------\n", "dataAccRw = pd.read_csv(pik(pathToRwFiles, \"acc\"))\n", "dataGyroRw = pd.read_csv(pik(pathToRwFiles, \"gyro\"))\n", "\n", "#EB data ---------------------------------------\n", "dataAccX = pd.read_csv(pik(pathToEBFiles, \"accx\"), header = None)\n", "dataAccY = pd.read_csv(pik(pathToEBFiles, \"accy\"), header = None)\n", "dataAccZ = pd.read_csv(pik(pathToEBFiles, \"accz\"), header = None)\n", "dataGyroX = pd.read_csv(pik(pathToEBFiles, \"gyrox\"), header = None)\n", "dataGyroY = pd.read_csv(pik(pathToEBFiles, \"gyroy\"), header = None)\n", "dataGyroZ = pd.read_csv(pik(pathToEBFiles, \"gyroz\"), header = None)\n", "\n", "#Some infos:\n", "lengthT = len(dataAccRw)/baseFreq\n", "print(\"Length of the signal, in seconds: {}s\".format(lengthT))\n", "print(\"--------------------------------------------------------------------\")\n", "avgFreqAccX = len(dataAccX)/lengthT\n", "print(\"Average frequency for the x accelerometer: {}Hz\".format(avgFreqAccX))\n", "avgFreqAccY = len(dataAccY)/lengthT\n", "print(\"Average frequency for the y accelerometer: {}Hz\".format(avgFreqAccY))\n", "avgFreqAccZ = len(dataAccZ)/lengthT\n", "print(\"Average frequency for the z accelerometer: {}Hz\".format(avgFreqAccZ))\n", "print(\"--------------------------------------------------------------------\")\n", "avgFreqGyroX = len(dataGyroX)/lengthT\n", "print(\"Average frequency for the x gyroscope: {}Hz\".format(avgFreqGyroX))\n", "avgFreqGyroY = len(dataGyroY)/lengthT\n", "print(\"Average frequency for the y gyroscope: {}Hz\".format(avgFreqGyroY))\n", "avgFreqGyroZ = len(dataGyroZ)/lengthT\n", "print(\"Average frequency for the z gyroscope: {}Hz\".format(avgFreqGyroZ))\n" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "ExecuteTime": { - "start_time": "2020-02-13T17:16:54.917Z" + "start_time": "2020-02-14T16:59:51.682Z" } }, "outputs": [], "source": [ - "accx = [list(dataAccX[0]),list(dataAccX[1]),list(dataAccX[2])]\n", - "accy = [list(dataAccY[0]),list(dataAccY[1]),list(dataAccY[2])]\n", - "accz = [list(dataAccZ[0]),list(dataAccZ[1]),list(dataAccZ[2])]\n", + "delta = []\n", + "for lbl in ['x','y','z'] :\n", + " ls = list(dataAccRw[lbl])\n", + " for i in range(1,len(ls)):\n", + " dif = abs(ls[i]-ls[i-1])\n", + " if dif != 0:\n", + " delta.append(dif)\n", + "print(\"Acc Min delta: \", min(delta))\n", + "print(\"Acc Max delta: \", max(delta))\n", + "print(\"Acc Avg delta: \", np.average(delta))\n", + "print()\n", + "minDvAcc = [min(delta),min(delta),min(delta)]\n", "\n", - "gyrox = [list(dataGyroX[0]),list(dataGyroX[1]),list(dataGyroX[2])]\n", - "gyroy = [list(dataGyroY[0]),list(dataGyroY[1]),list(dataGyroY[2])]\n", - "gyroz = [list(dataGyroZ[0]),list(dataGyroZ[1]),list(dataGyroZ[2])]\n", + "delta = []\n", + "for lbl in ['x','y','z'] :\n", + " ls = list(dataGyroRw['z'])\n", + " for i in range(1,len(ls)):\n", + " dif = abs(ls[i]-ls[i-1])\n", + " if dif != 0:\n", + " delta.append(dif)\n", + "print(\"Gyro Min delta: \", min(delta))\n", + "print(\"Gyro Max delta: \", max(delta))\n", + "print(\"Gyro Avg delta: \", np.average(delta))\n", + "print()\n", + "minDvGyro = [min(delta),min(delta),min(delta)]\n", "\n", - "dictSigNames = {'accx': 0,'accy': 1,'accz': 2,'gyrox': 3,'gyroy': 4,'gyroz': 5}\n", - "\n", - "stream = eb.VectorStreamer()\n", - "stream.load(dictSigNames.keys(),accx,accy,accz,gyrox,gyroy,gyroz)" + "minDvTot = minDvAcc\n", + "minDvTot.extend(minDvGyro)\n", + "print(minDvTot)" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "ExecuteTime": { - "start_time": "2020-02-13T17:16:54.919Z" + "start_time": "2020-02-14T16:59:51.683Z" } }, "outputs": [], "source": [ - "f1 = [1,1,1,1]\n", - "f2 = [1,2,3,4,5,6]\n", - "f3 = [-1,-2,-3,-2,-1,0,1,2,3,2,1,0,-1,-2,-3,-2,-1,0]\n", - "a = eb.inertial.inertial(filtersTaps=[f1,f2,f3])" + "accx = [list(dataAccX[0]),list(dataAccX[1]),list(dataAccX[2])]\n", + "accy = [list(dataAccY[0]),list(dataAccY[1]),list(dataAccY[2])]\n", + "accz = [list(dataAccZ[0]),list(dataAccZ[1]),list(dataAccZ[2])]\n", + "\n", + "gyrox = [list(dataGyroX[0]),list(dataGyroX[1]),list(dataGyroX[2])]\n", + "gyroy = [list(dataGyroY[0]),list(dataGyroY[1]),list(dataGyroY[2])]\n", + "gyroz = [list(dataGyroZ[0]),list(dataGyroZ[1]),list(dataGyroZ[2])]\n", + "\n", + "dictSigNames = {'accx': 0,'accy': 1,'accz': 2,'gyrox': 3,'gyroy': 4,'gyroz': 5}\n", + "\n", + "stream = eb.VectorStreamer()\n", + "stream.load(dictSigNames.keys(),accx,accy,accz,gyrox,gyroy,gyroz)" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "ExecuteTime": { - "start_time": "2020-02-13T17:16:54.920Z" + "start_time": "2020-02-14T16:59:51.685Z" } }, "outputs": [], "source": [ - "l = np.array([1,2,3,4,5])\n", - "q = [5,7,9]\n", - "q.extend(l)\n", - "q" + "a = eb.inertial.inertial(filtersTaps=[tapsLow,tapsHigh], dv = minDvTot)" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "ExecuteTime": { - "start_time": "2020-02-13T17:16:54.922Z" + "start_time": "2020-02-14T16:59:51.687Z" }, "scrolled": true }, "outputs": [], "source": [ "numWindow = 0\n", - "batchWindow = 10\n", + "batchWindow = 50\n", "state = 'train'\n", "windowPass = False\n", "\n", "for elem in stream.get():\n", " t,v,lbl,dim = elem[0],elem[1],elem[2],dictSigNames[elem[3]]\n", " if state == 'train':\n", " windowPass = a.buildTrain(t,v,lbl,dim)\n", " else:\n", - " res,orig = a.test(t,v,lbl,dim)\n", + " res,orig,lblsTot = a.test(t,v,lbl,dim)\n", " if type(orig) != type(None):\n", " print()\n", " print(\"/*\"*35)\n", - " print(\"Results: {} vs original: {}\".format(res,orig))\n", + " print(\"Results: {}({} --> {}) vs original: {}\".format(res,lblsTot,lblsTot[np.argmax(res)],orig))\n", " print(\"/*\"*35)\n", " print()\n", " state = 'train'\n", " if windowPass:\n", " windowPass = False\n", " numWindow +=1\n", " print(\"Window Number: \",numWindow)\n", " if numWindow % batchWindow == 0:\n", " success = a.fit()\n", "# print(\"/*\"*35)\n", "# print(success)\n", "# print(\"/*\"*35)\n", " if success:\n", " state = 'test'\n", " else:\n", " print(\"Error in fitting !\")\n", " print(\"----------------------------------------------------\")\n", " \n", " " ] }, { "cell_type": "code", "execution_count": null, "metadata": { "ExecuteTime": { - "start_time": "2020-02-13T17:16:54.924Z" + "start_time": "2020-02-14T16:59:51.687Z" } }, "outputs": [], "source": [ "from sklearn.svm import SVC\n", "svm = SVC(kernel='linear', probability = True, decision_function_shape = \"ovr\")\n", "vec = [[1,2,3,4],[7,6,8,9],[1,2,90,5],[4,6,78,0]]\n", "lbl = [4,4,3,3]" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "ExecuteTime": { - "start_time": "2020-02-13T17:16:54.926Z" + "start_time": "2020-02-14T16:59:51.689Z" } }, "outputs": [], "source": [ "a = svm.fit(vec,lbl)" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "ExecuteTime": { - "start_time": "2020-02-13T17:16:54.928Z" + "start_time": "2020-02-14T16:59:51.691Z" } }, "outputs": [], "source": [ "# vec.append([1,4,88,2])\n", "# lbl.append(0)\n", "svm.predict_proba(vec)" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": { "ExecuteTime": { - "start_time": "2020-02-13T17:16:54.930Z" + "start_time": "2020-02-14T16:59:51.692Z" } }, "outputs": [], "source": [ "eb = signalDescriptor(dataAccX[0],dataAccX[1])\n", "rw = signalDescriptor(np.arange(len(dataAccRw['x'])),dataAccRw['x'])" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "ExecuteTime": { - "start_time": "2020-02-13T17:16:54.932Z" + "start_time": "2020-02-14T16:59:51.694Z" } }, "outputs": [], "source": [ "fEB = eb.nudft(startT = start, endT = stop)\n", "fRW = rw.nudft(startT = start, endT = stop-1)\n", "print(\"Finished nuDft\")\n", "fEBres = eb.resampledDFT(startT = start, endT = stop)\n", "fRWres = rw.resampledDFT(startT = start, endT = stop)\n", "print(\"Finished reDft\")\n", "f = np.fft.fft(dataAccRw['x'][start:stop])\n", "print(\"Finished numpy fft\")" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "ExecuteTime": { - "start_time": "2020-02-13T17:16:54.933Z" + "start_time": "2020-02-14T16:59:51.695Z" } }, "outputs": [], "source": [ "fRWScale = generateFreqScale(dT,len(f))\n", "fEBScale = generateFreqScale(dT,len(fEB),len(f)/len(fEB))\n", "#--------------------------------------\n", "plt.figure()\n", "plt.title(\"data Raw\")\n", "plt.plot(dataAccRw['x'][start:stop])\n", "\n", "plt.figure()\n", "plt.title(\"data EB sampled\")\n", "t,v = eb.findInterval(start,stop)\n", "plt.plot(t,v)\n", "#--------------------------------------\n", "plt.figure()\n", "plt.title(\"Non uniform DFT on Event Based\")\n", "plt.plot(fEBScale,np.abs(eb.shift(fEB)))\n", "\n", "plt.figure()\n", "plt.title(\"Non uniform DFT on Raw\")\n", "plt.plot(fRWScale,np.abs(rw.shift(fRW)))\n", "#---------------------------------------\n", "plt.figure()\n", "plt.title(\"uniform DFT on re-sampled Event Based\")\n", "plt.plot(fRWScale,np.abs(eb.shift(fEBres)))\n", "\n", "plt.figure()\n", "plt.title(\"uniform DFT on raw signal\")\n", "plt.plot(fRWScale,np.abs(rw.shift(fRWres)))\n", "#---------------------------------------\n", "plt.figure()\n", "plt.title(\"Numpy FFT on the original signal\")\n", "plt.plot(fRWScale,np.abs(rw.shift(f)))" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": { "ExecuteTime": { - "start_time": "2020-02-13T17:16:54.947Z" + "start_time": "2020-02-14T16:59:51.710Z" } }, "outputs": [], "source": [ "import scipy.stats as stats" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "ExecuteTime": { - "start_time": "2020-02-13T17:16:54.947Z" + "start_time": "2020-02-14T16:59:51.711Z" } }, "outputs": [], "source": [ "stats.entropy([1/2, 1/2], base=2)" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "ExecuteTime": { - "start_time": "2020-02-13T17:16:54.949Z" + "start_time": "2020-02-14T16:59:51.711Z" } }, "outputs": [], "source": [ "a = [7,7,7,7,7,8,9,9,9,9,9,3,3,2,1,1,-4,-3,-9]\n", "sum(np.abs(a))" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "ExecuteTime": { - "start_time": "2020-02-13T17:16:54.950Z" + "start_time": "2020-02-14T16:59:51.713Z" } }, "outputs": [], "source": [ "a = [7,7,7,7,7,8,9,9,9,9,9,3,3,2,1,1]\n", "np.median(a)" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "ExecuteTime": { - "start_time": "2020-02-13T17:16:54.951Z" + "start_time": "2020-02-14T16:59:51.714Z" } }, "outputs": [], "source": [ "len(a)" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "ExecuteTime": { - "start_time": "2020-02-13T17:16:54.952Z" + "start_time": "2020-02-14T16:59:51.715Z" } }, "outputs": [], "source": [ "sorted(a)" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "ExecuteTime": { - "start_time": "2020-02-13T17:16:54.954Z" + "start_time": "2020-02-14T16:59:51.717Z" } }, "outputs": [], "source": [ "b = a[0:11]\n", "c = a[11:len(a)]\n", "\n", "mb = np.median(b)\n", "mc = np.median(c)\n", "q = " ] }, { "cell_type": "code", "execution_count": null, "metadata": { "ExecuteTime": { - "start_time": "2020-02-13T17:16:54.955Z" + "start_time": "2020-02-14T16:59:51.718Z" } }, "outputs": [], "source": [ "a = [2,3,5,6,8,4,345,567,3]\n", "sum([x**2 for x in a])" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, - "metadata": {}, + "metadata": { + "ExecuteTime": { + "start_time": "2020-02-14T16:59:51.723Z" + } + }, "outputs": [], - "source": [] + "source": [ + "taps = tapsLow\n", + "f = 19.5\n", + "t = np.arange(0.1,100,0.02)\n", + "y = np.sin(2*np.pi*f*t)\n", + "delay = (len(taps)-1)//2\n", + "filt = []\n", + "tFilt =[]\n", + "gain = []\n", + "\n", + "for i in range(delay,len(t)-delay):\n", + " res = np.dot(y[i-delay:i+delay+1],taps)\n", + " filt.append(res)\n", + " tFilt.append(i*0.02)\n", + "\n", + "filt = np.array(filt)\n", + "\n", + "plt.figure()\n", + "plt.title('filtered')\n", + "plt.plot(tFilt,filt)\n", + "#plt.figure()\n", + "#plt.title('Original')\n", + "plt.plot(t,y)" + ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.7.1" }, "varInspector": { "cols": { "lenName": 16, "lenType": 16, "lenVar": 40 }, "kernels_config": { "python": { "delete_cmd_postfix": "", "delete_cmd_prefix": "del ", "library": "var_list.py", "varRefreshCmd": "print(var_dic_list())" }, "r": { "delete_cmd_postfix": ") ", "delete_cmd_prefix": "rm(", "library": "var_list.r", "varRefreshCmd": "cat(var_dic_list()) " } }, "types_to_exclude": [ "module", "function", "builtin_function_or_method", "instance", "_Feature" ], "window_display": false } }, "nbformat": 4, "nbformat_minor": 2 } diff --git a/proc/highPass.fir b/proc/highPass.fir new file mode 100644 index 0000000..49acb96 --- /dev/null +++ b/proc/highPass.fir @@ -0,0 +1,341 @@ +-1.87744309766349e-05 +-1.32506328857496e-05 +-7.64591082430952e-06 +-1.93308634239090e-06 +3.91563251130230e-06 +9.92849977621546e-06 +1.61340752645808e-05 +2.25610638103370e-05 +2.92381505669683e-05 +3.61938328576359e-05 +4.34562490951040e-05 +5.10530053056054e-05 +5.90109998022162e-05 +6.73562465662804e-05 +7.61136979045506e-05 +8.53070669589730e-05 +9.49586506519924e-05 +0.000105089153656573 +0.000115717513982355 +0.000126860730772267 +0.000138533694902904 +0.000150749022981455 +0.000163516895326989 +0.000176844898520315 +0.000190737873098248 +0.000205197766960386 +0.000220223495045263 +0.000235810805821544 +0.000251952155124712 +0.000268636587855687 +0.000285849628039315 +0.000303573177723108 +0.000321785425174916 +0.000340460762818218 +0.000359569715318499 +0.000379078878211122 +0.000398950867434119 +0.000419144280102709 +0.000439613666832982 +0.000460309515894002 +0.000481178249435283 +0.000502162232006304 +0.000523199791551049 +0.000544225253028676 +0.000565168984775198 +0.000585957457688754 +0.000606513317283699 +0.000626755468624841 +0.000646599174114904 +0.000665956164074958 +0.000684734760017794 +0.000702840010480627 +0.000720173839244429 +0.000736635205733962 +0.000752120277352463 +0.000766522613473622 +0.000779733360774545 +0.000791641459561057 +0.000802133860700985 +0.000811095752749854 +0.000818410798817751 +0.000823961382698209 +0.000827628863746419 +0.000829293839967222 +0.000828836418741982 +0.000826136494601106 +0.000821074033418418 +0.000813529362383507 +0.000803383465083609 +0.000790518281006618 +0.000774817008755975 +0.000756164412253804 +0.000734447129189208 +0.000709553980957602 +0.000681376283322907 +0.000649808157027262 +0.000614746837560138 +0.000576092983298052 +0.000533750981217625 +0.000487629249386302 +0.000437640535430580 +0.000383702210190467 +0.000325736555765667 +0.000263671047171732 +0.000197438626827597 +0.000126977971112919 +5.22337482373576e-05 +-2.68431333115253e-05 +-0.000110295286466196 +-0.000198158618124377 +-0.000290462113224655 +-0.000387227631009619 +-0.000488469714115528 +-0.000594195411093534 +-0.000704404112946686 +-0.000819087404231046 +-0.000938228929244725 +-0.00106180427378789 +-0.00118978086295128 +-0.00132211787534665 +-0.00145876617416354 +-0.00159966825539141 +-0.00174475821351342 +-0.00189396172493222 +-0.00204719604935186 +-0.00220437004929397 +-0.00236538422788722 +-0.00253013078502138 +-0.00269849369191971 +-0.00287034878413194 +-0.00304556387291290 +-0.00322399887490026 +-0.00340550595996820 +-0.00358992971707959 +-0.00377710733792560 +-0.00396686881808716 +-0.00415903717541873 +-0.00435342868530286 +-0.00454985313238936 +-0.00474811407838460 +-0.00494800914542088 +-0.00514933031449211 +-0.00535186423840642 +-0.00555539256866420 +-0.00575969229564071 +-0.00596453610140895 +-0.00616969272451323 +-0.00637492733596383 +-0.00658000192570287 +-0.00678467569874974 +-0.00698870548022275 +-0.00719184612839416 +-0.00739385095492654 +-0.00759447215140356 +-0.00779346122126418 +-0.00799056941621786 +-0.00818554817621721 +-0.00837814957204084 +-0.00856812674953827 +-0.00875523437457423 +-0.00893922907770915 +-0.00911986989764374 +-0.00929691872246433 +-0.00947014072771391 +-0.00963930481033161 +-0.00980418401749797 +-0.00996455596944049 +-0.0101202032752595 +-0.0102709139408508 +-0.0104164817680126 +-0.0105567067438515 +-0.0106913954196074 +-0.0108203612780558 +-0.0109434250886544 +-0.0110604152496426 +-0.0111711681163152 +-0.0112755283147354 +-0.0113733490401726 +-0.0114644923395950 +-0.0115488293775723 +-0.0116262406849887 +-0.0116966163899972 +-0.0117598564306972 +-0.0118158707490443 +-0.0118645794655557 +-0.0119059130344115 +-0.0119398123785993 +-0.0119662290047922 +-0.0119851250977014 +-0.0119964735936836 +0.988021261219908 +-0.0119964735936836 +-0.0119851250977014 +-0.0119662290047922 +-0.0119398123785993 +-0.0119059130344115 +-0.0118645794655557 +-0.0118158707490443 +-0.0117598564306972 +-0.0116966163899972 +-0.0116262406849887 +-0.0115488293775723 +-0.0114644923395950 +-0.0113733490401726 +-0.0112755283147354 +-0.0111711681163152 +-0.0110604152496426 +-0.0109434250886544 +-0.0108203612780558 +-0.0106913954196074 +-0.0105567067438515 +-0.0104164817680126 +-0.0102709139408508 +-0.0101202032752595 +-0.00996455596944049 +-0.00980418401749797 +-0.00963930481033161 +-0.00947014072771391 +-0.00929691872246433 +-0.00911986989764374 +-0.00893922907770915 +-0.00875523437457423 +-0.00856812674953827 +-0.00837814957204084 +-0.00818554817621721 +-0.00799056941621786 +-0.00779346122126418 +-0.00759447215140356 +-0.00739385095492654 +-0.00719184612839416 +-0.00698870548022275 +-0.00678467569874974 +-0.00658000192570287 +-0.00637492733596383 +-0.00616969272451323 +-0.00596453610140895 +-0.00575969229564071 +-0.00555539256866420 +-0.00535186423840642 +-0.00514933031449211 +-0.00494800914542088 +-0.00474811407838460 +-0.00454985313238936 +-0.00435342868530286 +-0.00415903717541873 +-0.00396686881808716 +-0.00377710733792560 +-0.00358992971707959 +-0.00340550595996820 +-0.00322399887490026 +-0.00304556387291290 +-0.00287034878413194 +-0.00269849369191971 +-0.00253013078502138 +-0.00236538422788722 +-0.00220437004929397 +-0.00204719604935186 +-0.00189396172493222 +-0.00174475821351342 +-0.00159966825539141 +-0.00145876617416354 +-0.00132211787534665 +-0.00118978086295128 +-0.00106180427378789 +-0.000938228929244725 +-0.000819087404231046 +-0.000704404112946686 +-0.000594195411093534 +-0.000488469714115528 +-0.000387227631009619 +-0.000290462113224655 +-0.000198158618124377 +-0.000110295286466196 +-2.68431333115253e-05 +5.22337482373576e-05 +0.000126977971112919 +0.000197438626827597 +0.000263671047171732 +0.000325736555765667 +0.000383702210190467 +0.000437640535430580 +0.000487629249386302 +0.000533750981217625 +0.000576092983298052 +0.000614746837560138 +0.000649808157027262 +0.000681376283322907 +0.000709553980957602 +0.000734447129189208 +0.000756164412253804 +0.000774817008755975 +0.000790518281006618 +0.000803383465083609 +0.000813529362383507 +0.000821074033418418 +0.000826136494601106 +0.000828836418741982 +0.000829293839967222 +0.000827628863746419 +0.000823961382698209 +0.000818410798817751 +0.000811095752749854 +0.000802133860700985 +0.000791641459561057 +0.000779733360774545 +0.000766522613473622 +0.000752120277352463 +0.000736635205733962 +0.000720173839244429 +0.000702840010480627 +0.000684734760017794 +0.000665956164074958 +0.000646599174114904 +0.000626755468624841 +0.000606513317283699 +0.000585957457688754 +0.000565168984775198 +0.000544225253028676 +0.000523199791551049 +0.000502162232006304 +0.000481178249435283 +0.000460309515894002 +0.000439613666832982 +0.000419144280102709 +0.000398950867434119 +0.000379078878211122 +0.000359569715318499 +0.000340460762818218 +0.000321785425174916 +0.000303573177723108 +0.000285849628039315 +0.000268636587855687 +0.000251952155124712 +0.000235810805821544 +0.000220223495045263 +0.000205197766960386 +0.000190737873098248 +0.000176844898520315 +0.000163516895326989 +0.000150749022981455 +0.000138533694902904 +0.000126860730772267 +0.000115717513982355 +0.000105089153656573 +9.49586506519924e-05 +8.53070669589730e-05 +7.61136979045506e-05 +6.73562465662804e-05 +5.90109998022162e-05 +5.10530053056054e-05 +4.34562490951040e-05 +3.61938328576359e-05 +2.92381505669683e-05 +2.25610638103370e-05 +1.61340752645808e-05 +9.92849977621546e-06 +3.91563251130230e-06 +-1.93308634239090e-06 +-7.64591082430952e-06 +-1.32506328857496e-05 +-1.87744309766349e-05 diff --git a/proc/lowPass.fir b/proc/lowPass.fir new file mode 100644 index 0000000..80ca42c --- /dev/null +++ b/proc/lowPass.fir @@ -0,0 +1,69 @@ +-0.000440335526569633 +0.000752073665275446 +-0.000831122309556253 +0.000588957101340026 +-3.46436249221642e-18 +-0.000827459944534284 +0.00161026482296418 +-0.00193919666600433 +0.00143996079500368 +-7.17747082065639e-18 +-0.00204911840998437 +0.00392118417011142 +-0.00461031952047503 +0.00333139521022357 +-1.29174204891536e-17 +-0.00448625351020396 +0.00837009222232793 +-0.00961845095399045 +0.00681249875080015 +-1.94806966356424e-17 +-0.00890110319994063 +0.0164485237041207 +-0.0188000252487688 +0.0133068390500081 +-2.54911547218713e-17 +-0.0176724993310351 +0.0333062236509138 +-0.0392384317790799 +0.0290414086223355 +-2.96885614614232e-17 +-0.0453323576128337 +0.0991619338637596 +-0.150215709264938 +0.186775372818389 +0.800191309660683 +0.186775372818389 +-0.150215709264938 +0.0991619338637596 +-0.0453323576128337 +-2.96885614614232e-17 +0.0290414086223355 +-0.0392384317790799 +0.0333062236509138 +-0.0176724993310351 +-2.54911547218713e-17 +0.0133068390500081 +-0.0188000252487688 +0.0164485237041207 +-0.00890110319994063 +-1.94806966356424e-17 +0.00681249875080015 +-0.00961845095399045 +0.00837009222232793 +-0.00448625351020396 +-1.29174204891536e-17 +0.00333139521022357 +-0.00461031952047503 +0.00392118417011142 +-0.00204911840998437 +-7.17747082065639e-18 +0.00143996079500368 +-0.00193919666600433 +0.00161026482296418 +-0.000827459944534284 +-3.46436249221642e-18 +0.000588957101340026 +-0.000831122309556253 +0.000752073665275446 +-0.000440335526569633