Smys Umgebung und plot

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

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

TheFinger

Forum-Anfänger


	Beiträge: 28

	Anmeldedatum: 14.05.11

	Wohnort: ---

	Version: ---

Verfasst am: 30.05.2011, 18:06 Titel: Smys Umgebung und plot

Hallo jungs,
ich bin verzweifelt, ich muss eine Kraft anhand einer Funktion rechnen, die sowohl Integral als auch partielle ableitungen beinhaltet. Dann soll ich diese Funktion plotten.
Das problem ist ohne syms umgebung sagt mir Matlab, er kennt meine variablen nicht oder bringt eine andere Fehler meldung. Wenn ich das mit double lösen möchte, ist er 30 min lang busy und es kommt nichts raus.
Hier ist mein Code, vielleicht kann mir jemand helfen:

Code:

clear
close all

fpath='G:/Backup2/Simulations/Experiments_Lemoine/OneThirdRule/Versuch2_Experiment1_mitVerdampfung/'

%Import data
mass_data=importdata([fpath 'mass00.his']);
temp_data=importdata([fpath 'temper00.his']);

%Evaluation

rhol=0.78465;
Vol=mass_data(:,2)/rhol;
r_s=0.5*(6*Vol/pi).^(1/3); %Droplet radius
r_0=0.5* (6*Vol(1)/pi)^(1/3);
r_end=0.5* (6*Vol(end)/pi)^(1/3);

Ts=temp_data(:,4); %Surface temperature
Tc=516.2;
T_r=Ts/Tc;

A=131.38;
B=5.5437;
C=-8.4826;
D=4.3164;

sigma=A.*(1-T_r).^(B+C.*T_r+D.*(T_r).^2); %Surface tension

syms sig Tr y
sig=A.*(1-Tr).^(B+C.*Tr+D.*(Tr).^2); % Hier veruche sigma über T_r abzuleiten
dsdT=diff(sig,Tr);

%Marangoni force

F_marangoni=1./(pi*r_s)*dsdT*int(sin(y)*sin(y),y,0,2*pi);
F=double(F_marangoni);

q=1000;
figure
hold on
plot(temp_data(1:q:end,1),F_marangoni(1:q:end,1))

Funktion ohne Link?

P.S.: habe riesen Daten, das sind alle Vektoren der Größe (730000X1)

Thomas84

Forum-Meister


	Beiträge: 546

	Anmeldedatum: 10.02.10

	Wohnort: ---

	Version: ---

Verfasst am: 01.06.2011, 08:13 Titel:

Code:

A=131.38;
B=5.5437;
C=-8.4826;
D=4.3164;
syms sig Tr y
sig=A.*(1-Tr).^(B+C.*Tr+D.*(Tr).^2); % Hier veruche sigma über T_r abzuleiten
dsdT=diff(sig,Tr);
syms r_s
F_marangoni=1./(pi*r_s)*dsdT*int(sin(y)*sin(y),y,0,2*pi);

Funktion ohne Link?

nun noch Tr und r_s mit subs ersetzen und in double umwandeln.

Code:

F_v = subs(F_marangoni,{Tr,r_s},{Tr_data,R_s_data});

Funktion ohne Link?


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 - 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.