Reconstrucción de Sonido con Láser y Fotosensor: De la Teoría al Código

2026-03-15T00:00:00+00:00

El Problema: ¿Puede un Rayo de Luz Transportar Sonido?</h2>

La idea surgió de una pregunta aparentemente simple en clase de física: si el sonido es una onda mecánica que vibra el aire, ¿podría esa vibración modular un haz de luz y recuperarse después?

La respuesta es sí, y el principio es elegante. El sonido dentro de una caja insonorizada hace vibrar un panel de vidrio. Un láser de 650 nm apunta al vidrio desde afuera; el reflejo modulado por las vibraciones es capturado por un fotodiodo BPW34. La corriente fotoeléctrica resultante pasa por un amplificador de transimpedancia (TIA LM358) que la convierte en voltaje, luego por un MAX9814 con control de ganancia automático (AGC), y finalmente llega al jack de micrófono de la PC donde MATLAB la digitaliza y reconstruye el audio original:

$$V(t) = A \cdot \sin(2\pi f t) + \eta(t)$$

donde $\eta(t)$ es el ruido introducido por el sistema óptico-electrónico. El reto es precisamente ese: dado que $\eta(t)$ siempre existe, ¿cómo recuperamos la señal $A \cdot \sin(2\pi f t)$ con la mayor fidelidad posible?

Este post documenta la fase de simulación completa en MATLAB que valida los algoritmos antes de construir el hardware real.

Fundamentos Teóricos</h2>
Transformada de Fourier (FFT)</h3>
El teorema de Fourier establece que cualquier señal periódica puede descomponerse como suma de senos y cosenos. La Transformada de Fourier Discreta (DFT) hace esto numéricamente:
$$X[k] = \sum_{n=0}^{N-1} x[n] \cdot e^{-j2\pi kn/N}$$
La FFT (Fast Fourier Transform) es un algoritmo eficiente para calcular la DFT en $O(N \log N)$ en lugar de $O(N^2)$. Para este proyecto, la FFT cumple dos funciones:

Identificar qué frecuencias componen la señal capturada</li>
Filtrar el ruido eliminando componentes espectrales débiles</li> </ol>
Descomposición en Valores Singulares (SVD)</h3>
Dado que el proyecto conecta con el tema de espacios vectoriales del semestre, también apliqué SVD para analizar la estructura de cada señal. Dada una matriz de trayectoria $M$ construida a partir de la señal:
$$M = U \cdot \Sigma \cdot V^T$$
Los valores singulares $\sigma_i$ en $\Sigma$ revelan cuántos “modos independientes” necesita la señal para representarse. Una señal simple (tono puro) tiene muy pocos valores singulares grandes; una señal compleja los tiene distribuidos.

Conexión con el curso: El SVD es esencialmente un cambio de base ortogonal. Los vectores columna de $U$ son la nueva base en el espacio de la señal, y los valores singulares son las “coordenadas” de importancia de cada vector base. Exactamente el mismo concepto de base ortonormal que vimos en clase. </blockquote>

El Hardware</h2>
El sistema físico que se simulará y construirá está formado por los siguientes componentes, elegidos por su precisión y bajo ruido:
Componente</th> Función</th> Parámetro clave</th></tr></thead>

Láser D18×65mm 650nm 5mW</td> Fuente de luz, punto continuo</td> Enfocable, montado en soporte Neiko ~30°</td></tr>
BPW34 (fotodiodo PIN)</td> Luz → corriente</td> R = 0.3 A/W @ 650 nm, en tubo negro</td></tr>
LM358 (TIA)</td> Corriente → voltaje</td> Rf = 10 kΩ, Cf = 100 pF</td></tr>
Capacitor C_ac 10 µF</td> Bloqueo de offset DC</td> fc = 1.6 Hz</td></tr>
MAX9814 HiLetgo</td> Amplificador AGC</td> Av = 40 dB (GAIN flotante)</td></tr>
Capacitor C_sal 47 µF</td> Bloqueo offset salida</td> fc = 0.34 Hz</td></tr>
Módulo Jack ShillehTek 3.5mm</td> Interfaz PC</td> Entrada micrófono, breakout TRRS</td></tr>
ESP32</td> Fuente de alimentación 5V</td> Pin VIN → rail protoboard</td></tr>
Soporte Neiko</td> Posicionar láser en ángulo</td> Ajuste ~30° sobre el vidrio</td></tr> </tbody></table>
La ventaja principal del BPW34 sobre un LDR es su tiempo de respuesta: 20 ns frente a ~1 ms, permitiendo capturar audio hasta 20 kHz sin ninguna degradación por el sensor.
Diagrama del circuito</h3>

Montaje experimental</h3>

La escena 3D muestra el montaje completo: caja insonorizada con vidrio vibrante, láser industrial D18×65 mm sobre soporte Neiko apuntando al vidrio en ~30°, BPW34 en tubo negro en el ángulo de reflexión especular, y el circuito completo (ESP32 + LM358 + MAX9814 + módulo jack) en la protoboard. La animación sincroniza la vibración del vidrio con las señales en tiempo real en cada etapa del pipeline.

Metodología de Simulación</h2>
Para validar los algoritmos sin hardware, generé cuatro tipos de señal en MATLAB y simulé el proceso completo con el pipeline óptico-electrónico real: generación → modulación óptica → BPW34 → TIA → MAX9814 → filtrado FFT → reconstrucción → validación.
A diferencia de simplemente sumar ruido gaussiano, el modelo simula la cadena física completa:
function v_out = pipeline_sensor(audio, P_laser, mod_depth, ... R_bpw34, Rf_TIA, Av_MAX, ruido) % 1. Modulacion optica: el vidrio vibrante modula el reflejo P_luz = P_laser * (1 + mod_depth * audio); % 2. BPW34: fotocorriente = responsividad x potencia + shot noise I_foto = R_bpw34 * P_luz + ruido * 1e-4 * randn(size(audio)); % 3. TIA LM358: I → V, eliminar offset DC (como el capacitor C_ac) V_tia = I_foto * Rf_TIA; V_tia = V_tia - mean(V_tia); % 4. MAX9814: amplificar y saturar (±1.65V) V_max = Av_MAX * V_tia; v_out = max(-1.65, min(1.65, V_max)); end </code></pre> El ruido del BPW34 (shot noise) es de apenas 0.8% de la señal — notablemente menor que el 5% típico de un LDR, lo que explica los resultados de SNR superiores. Las Cuatro Señales de Prueba</h3> Señal</th> Frecuencias</th> Por qué es interesante</th></tr></thead> Tono puro</td> 1000 Hz</td> Caso base, señal más simple posible</td></tr> Acorde Do Mayor</td> 262 + 330 + 392 Hz</td> Tres frecuencias simultáneas (Do, Mi, Sol)</td></tr> Barrido (chirp)</td> 100 → 1000 Hz</td> Frecuencia variable en el tiempo</td></tr> Señal compleja</td> 262·k Hz, k=1..5</td> Serie armónica, simula voz humana</td></tr> </tbody></table> El Algoritmo de Reconstrucción</h3> El núcleo del sistema es el filtrado espectral: se aplica una máscara binaria al espectro FFT que conserva solo las componentes que superan el 10% del pico máximo, luego se reconstruye la señal con la IFFT. La clave es mantener la simetría hermitiana del espectro para que la IFFT devuelva valores reales: function senial_rec = reconstruir_por_fft(senial_ruidosa) N = length(senial_ruidosa); Y = fft(senial_ruidosa); % Magnitud del espectro unilateral (DC hasta Nyquist inclusive) mag_pos = abs(Y(1 : floor(N/2) + 1)); umbral = max(mag_pos) * 0.1; % Máscara binaria: 1 donde la señal supera el umbral mascara_pos = double(mag_pos > umbral); % Máscara completa con simetría hermitiana para señal real mascara = zeros(N, 1); mascara(1 : floor(N/2) + 1) = mascara_pos; mascara(floor(N/2)+2 : N) = flipud(mascara_pos(2 : floor(N/2))); senial_rec = real(ifft(Y .* mascara)); end </code></pre> Error que cometí: En versiones anteriores usé N/2</code> directamente para indexar, lo que fallaba cuando N</code> es impar. La solución correcta es floor(N/2)</code> en todos los índices, y construir la simetría conjugada explícitamente. Esto asegura que la máscara tenga exactamente N</code> elementos sin importar la paridad de la longitud de la señal. </blockquote> Resultados</h2> Gráficas Completas</h3> La siguiente figura muestra el pipeline completo para las cuatro señales: señal ideal → señal capturada por el sensor (BPW34→MAX9814) → señal reconstruida → espectro FFT con el SNR obtenido. Lo más notable visualmente: la columna “Reconstruido” (verde) recupera la forma de onda con una fidelidad que hace casi indistinguible el resultado del original, a pesar del ruido visible en la columna “Capturado”. Métricas Cuantitativas</h3> Para no quedarme solo con la inspección visual, calculé tres métricas de error para cada señal. La señal reconstruida se escala antes de comparar (producto interno mínimo cuadrados) para separar el error de fase/forma del error de amplitud: Señal</th> Error RMS</th> SNR (dB)</th> Correlación</th> Resultado</th></tr></thead> Tono puro</td> 2.79 × 10⁻⁶</td> 102.1 dB</td> 1.0000</td> Excelente</td></tr> Acorde</td> 2.87 × 10⁻⁴</td> 62.1 dB</td> 1.0000</td> Excelente</td></tr> Barrido</td> 1.06 × 10⁻²</td> 30.5 dB</td> 0.9996</td> Muy bueno</td></tr> Compleja</td> 3.00 × 10⁻⁴</td> 63.1 dB</td> 1.0000</td> Excelente</td></tr> </tbody></table> Referencia: En audio, un SNR por encima de 20 dB ya es considerado aceptable para reproducción de voz. Por encima de 40 dB es transparente para el oído humano en la mayoría de condiciones. Los resultados de este sistema están entre 30.5 y 102.1 dB. </blockquote> Análisis SVD</h3> El análisis de valores singulares confirma la intuición teórica sobre la complejidad de cada señal: Señal</th> Varianza acumulada en 5 modos</th> </th></tr></thead> Tono puro</td> 99.0%</td> casi toda la energía en 1–2 modos</td></tr> Acorde</td> 84.3%</td> ~3 modos significativos (uno por frecuencia)</td></tr> Barrido</td> 9.3%</td> energía muy distribuida</td></tr> Compleja</td> 88.6%</td> varios modos (armónicos)</td></tr> </tbody></table> El caso del barrido es el más revelador: como la frecuencia cambia continuamente de 100 a 1000 Hz durante 3 segundos, ningún modo fijo puede capturar bien la señal. Necesita cientos de vectores singulares para representarse fielmente. Esto explica matemáticamente por qué tiene el SNR más bajo: el filtro FFT de umbral fijo no es el mejor para señales variantes en el tiempo. El Código Completo</h2> El script completo está disponible en el repositorio del proyecto: Código fuente: script_prueba_demostracion.m</a> — MATLAB R2016b o superior, requiere Signal Processing Toolbox. </blockquote> Conclusiones</h2> Los resultados de la simulación son claros: El filtrado espectral FFT funciona. SNR de hasta 102.1 dB para el tono puro y correlación de 1.0000 en tres de las cuatro señales. El sistema recupera la forma de onda con fidelidad suficiente para reproducción de audio. </li> La señal más difícil es el barrido (chirp). Con 30.5 dB y correlación 0.9996, sigue siendo muy buena, pero revela la limitación del filtro de umbral fijo: no es el mejor para señales cuya frecuencia varía en el tiempo. La STFT (Short-Time Fourier Transform) o una transformada Wavelet serían mejores herramientas para ese caso. </li> El SVD confirma la teoría de espacios vectoriales. El rango efectivo de la señal (medido por los valores singulares) refleja directamente su complejidad: 1–2 modos para el tono puro, cientos para el barrido. </li> El BPW34 supera al LDR. El shot noise del BPW34 es de apenas 0.8% de la señal, frente al 5% típico de un LDR. Esto se traduce directamente en SNR más alto y reconstrucción más fiel. </li> La simulación valida la arquitectura antes del hardware. Ya sé que los algoritmos funcionan. Ahora puedo construir el sistema real con confianza. </li> </ol> Próximos Pasos</h2> El hardware está definido y en camino: Montaje físico: láser D18×65mm sobre soporte Neiko, BPW34 en tubo negro, circuito en protoboard 830 puntos</li> Primera captura real: comparar SNR experimental con el simulado (objetivo: >20 dB)</li> Señales de prueba: tono puro de generador, luego voz, luego música</li> Futuro: STFT para mejorar la reconstrucción del chirp y señales variantes en el tiempo</li> </ul> El siguiente post documentará la construcción del hardware y la primera captura real de audio. Este proyecto es parte del curso de Física II (Ondas, Oscilaciones y Espacios Vectoriales), segundo semestre. Todo el código es de acceso libre en GitHub</a>. Double Slit Experiment Simulation 2025-10-28T00:00:00+00:00 Double Slit Experiment Simulation Project</h1> Project Overview</h2> This project presents a comprehensive computational simulation of the famous double-slit experiment, one of the most profound demonstrations in quantum mechanics. The simulation is designed to run on a Raspberry Pi and includes four different visualization modes, each revealing different aspects of wave-particle duality and quantum behavior. The double-slit experiment, originally performed by Thomas Young in 1801, demonstrates that light and matter exhibit both wave-like and particle-like properties. This project brings this fundamental quantum phenomenon to life through interactive Python simulations. The Four Simulations</h2> 1. Basic Simulation - Static Interference Pattern</h3> The basic simulation displays the classic interference pattern produced when coherent waves pass through two slits. This visualization shows the fundamental wave behavior with bright fringes (constructive interference) and dark fringes (destructive interference). Key Features: Static interference pattern visualization</li> Heat map showing intensity distribution</li> Marked slit positions</li> Clear demonstration of wave superposition</li> </ul> Visual Result: </a> The image shows multiple bright and dark bands, demonstrating how waves from both slits interfere with each other. The cyan stars mark the positions of the two slits. 2. Animated Simulation - Waves in Motion</h3> This is the most visually striking simulation, showing waves propagating in real-time from both slits. The animation provides an intuitive understanding of how interference patterns form dynamically. Key Features: Real-time wave propagation animation</li> Two-panel display: wave amplitude and intensity</li> Color-coded visualization (red = wave crests, blue = wave troughs)</li> Continuous animation showing wave evolution</li> </ul> Visual Result: </a> The left panel shows the wave amplitude with red and blue colors representing crests and valleys. The right panel displays the observable intensity pattern that would be detected on a screen. 3. Interactive Simulation - Parameter Exploration</h3> The interactive simulation allows real-time adjustment of experimental parameters, providing hands-on exploration of how wavelength and slit separation affect the interference pattern. Key Features: Adjustable wavelength (λ) slider</li> Adjustable slit separation (d) slider</li> Toggle between 1 and 2 slits</li> Real-time intensity profile graph</li> Display of Young’s formula: Δy = λL/d</li> </ul> Visual Result: </a> The simulation shows both the 2D interference pattern and a 1D intensity profile, with interactive controls at the bottom for parameter adjustment. 4. Quantum Simulation - The Observer Effect</h3> This is the most impressive and conceptually profound simulation. It demonstrates the quantum mechanical observer effect: the phenomenon where the act of measurement fundamentally changes the behavior of quantum particles. Key Features: Individual particle detection simulation</li> “OBSERVE” button to toggle measurement</li> Side-by-side comparison of theoretical and experimental patterns</li> Accumulation of particle detections over time</li> Clear demonstration of wave function collapse</li> </ul> Visual Result: </a> The left panel shows the theoretical probability distribution, while the right panel accumulates individual particle detections. When observation is activated, the interference pattern disappears. The Observer Effect - The Heart of Quantum Mechanics</h2> The most striking result of this project is the demonstration of how observation changes quantum behavior: </a> Without Observation (Left - Blue):</h3> Multiple interference fringes appear</li> Each particle behaves as a wave</li> The particle passes through both slits simultaneously (quantum superposition)</li> Result: Interference pattern with 6-8 visible bands</li> </ul> With Observation (Right - Red):</h3> Only 2 bands appear</li> Each particle behaves as a classical particle</li> The particle “chooses” one specific slit</li> Result: No interference - just two separate distributions</li> </ul> This is the central mystery of quantum mechanics: the act of observing changes reality itself. Complete Comparison</h2> All four simulations side by side: </a> Top Left: Basic static pattern</li> Top Right: Animated wave propagation</li> Bottom Left: Interactive with adjustable parameters</li> Bottom Right: Quantum simulation with probability distribution</li> </ol> Physics Explained</h2> Wave Behavior (Classical Physics)</h3> When light passes through two slits, each slit acts as a new source of circular waves (Huygens’ Principle). These waves overlap and interfere: Constructive Interference: Wave crests align → bright fringe</li> Destructive Interference: Crest meets trough → dark fringe</li> </ul> The spacing between fringes follows Young’s formula: Δy = λL/d Where: λ = wavelength of light</li> L = distance to screen</li> d = separation between slits</li> </ul> Quantum Behavior</h3> The quantum version reveals three fundamental principles: Wave-Particle Duality: Quantum entities exhibit both wave and particle properties depending on how they are observed. </li> Quantum Superposition: An unobserved particle exists in multiple states simultaneously. In the double-slit experiment, each particle passes through both slits at once. </li> Wave Function Collapse: When we measure which slit a particle passes through, the wave function collapses. The particle is forced to “choose” one path, and the interference pattern disappears. </li> </ol> As physicist Richard Feynman stated: “The double-slit experiment contains the only mystery of quantum mechanics.” </blockquote> Technical Implementation</h2> Platform</h3> Hardware: Raspberry Pi (any model with Python 3)</li> Language: Python 3</li> Libraries: NumPy (numerical computation), Matplotlib (visualization)</li> </ul> Installation</h3> pip3 install numpy matplotlib </code></pre> Execution</h3> Each simulation can be run independently: python3 simulacion_basica.py # Basic simulation python3 simulacion_animada.py # Animated simulation python3 simulacion_interactiva.py # Interactive simulation python3 simulacion_cuantica.py # Quantum simulation </code></pre> Recommended Presentation Order</h2> For maximum impact in a physics competition: Start with Animated Simulation - Captures attention with visual movement</li> Continue with Interactive Simulation - Demonstrates understanding of parameters</li> Finish with Quantum Simulation - Delivers the “wow” moment with the observer effect</li> </ol> Project Significance</h2> This project successfully demonstrates: ✅ Classical wave interference through computational modeling ✅ Quantum superposition and wave function behavior ✅ The observer effect - measurement changing reality ✅ Interactive parameter exploration of physical phenomena ✅ Real-time visualization of abstract quantum concepts The simulations make abstract quantum mechanics concepts tangible and visually accessible, providing an excellent educational tool for understanding one of physics’ most profound experiments. Educational Value</h2> This project bridges the gap between theoretical quantum mechanics and practical visualization. It allows students and audiences to: See wave interference in action</li> Explore how parameters affect patterns</li> Understand the quantum measurement problem</li> Experience the strangeness of quantum behavior</li> </ul> The double-slit experiment has been called “the most beautiful experiment in physics” by Physics World readers. This simulation brings that beauty to life in an interactive, programmable format. Conclusion</h2> The double-slit experiment simulation project successfully demonstrates the fundamental principles of quantum mechanics through four complementary visualizations. From the basic interference pattern to the profound observer effect, each simulation reveals a different facet of wave-particle duality. The project proves that complex quantum phenomena can be effectively simulated and visualized using accessible technology like the Raspberry Pi, making cutting-edge physics education available to everyone. Key Takeaway: Quantum mechanics shows us that nature behaves fundamentally differently at the microscopic scale. The act of observation is not passive - it actively participates in creating the reality we measure. This project makes that profound truth visible and interactive. Project Created By: B.E. Alejandro Technology: Python 3, NumPy, Matplotlib Platform: Raspberry Pi Purpose: Physics Competition / Educational Demonstration Date: October 2025

Double Slit Experiment Simulation

Double Slit Experiment Simulation Project</h1>

The Four Simulations</h2>

The Observer Effect - The Heart of Quantum Mechanics</h2>
The most striking result of this project is the demonstration of how observation changes quantum behavior:</p> </a>

Physics Explained</h2>

Technical Implementation</h2>

Nezzontli - Physics

Reconstrucción de Sonido con Láser y Fotosensor: De la Teoría al Código

Fundamentos Teóricos</h2>

Diagrama del circuito</h3>
</p>

Resultados</h2>

Nezzontli - Physics

Reconstrucción de Sonido con Láser y Fotosensor: De la Teoría al Código

Fundamentos Teóricos</h2>

Diagrama del circuito</h3> </p>

Resultados</h2>

Double Slit Experiment Simulation

Double Slit Experiment Simulation Project</h1>

The Four Simulations</h2>

Physics Explained</h2>

Technical Implementation</h2>

Diagrama del circuito</h3>
</p>