Variational principles for ratios of linear neutron flux functionals are extended to treat boundary condition and boundary perturbations in nuclear systems which are described by the steady-state inhomogeneous and homogeneous Boltzmann equations. The relationship of these principles with first-order perturbation theory is established. Their applicability to perturbations of the external boundary of the system is investigated through numerical examples.