Discussion :

[oliv27400]

Bonjour, je réalise diverse opération sur des matrices (multiplication,division termes à termes).Voici le code utilisé, j'expliquerai mon problème par la suite.

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
Amp=[1 10];
eta=[1E-3 1E-2];
f=[5 10];
tmax=1;
 
omega=2*pi()*f;
fe=20*f;
omegaN=omega./sqrt(1-(eta.^2));
a=omegaN/(2*pi());
b=2*eta.*omegaN;
c=eta.*omegaN;
k1=(a.^2)./((a-b).*((c-a).^2+omega.^2));
k2=(b.^2)./((a-b).*((c-b).^2+omega.^2));
k3=sqrt(((c.^2-omega.^2).^2+4*c.^2.*omega.^2)./(((b-c).^2+omega.^2).*((a-c).^2+omega.^2)));
teta=atan((-2*c.*omega)./(c.^2-omega.^2))-atan(omega./(a-c))-atan(omega./(b-c));
 
t=0:1/fe:tmax
t=t';

Lorsque je veux calculer ma variable t (correspondant aux temps en fonction de la frequence fe) matlab m'indique une erreur de dimension de matrice.Or je ne vois pas comment lui faire comprendre que la variable t est différente selon ma fréquence f choisie. (l'objectif étant d'avoir une matrice t à deux dimensions.

L'idée de ce code est de calculé une variable W variant en fonction des paramètres Amp,eta,f sans utilisé de boucle for.

[shaiHulud]

Bonjour,

fe est un vecteur à deux éléments, et 1\fe n'existe pas. 1./fe a du sens, mais aura deux éléments aussi, donc 1:1./fe:tmax n'a aucun sens.

[oliv27400]

Merci, je me rend bien compte que cela n'a pas de sens, mais existe t il une solution sans passer par une boucle for?

[Jerome Briot]

Tu peux utiliser un tableau de cellules :

Code :

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

1
2
t{1} = 0:1/fe(1):tmax
t{2} = 0:1/fe(1):tmax

[shaiHulud]

En admettant que tu veuilles aller de 0 à tmax pour tous les pas contenus dans fe : les grilles résultants n'ont pas la même taille, donc il faut un conteneur de données non homogène (un cell array) comme proposé par Dut

Code :

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

1
2
3
 
t{1} = 0:1/fe(1):tmax
t{2} = 0:1/fe(2):tmax

ou (version boucle)

Code :

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

1
2
 
t= arrayfun(@(x) 0:1/x:tmax,fe,'uni',0);

Dans tous les cas, t n'est pas un matrice mais un ensemble de vecteurs de tailles différentes, donc t' ne donnera surement pas ce que tu attends, et tu risques de passer par une boucle pour faire tes plots.

Une autre solution serait de donner artificiellement à tes matrices la même taille, en répliquant des éléments ou en mettant des nan.

[tachmou]

Bonjour,

Euh, je n'ai pas tout lu dans les détails, mais tu peux faire ça pour te rapprocher de ce que tu cherches à faire, je pense :

Code :

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

1
2
3
4
5
6
7
8
9
10
t1 = 0:1/fe(1):tmax;
t2 = 0:1/fe(2):tmax;
(l1, c1) = size(t1);
(l2, c2) = size(t2);
l = 2;
c = max(c1, c2);
t = zeros(l,c);
t(1,1:c1) = t1;
t(2,1:c2) = t2;
t = t';

N'oublies pas qu'en procédant ainsi, tu n'auras pas autant de valeur dans chaque ti... Et effectivement, si tu dois utiliser des ti avec i très grand, tu peux faire un mix entre ma méthode et celle de Dut. Il suffit de remplacer mes t1 et t2 par ti{1} et ti{2}. Ainsi tu pourras faire une boucle permettant de ne pas être limité à seulement 2 séries de données.

Cordialement,

[EDIT] tiens, j'ai l'impression que la version que je t'ai donné soit résumée par celle de VV33D haha

[oliv27400]

Merci pour vos conseils!