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

Mein MATLAB Forum - goMatlab.de

Mein MATLAB Forum

 
Gast > Registrieren       Autologin?   

Partner:




Forum
      Option
[Erweitert]
  • Diese Seite per Mail weiterempfehlen
     


Gehe zu:  
Neues Thema eröffnen Neue Antwort erstellen

Ausführung plot dauert unglaublich lange und eps fehlerhaft

 

TheFinger
Forum-Anfänger

Forum-Anfänger


Beiträge: 28
Anmeldedatum: 14.05.11
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 26.05.2011, 22:10     Titel: Ausführung plot dauert unglaublich lange und eps fehlerhaft
  Antworten mit Zitat      
Hallo People,
folgendes Problem: ich rechne ein paar werte aus und dann will ich sie plotten, dannn hängt sich das Programm auf und mein Rechner ist (i3 cpu) überfordert. Nach ein paar Minuten kommt zwar eine figure aber leider in sehr schechter Qualität und die Legende ist falsch übernommen. Das Bild ist auch riesig im vergleich zu ähnlichen diagrammen (200KB): 1.7 MB
So ein Bild habe ich schon mal gemacht ohne Probleme zu bekommen, jetzt klappt es leider nicht. Ich wäre euch dankbar wenn jemdand mir helfen kann, ich muss nämlich meine Arbeit abgeben Sad
Code:

fpath='G:/Backup2/Simulations/DIFFUSIVE_case_depredurand_castanet_oldmodel/cdc2/c4.5/'
figbasename='cdc2_c4.5'

area_data=importdata([fpath 'area00.his']);
center_data=importdata([fpath 'center00.his']);
mass_data=importdata([fpath 'mass00.his']);
temper_data=importdata([fpath 'temper00.his']);

xkorr(:,1)=center_data(:,14);        % xkorr (=acceleration*delta_t)
a_proj(:,1)=area_data(:,2)/4;       % projected area, only valid for spheres
t=center_data(:,1);                 % time of simulation
m_l(:,1)=mass_data(:,2);         % total mass of liquid phase
rhog=5.421E-4;
rhol=0.78465;
mug=0.00032261;
C=4.5;
% Droplet velocity
for i=1:size(center_data)
v(i,1)=center_data(i,11)-1022.5;           % Droplet velocity in x-direction (-> uleft)    
end
vv=abs(v);
%COMPUTE DV/DT
dt=zeros(size(xkorr,1),1);
for i=2:size(xkorr,1)
    dt(i)=t(i)-t(i-1);
end
vdot=xkorr./dt;
%compute deceleration
n=size(v,1);
for i=1:n-1
    a(i)=(v(i+1)-v(i))/dt(i);
end
a=transpose(a);
%COMPUTE C_W
% pp=polyfit(t,v,2);
% vp=polyval(pp,t);
% vdot=abs(pp(1));

%FS3D-Computation a(1:end,1).
cw=2*m_l(1:end,1)*37.5E3/rhog./a_proj(1:end,1)./(v(1:end,1).^2);


%Reynolds number of the gaseous phase
Vol=mass_data(:,2)/rhol;
d= (6*Vol/pi).^(1/3);
d0= (6*Vol(1)/pi)^(1/3);
Re=abs(v(1:end,1)).*d.*rhog./mug;

%Virepinte Correlation
cd_V=53.2.*(Re.^(-1.524)*C^0.6);


% Cd_iso according to Mulholland et al.
cd_iso=24./Re.*(1+0.1*Re).^0.687;
rCd_iso=1-0.86*exp(-0.053*C);
cd_M=rCd_iso.*cd_iso;

%Cw according to von Brauer for Re<80 for isolated droplet
cw_B=24./Re+5.48./Re.^(0.573)+0.36;

%According to Sirignano Cd/Cdiso
Dv=4.3805E-001
Lv=5.3821E9;
cpl=2.137E7;
cpv=1.66121E7;
mdot=mass_data(2:end,13);
Tp=polyfit(temper_data(1:end,1),temper_data(1:end,2),3);
Tfit=polyval(Tp,temper_data(1:end,1));
dTmdt=(Tfit(2:end)-Tfit(1:end-1))./(temper_data(2:end,1)-temper_data(1:end-1,1));
Ts=temper_data(1:end,4); % Surface temperature
Tm=temper_data(1:end,2);% Averaged droplet temperature
N=size(dTmdt,1);
Ql=m_l(1:N-1).*cpl.*dTmdt(1:N-1)+cpl.*mdot(1:N-1).*(Ts(1:N-1)-Tm(1:N-1));
Tamb=temper_data(1:end,3);
BT=cpv*(Tamb(1:N-1)-Ts(1:N-1))./(Lv+Ql(1:N-1)./mdot(1:N-1));

s=size(BT);
Rcd_iso=0.549.*Re(1:s,1).^0.098.*(1+BT).^0.132*Dv^0.275;
Cd=Rcd_iso.*cw_B(1:s,1);

%Atthasit Correlation according to paper Atthasit et al
a= 3.43E-5*Re.^3-1.359E-3*Re.^2+1.8753E-2*Re+0.0484016;
b= -1.07E-5*Re.^3+2.106E-4*Re.^2+7.668E-4*Re+0.0100236;
rcdis=a*log(C)+b;
cdA_iso=24./Re(1:s,1).*(1+0.2*Re(1:s,1).^0.63).*(1+BT).^(-0.2);
cd_A=cdA_iso.*rcdis(1:s,1);

%Cd according to Choi in Virpinte dissertation
Cd_VD=2.3*28^(-0.37);
%Cd according to Chiang in Virpinte dissertation
RCd_iso_C=0.536.*Re(1:s,1).^0.118.*(1+BT).^(-0.144)*C^0.345;
Cd_C=RCd_iso_C.*cw_B(1:s,1); %passt nicht!!

%% Make Plots

figure(1)
hold on
set(gca, ...
  'Box'        , 'on'    , ...
  'TickDir'    , 'in'    , ...
  'TickLength'  , [0.01 0.01] , ...
  'XMinorTick'  , 'on'      , ...
  'XTick'      , (0:0.001:0.005) , ...
  'YMinorTick'  , 'on'      , ...
  'XColor'      , [.2 .2 .2], ...
  'YColor'      , [.2 .2 .2], ...
  'LineWidth'  , 1        ,...
   'FontSize'    , 12 ,...
   'XTickLabel', sprintf('%.3f|', 0:0.001:0.005))
plot(t(1:end-2,1),cd_A,'-c')
plot(t(1:end,1),cd_M,'-r')
plot(t(1:end,1),cd_V(1:end,1),'-b')
pl1=plot(t(1:end,1), cw)
set(pl1,'LineWidth',1,'LineStyle', '-','Color', [0 .8 0])
plot(t(1:end,1),Cd_VD(1:end,1),'-k')
plot(t(1:s,1),Cd(1:s,1),'-m')
set(gcf,'PaperUnits','centimeters','PaperSize',[15 12],'PaperPositionMode','auto')
legend('Attahsit','Mulholland','Virepinte','FS3D','Choi', 'Sirignano')
legend boxoff
ylim([0.5 3])
xlim([0 5E-3])
xlabel('Time [s]')
ylabel ('Drag coefficent [-]')
fname = strcat('cdc2_c45');
fig = gcf;
style = hgexport('factorystyle');
style.Bounds = 'tight';
hgexport(fig,fname,style,'applystyle', true);
drawnow;
print('-depsc', '-tiff',fname)
Private Nachricht senden Benutzer-Profile anzeigen


MaraBelle
Forum-Anfänger

Forum-Anfänger


Beiträge: 36
Anmeldedatum: 24.02.15
Wohnort: ---
Version: R2014b
     Beitrag Verfasst am: 23.07.2015, 10:31     Titel:
  Antworten mit Zitat      
Hallo!
Ich hoffe, du hast deine Arbeit noch rechtzeitig fertig bekommen? Hast du eine Lösung gefunden bzw. herausgefunden, wo das Problem liegt?
Ich habe nämlich ein ähnliches Problem. Ich habe ein Diagramm mit mehreren Plots, das auch schon problemlos funktioniert hat. Ich habe woanders im Skript was geändert, was an den eigentlich zu plottenden Werten aber nichts ändert (Diese wurden nur verwendet um einen weiteren Wert daraus zu berechnen). Jetzt dauert es ewig, bis die plots geplottet werden bzw. tut sich irgendwann scheinbar gar nichts mehr...
LG
Private Nachricht senden Benutzer-Profile anzeigen
 
Winkow
Moderator

Moderator



Beiträge: 3.842
Anmeldedatum: 04.11.11
Wohnort: Dresden
Version: R2014a 2015a
     Beitrag Verfasst am: 23.07.2015, 11:37     Titel:
  Antworten mit Zitat      
allgemein kann ich nur sagen das man mal mit dem profiler gucken kann wo die zeit verbraucht wird. wie das geht steht in der doc. oder auch zb in diesem buch
https://books.google.de/books?id=9o.....onepage&q&f=false
Zitat:
Ich habe woanders im Skript was geändert, was an den eigentlich zu plottenden Werten aber nichts ändert (Diese wurden nur verwendet um einen weiteren Wert daraus zu berechnen). Jetzt dauert es ewig, bis die plots geplottet werden bzw. tut sich irgendwann scheinbar gar nichts mehr...
LG

ohne den dazugehörigen code kann man unmöglich sagen woran es liegt. eröffne dazu aber am besten einen neuen thread.
_________________

richtig Fragen
Private Nachricht senden Benutzer-Profile anzeigen
 
MaraBelle
Forum-Anfänger

Forum-Anfänger


Beiträge: 36
Anmeldedatum: 24.02.15
Wohnort: ---
Version: R2014b
     Beitrag Verfasst am: 23.07.2015, 15:10     Titel:
  Antworten mit Zitat      
Vielen Dank für den Tipp mit dem Profiler. Das hätte ich ab und zu schon gut brauchen können.
Ich habe den Fehler bei mir gefunden... Durch meine Veränderung im Skript hat sich die Dimension meiner Variablen von 1440x1 auf 1440x1440 verändert. Im plot hat man den Unterschied nicht gesehen, da in Spalte 2-1440 nur Nullen drin standen. Im Nachhinein wunderts mich jetzt aber, dass MATLAB das überhaupt noch plotten konnte...
Private Nachricht senden Benutzer-Profile anzeigen
 
Neues Thema eröffnen Neue Antwort erstellen



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 | goMatlab RSS Button RSS

Hosted by:


Copyright © 2007 - 2025 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.