scilab hängt sich auf

WICHTIG: Der Betrieb von goMatlab.de wird privat finanziert fortgesetzt. - Mehr Infos...

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

ms22

Forum-Newbie


	Beiträge: 5

	Anmeldedatum: 07.08.08

	Wohnort: ---

	Version: ---

Verfasst am: 11.03.2009, 08:51 Titel: scilab hängt sich auf

Hallo,

habe das Problem, dass sich Scilab aufhängt, wenn ich den unten geschriebenen Code zum zweiten mal ausführe (also nachdem ich einmal "Execute in Scilab" gemacht hab und den code einmal komplett durchlaufen liess) . Meistens bringt mir Scilab dann die Fehlermeldung "Scilab may become unstable" und danach geht nichts mehr. woran könnte das liegen?

Code:

vb = 0.833;
nmmax = 3000.0;
ri = 0.25;
ra = 0.9;
r = ri:0.01:ra;
uemin = 30;
uemax = 150;
ue = uemin:(uemax-uemin)/(length(r)-1):uemax;
jg = 0.05;
jm = 0.23;
etag = 0.9;
a = 8.3;
mzug = 75.87;
stahllaenge = 0.27;
stahldichte = 7.85;

imax2 = (2*%pi*ri)*(nmmax/(vb*60));
jakt = 0.5*%pi*stahldichte*1000*stahllaenge*(r.^4-ri^4);

//mmg1 = (a.*ue/r).*(jm+jg);
//mj1 = (a.*ue/r)*(jakt ./ (ue^2*etag));
//m1 = mmg1 + mj1- mzug;

xgrid();
plot(r, jakt);
pause
mmg1 = (a.*ue./r).*(jm+jg);
mj1 = (a.*ue./r).*(jakt ./ (ue.^2*etag));
m1 = mmg1 + mj1- mzug;

length(m1)
length(mj1)

moment1 = (ue*a/ri)*(jm+jg) - mzug;
moment2 = (ue*a/ra)*(jm+jg) + (ue*a/ra).*(0.5*%pi*stahldichte*1000*stahllaenge*(ra^4-ri^4)./(ue.^2*etag))-mzug;

length(ue)
length(moment2)
length(moment1)

clf
subplot(141)
xgrid()
plot(ue, [moment1' moment2']);
legend('moment1', 'moment2');

moment10 = -(ue*a/ri)*(jm+jg)-mzug;
moment20 = -((ue*a/ra)*(jm+jg) + (ue*a/ra).*(0.5*%pi*stahldichte*1000*stahllaenge*(ra^4-ri^4)*etag./(ue.^2)))-mzug;

subplot(142)
xgrid()
plot(ue, [moment10' moment20']);
legend('moment10', 'moment20');

mmg2 = -(a.*ue/r)*(jm+jg);
mj2 = -(a.*ue/r).*(jakt.*etag./(ue^2));
m2 = mmg2+mj2-mzug;
//
//plot3d(ue,r,m1)

iopt = sqrt(jakt./(etag*(jg+jm)));
subplot(143)
xgrid()
xtitle('iopt = f(r)')
plot(r,iopt);

subplot(144)
xtitle('Momente = f(ue)')
xgrid()
plot(ue, [moment1' moment2' -moment10' -moment20']);
legend('moment1','moment2','moment10', 'moment20');

pause

clf
//plot3d(r',ue',m1')
xgrid
[R, UE] = meshgrid(r,ue);
M1 = (a*UE./R)*(jm+jg) + (a*UE./R)*(0.5*%pi*stahldichte*1000*stahllaenge*(R.^4-ri^4) ./ (UE.^2*etag)) - mzug;
//color = rand()*ones(1,length(M1));
surf(R,UE,M1)

maxbeschlg= max(m1)
maxBeschlg= max(M1)
minbeschlg = min(m1)
minBeschlg = min(M1)
pause;

M2 = -(a.*UE./R)*(jm+jg) - (a.*UE./R).*(0.5*%pi*stahldichte*1000*stahllaenge.*(R.^4-ri^4)*etag./(UE^2)) - mzug;
clf
xgrid
surf(R,UE,M2)

maxBrems= max(M2)
minBrems = min(M2)

Funktion ohne Link?

[/code]

danke schon mal

Verschoben: 11.03.2009, 09:55 Uhr von Martin
Von Grafiken und Plots nach Scilab / Xcos (ehemals Scicos) Forum


Einstellungen und Berechtigungen
Beiträge der letzten Zeit anzeigen:	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

Impressum | Nutzungsbedingungen | Datenschutz | FAQ |

RSS

Hosted by:

Copyright © 2007 - 2024 goMatlab.de | Dies ist keine offizielle Website der Firma The Mathworks

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.