Discussion :

[MaximeGeorges]

Bonjour,

Je ne sais pas si je suis au bon endroit pour cela, mais dans le cadre de mes études , je souhaite modifier le spectre d'un son enregistré , sur phyton.
Plus précisément, j'enregistre un son, avec audacity par exemple, j'en extrait le spectre sous forme de fichier texte, puis j'importe les listes Fréquences / Amplitudes sur python .
Ensuite, je les modifie , mais c'est là qu’apparaît mon problème : je souhaiterais reformer le signal temporel initial, j'ai lu que c’était possible sur python avec une transformée de Fourrier inverse, mais je ne sais pas du tout comment faire pour être honnête ..

Donc si quelqu'un à des éléments qui pourraient me faire avancer !

Merci d'avance !

[marco056]

Comme cela par exemple ?

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
import numpy as np
import tkinter as tk
import matplotlib
matplotlib.use('TkAgg')
from matplotlib.backends.backend_tkagg import FigureCanvasTkAgg, NavigationToolbar2TkAgg
from matplotlib.figure import Figure
from tkinter.constants import TOP, BOTH, BOTTOM, LEFT, RIGHT
 
Umax=1
omega=1
nb_val=500
t=np.linspace(0,4*np.pi/omega,nb_val)
 
def carre(nb_harm):
    tmp_carre=np.zeros(nb_val)
    for n in range(0,nb_harm):
        tmp_carre=tmp_carre+np.sin((2*n+1)*omega*t)/(2 *n +1)
    scale_carre=[4*Umax/np.pi]*nb_val
    tmp_carre=scale_carre*tmp_carre
    return tmp_carre
 
def triangle(nb_harm):
    tmp_triangle=np.zeros(nb_val)
    for n in range(0,nb_harm):
        tmp_triangle=tmp_triangle+np.cos((2*n+1)*omega*t)/((2*n+1)*(2*n+1))
    scale_triangle=[4*Umax/np.pi]*nb_val
    tmp_triangle=scale_triangle*tmp_triangle
    return tmp_triangle
 
def scie(nb_harm):
    tmp_scie=np.zeros(nb_val)
    for n in range(1,nb_harm):
        tmp_scie=tmp_scie+np.sin(2*n*omega*t)/(n)
    scale_scie=[-2*Umax/np.pi]*nb_val
    tmp_scie=scale_scie*tmp_scie
    return tmp_scie
 
def arche(nb_harm):
    tmp_arche=np.zeros(nb_val)
    for n in range(1,nb_harm):
        tmp_arche=tmp_arche+np.cos(2*n*omega*t)/(n)
    scale_arche=[-2*Umax/np.pi]*nb_val
    tmp_arche=1+scale_arche*tmp_arche
    return tmp_arche
 
def mono(nb_harm):
    tmp_mono=np.zeros(nb_val)
    for n in range(1,nb_harm):
        tmp_mono=tmp_mono+np.cos(2*n*omega*t)/(4*n**2-1)
    scale_mono=-(2/np.pi)
    tmp_mono=(1/np.pi)+(0.5*np.sin(omega*t))+scale_mono*tmp_mono
    return tmp_mono
 
def double(nb_harm):
    tmp_double=np.zeros(nb_val)
    for n in range(1,nb_harm):
        tmp_double=tmp_double+np.cos(2*n*omega*t)/(4*n**2-1)
    scale_double=-(4/np.pi)
    tmp_double=(2/np.pi)+scale_double*tmp_double
    return tmp_double
 
def essai(nb_harm):
    tmp_arche=np.zeros(nb_val)
    for n in range(1,nb_harm):
        tmp_arche=tmp_arche+np.cos(2*n*omega*t)/(n)
    scale_arche=[-2*Umax/np.pi]*nb_val
    tmp_arche=scale_arche*tmp_arche
    return tmp_arche
 
def figure_hide(figure):
    figure.canvas.get_tk_widget().pack_forget()
 
def figure_raz(figure):
    axe = figure.axes[0]
    axe.clear()    
    figure.canvas.draw()
 
def figure_show(figure):
    figure.canvas.get_tk_widget().pack(side=TOP, fill=BOTH, expand=1)
 
def plot_carre(nb_harm):
    nb_harm = val_harm.get()
    a.plot(carre(nb_harm))
    canvas.draw()
 
def plot_triangle(nb_harm):
    nb_harm = val_harm.get()
    a.plot(triangle(nb_harm))
    canvas.draw()
 
def plot_scie(nb_harm):
    nb_harm = val_harm.get()
    a.plot(scie(nb_harm))
    canvas.draw()
 
def plot_arche(nb_harm):
    nb_harm = val_harm.get()
    a.plot(arche(nb_harm))
    canvas.draw()
 
def plot_mono(nb_harm):
    nb_harm = val_harm.get()
    a.plot(mono(nb_harm))
    canvas.draw()
 
def plot_double(nb_harm):
    nb_harm = val_harm.get()
    a.plot(double(nb_harm))
    canvas.draw()
 
def plot_essai(nb_harm):
    nb_harm = val_harm.get()
    a.plot(essai(nb_harm))
    canvas.draw()
 
def quitter():
    root.quit()
    root.destroy()
 
def get_nb(val):
    global nb_harm
    nb_harm = val.get()
 
root = tk.Tk()
root.wm_title("Recomposition de Fourier")
 
frame = tk.Frame(root)
figure = Figure(figsize=(8,3), dpi=100)
a = figure.add_subplot(111)
 
canvas = FigureCanvasTkAgg(figure, master=root)
#zone = canvas.get_tk_widget()
#zone.pack(side=tk.TOP, fill=BOTH, expand=1)
toolbar = NavigationToolbar2TkAgg(canvas, root)
toolbar.update()
canvas._tkcanvas.pack(side=TOP, fill=BOTH, expand=1)
 
button1=tk.Button(frame,text='Carré',command=lambda: plot_carre(val_harm))
button1.pack(side=LEFT)
button2=tk.Button(frame,text='Triangle',command=lambda: plot_triangle(val_harm))
button2.pack(side=LEFT)
button3=tk.Button(frame,text='Scie',command=lambda: plot_scie(val_harm))
button3.pack(side=LEFT)
button4=tk.Button(frame,text='Arche',command=lambda: plot_arche(val_harm))
button4.pack(side=LEFT)
button5=tk.Button(frame,text='Mono',command=lambda: plot_mono(val_harm))
button5.pack(side=LEFT)
button6=tk.Button(frame,text='Double',command=lambda: plot_double(val_harm))
button6.pack(side=LEFT)
button7=tk.Button(frame,text='Essai',command=lambda: plot_essai(val_harm))
button7.pack(side=LEFT)
 
button10=tk.Button(frame, text='Raz', command=lambda f=figure: figure_raz(f))
button10.pack(side=LEFT)
button11=tk.Button(frame, text='Hide', command=lambda f=figure: figure_hide(f))
button11.pack(side=LEFT)
tk.Button(frame, text='show', command=lambda f=figure: figure_show(f)).pack(side=LEFT)
 
 
harmo=tk.Label(frame, text="Harmoniques")
harmo.pack(side=BOTTOM)
val_harm = tk.IntVar() # On definit val_harm
val_harm.set(1)        # On donne la valeur que prendra le curseur au départ
 
# Création d'un widget Curseur
echelle_harm = tk.Scale(frame,length=200, orient=tk.HORIZONTAL, troughcolor ='LightYellow2', \
sliderlength =30, showvalue=1,from_=1,to=30, resolution=1, tickinterval=6, \
variable = val_harm, command = get_nb(val_harm)) # width=10
echelle_harm.set(0),
echelle_harm.pack(side=BOTTOM)
 
button20=tk.Button(frame, text='Quit',command=quitter)
button20.pack(side=RIGHT)
 
frame.pack(side=BOTTOM)
tk.mainloop()