What is the advantage of linspace over the colon “:” operator?

by

It’s not just the usability. Though the documentation says:

The linspace function generates linearly spaced vectors. It is
similar to the colon operator :, but gives direct control over the
number of points.

it is the same, the main difference and advantage of linspace is that it generates a vector of integers with the desired length (or default 100) and scales it afterwards to the desired range. The : colon creates the vector directly by increments.

Imagine you need to define bin edges for a histogram. And especially you need the certain bin edge 0.35 to be exactly on it’s right place: 

edges = [0.05:0.10:.55];
X = edges == 0.35

edges =   0.0500    0.1500    0.2500    0.3500    0.4500    0.5500
X =  0     0     0     0     0     0

does not define the right bin edge, but:

edges = linspace(0.05,0.55,6);   %// 6 = (0.55-0.05)/0.1+1
X = edges == 0.35

edges =   0.0500    0.1500    0.2500    0.3500    0.4500    0.5500
X =  0     0     0     1     0     0

does.

Well, it’s basically a floating point issue. Which can be avoided by linspace, as a single division of an integer is not that delicate, like the cumulative sum of floting point numbers. But as Mark Dickinson pointed out in the comments:
You shouldn’t rely on any of the computed values being exactly what you expect. That is not what linspace is for. In my opinion it’s a matter of how likely you will get floating point issues and how much you can reduce the probabilty for them or how small can you set the tolerances. Using linspace can reduce the probability of occurance of these issues, it’s not a security.

That’s the code of linspace:

n1 = n-1
c = (d2 - d1).*(n1-1) % opposite signs may cause overflow
if isinf(c)
    y = d1 + (d2/n1).*(0:n1) - (d1/n1).*(0:n1)
else
    y = d1 + (0:n1).*(d2 - d1)/n1
end

To sum up: linspace and colon are reliable at doing different tasks. linspace tries to ensure (as the name suggests) linear spacing, whereas colon tries to ensure symmetry

In your special case, as you create a vector of integers, there is no advantage of linspace (apart from usability), but when it comes to floating point delicate tasks, there may is.

The answer of Sam Roberts provides some additional information and clarifies further things, including some statements of MathWorks regarding the colon operator.

More Related Contents:

  1. Generate a matrix containing all combinations of elements taken from n vectors
  2. Linear indexing, logical indexing, and all that
  3. Repeat copies of array elements: Run-length decoding in MATLAB
  4. Element-wise array replication in Matlab
  5. arrayfun can be significantly slower than an explicit loop in matlab. Why?
  6. Using transpose versus ctranspose in MATLAB
  7. How do I do multiple assignment in MATLAB?
  8. MATLAB adding array elements iteratively: time behavior
  9. Creating Indicator Matrix
  10. In MATLAB, when is it optimal to use bsxfun?
  11. Array of Matrices in MATLAB
  12. How do i define a structure in Matlab
  13. Element-wise array replication according to a count [duplicate]
  14. In an assignment A(:) = B, the number of elements in A and B must be the same
  15. How to use cell arrays in Matlab?
  16. how to replicate an array
  17. Is there anything like deal() for normal MATLAB arrays? [duplicate]
  18. convert string to number array in matlab
  19. How do I iterate through each element in an n-dimensional matrix in MATLAB?
  20. Is indexing vectors in MATLAB inefficient?
  21. How do I remove elements at a set of indices in a vector in MATLAB?
  22. Vector norm of an array of vectors in MATLAB
  23. Is A==0 really better than ~A?
  24. remove an index number from array
  25. Passing arrays as parameters in bash
  26. What is the ideal growth rate for a dynamically allocated array?
  27. Two-dimensional array in Swift
  28. How to iterate over array of objects in Handlebars?
  29. 1D Number Array Clustering
  30. Swift 3 2d array of Int

Leave a Comment