Large-signal (L-S) characterization of double-drift region (DDR) impact avalanche transit time (IMPATT) devices based on silicon designed to operate at different millimeter-wave (mm-wave) and terahertz (THz) frequencies up to 0.5 THz is carried out in this paper using an L-S simulation method developed by the authors based on non-sinusoidal voltage excitation (NSVE) model. L-S simulation results show that the device is capable of delivering peak RF power of 657.64 mW with 8.25% conversion efficiency at 94 GHz for 50% voltage modulation; whereas RF power output and efficiency reduce to 89.61 mW and 2.22% respectively at 0.5 THz for same voltage modulation. Effect of parasitic series resistance on the L-S properties of DDR Si IMPATTs is also investigated, which shows that the decrease in RF power output and conversion efficiency of the device due to series resistance is more pronounced at higher frequencies especially at the THz regime. The NSVE L-S simulation results are compared with well established double-iterative field maximum (DEFM) small-signal (S-S) simulation results and finally both are compared with the experimental results. The comparative study shows that the proposed NSVE L-S simulation results are in closer agreement with experimental results as compared to those of DEFM S-S simulation.