Computer implemented method for identifying alterations in a musculoskeletal system of a subject, wherein the method comprises the steps of: • associating, with said subject, non-invasive detecting means (Ri) capable of detecting physical quantities (Qi) representative of a state (St) of said subject, wherein said associating step comprises both the steps of: ∘ positioning said detecting means (Ri) on said subject; and ∘ applying said detecting means (Ri) on said subject. • detecting said physical quantities (Qi) through said detecting means (Ri); • identifying a biomechanical model (MBi) for said subject by an identifying unit; • determining combined movements (MCij) of muscles (Mi) and bones (Oj) as a function of said detected physical quantities (Qi) and said identified biomechanical model (MBi); • converting said combined movements (MCij) into a plurality of movement steps (Fkz) of segments (SEk) and joints (ATz), wherein, in said conversion step, said movement steps (Fkz) are computed as a function of: first joint parameters (Pz) of a first joint (ATz) and at least a second joint parameter (Pz+1) computed as a function of second joints (ATz+1) affected by the movement of the first joint (ATz); wherein, in said conversion step, said second joints (ATz+1) move relative to said first joint (ATz) in a condition between: a movement in synergy with said first joint (ATz); a movement compensating for the movement produced by said first joint (ATz); • providing a processing unit for processing alteration data; • computing parameters (PALi) representative of said alterations as a function of: ∘ said movement steps (Fkz) computed in said conversion step and ∘ said biomechanical model (MBi) identified by said identifying unit, and ∘ said joint parameters (Pz; Pz+1) wherein the computing is carried out by a computing module of said processing unit; • comparing said computed parameters (PALi) and predefined reference values (P REF i) representative of alteration thresholds (S AL) in said musculoskeletal system, wherein the comparing is carried out by a comparison module of said processing unit; • identifying said alterations as a function of a failed matching between said computed parameters (PALi) and said predefined reference values (P REF i), wherein the identifying is carried out by an identification module of said processing unit; • determining a belonging of an alteration to one or more groups of musculoskeletal pathologies as a function of said identified alteration, wherein the determining is carried out by a classification module of said processing unit. • associating, with said subject, non-invasive detecting means (Ri) capable of detecting physical quantities (Qi) representative of a state (St) of said subject, wherein said associating step comprises both the steps of: ∘ positioning said detecting means (Ri) on said subject; and ∘ applying said detecting means (Ri) on said subject. ∘ positioning said detecting means (Ri) on said subject; and ∘ applying said detecting means (Ri) on said subject. • detecting said physical quantities (Qi) through said detecting means (Ri); • identifying a biomechanical model (MBi) for said subject by an identifying unit; • determining combined movements (MCij) of muscles (Mi) and bones (Oj) as a function of said detected physical quantities (Qi) and said identified biomechanical model (MBi); • converting said combined movements (MCij) into a plurality of movement steps (Fkz) of segments (SEk) and joints (ATz), wherein, in said conversion step, said movement steps (Fkz) are computed as a function of: first joint parameters (Pz) of a first joint (ATz) and at least a second joint parameter (Pz+1) computed as a function of second joints (ATz+1) affected by the movement of the first joint (ATz); wherein, in said conversion step, said second joints (ATz+1) move relative to said first joint (ATz) in a condition between: a movement in synergy with said first joint (ATz); a movement compensating for the movement produced by said first joint (ATz); first joint parameters (Pz) of a first joint (ATz) and at least a second joint parameter (Pz+1) computed as a function of second joints (ATz+1) affected by the movement of the first joint (ATz); a movement in synergy with said first joint (ATz); a movement compensating for the movement produced by said first joint (ATz); • providing a processing unit for processing alteration data; • computing parameters (PALi) representative of said alterations as a function of: ∘ said movement steps (Fkz) computed in said conversion step and ∘ said biomechanical model (MBi) identified by said identifying unit, and ∘ said joint parameters (Pz; Pz+1) wherein the computing is carried out by a computing module of said processing unit; ∘ said movement steps (Fkz) computed in said conversion step and ∘ said biomechanical model (MBi) identified by said identifying unit, and ∘ said joint parameters (Pz; Pz+1) • comparing said computed parameters (PALi) and predefined reference values (P REF i) representative of alteration thresholds (S AL) in said musculoskeletal system, wherein the comparing is carried out by a comparison module of said processing unit; • identifying said alterations as a function of a failed matching between said computed parameters (PALi) and said predefined reference values (P REF i), wherein the identifying is carried out by an identification module of said processing unit; • determining a belonging of an alteration to one or more groups of musculoskeletal pathologies as a function of said identified alteration, wherein the determining is carried out by a classification module of said processing unit.