Skip to main content
Library homepage
 

Text Color

Text Size

 

Margin Size

 

Font Type

Enable Dyslexic Font
Mathematics LibreTexts

1.2: Chapter 1 Exercises

( \newcommand{\kernel}{\mathrm{null}\,}\)

Question 1.2.1

In order to refresh your matrix-vector multiply skills please calculate, by hand, the product ATGA in the 3 compartment case and write out the 4 equations in the vector equation we arrived at in step(S4): ATGAx=f

Feedback

The second equation should read

x1+2x2x3Ri+x2Rm=0


Question 1.2.2

We began our discussion with the 'hope' that a multicompartment model could indeed adequately capture the fiber's true potential and current profiles. In order to check this one should run fib1.m with increasing values of NN until one can no longer detect changes in the computed potentials.

  • (a) Please run fib1.m with N=8,16,32, and 64. Plot all of the potentials on the same (use hold) graph, using different line types for each. (You may wish to alter fib1.m so that it accepts NN as an argument).

Let us now interpret this convergence. The main observation is that the difference Equation ???, approaches a differential equation. We can see this by noting that

d(z)=lN

acts as a spatial 'step' size and that xkd(z) is approximately the value of the true potential at (k1)d(z). In a slight abuse of notation, we denote the latter

x((k1)d(z))

Applying these conventions to Equation ??? and recalling the definitions of Ri and Rm we see Equation ??? become

πa2ρix(0)+2x(d(z))x(2d(z))d(z)+2πad(z)ρmx(d(z))=0

or, after multiplying through by ρmπad(z)

aρmρix(0)+2x(d(z))x(2d(z))d(z2)+2x(d(z))=0

We note that a similar equation holds at each node (save the ends) and that as N and therefore d(z)0 we arrive at

d2dz2x(z)2ρiaρmx(z)=0

  • (b) With μ2ρiaρm show that

x(z)=αsinh(2μz)+βcosh(2μz)

satisfies Equation ??? regardless of α and β

We shall determine α and β by paying attention to the ends of the fiber. At the near end we find

πa2ρix(0)x(d(z))d(z)=i0

which, as d(z)0 becomes

ddzx(0)=ρii0πa2

At the far end, we interpret the condition that no axial current may leave the last node to mean

ddzx(l)=0

  • (c) Substitute Equation ??? into Equation ??? and solve for α and β and write out the final x(z).
  • (d) Substitute into x the l,a,ρi,ρm values used in fib1.m, plot the resulting function (using, e.g., ezplot) and compare this to the plot achieved in part (a).

This page titled 1.2: Chapter 1 Exercises is shared under a CC BY 1.0 license and was authored, remixed, and/or curated by Steve Cox via source content that was edited to the style and standards of the LibreTexts platform.

Support Center

How can we help?