C =======================================================================
C PROSEN (HACE LOS PROCESADOS DE LA SEQAL: FILTRADO PASO ALTO Y BAJO
C DERIVADAS,VALOR CUADRATICO Y INTEGRACION CON VENTANA MOVIL)
C AUTOR: PABLO LAGUNA/RAIMON JANE,
C DATA: 22-MAIG-87
C MODIFICAT: 19-FEBRER-88
C =======================================================================
C -----------------------------------------------------------------------
C Copyright (C) 2002 Pablo Laguna
C
C This program is free software; you can redistribute it and/or modify it
C under the terms of the GNU General Public License as published by the
C Free Software Foundation; either version 2 of the License, or (at your
C option) any later version.
C
C This program is distributed in the hope that it will be useful, but
C WITHOUT ANY WARRANTY; without even the implied warranty of
C MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
C General Public License for more details.
C
C You should have received a copy of the GNU General Public License along
C with this program; if not, write to the Free Software Foundation, Inc.,
C 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
C
C You may contact the author by e-mail (laguna@posta.unizar.es) or postal
C mail (Dpt. Ingenieria Electrónica y Comunicaciones, Grupo de Tecnología
C de las Comunicaciones (GTC), Centro Politécnico Superior. Universidad
C de Zaragoza, María de Luna, 3 (Pol. Actur), Edificio A, Despacho 3.16,
C 50015 Zaragoza. Spain). For updates to this software, please visit
C PhysioNet (http://www.physionet.org/).
C _______________________________________________________________________
c
c SE LE LLAMA DESDE ECGMAIN PARA REPRESENTAR POR PANTALLA,IMPRESORA
c Y PARA DETECTAR QRS (SU PROCESAMIENTO PREVIO)
c
C
C ---------------------------------------------------------------
C REALIZA UN FILTRADO PASO BAJO DE LA SEQAL
C ---------------------------------------------------------------
subroutine fpb(n,if,ifc,x,y,iretard)
dimension x(100000),y(100000)
c APLICA A LA SEQAL QUE ENTRA EN EL VECTOR x UN FILTRO PASO BAJO. lA
C SEQAL QUE ENTRA LA SUPONGO MUESTREADA A if hZ.
C EL FILTRO TIENE UN POLO DOBLE EN (0,1) Y ifc CEROS DOBLES EN
C ELA CIRCULO UNIDAD .LA GANANCIA DEL
C FILTRO SERIA DEl NUMERO de CEROS AL CUADRADO.
C EL RETARDO DEL FILTRO ES DE NUMERO de CEROS MENOS UNO
C LAS CONDICIONES INICIALES SUPONGO QUE LA SEQAL ES CONSTANTE ANTES DEL
C INSTANTE INICIAL CON VALOR 0
C LA ECUACION SERIA y(Nt)=2y(Nt-t)-y(Nt-2t)+x(Nt)-2x(Nt-CEROSt)
C +x(Nt-2*CEROSt).
C LA SEQAL FILTRADA LA DEVUELVE EN EL VECTOR y
C RM ES EL VALOR MEDIO DE LOS 10 PRIMEROS PUNTOS
c sum=x(1)
c do i=2,10
c sum=sum+x(i)
c end do
c rm=sum/10
c rm=x(1)
rm=0.
C l eS EL NUMERO DE CEROS
c ifc ES LA FRECUENCIA DE CORTE A CERO DEL FILTRO
c if ES LA FRECUENCIA DE MUESTREO DE LA SEQAL
c l=nint(1.*if/ifc)
c y(1)=x(1)+ l*rm
c do i=2,l
c y(i)=y(i-1)+x(i)-rm
c end do
c do i=l+1,n
c y(i)=y(i-1)+x(i)-x(i-l)
c end do
c do i=n,n+l
c y(i)=y(i-1)+x(n)-x(i-l)
c end do
c filtro de segundo orden (cuidado se puede inestabilizar
l=nint(1.*if/ifc)
y(1)=x(1)+(l**2-1)*rm
y(2)=(2*y(1)+x(2)-(l**2+1)*rm)
do i=3,l
y(i)=(2*y(i-1)-y(i-2)+x(i)-rm)
end do
do i=l+1,2*l
y(i)=(2*y(i-1)-y(i-2)+x(i)-2*x(i-l)+rm)
end do
do i=2*l+1,n
y(i)=(2*y(i-1)-y(i-2)+x(i)-2*x(i-l)+x(i-2*l))
end do
do i=n,n+l
y(i)=(2*y(i-1)-y(i-2)+x(n)-2*x(i-l)+x(i-2*l))
end do
C QUITO EL RETARDO DE LA SEQAL
iretard=l-1
if (n.le.(100000-l+1) ) then
do i=1,n
y(i)=y(i+l)
end do
else
do i=1,100000-l+1
y(i)=y(i+l)
end do
do i=100000-l+1,n
y(i)=0
end do
end if
C PONGO EL RESTO A CERO
do i=n+1,100000
y(i)=0
end do
return
end
C --------------------------------------------------------------
c REALIZA UN FILTRADO PASO BAJO DE UN TROZO DE SEQAL
C --------------------------------------------------------------
subroutine fpbc(ni,n,if,ifc,x,y,iretard)
dimension x(100000),y(400)
c APLICA A LA SEQAL QUE ENTRA EN EL VECTOR x UN FILTRO PASO BAJO. lA
C SEQAL QUE ENTRA LA SUPONGO MUESTREADA A if hZ.
C EL FILTRO TIENE UN POLO DOBLE EN (0,1) Y ifc CEROS DOBLES EN
C EL CIRCULO UNIDAD .LA GANANCIA DEL
C FILTRO SERIA DEl NUMERO de CEROS AL CUADRADO.
C EL RETARDO DEL FILTRO ES DE NUMERO de CEROS MENOS UNO
C LAS CONDICIONES INICIALES SUPONGO QUE LA SEQAL ES CONSTANTE ANTES DEL
C INSTANTE INICIAL CON VALOR 0
C lA ECUACION SERIA Y(nT)=2Y(nT-T)-Y(nT-2T)+X(nT)-2X(nT-cerosT)
c +X(nT-2*cerosT).
C lA SEQAL FILTRADA LA DEVUELVE EN EL VECTOR y
C RM ES EL VALOR MEDIO DE LOS 10 PRIMEROS PUNTOS
c sum=x(ni)
c do i=ni+2,ni+10
c sum=sum+x(i)
c end do
rm=x(ni)
c l ES EL NUMERO DE CEROS
C ifc ES LA FRECUENCIA A LA QUE EL FILTRO CAE A CERO
C if ES LA FRECUENCIA A LA QUE ESTA MUESTREADA LA SEQAL
C ni ES LA POSICION INICIAL A PARTIR DE LA CUAL QUIERO FILTRAR
c n ES EL NUMERO DE PUNTOS A FILTRAR
C SE DEBE CUMPLIR (n<1000-l) SI NO SE SALDRA DE LOS LIMITES DEL VECTOR
l=nint(1.*if/ifc)
y(1)=x(ni+1)+(l**2-1)*rm
y(2)=(2*y(1)+x(ni+2)-(l**2+1)*rm)
do i=3,l
y(i)=(2*y(i-1)-y(i-2)+x(ni+i)-rm)
end do
do i=l+1,2*l
y(i)=(2*y(i-1)-y(i-2)+x(ni+i)-2*x(ni+i-l)+rm)
end do
do i=2*l+1,n+l
y(i)=(2*y(i-1)-y(i-2)+x(ni+i)-2*x(ni+i-l)+x(ni+i-2*l))
end do
C QUITO EL RETARDO DE LA SEQAL
iretard=l
if (n.le.(400-l+1)) then
do i=1,n
y(i)=y(i+l)
end do
else
do i=1,400-l+1
y(i)=y(i+l)
end do
do i=400-l+1,n
y(i)=0
end do
end if
C PONGO EL RESTO A CERO
do i=n+1,400
y(i)=0
end do
return
end
C ---------------------------------------------------------------
C FILTRE PASSA ALT (primer ordre)
C ---------------------------------------------------------------
subroutine fpa(n,if,ifc,x,y,iretard)
dimension x(100000), y(100000)
C
c APLICA A LA SEQAL QUE ENTRA EN EL VECTOR x UN FILTRO PASO BAJO. lA
C SEQAL QUE ENTRA LA SUPONGO MUESTREADA A if hZ.
C EL FILTRO TIENE UN POLO DOBLE EN (0,1) Y ifc CEROS DOBLES EN
C ELA CIRCULO UNIDAD .LA GANANCIA DEL
C FILTRO SERIA DEl NUMERO de CEROS AL CUADRADO.
C EL RETARDO DEL FILTRO ES DE NUMERO de CEROS MENOS UNO
C LAS CONDICIONES INICIALES SUPONGO QUE LA SEQAL ES CONSTANTE ANTES DEL
C INSTANTE INICIAL CON VALOR 0
C LA ECUACION SERIA y(Nt)=2y(Nt-t)-y(Nt-2t)+x(Nt)-2x(Nt-CEROSt)
C +x(Nt-2*CEROSt).
C LA SEQAL FILTRADA LA DEVUELVE EN EL VECTOR y
C RM ES EL VALOR MEDIO DE LOS 10 PRIMEROS PUNTOS
c sum=x(1)
c do i=2,10
c sum=sum+x(i)
c end do
c rm=sum/10
c rm=x(1)
rm=0.
C l eS EL NUMERO DE CEROS
c ifc ES LA FRECUENCIA DE CORTE A CERO DEL FILTRO
c if ES LA FRECUENCIA DE MUESTREO DE LA SEQAL
c l=nint(1.*if/ifc)
c y(1)=x(1)+ l*rm
c do i=2,l
c y(i)=y(i-1)+x(i)-rm
c end do
c do i=l+1,n
c y(i)=y(i-1)+x(i)-x(i-l)
c end do
c do i=n,n+l
c y(i)=y(i-1)+x(n)-x(i-l)
c end do
c filtro de segundo orden (cuidado se puede desestabilizar)
l=nint(1.*if/ifc)
y(1)=x(1)+(l**2-1)*rm
y(2)=(2*y(1)+x(2)-(l**2+1)*rm)
do i=3,l
y(i)=(2*y(i-1)-y(i-2)+x(i)-rm)
end do
do i=l+1,2*l
y(i)=(2*y(i-1)-y(i-2)+x(i)-2*x(i-l)+rm)
end do
do i=2*l+1,n
y(i)=(2*y(i-1)-y(i-2)+x(i)-2*x(i-l)+x(i-2*l))
end do
do i=n,n+l
y(i)=(2*y(i-1)-y(i-2)+x(n)-2*x(i-l)+x(i-2*l))
end do
C PARTE PASO TODO MENOS PASO BAJO
iretard=l-1
gain=l*l
do i=1,n-iretard
y(i)=gain*x(i)-y(i+iretard)
end do
do i=n-iretard+1,n
y(i)=0.
end do
return
end
C ---------------------------------------------------------------
C REALIZA UN FILTRADO PASO ALTO DE LA SEQAL
C ---------------------------------------------------------------
C--- subroutine fpa(n,if,ifc,x,y,iretard)
c--- dimension x(100000),y(100000)
C
c aPLICA A LA SEQAL QUE ENTRA EN EL VECTOR x UN FILTRO PASO ALTO. lA
C SEQAL QUE ENTRA LA SUPONGO MUESTREADA A if hZ . eL FILTRO LO
c IMPLEMENTO RESTANDO DE UNO PASO TODO OTRO PASO BAJO DE FRECUENCIA
c DE CORTE MENOR QUE LA ANTERIOR .
C cON ESTO OBTENGO UNA FRECUENCIA DE CORTE DE FC HZ, Y EL
C lAS CONDICIONES INICIALES SUPONGOQUE LA SEQAL ES CONSTANTE ANTES DEL
C INSTANTE INICIAL CON VALOR 0.
C lA ECUACION SERIA y(Nt)=ceros*x(Nt-ceros/2*t)-y(Nt-t)-x(Nt)
c +x(Nt-ceros*t).
C lA SEQAL FILTRADA LA DEVUELVE EN EL VECTOR y
C
C RM ES EL VALOR DEL PRIMER PUNTO
C rm=x(1)
c-- rm=0.
C l eS EL NUMERO DE CEROS
c ifc ES LA FRECUENCIA DE CORTE A CERO DEL FILTRO
c if ES LA FRECUENCIA DE MUESTREO DE LA SEQAL
c-- l=nint(1.*if/ifc)
c-- y(1)=x(1)+(l**2-1)*rm
c-- y(2)=(2*y(1)+x(2)-(l**2+1)*rm)
c-- do i=3,l
c-- y(i)=(2*y(i-1)-y(i-2)+x(i)-rm)
c-- end do
c-- do i=l+1,2*l
c-- y(i)=(2*y(i-1)-y(i-2)+x(i)-2*x(i-l)+rm)
c-- end do
c-- do i=2*l+1,n
c-- y(i)=(2*y(i-1)-y(i-2)+x(i)-2*x(i-l)+x(i-2*l))
c-- end do
C do i=n,n+l
C y(i)=(2*y(i-1)-y(i-2)+x(n)-2*x(i-l)+x(i-2*l))
C end do
C PARTE PASO TODO MENOS PASO BAJO
c do i=1,(l-1)
c y(i)=l**2*rm-y(i)
c end do
c do i=l,n+l
c y(i)=l**2*x(i-(l-1))-y(i)
c end do
c-- do i=1,l-1
c-- y(i)=l**2*rm-y(i)
c-- end do
c-- do i=l,n
c-- y(i)=l**2*x(i-(l-1))-y(i)
c-- end do
C QUITO EL RETARDO DE LA SEQAL
c
c iretard=l-1
c if (n.le.(100000-l+1)) then
c do i=1,n
c y(i)=y(i+l-1)
c end do
c else
c do i=1,100000-l+1
c y(i)=y(i+l-1)
c end do
c do i=100000-l+1,n
c y(i)=0
c end do
c end if
c- if (n.le.(100000-l+1)) then
c- do i=1,n-L+1
c- y(i)=y(i+l-1)
c- end do
c- else
c- do i=1,100000-l+1
c- y(i)=y(i+l-1)
c- end do
c- do i=100000-l,n
c- y(i)=0
c- end do
c- end if
C PONGO EL RESTO A CERO
c- do i=n-L,100000
c- y(i)=0
c- end do
c-
c--- return
c--- end
C --------------------------------------------------------------
c REALIZA UN FILTRADO PASO AlTO DE UN TROZO DE SEQAL
C --------------------------------------------------------------
subroutine fpac(ni,n,if,ifc,x,y,iretard)
dimension x(100000),y(2000)
c APLICA A LA SEQAL QUE ENTRA EN EL VECTOR x UN FILTRO PASO ALTO. lA
C SEQAL QUE ENTRA LA SUPONGO MUESTREADA A if hZ.
C EL FILTRO TIENE UN POLO DOBLE EN (0,1) Y ifc CEROS DOBLES EN
C EL CIRCULO UNIDAD .LA GANANCIA DEL
C FILTRO SERIA DEl NUMERO de CEROS AL CUADRADO.
C EL RETARDO DEL FILTRO ES DE NUMERO de CEROS MENOS UNO
C LAS CONDICIONES INICIALES SUPONGO QUE LA SEQAL ES CONSTANTE ANTES DEL
C INSTANTE INICIAL CON VALOR X(NI)
C lA ECUACION SERIA Y(nT)=2Y(nT-T)-Y(nT-2T)+X(nT)-2X(nT-cerosT)
c +X(nT-2*cerosT).
C lA SEQAL FILTRADA LA DEVUELVE EN EL VECTOR y
C RM ES EL VALOR MEDIO DE LOS 10 PRIMEROS PUNTOS
c sum=x(ni)
c do i=ni+2,ni+10
c sum=sum+x(i)
c end do
rm=x(ni)
c l ES EL NUMERO DE CEROS
C ifc ES LA FRECUENCIA A LA QUE EL FILTRO CAE A CERO
C if ES LA FRECUENCIA A LA QUE ESTA MUESTREADA LA SEQAL
C ni ES LA POSICION INICIAL A PARTIR DE LA CUAL QUIERO FILTRAR
c n ES EL NUMERO DE PUNTOS A FILTRAR
C SE DEBE CUMPLIR (n<2000-l) SI NO SE SALDRA DE LOS LIMITES DEL VECTOR
l=nint(1.*if/ifc)
y(1)=x(ni+1)+(l**2-1)*rm
y(2)=(2*y(1)+x(ni+2)-(l**2+1)*rm)
do i=3,l
y(i)=(2*y(i-1)-y(i-2)+x(ni+i)-rm)
end do
do i=l+1,2*l
y(i)=(2*y(i-1)-y(i-2)+x(ni+i)-2*x(ni+i-l)+rm)
end do
do i=2*l+1,n+l
y(i)=(2*y(i-1)-y(i-2)+x(ni+i)-2*x(ni+i-l)+x(ni+i-2*l))
end do
C QUITO EL RETARDO DE LA SEQAL
iretard=l
if (n.le.(2000-l+1)) then
do i=1,n
y(i)=y(i+l)
end do
else
do i=1,2000-l+1
y(i)=y(i+l)
end do
do i=2000-l+1,n
y(i)=0
end do
end if
C PARTE PASO TODO MENOS PASO BAJO
do i=1,n
y(i)=l**2*x(ni+i)-y(i)
end do
C PONGO EL RESTO A CERO
do i=n+1,2000
y(i)=0
end do
return
end
C ---------------------------------------------------------------
C REALIZA UN TEXT SOBRE LA AMPLITUD POSIT. Y NEGAT. DE LA SEQAL
C ---------------------------------------------------------------
subroutine test(if,x,inv)
dimension x(100000)
logical inv
C inv dICE SI HE INVERTIDO LA SEQAL O NO
senp=0
senn=0
c bUSCO EN LOS PRIMEROS SEGUNDOS
do i=if,if*3
if (x(i).lt.senn) senn=x(i)
if (x(i).gt.senp) senp=x(i)
end do
inv=.false.
if (senp.lt.(-senn)) then
inv=.true.
end if
return
end
C ---------------------------------------------------------------
C INVIERTE LA SEQAL
C ---------------------------------------------------------------
subroutine invertir(sen,n)
dimension sen(100000)
do i=1,n
sen(i)=-sen(i)
end do
return
end
C ---------------------------------------------------------------
C REALIZA UN FILTRADO DERIVADA DE LA SEQAL
C ---------------------------------------------------------------
subroutine der(n,if,x,y)
dimension x(100000),y(100000)
c aPLICA A LA SEQAL QUE ENTRA EN EL VECTOR x UN FILTRO DERIVADA. lA
C SEQAL QUE ENTRA LA SUPONGO MUESTREADA A if hZ.
c ESTA DERIVADA SERA SERA TAL HASTA UNOS 30 hZ. DESPUES CAE LO CUAL
C ES BUENO PARA NO ENFATIZAR FRECUNCIAS ALTAS
C eL RETARDO DE LA DERIVADA ES DE l MUESTRAS.
C lAS CONDICIONES INICIALES SUPONGOQUE LA SEQAL ES CONSTANTE ANTES DEL
C INSTANTE INICIAL CON VALOR 0.
C lA ECUACION SERIA y(Nt)=25*(x(Nt)+2x(Nt-l/2*t)-x(Nt-(l/2+l)t)
c -x(Nt-(2*l)t)).
C lA SEQAL DERIVADA LA DEVUELVE EN EL VECTOR y
k=if/8
y(3)=if*(x(2)-x(1))
y(4)=(if/4)*(2*x(3)-2*x(1))
do i=5,n+2
y(i)=k*(x(i)+2*x(i-1)-2*x(i-3)-x(i-4))
end do
C QUITO EL RETARDO DE LA SEQAL
do i=1,n
y(i)=y(i+2)
end do
C PONGO EL RESTO A CERO
do i=n+1,100000
y(i)=0
end do
return
end
C ---------------------------------------------------------------
C REALIZA LA FUNCION CUADRADA DE LA SEQAL
C ---------------------------------------------------------------
subroutine quad(n,x,y)
dimension x(100000),y(100000)
C aPLICA A LA SEQAL QUE ENTRA EN EL VECTOR x EL OPERADOR ELEVER AL
C CUADRADO Y LA SACA EN EL VESTOR y
do i=1,n
y(i)=x(i)**2
end do
c do i=n+1,100000
c y(i)=0
c end do
return
end
C ---------------------------------------------------------------
C REALIZA LA INTEGRACION DE VENTANA MOVIL
C ---------------------------------------------------------------
subroutine mwint(n,if,x,y)
dimension x(100000),y(100000)
c aPLICA A LA SEQAL QUE ENTRA EN EL VECTOR x UNA INTEGRACION DE VENTANA MOVI
C MOVIL DE UNA ANCHURA DE 70 MUESTRAS QUE CORRESPONDE A 70 MS., MAS O
C MENOS LA ANCHURA DE UN qrs.
C SEQAL QUE ENTRA LA SUPONGO MUESTREADA A if
C lAS CONDICIONES INICIALES SUPONGOQUE LA SEQAL ES CONSTANTE ANTES DEL
C INSTANTE INICIAL CON VALOR x(1).
C lA SEQAL INTEGRADA LA DEVUELVE EN EL VECTOR y
l=nint(9.5*if/100.)
do i=1,n
if (i.lt.l) then
C La seguent instruccio tambe hi havia un sum=0 ???
C sum= (x(1)*(l-i))/l
sum=0
do j=1,i
sum=sum+x(j)/l
end do
end if
if (i.ge.l) then
sum=0
do j=i-l,i
sum=sum+x(j)/l
end do
end if
y(i)=sum
end do
c do i=n+1,100000
c y(i)=0
c end do
return
end
C ---------------------------------------------------------------
C REALIZA LA INTEGRACION DE VENTANA MOVIL
C ---------------------------------------------------------------
subroutine mwint_p(n,if,x,y,z,iqrs)
dimension x(100000),y(100000),z(100000)
c aPLICA A LA SEQAL QUE ENTRA EN EL VECTOR x UNA INTEGRACION DE VENTANA MOVI
C MOVIL DE UNA ANCHURA DE 70 MUESTRAS QUE CORRESPONDE A 70 MS., MAS O
C MENOS LA ANCHURA DE UN qrs.
C SEQAL QUE ENTRA LA SUPONGO MUESTREADA A if
C lAS CONDICIONES INICIALES SUPONGOQUE LA SEQAL ES CONSTANTE ANTES DEL
C INSTANTE INICIAL CON VALOR x(1).
C lA SEQAL INTEGRADA LA DEVUELVE EN EL VECTOR y
l=nint(30.5*if/100.)
l2=nint(30.5*if/100./2)
do i=1,n
if (i.lt.l2) then
C La seguent instruccio tambe hi havia un sum=0 ???
C sum= (x(1)*(l-i))/l
sum=0
do j=1,i+l2
sum=sum+x(j)*z(iqrs-i+j) /l
end do
end if
if ((i.ge.l2).and.(i.lt.(n-l2))) then
sum=0
do j=i-l2,i+l2
sum=sum+x(j)*z(iqrs-i+j)/l
end do
end if
if (i.ge.(n-l2)) then
sum=0
do j=i-l2,n
sum=sum+x(j)*z(iqrs-i+j)/l
end do
end if
y(i)=sum
end do
c do i=n+1,100000
c y(i)=0
c end do
return
end
C ---------------------------------------------------------------
C REALIZA UNA NORMALIZACION DE MODULO DE LA SEQAL
C ---------------------------------------------------------------
subroutine normaliz(if,max,y)
dimension y(100000)
C NORMALIZO AL MAXIMO VALOR max
rmax=y(1)
do i=2,2*if
if (ABS(y(i)).gt.rmax) then
rmax=ABS(y(i))
end if
end do
if (rmax.eq.0) rmax=1
C NORMALIZO
do i=1,100000
y(i)=max*y(i)/rmax
end do
return
end
C ---------------------------------------------------------------
C REALIZA LA NORMALIZACION INICIAL
C ---------------------------------------------------------------
subroutine normaliz_i(if,max,rmax,y)
dimension y(100000)
C rmax es el factor de normalitzacio (parametre de sortida)
C NORMALIZO AL MAXIMO VALOR max
rmax=y(1)
do i=2,2*if
if (ABS(y(i)).gt.rmax) then
rmax=ABS(y(i))
end if
end do
if (rmax.eq.0) rmax=1
C NORMALIZO
do i=1,100000
y(i)=max*y(i)/rmax
end do
return
end
C ---------------------------------------------------------------
C REALIZA LA NORMALIZACION "ciclica" con el factor RMAX
C ---------------------------------------------------------------
subroutine normaliz_c(if,max,rmax,y)
dimension y(100000)
C NORMALIZO
do i=1,100000
y(i)=max*y(i)/rmax
end do
return
end
C ---------------------------------------------------------------
C REALIZA UN CAMBIO EN LA FRECUENCIA DE MUESTREO DE LA SEQAL
C ---------------------------------------------------------------
subroutine nuevafrecuencia(if,sen,senam)
dimension sen(100000),senam(100000)
C IFN ES LA NUEVA FRECUENCIA A LA QUE SALE LA SEQAL MUESTREADA
C IFV ES LA FRECUENCIA A LA QUE ESTA MUESTREADA INICIALMENTE
ifv=if
5 write(6,10) ifv
10 format(1x,/,'$ nueva frecuencia_muestreo ? (<',i5,') if= ')
read(5,*,err=5) ifn
icon=ifv/ifn
if=ifv/icon
j=1
do i=1,100000,icon
sen(j)=sen(i)
senam(j)=senam(i)
j=j+1
end do
do i=j,100000
sen(i)=0
senam(j)=0
end do
return
end