ddefalco
Questo è un forum di discussione e non un sistema di messaggistica, in particolare con me.
Ciò vuol dire che ognuno che vi partecipa contribuisce alla discussione di un argomento (topic) cercando di risolvere i problemi che si pongono, con un intervento (post) che viene messo a disposizione di tutti: va pertanto inteso come un dibattito tra tutti nel quale gli interventi (e non le persone), hanno lo stesso peso.
Certamente non è una comunicazione ‘peer to peer’ con me (per fare questo esiste il ricevimento in stanza, l’email, la chat…) né tantomeno dev’essere un’esibizione del proprio sapere o non sapere (esame).

Pertanto bandiamo frasi del tipo “Ci scusiamo per eventuali errori”, “accettiamo consigli per eventuali correzioni”, peggio ancora “domanda per il professore”, …, e pensiamo sempre al “sodo”.
Concediamo alla forma, tutt’al più qualche volta, un saluto.
Domenico de Falco.

ddefalco

ING - IND 13 : DDEFALCO's Forum
 
IndiceCalendarioFAQCercaLista utentiGruppiRegistrarsiAccedi
Ultimi argomenti
Cerca
 
 

Risultati secondo:
 
Rechercher Ricerca avanzata

Condividere | 
 

 Esempio numerico pag.92

Vedere l'argomento precedente Vedere l'argomento seguente Andare in basso 
AutoreMessaggio
S.Campanile



Numero di messaggi : 7
Corso di appartenenza : Meccanica Applicata alle Macchine
NomeCognome : Silvio Campanile
Data d'iscrizione : 01.04.11

MessaggioOggetto: Esempio numerico pag.92   Mer Ott 30, 2013 3:26 pm

clear
clc
disp('Si consideri una particella che sia vincolata a muoversi su una retta')
disp('di un piano verticale xy, di equazione 3x+2y=0, sotto l''azione della')
disp('forza di gravità che ha la direzione dell''asse y positivo. Si vuole')
disp('trovare l''equazione del moto della particella')
disp('Condizioni iniziali')
dx=0
dy=0
m=1
disp(' La matrice dei coefficienti della equazione di vincolo è')
syms x y real
d=diff(sym('3*x+2*y'),x);
e=diff(sym('3*x+2*y'),y);
H=[d e];
A=sym2poly(H);
A
disp('-------------- premere un tasto per continuare -------------')
pause;
disp(' La matrice di massa del sistema è')
M=[m 0;0 m]
disp('-------------- premere un tasto per continuare -------------')
pause;
disp(' In virtù del sistema di riferimento adottato, il vettore forza è')
F=[0 9.81]
disp('-------------- premere un tasto per continuare -------------')
pause;
disp(' L''accelerazione del sistema in assenza di vincoli è')
a=[0 9.81]'
disp('-------------- premere un tasto per continuare -------------')
pause;
k=diff(d,x);
l=diff(e,y);
m=diff(d,y);
n=diff(e,x);
disp(' Il termine noto b è invece pari a')
c=-(k*(dx^2)+2*m*dx*dy+l*(dy^2));
b=sym2poly(c)
disp('-------------- premere un tasto per continuare -------------')
pause;
disp(' In questo caso il termine b è uno scalare')
disp(' Determiniamo ora la matrice dei vincoli B,sapendo che B=A*M^(-1/2)')
disp(' La matrice B è pari a')
D=M^(-1/2);
B=A*D
disp('-------------- premere un tasto per continuare -------------')
pause;
disp(' L''inversa generalizzata di Moore-Penrose della matrice B è')
MPB=pinv(B)
disp('-------------- premere un tasto per continuare -------------')
pause;
disp(' L''accelerazione effettiva del sistema è pari a')
ax=a+D*MPB*(b-A*a)


Ultima modifica di S.Campanile il Gio Nov 21, 2013 8:52 pm, modificato 1 volta
Tornare in alto Andare in basso
Vedere il profilo dell'utente
ddefalco
Admin


Numero di messaggi : 171
Località : Napoli
Corso di appartenenza : Meccanica Applicata alle Macchine
NomeCognome : Domenico de Falco
Data d'iscrizione : 29.12.08

MessaggioOggetto: Re: Esempio numerico pag.92   Gio Ott 31, 2013 10:09 pm

Interessante l'uso della matematica simbolica con Matlab di cui conoscevo l'esistenza ma che non mi è familiare.
Qualcuno ha voglia di confrontare i risultati ottenuti con il listato matlab proposto da campanile con quelli ottenuti, con gli setessi valori dei parametri ovviamente, con maple ?
Tornare in alto Andare in basso
Vedere il profilo dell'utente
Mario Guarnieri



Numero di messaggi : 4
Corso di appartenenza : Dinamica e Controllo dei Sistemi Meccanici
NomeCognome : Mario Guarnieri
Data d'iscrizione : 15.10.13

MessaggioOggetto: Re: Esempio numerico pag.92   Dom Nov 10, 2013 10:01 pm

non riesco a capire come mai il listato caricato mi dia il seguente errore:
Error using sym.svd
Too many input arguments.

Error in pinv (line 29)
[U,S,V] = svd(A,0);

Error in pinv (line 27)
X = pinv(A',varargin{:})';

Error in esercizio_nuovo (line 50)
MPB=pinv(B)
Tornare in alto Andare in basso
Vedere il profilo dell'utente
Contenuto sponsorizzato




MessaggioOggetto: Re: Esempio numerico pag.92   

Tornare in alto Andare in basso
 
Esempio numerico pag.92
Vedere l'argomento precedente Vedere l'argomento seguente Tornare in alto 
Pagina 1 di 1
 Argomenti simili
-
» Un esempio di allenamento anaerobico
» gare regolarità
» Esempio di acquisto globale per il completo detailing
» Catena cinetica chiusa o aperta?
» Come eliminare l'adipe senza restrizioni caloriche

Permesso di questo forum:Non puoi rispondere agli argomenti in questo forum
ddefalco :: CORSI :: Dinamica e Controllo di Sistemi Meccanici-
Andare verso: