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
177
178
179
180
181
182
183
| # -*- coding: utf-8 -*-
"""
Created on Mon Oct 16 11:19:15 2017
@author: marco3
"""
import numpy as np
import tkinter as tk
import matplotlib
matplotlib.use('TkAgg')
from matplotlib.backends.backend_tkagg import FigureCanvasTkAgg, NavigationToolbar2Tk
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 = NavigationToolbar2Tk(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() |
