simplyfing computations in arrays
1 Ansicht (letzte 30 Tage)
Ältere Kommentare anzeigen
int = zeros(1,length(x));
for n = 3:length(x)
y1 = x(2:n) ;
t1 = n*dt - (2*dt:dt:n*dt) ;
y2 = x(1:n-1) ;
t2 = n*dt - (dt:dt:(n-1)*dt) ;
int(n) = sum (t1.^(k-1)/factorial(k-1).*y1 ...
+ t2.^(k-1)/factorial(k-1).*y2)*dt/2 ;
end
I want to write t1= ((n-2)*dt:dt:0dt) is it okay? and also t2=((n-1)dt:dt:dt).
Is it possible? I sew also some errors...
0 Kommentare
Akzeptierte Antwort
Walter Roberson
am 4 Jul. 2012
Bearbeitet: Walter Roberson
am 4 Jul. 2012
Is t1= ((n-2)*dt:dt:0dt) intended to go to 0 ? If so then because n-2 is positive so you are going backwards from positive towards 0, you will need
(n-2)*dt: -dt : 0
2 Kommentare
Walter Roberson
am 5 Jul. 2012
You cannot replace the "for" loop with a different data structure. Your loop index, "n", is used as the endpoint for a colon index, so you are using expressions with different numbers of terms for different values of "n", and expressions with differing numbers of terms depending on the loop variable, cannot be vectorized.
To have any hope of vectorization, you will need to figure out the symbolic difference between the results for adjacent values of "n", and hope that that difference consists only of constants and individual array elements.
Weitere Antworten (1)
Siehe auch
Kategorien
Mehr zu Creating and Concatenating Matrices finden Sie in Help Center und File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!