Echantilloner un signal

Version imprimable

10/04/2008, 19h58
bi_turbo

Echantilloner un signal

Bonjour.
J'ai besoin d'une explication sur la démarche à suivre.
J'ai une courbe et j'aimerais des points sur cette courbe à intervalle régulier. Et à partir de ces points en tracer une nouvelle.
J'ai essayé une approximation non polynomiale pensant que c'est la solution, mais hélas non.
Merci d'avance pour vos conseils.
10/04/2008, 20h49
Caro-Line

Bonjour,

Un petit dessin pour illustrer.
Du code pour montrer comment est tracé la courbe.
Comment comptes-tu créer les points sur la courbe ?
10/04/2008, 21h53
bi_turbo
Voilà mon bout de code
Code:

1 2 3 4 Ts=(0:15:31560)'; %tous les 15secondes Y = [ones(size(Ts)) exp(-Ts) Ts.*exp(-Ts)]*a; plot(Ts,Y,'-',Tem,Ech,'p'), grid on
Tem et Ech sont des tabeaus contenant mes données que je veux tracer en fonction de Ts. et ensuite je veux placer un points sur la courbe tous les 15s.
Pour l'instant je peine à trouver la méthode pour placer ces points sur la courbe....
11/04/2008, 01h15
Pierre845

essayes interp
11/04/2008, 11h17
Jerome Briot

Citation:

Envoyé par Pierre845

essayes interp

INTERP est une fonction de la Signal Processing Toolbox

Si tu n'as pas cette toolbox, tu peux commencer par regarder la fonction INTERP1
11/04/2008, 15h36
bi_turbo

Rebonjour!
La commande interp1 aurait été une bonne solution.
en lisant les 'help ...' associées à ces commande, il est dit que mon vecteur ou matrice doit avoir des données de croissance monotone
au lieu de M=1:1:20
mon vecteur est de style mon_vec= [1 2 23 8 5 1 100 59 21 ....]
aucune des interpollations ne corresponds.

Je continue à chercher...

Essaie ceci :

Code:

1
2
3
4
5
6
7
8
9
10
11
12
V = [1 2 23 8 5 1 100 59 21];
nV = numel(V);
 
x = 1:nV;
xi = linspace(1,nV,30);
method = 'linear';
Vi = interp1(x,V,xi,method);
 
figure
p=plot(x,V,'b-s',xi,Vi,'ro');
legend(p,{'V' 'Vi'})
title(sprintf('%s interpolation',method))

Ou alors je n'ai pas bien compris ton problème... :?

Bonjour à tous.
Désolé de me manifester aussi tard
Citation:
Envoyé par Dut
Essaie ceci :
Code:

1 2 3 4 5 6 7 8 9 10 11 12 V = [1 2 23 8 5 1 100 59 21]; nV = numel(V); x = 1:nV; xi = linspace(1,nV,30); method = 'linear'; Vi = interp1(x,V,xi,method); figure p=plot(x,V,'b-s',xi,Vi,'ro'); legend(p,{'V' 'Vi'}) title(sprintf('%s interpolation',method))
Ou alors je n'ai pas bien compris ton problème... :?
Merci DUT, La ligne en gras rouge correspond à ce que je recherche. Toute fois celà m'amène à me poser une question. Peux -t-on s'affranchir du lien qu'il y a entre X et V? Je veux dire au lieu que x fasse de 1 à 9 par pas de 1,
je gradue x en temps (s) x=1:105:3600.
Du moins comment contourner

Citation:

Y must have length(X) rows.

?
Merci