You could create a local copy of the function.
Somes it makes sense to replace dirac with piecewise and use matlabFunction with File output (matlabFunction cannot convert piecewise for anonymous functions.) This does restrict the input to be scalar, but on the other hand it can permit simplifications, sometimes large simplifications.
Or sometimes it just makes sense to transform x*dirac(expression) to (expression == 0).*x provided that x can never be infinity. This can be a bit tricky to do right automatically, but the new R2019a https://www.mathworks.com/help/symbolic/sym.mapsymtype.html can help. Still, sometimes it is just easier to matlabFunction with file output and edit the result to stick in an anonymous function that maps dirac to @(x) x==0
x == 0 (producing 0 or 1) is of course not the real way that dirac is defined: dirac is defined as returning infinity at 0. But dirac as used in practice in integration has an area of 1 at the one point, so value * dirac(expression) inside integration expresses that the contrabution when the expression is exactly 0 is the value and 0 at any other location.
dirac with matlabFunction gets a bit questionable since when processed numerically it contributes either 0 or infinity and infinity times whatever it is multiplied by is infinity and those infinities tend to "poison" the calculations. If you are using dirac() with integral() then you need to understand the limitations: if integral() does not happen to evaluate at exactly the matching value then there will be no contribution, and if it does evaluate at exactly the matching value you will get an infinity. integral() does not know anything about the special properties of dirac.
