ejemplo de un reporte de practicas

Ingeniera Electrnica Laboratorio Electrnica DigitalMateria P.D.S Practica No Filtro Fir con Arduino Due Ttulo Procesamiento digital de seales Integrantes Lee Fonseca Francisco Rodrguez Caedo Pedro Delgado Caedo Jos Manuel Castro Viedas Juan Carlos Iribe Soto Vctor Alonso Flores Prez Hctor Alejandro Martnez Zavala Jess Alberto Duracin (hrs.) 2 hrs. Profesor M.C. Heraclio HerediaUreta.Materiales Arduino Due Amplificador Operacional LM324 Equipo 3 Generadores de seal 1 Osciloscopio 2 Fuentes triplesSeguridad e Higiene Seguir las medidas de seguridad e higiene indicadas en el reglamento interno del Laboratorio.Introduccin El modelo matemtico de los filtros FIR tambin se fundamenta la ecuacin de diferencia, pero con la particularidad de que todos los coeficientes ak son iguales a cero.Se tiene entonces que la ecuacin que los describe es funcin del conjunto de coeficientes bk Y de la secuencia de entrada x(n). EMBED PBrush Donde M1 corresponde a la longitud del filtro. Este sistema considera slo las ultimas M1 muestras de la seal de entrada y las pondera mediante los coeficientes bk. A este sistema se le denomina FIR, ya que su respuesta al impulso unitario (dada por loscoeficientes bk) es finita. Su diseo requiere la seleccin de la secuencia que mejor representa la respuesta al impulso de un filtro ideal. Los filtros FIR son siempre estables y son capaces de tener una respuesta lineal en fase. Frente a los filtros IIR presentan la desventaja de requerir un orden mucho mayor. Cdigos de MatLab para obtener los coeficientes y graficas.- Para el filtro pasa bajas.-clc close all passband_cutoff1100 stopband_cutoff1500 rbp0.01 rsb0.01 frecuencia_de_muestreo8000 n, Wn, beta, typkaiserord( passband_cutoff stopband_cutoff, 1, 0, rbp rsb, frecuencia_de_muestreo ) Bfir1(n, Wn, typ, kaiser(n1, beta), noscale) Fx00.00050.5 f1100 f21000 f32000 n1011000 fm2f3 T1/fm figure(1) Hfreqz(B,1,2pi.Fx) plot(Fxfm,abs(H)) X 2 cos((2pi(f1/fm)).n1) 0.5.cos((2pi(f2/fm)).n1)0.5.cos((2pi(f3/fm)).n1) figure(2) plot(n1,X) yfilter(B,1,X) figure(3) plot(n1,y) Para el fltro pasa bandas.- clc close all fcuts 200 900 1100 1900 devs 0.01 0.01 0.01 frecuencia_de_muestreo8000 n, Wn, beta, typkaiserord(fcuts, 0 1 0, devs, frecuencia_de_muestreo) Bfir1(n, Wn, typ, kaiser(n1, beta), noscale) Fx00.00050.5 f1100 f21000 f32000 n1011000 fm2f3 T1/fm figure(1) Hfreqz(B,1,2pi.Fx)plot(Fxfm,abs(H)) X 2 cos((2pi(f1/fm)).n1) 0.5.cos((2pi(f2/fm)).n1) 0.5.cos((2pi(f3/fm)).n1) figure(2) plot(n1,X) yfilter(B,1,X) figure(3) plot(n1,y) Para el filtro pasa altas.- clc close all passband_cutoff300 stopband_cutoff600 rbp0.01 rsb0.01 frecuencia_de_muestreo8000 n, Wn, beta, typkaiserord( passband_cutoff stopband_cutoff, 0, 1, rbp rsb, frecuencia_de_muestreo) Bfir1(n, Wn, typ,kaiser(n1, beta), noscale) Fx00.00050.5 f1100 f21000 f32000 n01200 fmf32 T1/fm Hfreqz(B,1,2pi.Fx) figure(1) plot(Fxfm,abs(H)) x 2 cos((2pi(f1/fm)).n) 0.5.cos((2pi(f2/fm)).n) 0.5.cos((2pi(f3/fm)).n) y filter(B,1,x) figure(2) subplot(211), plot(n.T,x) subplot(212), plot(n.T,y) Codigo para el arduino due.- define FILTERTAPS 31 // El nmero de coeficientes // Se declara la formacin de valores floatvaluesFILTERTAPS 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 float input 1 float output 0 // Contador de nuestras aplicaciones de la vida real no necesitar esta byte n 0 int interrupcion 0 void setup() startTimer(TC1, 0, TC3_IRQn, 4000) //TC1 channel 0, the IRQ for that channel and the desired frequency // 20,000 seems to give 100uS/divsignal, which is equivalent to 10,000Hz or 10Khz // put your setup code here, to run once analogReadResolution(12) analogWriteResolution(12) void loop() if (interrupcion 1) // put your main code here, to run repeatedly analogWrite(DAC0,fir((float)analogRead(A0))) //analogWrite(DAC0,analogRead(A0)) interrupcion 0 float fir(float in) static byte k // k almacena...