BJT schmitt trigger threshold voltage calculation
I'm sorry for the delay. I've been busy. I hope you won't mind if I take a moment now to expand on my earlier comment about an approach to the \$\uparrow\$ transition voltage calculation. You can decide if it is useful.

Preview

The turn-off of \$Q_2\$ starts when it goes out of saturation and ends right at the point when there is no collector or base current for \$Q_2\$. Right then, given there is no \$Q_2\$ base current there is also no voltage drop across \$R_{\text{B}_2}\$ and therefore \$Q_1\$'s \$V_{\text{CE}_1}\$ has exactly one diode drop across it -- as would be expected just before \$Q_1\$ starts moving into saturation.

Essentially, this means the just barely active mode collector current in \$Q_1\$ should be just \$I_{\text{C}_1}=\frac{V_\text{CC}-V_{\text{BE}_2}}{R_\text{E}+R_{\text{C}_1}}\$. Knowing \$\beta_1\$, that allows you to work out the voltage drop across \$R_{\text{B}_1}\$ with the only remaining problem being to find out the voltage drop across \$R_\text{E}\$.

So we expect to see a confirmation about the expected value for \$I_{\text{C}_1}\$ and enough additional information to find the voltage drop across \$R_\text{E}\$.

Currents

There are three significant currents (neglecting the base current of \$Q_1\$): \$I_{\text{C}_1}\$, \$I_{\text{C}_2}\$, and \$I_{\text{B}_2}\$. The sum of them will pass through \$R_\text{E}\$. \$I_{\text{C}_1}\$ merely diverts current around \$Q_2\$ but doesn't actually change the current in \$R_\text{E}\$, so this leads rapidly to the conclusion that while in this state and up until just before transition the voltage across \$R_\text{E}\$ will be constant.

In this context, \$I_{\text{C}_2}\$ and \$I_{\text{B}_2}\$ are a function of \$I_{\text{C}_1}\$.

Two KVL Equations

There are just two KVL equations to solve simultaneously for \$I_{\text{C}_2}\left(I_{\text{C}_1},\beta_2,V_{\text{CE}_2}\right)\$ and \$I_{\text{B}_2}\left(I_{\text{C}_1},\beta_2,V_{\text{CE}_2}\right)\$:

$$\begin{align*}
V_\text{CC}-I_{\text{C}_2}\:R_{\text{C}_2}-V_{\text{CE}_2}-\left(I_{\text{C}_1}+I_{\text{C}_2}+I_{\text{B}_2}\right)\:R_\text{E}&amp;=0\:\text{V}\\\\
V_\text{CC}-\left(I_{\text{C}_1}+I_{\text{B}_2}\right)\:R_{\text{C}_1}-I_{\text{B}_2}\:R_{\text{B}_2}-V_{\text{BE}_2}-\left( I_{\text{C}_1} + I_{\text{C}_2} + I_{\text{B}_2} \right)\:R_\text{E} &amp;=0\:\text{V}
\end{align*}$$

These can be solved into two parametric equations based on \$I_{\text{C}_1}\$.

Transition

\$\beta_2\$ is saturated and low to start, but gradually rises towards its active mode value. At that moment the system switches state until the collector current of \$Q_2\$ effectively goes to zero (implying a transition towards \$\beta_2=\frac{I_{\text{C}_2}=0\:\text{A}}{I_{\text{B}_2}}=0\$.) This is why I wrote my earlier comment.

Another way of saying this is that the system goes through a transition where the voltage drop across \$R_{\text{C}_2}\$ goes to zero and \$V_{\text{CE}_2}\$ occupies all of the voltage that remains after subtracting the emitter voltage.

Algorithm


First, solve the simultaneous KVL equations for \$I_{\text{C}_2}\left(I_{\text{C}_1},\beta_2,V_{\text{CE}_2}\right)\$ and \$I_{\text{B}_2}\left(I_{\text{C}_1},\beta_2,V_{\text{CE}_2}\right)\$.
Find \$I_{\text{C}_1}\left(V_{\text{CE}_2}\right)\$ where \$\frac{I_{\text{C}_2}\left(I_{\text{C}_1},\beta_2,V_{\text{CE}_2}\right)}{I_{\text{B}_2}\left(I_{\text{C}_1},\beta_2,V_{\text{CE}_2}\right)}=\beta_2\$.
Substitute \$I_{\text{C}_1}\$ in and solve \$I_{\text{C}_2}\left(I_{\text{C}_1},\beta_2\right)=0\$ in order to solve for \$V_{\text{CE}_2}\$.
Solve \$I_{\text{C}_1}\$.


The point here is to look for the transition where \$V_{\text{CE}_2}\$ goes from saturated to a fully open state. This is handled between step 2 and step 3 above, which finds the \$V_{\text{CE}_2}\$ needed at the end of the transition. (Step 2 finds the \$I_{\text{C}_1}\$ needed at the beginning of the transition.)

\$\uparrow\$ Threshold Voltage

Performing all of the above results in a very simple equation for the transition collector current of \$Q_1\$:

$$I_{\text{C}_1}=\frac{V_\text{CC}-V_{{BE}_2}}{R_\text{E}+R_{\text{C}_1}}$$

Note that this confirms the rough idea proposed at the outset. But now there is more information from which to work out the rest. The sum of the three currents can be computed and therefore the voltage across \$R_\text{E}\$. The base current can be estimated for \$Q_1\$ and therefore the voltage drop across \$R_{\text{B}_1}\$. Plug in \$V_{\text{BE}_1}\$ and you have the voltage.

Summary

The reasoning for the above algorithm is the following chain:


The circuit has three currents of significance, the collector currents for \$Q_1\$ and \$Q_2\$, plus the base current for \$Q_2\$.
The KVL equations provide only two equations. Clearly, we cannot solve for all three currents with only two equations.
So solve the KVL equations into two parametric equations, dependent upon the parameter \$I_{\text{C}_1}\$.
Without knowing \$I_{\text{C}_1}\$, we can't determine the value of \$V_\text{E}\$ since we need all three currents for that. But without knowing \$V_\text{E}\$, we can't know the value of \$I_{\text{C}_1}\$, either.
The critical time is when \$Q_2\$ is moving rapidly out of saturation and into active mode. In this narrow window of time, \$\frac{I_{\text{C}_2}}{I_{\text{B}_2}}=\beta_2\$ is always true. So solve for the active mode equation for \$I_{\text{C}_1}\$ as a function \$V_{\text{CE}_2}\$.
Choose a convenient value for \$V_{\text{CE}_2}\$. This would be best at the moment when there is \$I_{\text{C}_2}=0\:\text{A}\$ and \$V_{\text{CE}_2}\$ is maximum. However, we don't know \$V_\text{E}\$, yet.
However, we do have a parametric equation for it and we do have an active mode equation for \$I_{\text{C}_1}\$ and we can insert that into our parametric equation for \$I_{\text{C}_2}\$ to provide an equation for \$I_{\text{C}_2}\$ as a function of \$V_{\text{CE}_2}\$. Set this equation to \$0\:\text{A}\$ to solve for \$V_{\text{CE}_2}\$.
Apply this value of \$V_{\text{CE}_2}\$ to the equation developed in step 5 above to get the last transition value of \$I_{\text{CE}_1}\$.
Now, all three currents are known and the circuit is solved for this transition.


I hope that explains it well. Of course, all this is subject to testing. I will need to do some runs to verify it.

--- IMPLEMENTATION ---

Here is the sympy script I used to implement the above discussion, without having to use up a bunch of paper (and to avoid tl;dr texts here):

    var('ic1 ic2 ib2 beta1 beta2 rb1 rb2 re rc1 rc2 vbe1 vbe2 vce2 vcc')
    kcl = solve([Eq(vcc - ic2*rc2 - vce2 - (ic2+ic1+ib2)*re, 0),
                 Eq(vcc - (ic1+ib2)*rc1 - ib2*rb2 - vbe2 - (ic2+ic1+ib2)*re, 0)],
                 [ic2, ib2])
    ic1_active_q2 = solve(Eq(kcl[ic2]/kcl[ib2], beta2),ic1)[0]
    vce2_final = solve(Eq(kcl[ic2].subs({ic1:ic1_active_q2}),0), vce2)[0]
    ic1_final = ic1_active_q2.subs({vce2:vce2_final})


The result of that last line is the following equation for the final collector current of \$Q_1\$ at the end of the transition period:

$$I_{\text{C}_1}=\frac{V_\text{CC}-V_{\text{BE}_2}}{R_\text{E}+R_{\text{C}_1}}$$

Adding these lines, I can compute the value of the threshold voltage just at the point where \$Q_2\$ leaves saturation as:

    ic2_final = kcl[ic2].subs({ic1:ic1_final})
    ib2_final = kcl[ib2].subs({ic1:ic1_final})
    ire_final = ic1_final + ic2_final + ib2_final
    vth_high = ire_final*re + vbe1 + rb1*ic1_final/beta1
    vth_low = ic1_final*re + vbe1 + rb1*ic1_final/beta1


Suppose \$R_{\text{B}_1}=27\:\text{k}\Omega\$, \$R_{\text{B}_2}=100\:\Omega\$, \$R_{\text{C}_1}=1.5\:\text{k}\Omega\$, \$R_{\text{c}_2}=1\:\text{k}\Omega\$, and \$R_\text{E}=47\:\Omega\$. Let's assume \$V_\text{BE}=700\:\text{mV}\$ for both transistors and that both have \$\beta=200\$ when active mode. Assuming that just leaving saturation and entering into active mode occurs at \$V_{\text{CE}_2}=700\:\text{mV}\$, then I get: \$V_\text{H}\approx 1.388\:\text{V}\$ and \$V_\text{L}\approx 1.206\:\text{V}\$.

Let's see what spice shows:



Using .MEAS triggers in LTspice, it reports: \$V_\text{H}\approx 1.374\:\text{V}\$ and \$V_\text{L}\approx 1.209\:\text{V}\$. Hmm. Not bad, just doing this on the fly, so to speak.



Let set \$R_\text{E}=100\:\Omega\$. From this I get: \$V_\text{H}\approx 1.678\:\text{V}\$ and \$V_\text{L}\approx 1.332\:\text{V}\$.



Using .MEAS triggers in LTspice, it reports: \$V_\text{H}\approx 1.669\:\text{V}\$ and \$V_\text{L}\approx 1.340\:\text{V}\$. Also not bad.