Verfasst am: 12.12.2018, 10:44
Titel: Unterschiede Itegration quadgk und integral
Hallo,
ich habe ein kleines Probelm und weiß nicht, ob es an einen Syntax-Fehler liegt oder ob die Funktion einen anderen Wert berechnet.
Im Rahmen meiner Abschlussarbeit entwickle ich eine Fit-Funktion für ein physikalisches Thema dabei steht die Integration einer komplexen Stammfunktion an. Bevor ich eine Fit-Funktion erstellt habe, habe ich das Programm manuell getestet und dabei folgende Integration verwednet:
Minimalbeispiel
Code:
%A= Konstante
%B= Konstante
%C= Konstante
%D= Konstante
%E= Vektor Länge omega
%Fermi= Aufruf einer Funktion :
%function y=Fermi(x,Ef,T)
%kB = 8.617e-5;
%y=1./(1+exp((x-Ef)./(kB*T)));
%end
%T=Konstante
%EF= Konstante
omega=0:4000;
for j=1:length(omega)
siginter(j) = A,*omega.*quadgk(@(E)[(1+(B^2./E.^2))./((2*E).^2-(C.*omega(j)+1i*D).^2)].*[(Fermi(E,Ef,T)-Fermi(-E,Ef,T))], 0, Inf,'RelTol',0,'AbsTol',1e-12);
end
In dieser Berechnung habe ich keinen Fehler entdeckt und sie funktioniert so wie sie soll. Da ich jetzt eine Fit-Funktion durch fminsearch durchführe erweist sich die for-schleife als Flaschenhals und ich habe diese Berechnung durch die Funktion integral ersetzt.
Der Code ist wie folgt:
%Selben Konstanten wie zuvor
Code:
%Selben Konstanten wie zuvor
siginter =A.*omega.*integral(@(E)((1+((B^2)./(E.^2)))./((2*E).^2)-(C.*omega+1i.*D).^2).*[(Fermi(E,Ef,T)-Fermi(-E,Ef,T))], 0, 1000,'ArrayValued',true,'RelTol',1e-10,'AbsTol',1e-12);
Was nicht funktioniert ist eine Integrationsgrenze von Unendlich und auch die Funktion erscheint komplett verschieden berechnet zu werden. Erkennt jemand den Fehler.
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.
RSS