Discussion :

[Novae]

Bonjour.
Je débute sous MATLAB
J'aurais 2 questions :
En langage C par exemple , pour tracer une fonction y=x^2 par exemple , il fallait multiplier par la taille du pixel ( car en C on travaille sur les pixels ) pour pouvoir observer une figure .
et avec MATLAB ? quand on lui demande de tracer y=x^2 ( en imagesc ) il travaille comment lui ? Directement sur les x ou sur les pixels ? Dois-je multiplier par la taille du pixel ? Autrement dit , dois-je faire
y=x^2 ou y=taille_pixel*taille_pixel*x^2 ?

2ème question ; j'ai par exemple dans un fchier FONC.m

Code :

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

1
2
3
4
function Z = FONC(X,Y)
Z= X.^2 + Y.^2 
%% d’ailleurs petite question , il faut mettre .^ ? on considère les x,y comme des tableaux ou juste des variables ? 
end

Dans le programme :

Code :

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

1
2
3
z=FONC(x,y);
imagesc(x,y,z)
colormap(gray);

Ca marche comme commande imagesc avec des fonctions f(x,y) ?
voila
merci de vos futures réponses
cordialement
Novae

[le fab]

pour les pixel, j'en ai aucune idée
par contre .^ rend ta fonction plus flexible, lui permettant de travailler sur des vecteur et/ou des matrices
si tu ne met pas le point, tu devra utiliser une boucle for (ou une fonction du type arrayfun) pour calculer FONC pour toutes les valeurs de X et Y

[FLB]

Salut,
imagesc n'utilise pas de notion de pixel d'écran, regarde ce que fait ceci :

Code :

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

1
2
figure;imagesc(rand(10,10));
figure;imagesc(rand(100,100));

Tu peux gérer la résolution avec axis.

Ton utilisation des fonctions est bonne, mais je ne vois pas ce que tu veux dire avec

Ca marche comme commande imagesc avec des fonctions f(x,y) ?

[Novae]

Merci à vous deux.
En fait , dans TOUS les exemples que j'ai vu , imagesc c'est pour afficher des images qu'on charge par exemple. elles sont déjà existantes et quand on a des fonctions mathématiques , tous les exemples nous donnent des "plot". Ce que je veux faire en fait c'est un genre de DensityPlot comme sur mathematica en niveau de gris


J'ai écris dans le fichier fonc.m ( pas faire attention à a et b , j'ai mis ca pour l'exemple

Code :

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

1
2
 
Intensité = a*X.^2 + b*Y.^2

dans le programme principale

Code :

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

1
2
3
4
5
 
x=-1:0.1:1;
y=-1:0.1:1;
In= fonc(x,y)
imagesc(x,y,In)

Bon j'ai quelque chose en sortie mais c'est pas du tout çà mais au moins ça bug pas ^^

[le fab]

Ce que je veux faire en fait c'est un genre de DensityPlot comme sur mathematica en niveau de gris

dans ce cas regarde plutôt du coté de fonctions comme pcolor peut etre
(et colormap)

[Novae]

rha il me met Color data input must be a matrix .

Pourtant X et Y sont des vecteurs. L'intensité que je calcule est forcement une matrice non ?

[le fab]

si tu as écris :

Code :

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

1
2
3
4
x=-1:0.1:1;
y=-1:0.1:1;
In= fonc(x,y)
pcolor(x,y,In)

In n'est probablement pas une matrice !
va voir du coté de

Code :

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

doc meshgrid

[FLB]

Le problème est dans ton calcul de Z : si x et y sont des vecteurs de taille n et m, il faut que Z soit de taille n*m.
Le plus simple est de passer par meshgrid :

Code :

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

1
2
3
4
5
x = -10:10;
y = -1:0.1:1;
[xx,yy] = meshgrid(x,y);
pcolor(x,y,sin(xx).*sin(yy));
colormap('gray');

[Novae]

Merci , j'ai fais meshgrid.
En fait je crois comprendre ce qu'il ne va pas( mais j'arrive pas à la résoudre )
c'est dans ma fonction mathématique ;
en gros :
j'ai une fonction R qui dépend de x et y
size(R) = 201 201

Une fonction O qui dépend de x et y
size(O) = 201 201

et la le problème : J'ai une fonction Somme qui fait une somme discrète de fonction de x et y , genre

Code :

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

1
2
3
 
indice = 1:10
Somme = sum(A(indice).*X.^2 + B(indice).*Y.^2

et la size(Somme) = 1 201
le problème c'est que le multiplie O par Somme donc ça va pas au niveau des dimensions

EDIT : Bon apparemment la fonction sum fait le produit scalaire du tableau A(indice) par X donc ça foire tout. J'ai utilisé la boucle for et les dimensions sont correctes.

Waou c'est dur la manipulation vecteurs/matrices sous matlab !

[FLB]

J'ai utilisé la boucle for et les dimensions sont correctes.

La grande majorité des opérations peuvent se faire sous la forme vectorielle, matlab n'est pas fait pour traiter les boucles rapidement. Montre nous ta boucle si tu veux qu'on jette un coup d'oeil pour vectoriser l'opération.

Waou c'est dur la manipulation vecteurs/matrices sous matlab !

C'est vrai, on met un peu de temps à s'y faire, mais ça devient un outil puissant une fois maitrisé .

[Novae]

ahhh ca marche !
Bon je met la boucle for :

Code :

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

1
2
3
4
5
 
 
for indice = 1:10
SOMME=SOMME + a(indice).*K2xeq(indice).*K2yeq(indice).*(Nxeq(indice).*X.^2+Nyeq(indice).*Y.^2)).*(Mxeq(indice).*X.^2+ Myeq(indice).*Y.^2));
end

Il manque peut être des parenthèses , la formule est plus longue mais ce sont des constantes donc je l'ai ai supprimées sur le forum , j'ai peut être oublié des parenthèses.


Pour trouver mon erreur j'avais fais un petit test dans la console :

Code :

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

1
2
3
4
5
 
A=[1 2 3 4];
indice=1:4;
Somme = sum(A(indice).*A)
Somme=30 ( résultat )

Donc il me donne un scalaire alors que je voulais un vecteur , c'est pour ça que j'ai utilisé for. Sum fait le produit scalaire non ?

De plus je pense que A(indice) était un scalaire ( car A(1) =1 ) mais A(indice)*A erreur de dimension.

[le fab]

non, sum fait la somme de tous les éléments de ton vecteurs ou de tous les éléments de chaque colonnes d'une matrice (dans ce cas le retour est un vecteur)
A(indice) est une matrice :

Code :

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

1
2
3
4
5
A(indice)
 
ans =
 
     1     2     3     4

en fait cette technique permet de trier/filtrer une matrice/un vecteur
pour avoir un scalaire il faut que indice soit aussi un scalaire

edit :
peut être voulais tu utiliser cumsum ??

[Novae]

Arf avec la commande pcolor , c'est pas trop la même chose que densityplot , y'a comme une espèce de quadrillage

[le fab]

ajoutes la ligne

Code :

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

shading interp

ou

Code :

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

set(h, 'LineStyle', 'none');

après avoir écrit

Code :

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

h = pcolor(...);

[Novae]

Ahhhh merci lefab c'est beaucoup mieux !

J'ai regardé un peu cumsum , je vais tester ça !
merci

EDIT : Ah ouai , cumsum ça donne pareil , niquel
Merci beaucoup

[Novae]

J'aurais 2 petites question svp
Une fois ma figure faite ( grâce a pcolor ) , les commandes imagesc,imread etc ne servent à rien non ? ces commandes sont utilisées quand on a des images extérieures non ?

Aussi , j'ai ma belle figure et je voudrais la trasnformée un peu ( les niveaux de gris ) :
J'ai :

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
 
figure(1);
h=pcolor(X,Y,I);
set(h, 'LineStyle', 'none');
xlabel('x(m)');
ylabel('y(m)');
colormap(gray);
 
 
figure(2);
size(I);
I(200,200)=0;
hh=pcolor(X,Y,I);
set(hh, 'LineStyle', 'none');
xlabel('salut');
ylabel('salut');
colormap(gray);

Je veux imposer une certaine intensité sur les bords de l'image par exemple
I(100:200,100:200)=0;
Mais cela ne fait rien , il ne me trace pas la 2eme figure enfin c'est tout blanc mais pas de message d'erreur. Je me demandais si on pouvait transformer directement comme ca ?
merci

[FLB]

Salut,
Si tu souhaites ajouter un bord noir :

Code :

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

1
2
3
4
5
6
7
8
figure(2);
[Im, In] = size(I);
noir = min(I);%pour du blanc : max
I(:,1)=noir;
I(1,:)=noir;
I(:,Im)=noir;
I(In,:)=noir;
hh=pcolor(X,Y,I);

imread sert à lire des fichiers images, imagesc à afficher un tableau contenant des données de couleurs.

[Novae]

Merci de ta réponse. En fait par bord noir je me suis mal exprimé je voulais plutôt dire ajuster la largeur de ma fonction. Car par exemple en fait mon intensité c'est des interférences lumineuses dons des sinusoïdes qui prennent toute la figure ( normale il y a une intensité pour chaque valeur de X et Y ) en fait ce que je veux faire c'est rétrécir un peu la figure d’interférence , qu'elle soit juste au centre et qu'on voit juste quelque sinusoïdes. C'est un peu comme si on ajustant la largeur a mi hauteur d'une gaussienne. Je voudrais avoir mes sinusoïdes et sur le reste de la figure du gris

Je vais m’inspirer de ton code. Au départ , j'avais tenté une condition if genre

Code :

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

1
2
3
4
if (X<200) 
    I=0 
else 
    I= f(x,y)

Mais apparemment il aime pas trop le if (X<200) car X est un vecteur

[le fab]

met NaN au lieu de 0

Code :

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

I(100:200,100:200)=NaN;

sinon tu fais quelque chose comme :

Code :

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

pcolor(X(100:200),Y(100:200),I(100:200,100:200))

ou

Code :

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

pcolor(X(100:200,100:200),Y(100:200,100:200),I(100:200,100:200))

en fonction de comment tu as écrit X et Y

[Novae]

Merci
Mais alors le pcolor , il va juste tracer pour l'intervalle de X voulu et l'intensité correspondante. Il ne va pas tracer la vrai fonction pour les autre X

Code :

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

1
2
 
pcolor(X(100:200),Y(100:200),I(100:200,100:200))

En fait je pense qu'on se complique non ? matlab doit avoir avoir un moyen de faire par exemple :

I = sin(x+y) SI X ∈ [100;200] et Y ∈[100;200] 
  = 10 ( chiffre au hasard ) sinon