Discussion :

[MikaFSL]

Bonjour,

Je cherche à re-echantillonner une courbe (y=f(x) avec y: parametres S et x: frequences) et utilise donc la commande "spline". Suite à une analyse approfondie de mon programme (taille des vecteurs, type de vecteur: ligne ou colonne, etc.), je ne comprends toujours pas pourquoi celle-ci ne fonctionne pas. Peut etre que l'un d'entre vous saurait resoudre ce probleme?!

Mon programme est le suivant:

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
% Import fichier de mesures (Lecture des fichiers touchstone du systeme d'injection et de la charge)
 
[S_inj,freq1] = read_snp(6); %read_snp est une fonction speciale permettant de lire les fichiers touchstone
[S_load,freq2] = read_snp(6);
 
% Re-echantillonage (but: vecteur frequence identique pour les 2 fichiers)
 
x_inj=freq1;
x_load=freq2;
min_freq = max(min(freq1), min(freq2));
max_freq = min(max(freq1), max(freq2));
freq = min_freq:1000:max_freq;
 
Size_matrix=size(S_inj);
y_inj=[];
S_inj_interp=[];
S_inj_int=[];
 
for j=1:1:Size_matrix(1)
    for i=1:1:Size_matrix(2)
        for tt=1:1:size(freq1')
 
            Sinj=100*abs(S_inj(i,j,tt));
            y_inj=[y_inj,Sinj];
 
        end
 
        S_inj_interp=spline(x_inj, y_inj,freq)';
        y_inj=[];
 
    end
end

[Jerome Briot]

pourquoi celle-ci ne fonctionne pas.

Et quel est le problème exactement ?
Les résultats donnés ne sont pas bons ?
MATLAB renvoie un message d'erreur ?

[MikaFSL]

En fait Matlab me renvoi un message d'erreur qui est le suivant:

Code :

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

1
2
3
4
NaN/Inf breakpoint hit for ppval on line 71.
Stopping at next line.
71  % adjust for troubles, like evaluation sites that are NaN or +-inf
72  infxs = find(xs==inf); if ~isempty(infxs), index(infxs) = l; end

Le probleme vient du nouveau vecteur frequence a prioris (jme trompe peut etre). "ppval" est utilisé par la commande "spline". Apparement une valeur des valeur du vecteur "freq" lui poserait probleme à moins le probleme vienne d'ailleur....