Simulated annealing with Discrete searching space

Hi, I want to implement a discrete SA. I am able to use the builded-in function

[x,fval] = simulannealbnd(ObjectiveFunction,x0,lb,ub)

However, it is a continous search for the variables x and with lb and ub constrains.

Can I discrete the lb bound and ub for a meshgrid and feed into some disrete SA agorithm?

If yes, where can i find the related documetary or Matlab builded-in function?

If no, what kind of method will be recommended for implementing discrete variable space SA.

Becuase I want to compare the results with exhaustive search.