We present results for the matrix elements of a variety of four-fermion operators calculated using quenched Wilson fermions. Our simulations are done on 170 lattices of size 32 3 ϫ64 at ϭ6.0. We find B K ϭ0.74Ϯ0.04Ϯ0.05, B D ϭ0.78Ϯ0.01, B 7 3/2 ϭ0.58Ϯ0.02 Ϫ0.03 ϩ0.07 , B 8 3/2 ϭ0.81Ϯ0.03 Ϫ0.02 ϩ0.03 , with all results being in the NDR scheme at ϭ2 GeV. We also calculate the B parameter for the operator Q s , which is needed in the study of the difference of B-meson lifetimes. Our best estimate is B S (NDR,ϭ1/aϭ2.33 GeV) ϭ0.80Ϯ0.01. This is given at the lattice scale since the required two-loop anomalous dimension matrix is not known. In all these estimates, the first error is statistical, while the second is due to the use of truncated perturbation theory to match continuum and lattice operators. Errors due to quenching and lattice discretization are not included. We also present new results for the perturbative matching coefficients, extending the calculation to all Lorentz scalar four-fermion operators, and using NDR as the continuum scheme. ͓S0556-2821͑97͒05307-1͔