Egalité de parseval et Matlab

[invite55481d20]

Bonjour,

Je cherche à vérifier l'égalité de Parseval en Matlab sur un bruit blanc.
Voilà mon code:

N=2^10;
X=rand(N,1);  //def de mon bruit blanc

Fs=100;
Ts=1/Fs;

// 1. Calcul de la DSP
TF_X=fft(X);
DSP_X=abs(TF_X)^2;
puiss_X=sum(DSP_X)*(Fs/Nf)

// 2. Calcul de la puissance du signal en temporel
A=Ts*abs(X)^2;
puissance_X=(1/N)*sum(A)*Ts

Bien évidemment je ne trouve pas la même chose...
Malheureusement je ne maitrise pas trop ces histoires de puissances, Parseval etc...
Est-ce que quelqu'un pourrait m'aider ?
merci beaucoup !
-----

[invite55481d20]

Bonjour,

Je cherche à vérifier l'égalité de Parseval en Matlab sur un bruit blanc.
Voilà mon code:

N=2^10;
X=rand(N,1);  //def de mon bruit blanc

Fs=100;
Ts=1/Fs;

// 1. Calcul de la DSP
TF_X=fft(X);
DSP_X=abs(TF_X)^2;
puiss_X=sum(DSP_X)*(Fs/Nf)

// 2. Calcul de la puissance du signal en temporel
A=Ts*abs(X)^2;
puissance_X=(1/N)*sum(A)*Ts

Bien évidemment je ne trouve pas la même chose...
Malheureusement je ne maitrise pas trop ces histoires de puissances, Parseval etc...
Est-ce que quelqu'un pourrait m'aider ?
merci beaucoup !

Personne pour m'aider?

[invite55481d20]

J'ai encore essayer mais je n'y arrive pas. Est-ce que qqun pourrait m'aider ?

[phryte]

Bonjour.
AS-tu mis le point avant la puissance :
DSP_X=abs(TF_X).^2;
A=Ts*abs(X).^2;
Et que vaut Nf ?

[A voir en vidéo sur Futura]

[invite55481d20]

Merci de ta réponse. Que signifie le point ?
En attendant des réponse j'ai codé ceci :

N=2^10;
X=rand(N,1);

      puiss_temp=sum(abs(X)^2);

tf_X=fft(X);
DSP_X=abs(tf_X)^2;

      puiss_freq=(1/N)*sum(DSP_X);

Et là je trouve bien que puiss_temp = puiss_freq.

Mais que signifie le .^2 ?

[phryte]

Bonjour.
Le point fait faire le calcul élément par élément.

[invite55481d20]

Bonjour.
Le point fait faire le calcul élément par élément.

ok merci.

J'ai vraiment du mal avec ces histoires de Parseval etc...

N=2^10;
X=rand(N,1);

      puiss_temp=sum(abs(X).^2);

tf_X=fft(X);
DSP_X=abs(tf_X).^2;

      puiss_freq=(1/N)*sum(DSP_X);

ce code me donne puiss_temp = puiss_freq

mais il n'y a pas un souci avec les (1/N) ???
apparemment la puiss_temp aurait u facteur 1/N normalement, non ?