Si hay algun iluminado por ahi, o un buen samaritano(perdonad mis explicaciones pero es complicado explicar matematicas cuando siempre las has dado en ingles).
EStoy con una parte del ejercicio que me esta tocando las narices.
dt,N1,N2,N3, h son cifras que te dan.
tengo que calcular los valores de u en ese grafico. Resultado de un vector t.
Pongo aqui el codigo incluida la parte donde calculo la transfer function.
function [G1,G2,G3,time,u,x1,x2,f ] = m286_r7_suspension_200740111(m1,m2,c1,c2,k1,k2,dt, N1,N2,N3,h )
D1=[m2, c1+c2, k1+k2];
D2=[m1, c1, k1];
D3=[0, c1, k1];
Delta1=conv([D1],[D2]);
Delta2=conv([D3],[D3])
Delta= Delta1+Delta2
Num1=conv([c1,k1],[c2,k2]);
Num2=conv([m1,c1,k1],[c2,k2]);
Num3a=conv([c1,k1],[c2,k2]);
Num3=conv([Num3a],[m1,0,0]);
G1=tf(Num1,Delta);
G2=tf(Num2,Delta);
G3=tf(Num3,Delta);
tmax=N3*dt; %setting up tmax
time= 0:dt:tmax;
LT=length(time);
u=zeros(LT,1);
N1dt=N1*dt;
N2dt=N2*dt;
Grad=(h)/(dt*(N2-N1));
for z=1:1:LT
t=time(z);
if t>N1dt & t<N2dt %SEGMENT NEEDS CHECKING
u(z)=t*Grad
if t>=N2dt
u(z)=h;
else u(z)=0;
end
end
end
a lo de X1, x2,f aun no he llegado...
la parte en rojo es la que me esta dando problemas, me calcula lo que quiero, pero en vez de almacenarlo lo va borrando poco a poco. Seguramente habran mas cosas que no estaran muy bien hechas, una vez aprueba esto nunca mais. Cualquier buen samaritano sera recompensado con eterna gratitud y muchos me gusta
upeo
Permisos de publicación
Citar