4.10: Numerical Experiment (Circuit Design)
- Page ID
- 10118
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\( \newcommand{\dsum}{\displaystyle\sum\limits} \)
\( \newcommand{\dint}{\displaystyle\int\limits} \)
\( \newcommand{\dlim}{\displaystyle\lim\limits} \)
\( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)
( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\id}{\mathrm{id}}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\kernel}{\mathrm{null}\,}\)
\( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\)
\( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\)
\( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)
\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)
\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)
\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vectorC}[1]{\textbf{#1}} \)
\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)
\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)
\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\(\newcommand{\longvect}{\overrightarrow}\)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\(\newcommand{\ket}[1]{\left| #1 \right>}\)
\(\newcommand{\bra}[1]{\left< #1 \right|}\)
\(\newcommand{\braket}[2]{\left< #1 \vphantom{#2} \right| \left. #2 \vphantom{#1} \right>}\)
\(\newcommand{\braopket}[3]{\left< #1 \vphantom{#2}\vphantom{#3} \right| #2 \vphantom{#1}\vphantom{#3} \left| #3 \vphantom{#1}\vphantom{#2} \right>}\)
\(\newcommand{\qmvec}[1]{\mathbf{\vec{#1}}}\)
\(\newcommand{\op}[1]{\hat{\mathbf{#1}}}\)
\(\newcommand{\expect}[1]{\langle #1 \rangle}\)
\(\newcommand{\dfn}[1]{\emph{\textbf{#1}}}\)
Our analysis in Example 1 from "Linear Algebra: Circuit Analysis" and Problem 1 from "Linear Algebra: Circuit Analysis" indicates that not enough current will flow through the lamp to make it glow. We now wish to change the resistance of the 100 ohm resistor to a new value so that the lamp will glow. We replace 100 in the equations with an unknown resistance \(R\). Equation 2 is unchanged, but Equation 6 becomes
\[\begin{align}
& \tfrac{\left(v_{2}-v_{1}\right)}{50}+\tfrac{v_{2}}{300}+\tfrac{\left(v_{2}-v_{3}\right)}{R}=0 \nonumber \\
\Rightarrow & 6 R\left(v_{2}-v_{1}\right)+R v_{2}+300\left(v_{2}-v_{3}\right)=0 \nonumber \\
\Rightarrow &-6 R v_{1}+(7 R+300) v_{2}-300 v_{3}=0 \label{}.
\end{align} \nonumber \]
Equation 7 from "Linear Algebra: Circuit Analysis" becomes
\[\begin{align}
\tfrac{\left(v_{3}-v_{2}\right)}{R}+\tfrac{v_{3}}{2}=0 \nonumber \\
\Rightarrow \quad 2\left(v_{3}-v_{2}\right)+R v_{3}=0 \nonumber \\
\Rightarrow \quad 0 v_{1}-2 v_{2}+(R+2) v_{3}=0 \label{}
\end{align} \nonumber \]
The matrix form of these equations is
\[\left[\begin{array}{lll}
1 & & 0 & 0 \\
-6 R & 7 R & +300 & -300 \\
0 & & -2 & R+2
\end{array}\right]\left[\begin{array}{l}
v_{1} \\
v_{2} \\
v_{3}
\end{array}\right]=\left[\begin{array}{l}
5 \\
0 \\
0
\end{array}\right] \nonumber \]
Write a MATLAB function file called builda to accept \(R\) as an input and return the matrix \(A\) in Equation 3 as an output. The first line of your function file should be
>> function A = builda(R);
Now choose several values for \(R\). For each choice, use your function builda and solve the resulting matrix equation \(A\nu=b\) for the voltages. Each time you choose a different \(R\) to build a different matrix \(A\), check the determinant of A to be sure the equations have a unique solution:
>> det(A)
Make a table of \(R\) and the corresponding values of \(\nu _3\):
Now add a column to your table for the current through the lamp \(i=\nu _3/2\). Add rows to your table until you have found a value of \(R\) for which the lamp will glow. (\(i\) needs to be between 0.05 and 0.075 ampere.)