reconstruction d'un signal à partir de spectres
Bonjour,
je vous explique comment j'ai procédé pour reconstruire un signal pouvez vous me dire si c'est bien correct et si mon programme de reconstruction est bon:
1°) un logiciel de simulation me donne pour plusieurs points des efforts dans des directions x et y sous forme complexe pour plusieurs frequence
J'ai donc des données d'entree du type:
Fx=[5+2i 3+3i 2i 2]
Fy=[3i 3 2i 6i+9i]
2°) se qui m'interesse c'est la force dans une autre direction (radial)
j'ai donc fait en gros pour toutes les frequences et tous les points
Fr(n)=Fx(n)*cos(n)+Fy(n)*sin(n)
c'est à dire que j'ai projeté les vecteurs Fx et Fy dans une direction
3°) ensuite la grandeur qui m'interesse est une moyenne pour chaque frequence de Fr j'ai donc classé Fr dans un ordre pour pouvoir integrer pour chaques frequences.
4°) j'ai donc a present un nombre complexe pour chaques frequences:
=>j'en ai pris le module et l'argument
=>j'ai rentré cela dans la formule:
P=module(n)*cos(2*pi*n*t+argument(n)*+module(n)*i*sin(2*pi*n*t+argument(n);
Avec n les fréquences de mon signal qui m'interesse allant de 1 à 10Hz car seul les 10premiers harmoniques m'interessent.
Le signal d'entree etait de 1Hz
Code:
1 2 3 4 5 6 7 8 9 10
| z = 1; %initialisation d'un compteur pour les indices de P
P = zeros(1,length(0:0.0005:2)); % on initialise P cela donne P=0, utilie pour 1ere boucle
for t=0:0.0005:2;
for n=1:1:n_max;
P(z)=P(z)+module(n)*cos(2*pi*n*t+argument(n)*+module(n)*i*sin(2*pi*n*t+argument(n);
end
z = z+1;
end
P_physique=real(P);
T=[0:0.0005:2]; |