SISTEMA DINAMICO ARX IN MATLAB

[Benjamin]

Vorrei graficarmi col matlab il seguente sistema dinamico:

y(k+1)=0.3cos(y(k))+u(k)

E' un sistema ARX cioè autoregressivo esogeno.
Autoregressivo in quanto l'uscita dipende dall' uscita negli istanti precedenti (y(k+1)=......y(k)....)
Esogeno in quanto è presente un ingresso (u(k)).

In particolare come u(k) ho messo un rumore bianco gaussiano.

k rappresenta gli istanti, dunque un vettore di 100 elementi.

Essendo un sistema più complesso rispetto a una semplice funzione a una variabile, non riesco a capire cosa potrei usare del matlab per poterlo graficare. ES. come implemento il rumore bianco gaussiano??

SPERO DI ESSERE STATO ABBASTANZA CHIARO NELL' ESPOSIZIONE E CHE QUALCUNO SIA IN GRADO DI DARMI UNA MANO

[ext33n]

Una sola parola: simulink

[Benjamin]

non occorre simulink.

Occorre la riga di comando

[rini]

puoi generare un vettore di numero causali che seguono tale distribuzione. Una funzione che potresti utilizzare e la funzione randn del matlab da riga di comando che secondo la onorevole fonte rappresenta la distribuzione normale.
Link: http://www.mathworks.com/access/helpdes ... randn.html

Codice: Seleziona tutto

u = randn(1,100)

[Benjamin]

Ok, fin qua ci siamo.

Ma per il resto....

[rini]

non ho il matlab sotto mano ma potresti fare qualcosa del genre:

Codice: Seleziona tutto

u = randn(1,100);
y=1:1:100;

for i=1:100
  y(i+1)=0.3*cos(y(i)+u(i));
end

plot(y)


salvo errori di sintassi vedi se era quello che ti aspettavi!

[Benjamin]

Innanzitutto ti ringrazio dell' aiuto che mi stai dando, visto che sei stato l' unico a rispondere da quasi un mese a questa parte, GRAZIE.

Forse hai sbagliato perché hai messo u(i) dentro il coseno, invece dovrebbe essere fuori...

La scrittura contratta 1:1:100 cosa vuol dire, non capisco

[rini]

Si ho sbagliato il u(i)

Codice: Seleziona tutto

y=1:1:100

vuol dire che ti crea un vettore di cento elementi da 1 a 100 con distanza tra un elemento è l'altro di 1 unita
in pratica equivale a questo:

Codice: Seleziona tutto

y = [1 2 3 4 ...... 20 21 22 23 ......  98 99 100]

quando non capisci qualcosa la puoi provare direttamente sul prompt del matlab.
prova a scrivere queste due comandi

Codice: Seleziona tutto

y1 = 1:1:10
y2 = 1:0.1:10


i due elementi equivalgono a

Codice: Seleziona tutto

y1 = [1 2 3 4 5 6 7 8 9 10]
y2 = [1 1.1 1.2 1.3 ....... 9.7 9.8 9.9 10]


ho scritto y = 1:1:100 per prendere il posto del tuo vettore che hai già a disposizione. Io l'ho scritto per poter eseguire le 4 righe di codice di sopra

[Benjamin]

Ok, adesso lo provo direttamente su MATLAB e vedo se funziona.

Mi raccomando non abbandonarmi