Gradient optimization of RF amplifiers for digital communications

A gradient optimization methodology for use with nonlinear envelope simulation is presented. Emphasis is placed on efficient evaluation of cost function gradients. To this end, an envelope sensitivity equation is derived and methods for its efficient solution are proposed. In the case of circuits of moderate size, the solution of the sensitivity equation is shown to be simple and inexpensive. The more troublesome case of larger circuits is treated by a novel application of a recently developed iterative linear equation solver. The result is a general-purpose, rigorous optimization methodology for the fine tuning of gain, adjacent channel power, power efficiency, and related performance measures in radio frequency (RF) amplifiers for digital communications. As an illustrative example, the optimization of a feedforward-linearized power amplifier is presented. © 2000 John Wiley & Sons, Inc. Int J RF and Microwave CAE 10, 353–365, 2000.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/35228/1/4_ftp.pd

Sensitivity Analysis and Distortion Decomposition of Mildly Nonlinear Circuits

Volterra Series (VS) is often used in the analysis of mildly nonlinear circuits. In this approach, nonlinear circuit analysis is converted into the analysis of a series of linear circuits. The main benefit of this approach is that linear circuit analysis is well established and direct frequency domain analysis of a nonlinear circuit becomes possible. Sensitivity analysis is useful in comparing the quality of two designs and the evaluation of gradient, Jacobian or Hessian matrices, in analog Computer Aided Design. This thesis presents, for the first time, the sensitivity analysis of mildly nonlinear circuits in the frequency domain as an extension of the VS approach. To overcome efficiency limitation due to multiple mixing effects, Nonlinear Transfer Matrix (NTM) is introduced. It is the first explicit analytical representation of the complicated multiple mixing effects. The application of NTM in sensitivity analysis is capable of two orders of magnitude speedup. Per-element distortion decomposition determines the contribution towards the total distortion from an individual nonlinearity. It is useful in design optimization, symbolic simplification and nonlinear model reduction. In this thesis, a numerical distortion decomposition technique is introduced which combines the insight of traditional symbolic analysis with the numerical advantages of SPICE like simulators. The use of NTM leads to an efficient implementation. The proposed method greatly extends the size of the circuit and the complexity of the transistor model over what previous approaches could handle. For example, industry standard compact model, such as BSIM3V3 [35] was used for the first time in distortion analysis. The decomposition can be achieved at device, transistor and block level, all with device level accuracy. The theories have been implemented in a computer program and validated on examples. The proposed methods will leverage the performance of present VS based distortion analysis to the next level

Time-Varying Volterra Analysis of Nonlinear Circuits

Today’s advances in communication systems and VLSI circuits increases the performance requirements and complexity of circuits. The performance of RF and mixed-signal circuits is normally limited by the nonlinear behavior of the transistors used in the design. This makes simulation of nonlinear circuits more important. Volterra series is a method used for simulation of mildly nonlinear circuits. Using Volterra series the response of the nonlinear circuit is converted into a sum of multiple linear circuit responses. Thus, using Volterra series, simulation of nonlinear circuits in frequency-domain analysis becomes possible. However, Volterra series is not able to simulate strongly nonlinear circuits such as saturated Power Amplifiers. In this thesis, a new time-varying Volterra analysis is presented. The time-varying Volterra analysis is the generalization of conventional Volterra analysis where instead of using a DC expansion point a time-varying waveform has been used. Employing a time-varying expansion waveform for Volterra analysis, time-varying Volterra achieves better accuracy than conventional Volterra. The time-varying expansion waveforms are derived using a fast pre-analysis of the circuit. Using numerical examples, it has been shown that the time-varying Volterra is capable of simulating nonlinear circuits with better accuracy than conventional Volterra analysis. The time-varying Volterra analysis in both time and frequency domains are discussed in this thesis. The time-varying Volterra analysis has been used to simulate a saturated Class-F Power Amplifier in frequency-domain. The simulation results show good agreement with ELDO® steady-state and Harmonic Balance simulation results. The proposed method manages to simulate nonlinear circuits, such as saturated Power Amplifier, mixers and nonlinear microwave circuits, with good accuracy. Also, this method can be used to simulate circuit with large number of nonlinear elements without the convergence issues of Harmonic Balance

Spectral methods for circuit analysis

Thesis (Ph.D.)--Massachusetts Institute of Technology, Dept. of Electrical Engineering and Computer Science, 1999.Includes bibliographical references (p. 119-124).This electronic version was submitted by the student author. The certified thesis is available in the Institute Archives and Special Collections.Harmonic balance (HB) methods are frequency-domain algorithms used for high accuracy computation of the periodic steady-state of circuits. Matrix-implicit Krylov-subspace techniques have made it possible for these methods to simulate large circuits more efficiently. However, the harmonic balance methods are not so efficient in computing steady-state solutions of strongly nonlinear circuits with rapid transitions. While the time-domain shooting-Newton methods can handle these problems, the low-order integration methods typically used with shooting-Newton methods are inefficient when high solution accuracy is required. We first examine possible enhancements to the standard state-of-the-art preconditioned matrix-implicit Krylovsubspace HB method. We formulate the BDF time-domain preconditioners and show that they can be quite effective for strongly nonlinear circuits, speeding up the HB runtimes by several times compared to using the frequency-domain block-diagonal preconditioner. Also, an approximate Galerkin HB formulation is derived, yielding a small improvement in accuracy over the standard pseudospectral HB formulation, and about a factor of 1.5 runtime speedup in runs reaching identical solution error. Next, we introduce and develop the Time-Mapped Harmonic Balance method (TMHB) as a fast Krylov-subspace spectral method that overcomes the inefficiency of standard harmonic balance for circuits with rapid transitions. TMHB features a non-uniform grid and a time-map function to resolve the sharp features in the signals. At the core of the TMHB method is the notion of pseudo Fourier approximations. The rapid transitions in the solution waveforms are well approximated with pseudo Fourier interpolants, whose building blocks are complex exponential basis functions with smoothly varying frequencies. The TMHB features a matrix-implicit Krylov-subspace solution approach of same complexity as the standard harmonic balance method. As the TMHB solution is computed in a pseudo domain, we give a procedure for computing the real Fourier coefficients of the solution, and we also detail the construction of the time-map function. The convergence properties of TMHB are analyzed and demonstrated on analytic waveforms. The success of TMHB is critically dependent on the selection of a non-uniform grid. Two grid selection strategies, direct and iterative, are introduced and studied. Both strategies are a priori schemes, and are designed to obey accuracy and stability requirements. Practical issues associated with their use are also addressed. Results of applying the TMHB method on several circuit examples demonstrate that the TMHB method achieves up to five orders of magnitude improvement in accuracy compared to the standard harmonic balance method. The solution error in TMHB decays exponentially faster than the standard HB method when the size of the Fourier basis increases linearly. The TMHB method is also up to six times faster than the standard harmonic balance method in reaching identical solution accuracy, and uses up to five times less computer memory. The TMHB runtime speedup factor and storage savings favorably increase for stricter accuracy requirements, making TMHB well suited for high accuracy simulations of large strongly nonlinear circuits with rapid transitions.by Ognen J. Nastov.Ph.D

Multilevel harmonic balance analysis of large-scale nonlinear RF circuits via Newton-Krylov and tensor-Krylov methods

Orientador: Hugo Enrique Hernandez FigueroaTese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Eletrica e de ComputaçãoResumo: Este trabalho, tem como objetivo o desenvolvimento de novas técnicas, para análise de regime permanente não-autonoma de circuitos de alta-velocidade não-lineares em grande-escala. Para tal, é proposto um novo método do balanço harmônico (BH) fundamentado em uma eficiente metodologia de decomposição multi-níveis, que subdivide um circuito não-linear em grande escala em uma estrutura hierarquica de super-redes (SuRs) esparsamente interconectadas. Mais precisamente, em cada nível de hierarquia, o circuito é composto por SuRs intermediárias, SuRs de fundo, e redes de conexão (RCs). As SuRs de fundo são decompostas em um aglomerado de subredes não-lineares (SRNs) correspondendo a dispositivos semicondutores, que por sua vez, estão envolvidos por uma sub-rede linear (SRL). A equação de estado e de sonda das SuRs de fundo foram obtidas utilizando uma nova metodologia que combina a formulação de espaço de estado (FEE) para as SRNs com a formulação nodal modificada (FNM) para a SRL. Esta metodologia FEE/FNM produz um sistema quadrado de equações com menor tamanho possível. Para realização das conversões do sinal entre os domínios do tempo e da frequência, foram discutidas e implementadas diferentes transformadas de Fourier discreta (TFDs), para operação em regime multi-tons, incluindo sinais com modulação digital. A equação determinante do BH multi-níveis do circuito assume uma estrutura hierarquica do tipo bloco diagonal com borda , que pode ser eficientemente resolvida utilizando técnicas de processamento paralelo. A matriz jacobiana de cada SuR de fundo é processada utilizando eficientes técnicas de matrizes esparsas, junto com o conceito de espectro de derivada. Para a solução da equação determinante, foram utilizados os métodos de Newton e do tensor para problemas de pequena- e média-escala, e os métodos de Newton inexato e do tensor inexato para problemas em grande-escala. A globalização via pesquisa-em-linha com retrocedimento, foi adotada para nestes solucionadores não-lineares. Entretanto, para o método do tensor e do tensor inexato, também foi adotada a técnica de pesquisa-em-linha curvilinear. Nos métodos inexatos, técnicas de pré-condicionamento foram utilizadas, para aumentar a eficiência e a robustez do solucionador linear iterativo em subespaço de Krylov (GMRES, GMRES-Bt e TGMRES-Bt). Finalmente, a formulação proposta foi validada e a eficiência do método do tensor e do tensor inexato comparada com o método de Newton e de Newton inexato, para diferentes topologias de circuitos utilizando diodos, FETs e HBTs, e operando sob diferentes regimes de excitação multi-tons.Abstract: This work deals with the development of new techniques for nonautonomous nonlinear steady-state analysis of high-speed large-scale integrated circuits. To this end, it is proposed a novel harmonic balance (HB) method fundamented on a efficient multi-level decomposition methodology, that divides a large-scale circuit into hierarchical structure of sparsely interconnected supernetworks (SuNs). More precisely, the circuit is composed by intermediary SuRs, bottom SuRs and connection networks (CNs). The bottom SuNs are decomposed into a cluster of nonlinear subnetworks (NSNs) corresponding to the opto-electronic semiconductor devices, which in turn, are embedded by a linear subnetwork (LSN). Multi-port elements can be included in the LSN, in order to use measured data or results from electromagnetic analysis of structures with complex geometries. The formulation of the bottom SuN state and probe equations uses an improved table-oriented statespace formulation (SSF), that produces a square system with the lowest possible size, which is equal to the number of nonlinear state-variables (branch voltages and currents) that act as argument of the fuctions representing the semiconductor devices nonlinearities. The SSF is compared with the classical modified nodal formulation (MNF). For dealing with signal timefrequency conversions, discrete Fourier transform (DFT) techniques for different multi-tone regimes are discussed, including complex digitally modulated signals. The multi-level HB determining equation of the circuit assumes a hierarchical block bordered structure that can be efficiently tackled by parallel processing techniques. The HB jacobian matrix is handled using efficient sparse matrix techniques with a proper definition of the derivatives spectra. For the solution of a large-size HB problem, we investigated the applications of inexact tensor method based on Krylov-subspace techniques. Preconditioning are used to improve the robustness of the iterative tensor solver. To determine the circuit DC regime, we employ the tensor method. We adopted the backtracking linesearch technique as a globalisation strategy. However, for the tensor method, in particular, a curvilinear linesearch was also implemented. Finally, the formulation was validated and, the tensor and inexact tensor method efficiency was compared with the Newton and inexact Newton method, respectively, for several different circuits using diodos, FETs and HBTs, and operating under different multi-tone regimes.DoutoradoEngenharia de TelecomunicaçõesDoutor em Engenharia Elétric