# -*- coding: utf-8 -*-
"""
Created on Sat Jul 22 18:32:24 2017

@author: skate
"""

#%%

import random as rand
from scipy.signal import *
from cmath import *
import wave, struct, math
import numpy as np
import matplotlib.pyplot as plt
import scipy.io.wavfile as wavfile

#
N=10 ## complexité du modèle LPC
# 

## on charge le signal de la voyelle

signal = np.load("voyelle_a.npy")


#%% Question 1: analyse de la transformée de Fourier
y = np.fft.fft(signal[:10000])
plt.plot(np.abs(y)[:300])

#%% Question 2: analyse du pitch
analyse_pitch=[]


#%%
#### Question 3: equations de Yulle Walker.
m=40000
signal_final=[]
Rx=[]


## Calculer Rx

## Résoudre les équations de Yulle Walker

#%%

## Questions bonus: synthèse vocale

   
signal_reconstruit=[]

## compléter en remplissant le tableau du signal reconstruit avec un signal de synthèse constitué 
## d'un train d'impulsions passant à travers le filtre autoregressif calculé précédemment

wavfile.write("a_synthetise.wav",44100,np.array(signal_reconstruit))


#%%