ECSE343 Project Solved

In this project your task is to compute the DC solution of the circuit shown in Figure 1.  

  

Figure 1: Full wave rectifier

We can write the nodal current equations of the above circuit as,

                                                 𝑭(𝑿) = 𝑮 𝑿 + 𝒈(𝑿) − 𝑼 = 𝟎                                          (1)

 

where  𝑿 is the unknown vector, it consists of nodal voltages, and the current.   

 𝑼  is the source vector, it contains the contribution of input sources. The matrix  𝑮 mainly contains the contribution of the resistor. The vector 𝒈(𝑿)  consists of the diode currents and it is dependent on the vector 𝑿.

The equation (1) can be further elaborated as,

                                                                     𝐼𝑆 (𝑒𝑉1𝑉−𝑡𝑉2 −1)−𝐼𝑆 (𝑒−𝑉𝑉𝑡1 −1)                                            

0 0 0 1

0 0 0 −1

 𝑭(𝑿)=⏟[10   0 .002   𝑮 −0 1    00 ][⏟𝑉𝑉𝑉𝐼𝑿𝐸132]+        −𝐼𝐼𝑆𝑆  ((𝑒𝑒𝑉𝑉31𝑉𝑉−−𝑡𝑡𝑉𝑉22 −−11))−−𝐼𝐼𝑆𝑆  ((𝑒𝑒−𝑉𝑉3𝑉𝑡𝑉−3𝑡𝑉−2 −1)1)      −⏟[𝑉 𝑼000𝑖𝑛  ]=0                                                                       (  2)

[⏟            0            ]

𝒈(𝑿)

 

The system, 𝑭(𝑿), described in (2) contains nonlinear nodal equations of the circuit shown in Figure 2. Your task in this project is to implement a nonlinear solver to compute 𝑿. 

 

 

For your convenience you are provided with a Matlab function named nonlinearFunc() it takes 𝑿 as an input and provides 𝑭(𝑿) and 𝝏𝑭(𝑿). Note that the current function is hard-coded to compute

𝝏𝑿

𝑭(𝑿) and the jacobian for Vin  = 10V. You may edit this function as needed in order to achieve your objectives.

PART 1 
Design and implement a Matlab program that computes the solution for the output voltage V2. For this part you may use built-in Matlab functions of your choice. Verify if the solutions obtained are correct, you may either use SPICE simulator or verify the solution analytically. Justify the choice of functions used as well as the choice of parameters and options used when running these functions.

 

PART 2 
For this part implement a solver using an appropriate algorithm or combination of algorithms. Your algorithm should function without any meaningful initial guess provided by the user (i.e. it should determine its own initial guess or simply work for a trivial initial guess of all zeros). You may use the provided nonlinearFunc(). Do not hesitate to modify the nonlinearFunc method as per your requirements.  

Note: You do not need to implement every algorithm use. If any of the following algorithms are needed you may use matlab functions/operators: ‘\’ (backslash),  lu() (LU decomposition), qr() (QR decomposition), chol() (Cholesky decomposition), norm(), and Forward/ Backward Substitution (provided in assignment 2).

PART 3 
Design and implement a program that computes and plots the output V2 if Vin is a 60Hz sine wave with amplitude 5V. You may use any Matlab function of your choice including ones that you may have developed in this course.