logo

Crowdly

Un oscil·lador  x'' + \omega^2 x = F(t) està sotmès a la força periò...

✅ The verified answer to this question is available below. Our community-reviewed solutions help you understand the material better.

Un oscil·lador 

 x'' + \omega^2 x = F(t) x'' + \omega^2 x = F(t)

està sotmès a la força periòdica F(t) donada per

n=0, 1, 2, 3...

a) Defineix adequadament les variables y1, y2, i escriu el sistema d'edos de primer ordre equivalent a l'equació de segon ordre (les respostes són cadenes de caràcters).

y1 ≡ y1' =  
y2 ≡ ;   y2' = 

 

b) Integra el sistema numèricament per a   ω=1, amb les condicions inicials x(0)=x'(0) =0. Emplena amb 4 xifres significatives, els valors de la taula següent. 

Utilitza RK2 amb un pas d'integració de 10-3.

 t=1t=5t=10t=20t=100
x(t)
x'(t)

 

És important, per definir la força F(t), referir el temps t al primer periode (n=0), ja que els valors de F(t) es van repetint a intervals de 1 s. 

Si no, caldria fer un if else molt llarg, amb un codi llarguíssim i inútil. Per contra, tot això es pot fer d'una forma eficient amb un únic IF / ELSE, si utilitzes la funció rem de MATLAB:  rem(t,1) dóna el residu de dividir t entre 1.  

Per exemple, si t=3.26: rem(t,1)=0.26. 

I si t=3.999: rem(t,1)=0.999. 

Prova-ho. Aleshores, rem(t,1) és un valor que queda sempre entre 0 i 1, i amb això pots reduir el valor de F(t) en qualsevol interval n al valor que pren al primer periode (n=0). Avalua amb un if la condició de quant val rem(t,1), si menys de 0.5 , o més de 0.5, i assigna-li el valor de F en cada cas.

More questions like this

Want instant access to all verified answers on campusvirtual.urv.cat?

Get Unlimited Answers To Exam Questions - Install Crowdly Extension Now!