Programmare in matlab

[enricobozzetti]

Buongiorno,
sono alle prime armi con matlab e mi serve un aiuto con il programma. Devo eseguire questo programma:
M*(D2x)+C*(Dx)+K*x(t) = b*H(t)*[x(t-tau)-x(t)]+u(t)
M = m*I
C = 2*m*zeta*o*I
K = m*(o^2)*I
I matrice identica
H(t) = [Hxx(t),Hxy(t);Hyx(t),Hyy(t)]
Hxx(t)= sommatoria di j da 1 a 2(s[phij(t)]*[Kt*cj(t)+Kn*sj(t)]*sj(t)
Hxy(t)= sommatoria di j da 1 a 2(s[phij(t)]*[Kt*cj(t)+Kn*sj(t)]*cj(t)
Hyx(t)= sommatoria di j da 1 a 2(s[phij(t)]*[-Kt*sj(t)+Kn*cj(t)]*sj(t)
Hyy(t)= sommatoria di j da 1 a 2(s[phij(t)]*[-Kt*sj(t)+Kn*cj(t)]*cj(t)
cj(t)= cos(phij(t))
sj(t)= sin(phij(t))
s[phij(t)]= 1 se phie 0 altrimenti
phie=0 phia=1-arccos(1-(a/R))
phij(t)=((2*3.14*OMEGA)/60)*t+((2*3.14*(j-1))/n)
u(t)= -[muxT*g;muyT*g]
mux=dx*Ix*g
muy=dy*Iy*g
e mi vengono dati questi valori:
n Number of teeth 2
m Mass 0.013 kg
Kt Linear tangential cutting coefficients 6  108kg=ðms2Þ
Kn Linear normal cutting coefficients 2  108kg=ðms2Þ
z Damping ratio 0.05
xn Natural frequency 2p  778rad=s
a Radial depth of cut 2.5mm
r Radius of cutter 2.5mm
fz Feed per tooth 0.1mm
b Axial depth of cut Variable
OMEGA Spindle speed in rpm Variable
per avere un'idea più chiara mando un link dove si illustra il problema:
https://www.dropbox.com/s/y8ukgr0d1trip ... 7.pdf?dl=0.
Mi serve per la tesi e se qualcuno può aiutarmi ne sarei grato.
Uso matlab 2011.
Ringrazio in anticipo.

[Linux1987]

Ciao Enrico, non ho capito devi scrivere o eseguire un programma?