python

import matplotlib.pyplot as plt
import numpy as np

# Configuració
vel_so = 340  # m/s
distancies = np.linspace(0, 15, 500)  # metres
retards_ms = (2 * distancies / vel_so) * 1000  # ms

# Crear figura
fig, ax = plt.subplots(figsize=(14, 8))

# Dibuixar zones de color
# Zona 1: Fusió total (<1 ms)
ax.fill_between(distancies, 0, retards_ms, where=(retards_ms <= 1), 
color='#006400', alpha=0.3, label='Fusió total (<1 ms)')

# Zona 2: Precedència per a clics (1-5 ms)
ax.fill_between(distancies, 0, retards_ms, where=((retards_ms > 1) & (retards_ms <= 5)), 
color='#FFD700', alpha=0.4, label='Precedència clics (1-5 ms)')

# Zona 3: Intel·ligible per a parla (5-40 ms)
ax.fill_between(distancies, 0, retards_ms, where=((retards_ms > 5) & (retards_ms <= 40)), 
color='#90EE90', alpha=0.4, label='Intel·ligible parla (5-40 ms)')

# Zona 4: Límit crític (40-50 ms)
ax.fill_between(distancies, 0, retards_ms, where=((retards_ms > 40) & (retards_ms <= 50)), 
color='#FFA500', alpha=0.5, label='Límit crític (40-50 ms)')

# Zona 5: Eco + inintel·ligible (>50 ms)
ax.fill_between(distancies, 0, retards_ms, where=(retards_ms > 50), 
color='#8B0000', alpha=0.4, label='Eco perceptible (>50 ms)')

# Dibuixar línia principal
ax.plot(distancies, retards_ms, 'b-', linewidth=2.5, label='Retard = 2d/340 × 1000')

# Línies verticals de referència
ax.axvline(x=0.17, color='gray', linestyle='--', alpha=0.7, linewidth=1)
ax.axvline(x=0.85, color='gray', linestyle='--', alpha=0.7, linewidth=1)
ax.axvline(x=6.8, color='gray', linestyle='--', alpha=0.7, linewidth=1)
ax.axvline(x=8.5, color='red', linestyle='--', linewidth=2, label='Llindar eco (8.5 m / 50 ms)')

# Línies horitzontals de referència
ax.axhline(y=1, color='gray', linestyle=':', alpha=0.7, linewidth=1)
ax.axhline(y=5, color='gray', linestyle=':', alpha=0.7, linewidth=1)
ax.axhline(y=40, color='orange', linestyle=':', linewidth=1.5)
ax.axhline(y=50, color='red', linestyle=':', linewidth=1.5)

# Punts clau
punts_dist = [0.17, 0.85, 1.7, 3.4, 5.1, 6.8, 8.5, 10, 12]
punts_ms = [1, 5, 10, 20, 30, 40, 50, 58.8, 70.6]
ax.scatter(punts_dist, punts_ms, color='red', zorder=5, s=50)

# Etiquetes pels punts
for d, ms in zip(punts_dist, punts_ms):
ax.annotate(f'{d:.1f}mn{ms:.0f}ms', 
xy=(d, ms), 
xytext=(5, 10), 
textcoords='offset points',
fontsize=8,
bbox=dict(boxstyle='round,pad=0.3', facecolor='white', alpha=0.8))

# Configurar eixos
ax.set_xlabel('Distància a la paret (metres)', fontsize=12, fontweight='bold')
ax.set_ylabel('Retard del primer eco (milisegons)', fontsize=12, fontweight='bold')
ax.set_title('Efecte de precedència (Wallach et al. 1949)nRelació distància-paret vs. intel·ligibilitat de la parla', 
fontsize=14, fontweight='bold', pad=20)

ax.set_xlim(0, 15)
ax.set_ylim(0, 90)
ax.grid(True, alpha=0.3, linestyle='--')

# Llegenda
ax.legend(loc='upper left', fontsize=10, framealpha=0.9)

# Text explicatiu
text_explicatiu = """Zones d'intel·ligibilitat:
• <1 ms: Fusió total (no rellevant en sales reals)
• 1-5 ms: Efecte de precedència màxim per a clics
• 5-40 ms: Zona segura per a parla intel·ligible
• 40-50 ms: Límit crític, comença a afectar
• >50 ms: Eco perceptible (8.5m) → inintel·ligible

Velocitat del so: 340 m/s
Retard (ms) = (2 × distància / 340) × 1000"""

ax.text(0.98, 0.02, text_explicatiu, transform=ax.transAxes,
fontsize=9, verticalalignment='bottom', horizontalalignment='right',
bbox=dict(boxstyle='round', facecolor='wheat', alpha=0.8))

plt.tight_layout()
plt.savefig('reverberacio_precedencia.png', dpi=150, bbox_inches='tight')
plt.show()

# Taula addicional
print("n" + "="*70)
print("TAULA DE CONVERSIÓ DISTÀNCIA → RETARD")
print("="*70)
print(f"{'Distància (m)':<15} {'Retard (ms)':<15} {'Percepció':<30}")
print("-"*70)

valors_mostra = [0.17, 0.5, 0.85, 1.0, 1.7, 2.0, 3.0, 3.4, 4.0, 5.0, 5.1, 6.0, 6.8, 7.0, 8.0, 8.5, 9.0, 10.0, 12.0, 15.0]

for d in valors_mostra:
retard = (2 * d / 340) * 1000
if retard <= 1:
percepcio = "Fusió total"
elif retard <= 5:
percepcio = "Precedència clics"
elif retard <= 40:
percepcio = "Intel·ligible parla"
elif retard <= 50:
percepcio = "Límit crític"
else:
percepcio = "ECO + inintel·ligible"

print(f"{d:<15.2f} {retard:<15.1f} {percepcio:<30}")

print("="*70)
print("n Gràfic desat com a 'reverberacio_precedencia.png'")

---------