I am trying to derive an optimal hedging interval for a delta-replicated European call option from the Black-Scholes model. To do that, I would like to compare the resulting hedging error with my trading costs.
I am currently hedging (i.e. adjusting the delta of) the position every day and am trying to find out whether more frequent intraday hedging makes sense.
Unfortunately, I am making a mistake in my assumptions somewhere which I am hoping you can point out.
What I tried so far:
Lets assume volatility and time to maturity stay the same and risk-free-rate = dividends = 0.
I'll also use $\Delta \delta \approx \gamma \cdot \Delta S$ to describe the change in Delta for a small change in price.
I can roughly calculate the hedging error for a small price change using a Taylor expansion of the option price. $C(S) \approx C(S_0) + \delta(S_0) \cdot (S - S_0) + \frac{\gamma(S_0)}{2} \cdot (S - S_0)^2+... $
Since I hold $\delta(S_0)$ in my portfolio, the error for a small change in price is roughly $\frac{1}{2} \gamma (\Delta S)^²$
Additionally, I can calculate the trading costs as: $c \cdot \left| (\delta_t - \delta_0) \cdot S_{t} \right| \approx c \cdot \left| \gamma \cdot \Delta S \cdot (S_0 + \Delta S) \right|$
c is a constant representing my trading costs per trading volume.
Now, using geometric brownian motion $dS_t = \mu S_t \, dt + \sigma S_t \, dW_t$, I find the expected value for my hedging error and trading costs.
$\mathbb{E}[\frac{1}{2} \gamma (\Delta S)^²] \approx \frac{1}{2} \gamma S^2 \sigma^2 \Delta t$
$\mathbb{E}[c \left| \gamma \Delta S (S_0 + \Delta S) \right|] \approx c \gamma S^2 \sigma \sqrt{\Delta t} \sqrt{\frac{2}{\pi}}$
This seems to be correct in my simulations.
However, when I try to calculate the total costs over a greater timespan T (one day for example), the hedging error seems to be the same no matter the hedging frequency. Let $n = \frac{T}{\Delta t}$ the number of intervals/hedges, then $\sum_{i=1}^{n} \frac{1}{2} \gamma S^2 \sigma^2 \Delta t = n \cdot \frac{1}{2} \gamma S^2 \sigma^2 \Delta t = \frac{1}{2} \gamma S^2 \sigma^2 T$
Adding this error to my costs (which decrease by increasing $\Delta t$) doesn't make much sense. Are my calculations wrong or my assumptions too simple? Do I have to calculate it completely differently? What am I doing wrong?
Please only use Black-Scholes and GBM in your answer and please also try to simplify it as much as possible. Thank you very much!