Suppose the Fresnel equations give us complex reflexion co-effcients $R_p$ and $R_s$ for $p$- and $s$-polarized light, respectively. Then the intensity reflexion co-efficient (power reflexion coefficient) for depolarized light is (in most cases):
$\frac{1}{2} (|R_s|^2 + |R_p|^2)$
You do likewise for the transmission co-efficients, so that the transmitted power ratio is:
$\frac{1}{2} (|T_s|^2 + |T_p|^2) = 1- \frac{1}{2} (|R_s|^2 + |R_p|^2)$
where $T_p$ and $T_s$ are the Fresnel equation-derived complex transmission co-efficients for $p$- and $s$-polarized light. Forming average square magnitudes like this is often called "incoherent summing".
To understand fully how to do your calculation, you need to understand exactly what depolarized light is, and it has quite a complicated description: it is bound up with decoherence and partially coherent light, a topic which Born and Wolf in "Principles of Optics" give a whole chapter to. A classical description, roughly analogous to Born and Wolf's is as follows: if the transverse (normal to propagation) plane is the $x,y$ plane, then we represent the electric field at a point as:
$\mathbf{E} = \left(\begin{array}{cc}E_x(t) \cos(\omega t + \phi_x(t))\\E_y(t) \cos(\omega t + \phi_y(t))\end{array}\right)$
where $\omega$ is the centre frequency and the phases $\phi_x(t)$, $\phi_y(t)$ and envelopes $E_x(t)$, $E_y(t)$ are stochastic processes, which can be as complicated as you like. The formulas I cite above just assume that:
- $E_x$, $E_y$ and $\phi$ behave like independent random variables, and
- They vary with time swiftly compared to your observation interval (the time interval whereover you gather light in a sensor to come up with an "intensity" measurement) but not so swiftly that the light's spectrum broadened so much that we cannot still think of the light as roughly monochromatic.
A simple quantum description is actually conceptually clearer than Born and Wolf's classical one, as long as light states do not become entangled. Each photon can be thought of as a perfectly coherent wave propagating following Maxwell's equations. The Fresnel equations thus apply to each photon as they would to a perfectly coherent wave. For each photon, therefore, you calculate the intensity of reflexion and transmission, and then average this intensity over all photon polarization states - we assume the source is creating "random" pure states. Thus, suppose the Fresnel equations give us complex reflexion co-effcients $R_p$ and $R_s$ for $p$- and $s$-polarized light: the complex amplitude reflexion co-efficient for a general polarization state is then:
$R(\alpha, \phi) = \alpha R_p e^{i \frac{\phi}{2}} + \sqrt{1-\alpha^2} R_s e^{-i \frac{\phi}{2}}$
where $\alpha \in [0, 1]$ and $\phi \in [0, 2\,\pi)$. Summing intensities over all values of $\phi$ (assuming all phases equally likely) yields:
$\frac{1}{2\pi}\int\limits_0^{2\pi} \left(\alpha^2 |R_p|^2 + (1-\alpha^2) |R_s|^2 + 2 \alpha\sqrt{1-\alpha^2} |R_p| |R_s| \cos\phi\right)\mathrm{d}\phi = \alpha^2 |R_p|^2 + (1-\alpha^2) |R_s|^2$
and then summing intensities over all values of $\alpha^2$ (assuming the $\alpha^2$ is uniformly distributed in $[0,1]$) leaves the formulas above.
This will not give a full picture for general entangled polarization states, when you have to resort to more general coherent and cross correlation functions to describe what is going on. Likewise for Born and Wolf's classical description. But it is an excellent first approximation and it is probably true to say that it is hard to arrange for it not to hold in the laboratory. Deviations from it are likely to be seen if you sample the light intensities over very short sampling intervals, when you will see complicated, extremely swift fluctuations in scattered and transmitted intensities, often following white noise processes.