Qcno conosce MATLAB?

[Uzi[WNCT]]

Chi ha voglia/tempo di aiutare un povero studente di informatica che deve fare l'esame di calcolo numerico ma il prof nn ha tempo di spiegare matlab??

Se ho risposte affermative () dopo posto delle domandine

[soalle]

http://www.frattali.it/matlab/intromatlab.html

Qualcosa so...

[ottobrerosso]

Originariamente inviato da Uzi[WNCT]
Chi ha voglia/tempo di aiutare un povero studente di informatica che deve fare l'esame di calcolo numerico ma il prof nn ha tempo di spiegare matlab??

Se ho risposte affermative () dopo posto delle domandine

il prof 'non ha tempo' ? ma non c'è un esercitatore o meglio una esercitatrice ?

[Urizen]

scusa ... ma dove studi informatica ????? Venezia ???? Padova ???

[Uzi[WNCT]]

C'è l'esercitatore ma lui "confida nel fatto che siamo informatici e nn abbiam difficoltà ad imparare un nuovo linguaggio ".

Cmq, chi ha idea di come risolvere sta roba?

date due matrici A e B di dimensione n = 2500, calcoli il prodotto "per componenti" delle prime n-1 righe di A e di B, ne calcoli il tempo di esecuzione, e poi faccia lo stesso con le prime n-1 colonne

lasciam perdere il 2500 di dimensione che ci mette na cifra anche col mio xp1800, tanto mi basta sapere come si implementa la cosa.

Ho fatto così, ma nn sono sicuro che sia giusto:

max_n = 2500;

A = rand(max_n,max_n);
B = rand(max_n,max_n);


% per colonne
C=zeros(max_n,max_n);
tic
for i=1:max_n
  for j=1:max_n
    for k=1:max_n
      C(i,j)=C(i,j)+A(i,k)*B(k,j);
    end
  end
end
tempo1 = toc

% per righe
C=zeros(max_n,max_n);
tic
for i=1:max_n
  for j=1:max_n
    for k=1:max_n
      C(i,j)=C(i,j)+A(j,k)*B(k,i);
    end
  end
end
tempo2 = toc

[Sephiroth83]

Ahhhhhhh, anch'io ho l'esame di calcolo numerico (lunedì) e domani l'esonero! come ti capisco!!!
Ecco il sito del corso che frequento, c'è un manuale in pdf e qualche esercitazione:
http://calvino.polito.it/~scuderi/INF03.html
Ho una domanda: cosa intendi precisamente per "prodotto per componenti".Centra qualcosa il prodotto puntuale (quello fatto tra elementi con le stesse 'coordinate')?

[Uzi[WNCT]]

ecco...... me lo domando anche io... il prof ha scritto così.... "per componenti". (cmq credo sia quello)

[Sephiroth83]

Mi sembra che tu abbia implementato il prodotto normale tra matrici. Io farei qualcosa del genere:

% per righe
C=zeros(max_n - 1,max_n);
tic
for i=1 : max_n - 1
for j=1:max_n
C(i, j) = A(i, j) * B(i, j);
end
end

In più aggiungi tic e toc come hai fatto tu, io non mi sono dilungato. Conta che questo codice l'ho scrito in due minuti senza testarlo, quidni prob non funge. Domani (dopo l'esonero) lo provo e ti faccio sapere.

[Uzi[WNCT]]

in effetti è vero...... è il prodotto normale.... mo guarderò....

[cipi]

A parte il fatto che, secondo me, il tuo prof. non ha le idee molto chiare sui prodotti matriciali... io risolverei così:

clear
clc
n_max=3;
A=rand(n_max,n_max);
B=rand(n_max,n_max);
tic
for i=1: (n_max-1),
A( :,i ).*B( :,i ); % per colonne
end
tempo1=toc
for i=1: (n_max-1),
A( i,: ).*B( 1,: ); % per righe
end
tempo2=toc

ciao

[Uzi[WNCT]]

Sinceramente credo che il programma nn fosse per niente riguardo alle matrici, ma solo moltiplicazione "numero per numero" solo che nn è chiaro.......

[Sephiroth83]

Non basterebbe semplicimente fare così?:

n_max=3;
A=rand(n_max,n_max);
B=rand(n_max,n_max);

tic
A( :, 1 : n_max-1 ).*B( :, 1 : n_max-1); % per colonne

tempo1=toc
A(1 : n_max-1,.*B(1 : n_max-1,: ); % per righe
tempo2=toc

[Uzi[WNCT]]

(1 : n_max-1,: )

si può fare una cosa del genere?!?! FIGATAAAA!!!

Cmq, bene o male quell'esercizio è stato risolto, adesso ce ne sono altri, assolutamente improponibili

[cipi]

Originariamente inviato da Uzi[WNCT]
(1 : n_max-1,: )

si può fare una cosa del genere?!?! FIGATAAAA!!!

Cmq, bene o male quell'esercizio è stato risolto, adesso ce ne sono altri, assolutamente improponibili

proponi, proponi...
cmq puoi ridurlo a:
n_max=3;
A=rand(n_max); % crea gia' matrici quadrate...
B=rand(n_max);

tic
A( :, 1 : n_max-1 ) .* B( :, 1 : n_max-1); % per colonne

tempo1=toc
A(1 : n_max-1, : ) .* B(1 : n_max-1, : ); % per righe
tempo2=toc

.... ma è meno leggibile per un principiante!

[bsummer]

Originariamente inviato da Uzi[WNCT]
(1 : n_max-1,: )

si può fare una cosa del genere?!?! FIGATAAAA!!!

Si, ed è una buona regola in matlab cercare di evitare i cicli for allo scopo di non invecchiare nell'attesa che vengano effettuali i calcoli...
Solitamente per i calcoli sulle matrici, matlab fornisce soluzioni che non utilizzano cicli for (dipende dai casi, cmq).

Cmq, bene o male quell'esercizio è stato risolto, adesso ce ne sono altri, assolutamente improponibili

Beh, tu postali pure che non hai nulla da perdere, anzi come in questo caso c'è chi potrebbe aiutarti

Aloha!

[cipi]

Originariamente inviato da bsummer
Si, ed è una buona regola in matlab cercare di evitare i cicli for allo scopo di non invecchiare nell'attesa che vengano effettuali i calcoli...
Solitamente per i calcoli sulle matrici, matlab fornisce soluzioni che non utilizzano cicli for (dipende dai casi, cmq).
Aloha!

Esatto bsummer,
fai conto che un o script di un algoritmo lanciato con matlab con un sacco di for su matrici ci stava circa 40 minuti a dare un risultato... Fatto in C lo stesso algoritmo.... 30 secondi
quindi... okkio alle lungaggini Uzi[WNCT]

[Uzi[WNCT]]

Ragazzi...... ve la siete cercata

Cmq, nn posso postarvi solo il test dell'esercizio altrimenti nn capireste nulla, quindi vi devo dare una intera pagina!

Per quanto riguarda il for, dopo sistemo il mio programma e metto come avete fatto voi, che mi pare moooooooolto + "profeeeescional"

La pagina? Eccola qua!!!

http://www.math.unipd.it/~marcuzzi/D...rmatica_3.html

Alla fine in basso c'è il testo dell' esercizio

[Uzi[WNCT]]

Ritorno a scocciarvi col Matlab

Come posso plottare in 3D una matrice del genere?

risultato =

    0.3333    0.0005    0.0001    0.0000
    0.3333    0.0033    0.0008    0.0003
    0.3333    0.0092    0.0023    0.0009
    0.3333    0.0173    0.0046    0.0019
    0.3333    0.0267    0.0075    0.0032
    0.3333    0.0365    0.0108    0.0047
    0.3333    0.0463    0.0143    0.0064
    0.3333    0.0558    0.0178    0.0082
    0.3333    0.0651    0.0214    0.0100
    0.3333    0.0741    0.0250    0.0118

Vorrei che:
- in X ci fossero valori 1:10 (cioè gli indici delle righe)
- in Y ci fossero valori 1:4 (cioè gli indici delle colonne)
- in Z ci fossero i valori che si vedono

C'è da usare plot3? Mesh?

[Uzi[WNCT]]

UP!

[Uzi[WNCT]]

Lasciate pure perdere le altre domande, che più o meno ho sistemato.



Ore però ne ho una "semplice"

Se io ho una funzione 3D "complicata", come faccio ad avere le coordinate di un punto Z, conoscendo X e Y??