#!/usr/bin/env python3
# -*- coding: utf-8 -*-
"""
Created on Thu Mar 23 15:02:00 2023

@author: julie
"""
import matplotlib.pyplot as plt
import numpy as np


S0 = 1  # Amplitude 
f0=1 # fréqence du signal moyen en Hz
omega0 = 2*np.pi*f0 # Pulsation de fréquence centrale ou moyenne
T0 = 1/f0  # Période moyenne du signal
print(T0)

""" Tracé de s(t) signal sinusoïdal pur """

def ssin(S,omega,t) :
    y = S*np.cos(omega*t)
    return y

NT = 5001 # nombre de dates
'Axe des temps'
t = np.linspace(0,25*T0,NT) # axe des temps

sin1 = ssin(S0,omega0,t)

""" Somme de deux signaux à omega0+- Deltaomega/2 """

Deltaomega=omega0/10
sin21 = ssin(S0,omega0-Deltaomega/2,t)
sin22 = ssin(S0,omega0+Deltaomega/2,t)
sin2 = (sin21 + sin22)/2
plt.plot(t,sin21 + 5,label = '$s_1$')
plt.plot(t,sin22 + 2.5, label = '$s_2$')
plt.plot(t,sin2, label = '$\dfrac{s_1+s_2}{2}$')
plt.xlabel('temps')
plt.title('Somme de deux signaux sinusoïdaux de fréquences proches')
plt.grid()
plt.legend()
plt.show()