I am new with Matlab, and I have the a code with two loops and I need vectorize, Can someone help me?
function f = eval_Funcion_RHE(x)
n = length(x);
f = 0;
for i = 1:n
sumaInt = 0;
for j = 1: i
sumaInt = sumaInt + (x(j)^2);
end
f = f + sumaInt;
end
end

 Akzeptierte Antwort

Bruno Luong
Bruno Luong am 7 Okt. 2018
Bearbeitet: Bruno Luong am 7 Okt. 2018

1 Stimme

x = x(:).';
f = sum(x.^2.*(length(x):-1:1))

7 Kommentare

Iori Yagami Kusanagi
Iori Yagami Kusanagi am 7 Okt. 2018
Thank you so much... I will try!
Iori Yagami Kusanagi
Iori Yagami Kusanagi am 7 Okt. 2018
Bearbeitet: Iori Yagami Kusanagi am 7 Okt. 2018
@Bruno Luong: Now, if the function is more complex, how can you optimize it?
for i = 1:n-1
a = max (i-1,1);
b = min (i+1,n);
prodc = 1;
for j = a: b
prodc = prodc * (x(j)^3);
end
f = f + prodc;
end
Bruno Luong
Bruno Luong am 7 Okt. 2018
Bearbeitet: Bruno Luong am 8 Okt. 2018
fp = exp(3*conv(log(x(:)),ones(3,1)));
f = real(sum(fp(2:n)))
Iori Yagami Kusanagi
Iori Yagami Kusanagi am 7 Okt. 2018
I tried, but i get this message: Error using conv2
Error using conv2 Invalid data type. First and second arguments must be numeric or logical.
Error in conv (line 43) c = conv2(a(:),b(:),shape);
Error in Funcion_EN (line 9) fp = exp(3*conv(log(x(1:n)),ones(3,1)));
Bruno Luong
Bruno Luong am 8 Okt. 2018
make sure x contain numbers, make sure you don't name something else ONES.
Put a breakpoint and check.
Bruno Luong
Bruno Luong am 8 Okt. 2018
Bearbeitet: Bruno Luong am 8 Okt. 2018
I run this:
% random input
x= rand(1,10);
% for-loop method
n = length(x);
f = 0;
for i = 1:n-1
a = max (i-1,1);
b = min (i+1,n);
prodc = 1;
for j = a: b
prodc = prodc * (x(j)^3);
end
f = f + prodc;
end
f
% vectorize method with conv
n = length(x);
fp = exp(3*conv(log(x(:)),ones(3,1)));
f = real(sum(fp(2:n)))
I get this
>> test
f =
0.0199
f =
0.0199
>>
Iori Yagami Kusanagi
Iori Yagami Kusanagi am 8 Okt. 2018
Thank you, again!!!

Melden Sie sich an, um zu kommentieren.

Weitere Antworten (0)

Produkte

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by