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

Sorry, die Kommentarfunktion ist geschlossen.