Discussion :

[ssarah.z]

salut j'ai un programme qui traite un signal sweep mais je ne maitrise pas le matlab,expliquez moi ce que dit le programme et comment l'écrire en langage C merci de m'aider;voici la bête lol

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
clear all
clc
 
t=[0:0.0001:2];
F1=10;F2=30;T=2; 
 
if t>(T/2)
At=0;
else 
At=1;
end 
a=(F2-F1)/(2*T);
ft=a*(t.^2)+F1*t;
Ot=2*pi*ft;
St=At*sin(Ot);
%%%%%%%représentation d'un up sweep%%%%%%%%%%%%%%%%%%%%%%%%%
plot(t,St);title('b');grid on
 
%%%%%%%représentation de la valeur absolu d'un up sweep%%%%%%%%%%%%%%%%%%%%%%%%%
b=abs(St);
% plot(t,b);title('b');grid on

[mr_samurai]

ton code ne fait pas traitements , et ne prend pas d'entrées, il généere juste un vecteur ( ou function ). C'est le up.

PS :

Code :

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

1
2
3
4
5
6
 
if t>(T/2)
At=0;
else 
At=1;
end

Dans ton code At est un vecteur je suppose ?
Il faut boucler sur t :

Code :

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

1
2
3
4
5
6
7
8
9
 
At = zeros(size(t));
for i=1:length(t)
     if t(i)>(T/2)
       At(i)=0;
     else 
       At(i)=1;
     end 
end

++

[ssarah.z]

merci pour la precision mais j ai pas trops compris qu est ce qu un up?et comment je fais si je veus avoir le spectre du sweep tout en utilisant la fft,merci

[mr_samurai]

Ce que j'ai compris est que le up est la fonction :

up(t) = X[0,T/2](t)*sin(2*pi*(a*t^2+ F1*t))

avec X[0,T/2] l'indicatrice de [0,T/2], Forcement up est de periode T.....

je pense qu'il faut faire un fft(up) ...., mais c'est pas ton boulot cela ?


++ bonne courage

[ssarah.z]

oui tu as raison ce n'est pas trop mon domaine. Moi je fais géophysique spécialisée en diagraphie mais j'en ai besoin pour résoudre un problème. J'ai besoin d injecter le sweep à la fft pour avoir son spectre;alors je pensais que ce programme en matlab le faisait et qu'il me restait qu'à le traduire,merci de ton aide