Forum / Grundlagen der Mehrkörpersystemdynamik / Aufgabe 7
Aufgabe 7
Wer kennt sich mit Matlab aus? Habe für den Kreiselkompass eine DGL 2.Ordnung, nichtlinear und habe folgenden Code geschrieben: was ist da falsch? freue mich über jeden Tipp:) function [ dzdt] = myode(t,z)
%% Definitionen
global g m G d r omegaKF I1 I2 I3
g = 9.81; %%m/s^2
m = 0.120; %%kg
G = mg; %%N
d = 0.01; %%m
r = 0.025; %%m
omegaKF = 7000; %%1/min
omegaE = 7.2910^(-5); %%1/s
alpha;
if alpha < 1
end
%%DGL System
I1 = (m/2)(r^2); %%kgm^2 I2 = (m/4r^2)+ (m/12*d^2); I3 = I2;
z1 = z(1,1); z2 = z(2,1);
dz1dt(1)=z2; dz2dt(2)=omegaE^2cos(alpha)^2sin(z1)cos(z1)(I2-I1)-I1omegaEcos(alpha)*sin(z1)omegaKF-kz2; dzdt = [dz1dt ; dz2dt];
%Integrationszeitraum t_span = [0 6];
% Anfangsbedingungen z_0 = zeros(2,1);
% Integration options = odeset('MaxStep',1e-3); odeset; [t,z] = ode45(@myode,t_span,y_0,options);
end