Nom : Capture d'écran 2023-11-24 104019.png
Affichages : 342
Taille : 81,9 Ko

Voici ce que j'ai fait mais je n'obtiens pas un résultat pou mon intégrale.
Pourriez-vous m'aider, s'il vous plaît ?

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
 
import pywt
import numpy as np
import matplotlib.pyplot as plt
import scipy.integrate as integrate
 
k=0.1
 
# Paramètres
wavelet_name = 'dmey'  # Ondelette de Meyer discrète
scale = 10  # Échelle de l'ondelette de Meyer
 
# Générer la fonction d'échelle
wavelet = pywt.Wavelet(wavelet_name)
phi, _, x = wavelet.wavefun(level=scale)
 
# Calculer la transformée de Fourier de la fonction d'échelle
fourier_transform_phi = np.fft.fftshift(phi)
modsquare = np.abs(fourier_transform_phi) ** 2
freq = np.fft.fftfreq(len(x), d=x[1]-x[0])
 
def integrand(t,n):
    intermediaire = np.exp(1j * t * n)
    moninter = np.abs(np.sinc(t/ 2)) ** (2 * k) * modsquare
    return intermediaire * moninter
 
 
n = 2
# Intégration numérique de la fonction intégrande
res = integrate.quad(integrand, -4 * np.pi / 3, 4 * np.pi / 3 , args=(n,))
 
 
print(res)