|
|
|
Trapezregel numerische Integration |
|
| K4ndo |
Forum-Newbie
|
 |
Beiträge: 4
|
|
|
|
Anmeldedatum: 29.07.14
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
 |
|
Verfasst am: 02.08.2014, 13:33
Titel: Trapezregel numerische Integration
|
 |
| |
 |
 |
Hallo zusammen,
ich versuche die Funktion k(x) =(0.75+0.1*cos(x))^0.5 numerisch zu approxieren. Dies versuch ich mit der Funktion integrate. Dabei ist x der Wert bis zu dem ich integriere und k die Funktion.Meine Grenzwerte sind dabei immer 0 bis x, wobei x aus einer anderen Funktion übergeben wird. Mein programm sieht wie folgt aus:
function[] = plotkoeffsmat()
eps =0.1;
a=0;
b=1;
lam = 0.75;
x= 0:0.1:100;
k = @(x)(lam+eps*cos(x)).^0.5;
function[y]= integrate(x,k)
p =0:0.01:x;
F =k(p);
y=trapz(p,F);
end
z = zeros(1,length(x));
for i = 1:length(x)
z(i) = 1/sqrt(k(x(i)))*(b/sqrt(k(0))*sin(integrate(x(i),k))+a*sqrt(k(0))*sin(integrate(x(i),k)));
end
plot(x,z,'k:')
legend('Multiskalen');
xlabel('Zeit t')
ylabel('Lösung')
titleString = ['\epsilon = ' num2str(eps)];
title(titleString);
figure(1)
end
Jedoch erhalte ich diese Fehlermeldung
Subscript indices must either be real positive integers or logicals.
Error in ipermute (line 23)
inverseorder(order) = 1:numel(order); % Inverse permutation order
Error in trapz (line 73)
if ~isempty(perm), z = ipermute(z,perm); end
Error in plotkoeffsmat/integrate (line 11)
y=trapz(p,F);
Error in plotkoeffsmat (line 15)
z(i) = 1/sqrt(k(x(i)))*(b/sqrt(k(0))*sin(integrate(x(i),k))+a*sqrt(k(0))*sin(integrate(x(i),k)));
Kann mir jemand sagen was ich falsch gemacht habe? Vielen Dank schonmal!
beste Grüße,
K4ndo
|
|
|
|
|
|
| Harald |
Forum-Meister
|
|
Beiträge: 24.448
|
|
|
|
Anmeldedatum: 26.03.09
|
|
|
|
Wohnort: Nähe München
|
|
|
|
Version: ab 2017b
|
|
|
|
|
|
|
Verfasst am: 04.08.2014, 22:18
Titel:
|
|
Hallo,
mit dem Debugger kannst du leicht rausfinden, wo das Problem liegt:
Deine erste obere Grenze ist 0, und von 0 bis 0 kann man schlecht ein Trapez bilden.
Grüße,
Harald
P.S.: das einzige, was ich grusliger als eine geschachtelte Funktion finde, ist eine geschachtelte Funktion mitten im Code...
|
|
|
|
| K4ndo |
Themenstarter
Forum-Newbie
|
|
Beiträge: 4
|
|
|
|
Anmeldedatum: 29.07.14
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
|
|
|
Verfasst am: 04.08.2014, 22:38
Titel:
|
|
das Problem hab ich inzwischen gelöst 
Mit der Funktion hast du wohl Recht  wollte diese eben mal testen und habs deshalb nich ausgelagert =)
Das Trapez von 0 bis 0 war nicht das Problem, das sollte nach Definition sofort 0 ergeben
Trotzdem vielen dank !
|
|
|
|
| Harald |
Forum-Meister
|
|
Beiträge: 24.448
|
|
|
|
Anmeldedatum: 26.03.09
|
|
|
|
Wohnort: Nähe München
|
|
|
|
Version: ab 2017b
|
|
|
|
|
|
|
Verfasst am: 04.08.2014, 23:15
Titel:
|
|
Hallo,
| Zitat: |
| Das Trapez von 0 bis 0 war nicht das Problem, das sollte nach Definition sofort 0 ergeben |
Ich würde da mal widersprechen. Wenn du die Schleife bei 2 anfängen lässt, läuft der Code nämlich wunderbar durch.
Das Integral von 0 bis 0 mag zwar definiert sein, aber trapz ist eben nicht darauf ausgelegt, diesen Sonderfall zu berücksichtigen.
Grüße,
Harald
|
|
|
|
| K4ndo |
Themenstarter
Forum-Newbie
|
|
Beiträge: 4
|
|
|
|
Anmeldedatum: 29.07.14
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
|
|
|
Verfasst am: 04.08.2014, 23:33
Titel:
|
|
ouh ok vielen dank für die info !
Jetzt läuft auch das mit dem trapz
|
|
|
|
|
|
|
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
|
|
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.
|
|
