Sistema oscilatório no scilab

Universidade Federal do Rio Grande do Norte
Escola de Ciências e Tecnologia
Computação Numerica
Aluno: Felipe Diogo Silva de Medeiros

Introdução

Neste laboratorio tivemos por objetivo aprender a utilizar ferramentas como o software silab para resoluções de problemas fisicos, no caso, um sistema oscilatorio, suspenção e amortecedor de um veiculo. Com uma modelagem pre estabelica e disponibilaza.

DESENVOLVIMENTO

Para a resolução do problema foi utilizado o software livre scilab: 01) descrevemos a equeção x(t) em função de k e t, dado m=1,2x10^3kg, t=0,05s, c=10^4kg/s . x(t) = Aexp(−λt) cos(μt) + Bexp(-λt) sin(μt), λ= c / 2m, μ =√|c² − 4mk| / 2m logo x(t) = A*exp(−4,16t)*cos((√|(10^4)² − 4*(1,2*10^3)k| / 2*(1,2*10^3))*t) + B*exp(-4,16t)*sin((√|(10^4)² − 4*(1,2*10^3)k| / 2*(1,2*10^3))*t). 02) tendo as informações montamos a equação f(k)=0. x(t) = A*exp(−0,208)*cos((√|(10^4)² − 4*(1,2*10^3)k| / 2*(1,2*10^3))*0,05) + B*exp(-0,208)*sin((√|(10^4)² − 4(1,2*10^3)k| / 2*(1,2*10^3))*0,05)=0. 03) foi plotado um grafico o grafico f(k) no silab. Linha de comandos: t=0.05;m=1200;c=10^4 c y=0.1;A=y;a=c/(2*m); for i=1:100 k=i*10^6; x(i)=A*(%e^(-la*t))*cos(t*(sqrt(abs(c^2-4*m*k))/(2*m)))+((-y*c)/(sqrt(abs(c^2-4*m*k))))*(%e^(-la*t))*sin(t*(sqrt(abs(c^2-4*m*k))/(2*m))); end k=[10^6:10^6:10^8]; plot2d(k,x) 4) Achando os valores de K. Linha de comandos: t=0.05; m=1200;c=10^4;y=0.1;A=y;la=c/(2*m);

function x=f(k) t=0.05; x=A*(%e^(-la*t))*cos(t*(sqrt(abs(c^2-4*m*k))/(2*m)))+((-y*c)/(sqrt(abs(c^2-4*m*k))))*(%e^(-la*t))*sin(t*(sqrt(abs(c^2-4*m*k))/(2*m))); endfunction a=10^6;b=5*10^6;z=(a+b)/2; iter=0; E=1.00; while E>=10^(-2) z0=z; z=(a+b)/2; E=abs((z-z0)/z); if