EDIT: Reread it some hours later and found my error. I figured I was doing something wrong. I was applying operations out of order when calculating the conditional probability. It is 1/2 in each case. I'll leave the commentary untouched.
I think the answer is Yes, or at least I'm not entirely convinced the answer is no.
I will provide an example below, but I don't find it very convincing since I just ad-hoc approached it, and don't have a nice "overarching" principle to take away from this. Basically, consider this more as a comment to get discussion going, than a full fledged answer.
The other answers show that the expectation value of measuring the system to be in a particular state is the same. Basically the density matrix of the ensemble is the same, but the density matrix of the first machine only has two possible outputs while the second has an infinite number. Focusing immediately on the ensemble average seems to be throwing away any possibility we have of distinguishing them.
Here's an attempt at distinguishing them:
Machine 1 possible output, only pure states
$|0\rangle$
$|1\rangle$
Machine 2 possible output, any state
$\frac{1}{\sqrt{2}}(|0\rangle + p |1\rangle)$
where $p = e^{i\theta}$ with $0 \le \theta < 2\pi$
Now take some other qubit B (it doesn't matter here physically what it is) of prepared state $\frac{1}{\sqrt{2}}(|0\rangle+|1\rangle)$ to get the product states:
machine 1
$\frac{1}{\sqrt{2}}(|0\rangle+|1\rangle)|0\rangle = \frac{1}{\sqrt{2}}(|00\rangle +|10\rangle)$
$\frac{1}{\sqrt{2}}(|0\rangle+|1\rangle)|1\rangle = \frac{1}{\sqrt{2}}(|01\rangle +|11\rangle)$
machine 2
$\frac{1}{2}(|0\rangle+|1\rangle)(|0\rangle + p |1\rangle) = \frac{1}{2}(|00\rangle+p|01\rangle + |10\rangle + p |11\rangle)$
Now let's introduce an interaction which can cause some interference:
$|00\rangle \rightarrow |00\rangle$
$|01\rangle \rightarrow \frac{1}{\sqrt{2}}(|01\rangle + |10\rangle)$
$|10\rangle \rightarrow \frac{1}{\sqrt{2}}(|01\rangle - |10\rangle)$
$|11\rangle \rightarrow |11\rangle$
now we have
machine 1
$\frac{1}{\sqrt{2}}(|00\rangle +|10\rangle) \rightarrow \frac{1}{\sqrt{2}}|00\rangle + \frac{1}{2}(|01\rangle-|10\rangle)$
$\frac{1}{\sqrt{2}}(|01\rangle +|11\rangle) \rightarrow \frac{1}{\sqrt{2}}|11\rangle + \frac{1}{2}(|01\rangle+|10\rangle)$
machine 2
$\frac{1}{2}(|00\rangle+p|01\rangle + |10\rangle + p |11\rangle) \rightarrow \frac{1}{2}(|00\rangle + p\frac{1}{\sqrt{2}}(|01\rangle+|10\rangle) + \frac{1}{\sqrt{2}}(|01\rangle-|10\rangle) + p |11\rangle)$
$ \ \ \ = \frac{1}{2}(|00\rangle + (p+1)\frac{1}{\sqrt{2}}|01\rangle+(p-1)\frac{1}{\sqrt{2}}|10\rangle + p |11\rangle)$
Now let's do two measurements. First measure the state of B to be 0 or 1, then measure the sate of the atom to be 0 or 1.
Conditional probability on the ensemble:
Given that we find B in state 1, what is the probability of finding the atom in state 0?
machine 1
(1/2) x 1 + (1/2) x (1/3) = 4/6
machine 2
$\frac{\frac{1}{2}(p-1)^2}{\frac{1}{2}(p-1)^2 + p^2} = \frac{\frac{1}{2}(2 - 2\cos\theta)}{\frac{1}{2}(2 - 2\cos\theta) + 1} = \frac{1 - \cos\theta}{2 - \cos\theta}$
Now averaging over $\theta$
$\mathrm{Prob} = \frac{1}{2\pi}\int_0^{2\pi}\frac{1 - \cos\theta}{2 - \cos\theta} d\theta = 1 - \frac{1}{\sqrt{3}}$
Now, it is quite possible I've made a mistake here. But my main point is that the other answers seem to be throwing away the useful information to obtain solely an average of the initial output states. As the answers stand now, they do not mathematically convince me that we can never obtain an effect by adding interactions and multiple measurements with conditional probability or maybe 'weak' measurements, since individually the states have much different density matrices. Hopefully I didn't make a mistake above, but even if I did, I'd still very much like to hear more in the other answers beyond what is currently written. This is a fascinating question, so I'm quite interested in discussing this further.