Verfasst am: 18.06.2014, 18:33
Titel: Vorliegender Quellcode in Octave zum Laufen bringen
Guten Abend an alle.
Folgendes "Problem":
Ich habe einen Matlab/Octave Quellcode vorliegen, welchen ich jedoch nicht in Octave zum laufen bekomme. Ich habe bislang keinen Kontakt zu beiden Programmen gehabt, was das ganze wohl erklärt....
Den nachfolgenden Quellcode habe ich bereits mehrmals in Octave eingebunden (Copy and Paste), jedoch scheine ich irgend etwas grundlegendes zu übersehen oder falsch zu machen, da bei drücken auf "Speichern und Ausführen" nichts passiert. Ich vermute das es mit der fehlenden Eingabedatei zu tun hat. Wie man diese erstellt, wenn dem so ist, weiß ich leider auch nicht.
Das ganze ist für eine Prüfungsleistung, worin es lediglich um das Einbinden in Octave und das Nachvollziehen der Kommandos geht. Das wäre dann auch schon das nächste Problem.
% Durchfluss bei n = 0 if dpzulVariante == 1
qp= sqrt(pin/(cvzul+rho/8/A^2));
else
c = 128*R*A^2*kinVis/(dzul^2*Azul);
qp= -c/2+sqrt((c/2)^2+8/rho*pin*A^2);
end
va = qp/2/A;
vr = va;
vrt= va*cosb*cosa;
vat= vrt;
m = rho*qp*R*vat;
if m>mt % Haftreibung überschritten?
% Iteration bis zum Momentengleichgewicht whileabs(mdiff) > MaxFehler*abs(m)
n = omega/2/pi;
msoll= mt+n*mf;
varm = omega*R;
% Verlustleistung in Druck umgerechnet
deltap= abs(msoll*omega)/qp;
if dpzulVariante == 1
dpzul= cvzul*qp^2;
else
vzul = qp/2/Azul; % qp/2 durch einen Arm
Re = abs(dzul*vzul/kinVis);
dpzul= 64/Re*rho/2*R/dzul*vzul^2;
end
it=it+1;
if it > 10000 fprintf('Fehler: it > 10000\n');
it,msoll,mdiff,ominc,alpha,beta
A,d,mt,mf,vr,va,vrt,varm,vat,omega
break;
end if(omega < 0.0062) fprintf('Fehler: omega < 0.0062\n');
it,omega
n=0;
break;
end end else
omega=0; n=0;
end
% Flugbahn
dtropf = sqrt(4*A/pi);
etaluft= 1.82e-5;
nyluft = etaluft/1.25;
v = va;
z = 1.0e-3;
sfl = 0.0;
vh = va*cosa;
vv = va*sina;
deltat = 0.005;
mtr = pi/6*dtropf^3*rho;
while z > 0 if(va<0.01) break;
end
Re = va*dtropf/nyluft;
% Abraham, The Physics of Fluids 13, S.2194
zeta= 24/Re*(1+0.11*sqrt(Re))^2;
Fwid= 1.25/2*va^2*pi/4*dtropf^2*zeta;
atr = Fwid/mtr;
sfl = sfl+vh*deltat;
z = z+vv*deltat;
vh = vh-atr*cosa*deltat;
vv = vv-g*deltat-atr*sina*deltat;
va = sqrt(vh^2+vv^2);
cosa= vh/va;
sina= vv/va;
end;
if sflVariante == 1
sfl=sfl;
else
sfl=sqrt((R+sinb*sfl)^2+(cosb*sfl)^2);
end
qp = 2*vr*A*60000;
pverh= deltap/(rho*g*h);
% Hilfsfunktionen function Value=Norm0(MinVal,MaxVal,data,row,col)% ohne Kodierung ifsize(data,2) < col || size(data,1) < row
Value = (MinVal+MaxVal)/2;
else
Value = data(row,col);
end function Value=Norm1(MinVal,MaxVal,data,row,col)% [-1;1] ifsize(data,2) < col || size(data,1) < row
Value = (MinVal+MaxVal)/2;
else
Value = MinVal+(MaxVal-MinVal)*(data(row,col)+1)/2;
end function Value=Norm2(MinVal,MaxVal,data,row,col)% [0;1] ifsize(data,2) < col || size(data,1) < row
Value = (MinVal+MaxVal)/2;
else
Value = MinVal+(MaxVal-MinVal)*data(row,col);
end function Value=Norm3(MinVal,MaxVal,data,row,col)% [1,2,...,ns] ifsize(data,2) < col || size(data,1) < row
Value = (MinVal+MaxVal)/2;
else
minStufe=min(data(:,col));data(:,col)=data(:,col)-minStufe+1;
maxStufe=max(data(:,col));diffStufe=maxStufe-1;
if diffStufe == 0
Value=(MinVal+MaxVal)/2;
else
Value=MinVal+(MaxVal-MinVal)*(data(row,col)-1)/diffStufe;
end end function[spalte,minus,plus]=SetD(Spalte,Minimal,Maximal)
spalte=Spalte;minus=Minimal;plus=Maximal; % Set Factor Data[/size]
Anhang A - Berechnungsmodell zum Fallbeispiel Rasensprenger.pdf
Dateigröße:
203.53 KB
Heruntergeladen:
4132 mal
Einstellungen und Berechtigungen
Du kannst Beiträge in dieses Forum schreiben. Du kannst auf Beiträge in diesem Forum antworten. Du kannst deine Beiträge in diesem Forum nicht bearbeiten. Du kannst deine Beiträge in diesem Forum nicht löschen. Du kannst an Umfragen in diesem Forum nicht mitmachen. Du kannst Dateien in diesem Forum posten Du kannst Dateien in diesem Forum herunterladen
MATLAB, Simulink, Stateflow, Handle Graphics, Real-Time Workshop, SimBiology, SimHydraulics, SimEvents, and xPC TargetBox are registered trademarks and The MathWorks, the L-shaped membrane logo, and Embedded MATLAB are trademarks of The MathWorks, Inc.