Beginner's question - why doesn't a variable overwrite when vectorizing a loop?

Question

B.M. 2014-2-5

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/115011-beginner-s-question-why-doesn-t-a-variable-overwrite-when-vectorizing-a-loop

编辑： B.M. 2014-2-7

I am learning using the book "Matlab for Engineers and Scientists." Unfortunately this book does not provide answers to the exercises.

I've never been good at writing for loops, or really understanding how to index well, so I am running into problems.

Here's the one at hand. I need to convert a for loop into a vectorization. For simplicity I've written a short example:

s=0;
for n=1:10;
    s=s+n
end

This gives the answers 1, 3, 6, 10, 15, 21, 28, 36, 45, 55. When I try to vectorize it like this, s does not overwrite itself after each iteration, but stays at zero:

s=0;
n=1:10;
s=s+n

This gives s=1:10.

Now if s is a simple variable to explain I could just make it a vector as well, but when I don't know what its values will be I can't specify it beforehand. Any help would be appreciated. I know this is simple, but going through this book, the earlier examples aren't relevant to this type of problem.

EDIT: Here is the actual question:

====================

Work out by hand the output of the following script for n = 4:

n = input( ’Number of terms? ’ );

s = 0;

for k = 1:n

s = s + 1 / (k ˆ 2);

end;

disp(sqrt(6 * s))

If you run this script for larger and larger values of n, you will find that the output approaches a well-known limit. Can you figure out what it is? Now rewrite the script using vectors and array operations. ====================

The limit is pi. When I try to do this with vectors I am entering the following:

====================

s=0;

n=4; (just to start)

k=1:1:n;

s=s+1./k.^2;

a=[sqrt(6.*s)];

====================

Of course, this doesn't work. What am I not getting? (Other that you can't overwrite variables if you aren't doing a for loop...but the exercise asks specifically to do this not using a for loop)

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Walter Roberson 2014-2-5

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/115011-beginner-s-question-why-doesn-t-a-variable-overwrite-when-vectorizing-a-loop#answer_123418

在 MATLAB Online 中打开

cumsum(1:10)

3 个评论
显示 1更早的评论隐藏 1更早的评论

Matt J 2014-2-6

编辑：Matt J 2014-2-6

Vectorization is all about doing the same thing to every value

@Walter

Actually the way I've seen "vectorization" used, it just means "feeding vector/matrix arguments to optimized MATLAB functions" instead of working on their indiviudal elements manually. Your use of CUMSUM qualifies perfectly as an example of vectorization, even though its not doing the same thing to every value.

B.M. 2014-2-7

Ok, I think I get what you're saying, thanks. I'm just going to post the actual question (to the original post) so you can see what the exercise asks.

请先登录，再进行评论。

Answer 2

B.M. 2014-2-5

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/115011-beginner-s-question-why-doesn-t-a-variable-overwrite-when-vectorizing-a-loop#answer_123422

The above was just an example. I just need to know how to get a variable to overwrite itself outside of a for loop.

7 个评论
显示 5更早的评论隐藏 5更早的评论

Matt J 2014-2-7

I have seen your EDIT and fail to see how it is not a sum. The only difference I see relative to your original example is that the formula for the terms in the series being summed is different 1./k.^2

B.M. 2014-2-7

编辑：B.M. 2014-2-7

Ah, you're right. Sorry. I've been doing other exercises since I posted this and didn't remember accurately what the issue was here. I'll go back and look again at the earlier examples.

EDIT: I got it. Jeez. I think I was fatigued or something after doing all those exercises, now it seems very clear. Thanks to all of you.

请先登录，再进行评论。