Über die Verwendung von Quad in MATLAB

Question

Ich stoße auf ein seltsames Problem mit quad Funktion. Ich war mit quad einfachen Integral zu berechnen, und es funktionierte für 10 bis 20-mal, dann gibt Matlab die folgenden Fehler:

Error using quad (line 75) 
The integrand function must return an output vector of the same length as the input vector. 
yteor(k) = quad(@(q)(exp(-(q).^2).*q.^2/(k.^2+1)), 0, 1); 

Hier q und k Skalare sind. Ich kann nicht bekommen, was falsch ist und warum es vor einigen Stunden funktionierte.

bearbeiten

Hier ist mein Code

for k=1:100, 
    xteor(k)=step*k; 
    yteor(k)=quad(@(q)(exp(-(q).^2).*q.^2/((step.*k+1).^2)),0,1); 
end plot(xteor,yteor,'r'); 

Answer 1

Der folgende Ausschnitt funktioniert für mich auf Octave (Matlab GNU-Klon)

step = 1; 
xteor = zeros(100,1); 
yteor = zeros(100,1); 

for k=1:100, 
    xteor(k)=step*k; 
    yteor(k)=quad(@(q)(exp(-(q).^2).*q.^2/((step.*k+1).^2)),0,1); 
end 
plot(xteor,yteor,'r'); 
pause 

Meine Hypothese ist, dass Ihr Fehler das ist Konsequenz von etwas anderes passiert früher in Ihrem Code (vielleicht bezogen auf step kein Skalar?). Anstatt sich auf diese Linie zu konzentrieren, wo der Fehler entsteht. Versuchen Sie zu suchen, was Sie geändert haben, bevor der Fehler angezeigt wird.