We propose an efficient wideband theory for adjoint variable sensitivity analysis of problems with general anisotropic materials. The method is formulated based on the transmission line numerical modeling technique. The anisotropic material properties of potential interest are the full tensors of permittivity, permeability, electrical conductivity, magnetic resistivity, magnetoelectric coupling, and electromagnetic coupling. The tensors may contain non-diagonal elements. Our method estimates the gradients of the desired response with respect to all designable parameters using at most one extra simulation, regardless of their number. In contrast, in the conventional sensitivity analysis method using central finite differences, the number of the required simulations scales linearly with the number of designable parameters. The theory has been implemented for sensitivity analysis of the two and three-dimensional structures. The available adjoint variable method (AVM) sensitivities enable the optimization-based design of anisotropic and dispersive anisotropic structures. We apply our AVM technique to optimization-based wideband invisibility cloak design of arbitrary-shape objects. Our method optimizes the voxel-by-voxel constitutive parameters of an anisotropic cloak. This results in a large number of optimizable parameters. The associated sensitivities of a wideband cloaking objective function are efficiently estimated using our anisotropic adjoint variable method technique. A gradient-based optimization algorithm utilizes the available sensitivity information to iteratively minimize the visibility objective function and to determine the constitutive parameters of the optimal cloak. / Thesis / Doctor of Philosophy (PhD)