[Matlab] Aspect ratio grafici

[PepePepe]

Buongiorno, sono un nuovo utente.

Sto lavorando in ambiente matlab ad un codice che fra le altre cose genera un abaco (credo non sia importante specificarlo, ma il contesto è quella della dinamica strutturale, in particolare gli spettri di risposta). Il grafico che devo produrre contiene delle rette inclinate a 45° rispetto l'asse delle ascisse. Entrambi gli assi cartesiani sono in scala logaritmica. Il problema che sto tendando di risolvere è quello di generare un grafico in cui le proporzioni fra i due assi siano tali da visualizzare le suddette rette esattamente a 45° (in altre parole, la singola "decade" deve avere la stessa lunghezza in entrambi gli assi). Al momento infatti il grafico viene generato con un aspect ratio che le fa sembrare inclinate di 60° circa. L'unica soluzione provvisoria che ho trovato è quella di specificare pbaspect a [sqrt(2) 1 1] ma è un valore approssimato trovato per tentativi. Vorrei capire invece se c'è modo di risolvere il problema in modo esatto. Spero di essere stato chiaro.

Grazie anticipatamente

P.S. Porzione di codice che genera il grafico:

Codice: Seleziona tutto

loglog(T0,D,'Linewidth', 0.5,'Color','black','linestyle','--')
grid on
xlabel 'T_0'
xticks([0.02 0.05 0.1 0.2 0.5 1 2 5 10 20 50])
xticklabels({'0.02','pi','-\pi','0','\pi','2\pi','3\pi'})
ylabel 'S_D'
xlim([0.02 50])
ylim([0.2 50])
pbaspect([sqrt(2) 1 1])

[RenzoDF]

Beh, basta che tu vada ad usare per l'asse x, ovvero per il primo parametro numerico di pbaspect il rapporto fra i logaritmi dei rapporti fra i limiti dei due intervalli, ovvero

Codice: Seleziona tutto

r=log10(50/0.02)/log10(50/0.2);
pbaspect([r 1 1])

[PepePepe]

Ottimo, grazie!

[rugweri]

PepePepe, intervengo solo per una piccola proposta: hai avuto la risposta che ti serviva, e va bene... ma hai voglia di ragionare un po' e magari spiegare qui sul forum perché quella risposta è giusta?

Non è un obbligo, nessuno ti impedisce di andartene via come fanno quasi tutti quelli che ottengono la risposta ai propri quesiti... ma se volessi farlo secondo me potrebbe esserti utile

[PepePepe]

Con piacere.

Il comando pbaspect consente di fissare il rapporto fra la lunghezza degli assi in matlab (diciamo ed ), ovvero la lunghezza dell'asse rispetto all'asse secondo la sintassi:

pbascpetc( [r 1 1] )

Volendo ottenere rette visualizzate a 45 gradi, serve che la lunghezza dell'intervallo tra due potenze del dieci consecutive sia uguale per entrambi gli assi o, in altre parole, è necessario che gli assi siano monometrici di modo che la lunghezza degli assi sia rappresentata con il medesimo fattore di scala . Esso è dato dal rapporto fra la lunghezza numerica ("reale") dell'asse e la lunghezza che l'asse assume in visualizzazione grafico di matlab: .

Allora deve essere:



da cui:



That's all!

P.S. Sicuramente c'è un modo più semplice di spiegarlo, ma per deformazione professionale, questo è come lo spiegherei io a qualcuno.