Encontré esta explicación en algún lugar y la anoté en mis notas personales. Lo explicaré con un ejemplo que creo que ejemplifica por qué Riemann-Stieltjes proporcionará la respuesta incorrecta.
En primer lugar, recordemos cómo podemos definir la integral de Riemann-Stieltjes a continuación
\begin{equation}\int_0^t Z(x)dZ(x)=\lim_{n\rightarrow\infty}\sum_{k=1}^n(Z(t_k)-Z(t_{k-1}))Z(t_{k-1})\tag{*}\end{equation}
donde $0=t_1<t_2<...<t_n=t$ cuando $n\rightarrow\infty$ , $(t_k-t_{k-1})$ va a cero y $Z(x)$ es continua y tiene una variación acotada. La integral de Ito se puede definir de la misma manera (suponiendo que $Z(t)$ para ser cualquier trayectoria browniana). Por lo tanto, en esta definición elemental no hay realmente ninguna diferencia, es sólo que cada uno está tratando con diferentes tipos de funciones. Pero las reglas de integración serán diferentes.
Podemos reescribir los términos de la suma como
$$\sum_{k=1}^n(Z(t_k)-Z(t_{k-1}))Z(t_{k-1})=\frac{1}{2}Z(t)^2-\frac{1}{2}\sum_{k=1}^n(Z(t_k)-Z(t_{k-1}))^2$$
Para ejemplificar, supongamos ahora que $Z(x)=x$ . Sabemos cómo resolver la integral en este caso, esto es $0.5t^2$ como en el primer término de la RHS anterior. Lo que ocurre con esta función es que podemos ignorar el segundo término de arriba, ya que irá a cero (prueba a parear $[0,t]$ con $t_k=kt/n$ por ejemplo). Esto se debe a que $Z(x)=x$ tiene una variación limitada.
¿Por qué no podemos simplemente hacer esto si $Z(t)$ ¿es una trayectoria browniana? El problema aquí es que $\frac{1}{2}\sum_{k=1}^n(Z(t_k)-Z(t_{k-1}))^2$ no se pondrá a cero cuando $n\rightarrow \infty$ . De hecho,
$$\lim_{n\rightarrow\infty}\frac{1}{2}\sum_{k=1}^n(Z(t_k)-Z(t_{k-1}))^2=t$$
Esto es así porque estamos considerando infinitas realizaciones de una variable normal con varianza $t_k-t_{k-1}$ como $n$ va al infinito. El límite de la trayectoria browniana anterior no llega a cero porque no satisface la variación acotada.
Por eso, aunque la integración de Ito y Riemann-Stieltjes parten de la misma definición (*) los resultados son muy diferentes. Si $Z(x)$ es un movimiento browniano obtenemos:
$$\int_0^t Z(x)dZ(x)=\frac{1}{2}Z(t)^2-\frac{1}{2}t$$
Mientras que si $Z(x)$ es una función continua con variación acotada obtenemos
$$\int_0^t Z(x)dZ(x)=\frac{1}{2}Z(t)^2$$