Hej!
Jag behöver hjälp, har räknat ut en uppgift i mekanik för hand. Nu vill jag köra ett program i matlab och få fram så att om jag stoppar in mina värden i matlab så ska jag få fram samma resultat.
Jag laddar upp bilder och min matlab kod här så får ni gärna kolla och hjälpa mig
http://www.ladda-upp.se/bilder/nznsjtzavtarvz/
http://www.ladda-upp.se/bilder/xrnytzryjekqaf/
här är matlab koden:
function yyy=uppgA21pde(beta0,v0,tend,tstep);
tv=[0:tstep:tend];
lv=length(tv)
betanu=beta0;
vnu=v0;
betav(1)=beta0;
vv(1)=v0;
dt=tstep;
for ii=1:lv-1,
tnu=ii*dt;
kk=-41.244*sin(betanu) + sin(betanu+0.1)*cos(betanu);
dv=kk*dt;
vny=vnu+dv;
vv(ii+1)=vny;
betany=betanu+(vny+vnu)/2*dt;
betav(ii+1) = betany;
vnu=vny;
betanu=betany;
end;
figure(1);
plot(tv,vv);
title('Hastighet som funktion av tid');
figure(2);
plot(tv,betav);
title('Vinkel som funktion av tid');
figure(3);
plot(betav,vv);
title('Hastighet som funktion av vinkeln');
yyy=vv;
och som koden säger har jag plottat vissa grafer här kommer dem:
http://www.ladda-upp.se/bilder/clmmbwqdhvmdv/
Där finns 3 frågor vi skulle svara på och här är ytterligare 2 frågor jag vill svara på:
1. Kan ni ge exempel på värden så man ser att formeln stämmer med det numeriska ?
2. Hur kommer det sig att både vinkelhastigheten och vinkeln beta ökar ?
Tack för all hjälp, uppskattas!