Matlab: so berechnen Sie das bestimmte integral einer Funktion über mehrere Grenzen?

Übernehmen:

z = [0.4 0.5 0.75]'  
function y = myfunct(x)  
y = quad(@sin, 0, x)

Möchte ich berechnen das bestimmte integral von sin(x) von 0 bis 0.4, 0.5 und 0.75, Verwendung: 

myfunct(z)

Doch Matlab gibt:

??? Error using ==> quad at 70  
The limits of integration must be scalars.

Ich wäre dankbar für alle konstruktiven Vorschläge.

InformationsquelleAutor skip | 2011-10-16
  1. 1

    Können Sie auch quadv um dies zu tun. ABER, anstatt einen Vektor der integration Grenzen gesetzt, machen Sie eine array-valued function, so dass beim integrieren Sie jedes element, der Bereich der integration wird 0-1.

    Um genauer zu sein, die Sie integrieren möchten von sin(x) von x = 0 bis z. Dies ist die gleiche wie die Integration von sin(u z)*z aus u = 0 1 (u-substitution). Eine array-Funktion 

    F = @(u) [sin( .4 * u) * .4, sin( .5 * u ) * .5, sin( .75 * u ) * .75 ];

    Dann tun 

    quadv(F,0,1)
    • Ich denke mal, deine Art und Weise vielleicht ein bisschen schneller als meine. Aber ich würde umstrukturieren es ein bisschen, da Sie wiederholen sich einige Informationen (z.B., wenn Sie ändern möchten, 1 limit, ändern Sie die 2 zahlen, wenn Sie wollen, um die Funktion zu ändern, müssen Sie ändern 3 IDS). Der folgende code ist äquivalent zu Ihnen: lims = [0.4 0.5 0.75]'; F = @(u)(sin(u.*lims).*lims); quadv(F,0,1)
    • vereinbart. das ist eine viel sauberere Weg, es zu tun, und verallgemeinert besser. Ich denke, ich habe versucht, zu explizit.
    • Danke Euch beiden für die hilfreichen Kommentare.
    InformationsquelleAutor MarkV
  2. 2

    Können Sie die arrayfun Funktion in neueren Versionen von MATLAB:

    z = [0.4 0.5 0.75]';
y = arrayfun(@(upperLim)(quad(@sin,0,upperLim)),z);
    • Das löst das problem. Vielen Dank!
    InformationsquelleAutor Egon
Schreibe einen Kommentar