Simulación De Un Biopotencial

clc
clear all
n=[1 15 30 55 75 100 115 135 165 170];
s1=[0 0 1.5 -0.28 -0.2 -0.13 -0.07 0 0 0];
d=(170/165);
p=[1:d:170];
t=[1/165:1/165:1];
y1=pchip(n,s1,p);
c=y1+0.28;o=c*143;
e=round (o);
e
plot(t,e)
grid on

S1 son los puntos de la imagen, De ahí interpolo los tres vectores con pchip.
C=y1 +0.28 porque mi punto más bajo era -0.28, despuésde sumar ese vector más 0.28 mi punto más bajo está en cero, después sólo multiplico para que mi puntos más alto se acerque lo más posible a 255 y lo redondeo, así ya tendrévalores enteros y positivos de 0 a 255 para ocupar 8 bits.
E són los puntos que meteras a proton y van a aparecerte en el workspace.

Device = 16F877A
XTAL = 20

Dim salidas AsPORTB
Dim BIO1[64] As Byte
Dim BIO2[64] As Byte
Dim x As Byte
ALL_DIGITAL=TRUE
Output salidas
StrBIO1=40,40,40,40,40,40,40,40,40,40,40,40,40,40,41,46,56,70,87,107,128,150,172,193,212,229,242,251,254,254,250,245,237,228,217,205,192,178,163,147,132,116,100,85,71,57,44,32,22,14,7,2,0,0,0,0,1,1,2,2,3,4,4,5;
StrBIO2=6,7,8,8,9,10,10,11,12,12,12,13,13,14,14,15,15,15,16,16,16,17,17,18,18,18,19,19,20,20,20,21,21,22,22,23,24,24,25,25,26,27,27,28,28,29,30,30,31,31,32,33,33,34,35,35,36,37,37,38,38,39,39,40;
While 1=1For x=0 To 63
salidas=BIO1[x]
DelayUS 6
Next x
For x=0 To 63
salidas=BIO2[x]
DelayUS 6
Next x
For x=0 To 35
salidas= 40
DelayMS 4
Next x
Wend
End

Lo del protoncreo que está muy sencillo, lo único que más o menos te explico es que el delay va de acuerdo con tus puntos yo tengo 163 puntos y mi señal debe durar un milisegundo, entonces paraobtener mi periodo divido 0.001s/163 puntos =6microsegundos por eso ese delay y los 40 que están en el último for son como 0 para dejar que haya un espacio entre las dos señales