例: 1 次のベッセル関数

関数 > 特殊関数 > ベッセル関数 > 例: 1 次のベッセル関数

例: 1 次のベッセル関数

関数J0、J1およびJnの関係を表します。関数とその尺度化された型との関係も表示します。

1. 2 つのレンジ変数を次のように定義します。

2. 関数 J0 と J1 をプロットします。2 次関数関数 Jn をプロットに次のように追加します。

<region id="ID0EF3E6" actualWidth="591.66857242139008" actualHeight="316.8" top="134.40000000000003" left="0" height="316.8" width="591.66857242139008" xmlns="http://schemas.mathsoft.com/worksheet50">
  <plot background-type="white" origin-positioning="true">
    <xyPlot>
      <title class="- topic/title " wwtype:type="Paragraph" xmlns:wwtype="urn:WebWorks-Type-Schema" />
      <legend />
      <traces>
        <trace resultRef="54">
          <traceStyle color="#FFED1D2F" symbol="none" line-weight="1" line-style="Solid">lines</traceStyle>
        </trace>
        <trace resultRef="55">
          <traceStyle color="#FF068149" symbol="none" line-weight="1" line-style="Dash">lines</traceStyle>
        </trace>
        <trace resultRef="56">
          <traceStyle color="#FF2E3192" symbol="none" line-weight="1" line-style="Dash">lines</traceStyle>
        </trace>
      </traces>
      <graph-size width="438.06857242139" height="240" />
      <axes>
        <xAxis rank="1" legend-position="PlotBoundaryBottom" start="-2" end="20">
          <axisLine position="origin" positionticmark="2" legendWidth="60.93" />
          <axisGrid>
            <gridFrequency>12</gridFrequency>
            <gridLabels display="true" />
            <gridLines />
            <tickMarks display="true" />
          </axisGrid>
          <axisLabel />
          <markers />
          <plotEquations>
            <plotEquation>
              <math>
                <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true" xmlns="http://schemas.mathsoft.com/math50">x</id>
              </math>
              <math>
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
          </plotEquations>
          <xyDomain scale-type="linear" auto-scale="true">
            <startValue resultRef="57">
              <real xmlns="http://schemas.mathsoft.com/math50">-2</real>
            </startValue>
            <secondTickValue resultRef="58">
              <real xmlns="http://schemas.mathsoft.com/math50">0</real>
            </secondTickValue>
            <endValue>
              <placeholder xmlns="http://schemas.mathsoft.com/math50" />
            </endValue>
          </xyDomain>
        </xAxis>
        <yAxis rank="1" legend-position="PlotBoundaryLeft" start="-0.5" end="1">
          <axisLine position="origin" positionticmark="1" legendWidth="113.39" />
          <axisGrid>
            <gridFrequency>7</gridFrequency>
            <gridLabels display="true" />
            <gridLines />
            <tickMarks display="true" />
          </axisGrid>
          <axisLabel />
          <markers />
          <plotEquations>
            <plotEquation>
              <math>
                <apply xmlns="http://schemas.mathsoft.com/math50">
                  <id labels="FUNCTION" xml:space="preserve">J0</id>
                  <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">x</id>
                </apply>
              </math>
              <math>
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
            <plotEquation>
              <math>
                <apply xmlns="http://schemas.mathsoft.com/math50">
                  <id labels="FUNCTION" xml:space="preserve">J1</id>
                  <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">x</id>
                </apply>
              </math>
              <math>
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
            <plotEquation>
              <math>
                <apply xmlns="http://schemas.mathsoft.com/math50">
                  <id labels="FUNCTION" xml:space="preserve">Jn</id>
                  <sequence>
                    <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">m</id>
                    <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">x</id>
                  </sequence>
                </apply>
              </math>
              <math>
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
          </plotEquations>
          <xyDomain scale-type="linear" auto-scale="true">
            <startValue resultRef="59">
              <real xmlns="http://schemas.mathsoft.com/math50">-0.5</real>
            </startValue>
            <secondTickValue resultRef="60">
              <real xmlns="http://schemas.mathsoft.com/math50">-0.25</real>
            </secondTickValue>
            <endValue>
              <placeholder xmlns="http://schemas.mathsoft.com/math50" />
            </endValue>
          </xyDomain>
        </yAxis>
      </axes>
    </xyPlot>
  </plot>
</region>

3. 次を示すプロットを作成します。J0(y)=Jn(0,y)目盛値をリセットして、X 軸を拡大してより詳細に表示します。

<region id="ID0EZ4E6" actualWidth="591.66857242139" actualHeight="297.6" top="460.80000000000007" left="0" height="297.6" width="591.66857242139" xmlns="http://schemas.mathsoft.com/worksheet50">
  <plot background-type="white" origin-positioning="true">
    <xyPlot>
      <title class="- topic/title " wwtype:type="Paragraph" xmlns:wwtype="urn:WebWorks-Type-Schema" />
      <legend />
      <traces>
        <trace resultRef="61">
          <traceStyle color="#FFED1D2F" symbol="none" line-weight="1" line-style="Solid">lines</traceStyle>
        </trace>
        <trace resultRef="62">
          <traceStyle color="#FF068149" symbol="circle" line-weight="1" line-style="None">lines</traceStyle>
        </trace>
      </traces>
      <graph-size width="438.06857242139" height="240" />
      <axes>
        <xAxis rank="1" legend-position="PlotBoundaryBottom" start="0" end="5">
          <axisLine position="origin" positionticmark="2" legendWidth="60.74" />
          <axisGrid>
            <gridFrequency>11</gridFrequency>
            <gridLabels display="true" />
            <gridLines />
            <tickMarks display="true" />
          </axisGrid>
          <axisLabel />
          <markers />
          <plotEquations>
            <plotEquation>
              <math>
                <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true" xmlns="http://schemas.mathsoft.com/math50">y</id>
              </math>
              <math>
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
          </plotEquations>
          <xyDomain scale-type="linear" auto-scale="true">
            <startValue>
              <placeholder xmlns="http://schemas.mathsoft.com/math50" />
            </startValue>
            <secondTickValue resultRef="63">
              <real xmlns="http://schemas.mathsoft.com/math50">0.5</real>
            </secondTickValue>
            <endValue>
              <placeholder xmlns="http://schemas.mathsoft.com/math50" />
            </endValue>
          </xyDomain>
        </xAxis>
        <yAxis rank="1" legend-position="PlotBoundaryLeft" start="-0.5" end="1">
          <axisLine position="origin" positionticmark="0" legendWidth="107.526666666667" />
          <axisGrid>
            <gridFrequency>7</gridFrequency>
            <gridLabels display="true" />
            <gridLines />
            <tickMarks display="true" />
          </axisGrid>
          <axisLabel />
          <markers />
          <plotEquations>
            <plotEquation>
              <math>
                <apply xmlns="http://schemas.mathsoft.com/math50">
                  <id labels="FUNCTION" xml:space="preserve">J0</id>
                  <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">y</id>
                </apply>
              </math>
              <math>
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
            <plotEquation>
              <math>
                <apply xmlns="http://schemas.mathsoft.com/math50">
                  <id labels="FUNCTION" xml:space="preserve">Jn</id>
                  <sequence>
                    <real>0</real>
                    <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">y</id>
                  </sequence>
                </apply>
              </math>
              <math>
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
          </plotEquations>
          <xyDomain scale-type="linear" auto-scale="true">
            <startValue resultRef="64">
              <real xmlns="http://schemas.mathsoft.com/math50">-0.5</real>
            </startValue>
            <secondTickValue resultRef="65">
              <real xmlns="http://schemas.mathsoft.com/math50">-0.25</real>
            </secondTickValue>
            <endValue>
              <placeholder xmlns="http://schemas.mathsoft.com/math50" />
            </endValue>
          </xyDomain>
        </yAxis>
      </axes>
    </xyPlot>
  </plot>
</region>

4. 次を示すプロットを作成します。J1(y)=Jn(1,y)目盛値をリセットして、X 軸を拡大してより詳細に表示します。

<region id="ID0EN6E6" actualWidth="591.66857242139008" actualHeight="297.6" top="883.20000000000016" left="0" height="297.6" width="591.66857242139008" xmlns="http://schemas.mathsoft.com/worksheet50">
  <plot background-type="white" origin-positioning="true">
    <xyPlot>
      <title class="- topic/title " wwtype:type="Paragraph" xmlns:wwtype="urn:WebWorks-Type-Schema" />
      <legend />
      <traces>
        <trace resultRef="66">
          <traceStyle color="#FFED1D2F" symbol="none" line-weight="1" line-style="Solid">lines</traceStyle>
        </trace>
        <trace resultRef="67">
          <traceStyle color="#FF2E3192" symbol="circle" line-weight="1" line-style="None">lines</traceStyle>
        </trace>
      </traces>
      <graph-size width="438.06857242139" height="240" />
      <axes>
        <xAxis rank="1" legend-position="PlotBoundaryBottom" start="0" end="5">
          <axisLine position="origin" positionticmark="2" legendWidth="60.74" />
          <axisGrid>
            <gridFrequency>11</gridFrequency>
            <gridLabels display="true" />
            <gridLines />
            <tickMarks display="true" />
          </axisGrid>
          <axisLabel />
          <markers />
          <plotEquations>
            <plotEquation>
              <math>
                <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true" xmlns="http://schemas.mathsoft.com/math50">y</id>
              </math>
              <math>
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
          </plotEquations>
          <xyDomain scale-type="linear" auto-scale="true">
            <startValue>
              <placeholder xmlns="http://schemas.mathsoft.com/math50" />
            </startValue>
            <secondTickValue resultRef="68">
              <real xmlns="http://schemas.mathsoft.com/math50">0.5</real>
            </secondTickValue>
            <endValue>
              <placeholder xmlns="http://schemas.mathsoft.com/math50" />
            </endValue>
          </xyDomain>
        </xAxis>
        <yAxis rank="1" legend-position="PlotBoundaryLeft" start="-0.5" end="1">
          <axisLine position="origin" positionticmark="0" legendWidth="107.526666666667" />
          <axisGrid>
            <gridFrequency>7</gridFrequency>
            <gridLabels display="true" />
            <gridLines />
            <tickMarks display="true" />
          </axisGrid>
          <axisLabel />
          <markers />
          <plotEquations>
            <plotEquation>
              <math>
                <apply xmlns="http://schemas.mathsoft.com/math50">
                  <id labels="FUNCTION" xml:space="preserve">J1</id>
                  <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">y</id>
                </apply>
              </math>
              <math>
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
            <plotEquation>
              <math>
                <apply xmlns="http://schemas.mathsoft.com/math50">
                  <id labels="FUNCTION" xml:space="preserve">Jn</id>
                  <sequence>
                    <real>1</real>
                    <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">y</id>
                  </sequence>
                </apply>
              </math>
              <math>
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
          </plotEquations>
          <xyDomain scale-type="linear" auto-scale="true">
            <startValue resultRef="69">
              <real xmlns="http://schemas.mathsoft.com/math50">-0.5</real>
            </startValue>
            <secondTickValue resultRef="70">
              <real xmlns="http://schemas.mathsoft.com/math50">-0.25</real>
            </secondTickValue>
            <endValue resultRef="71">
              <real xmlns="http://schemas.mathsoft.com/math50">1</real>
            </endValue>
          </xyDomain>
        </yAxis>
      </axes>
    </xyPlot>
  </plot>
</region>

5. シンボリック評価を使用して、1 次のベッセル関数とその尺度化された型の関係を次のように表示します。

6. 次を示すプロットを作成します。

J0.sc には複素数の成分が含まれているので、最初にRe関数を適用して数字の実部のみを表示します。

<region id="ID0EOEF6" actualWidth="534.66857242138985" actualHeight="297.6" top="1401.6000000000001" left="0" height="297.6" width="534.66857242138985" xmlns="http://schemas.mathsoft.com/worksheet50">
  <plot background-type="white" origin-positioning="true">
    <xyPlot>
      <title class="- topic/title " wwtype:type="Paragraph" xmlns:wwtype="urn:WebWorks-Type-Schema" />
      <legend />
      <traces>
        <trace resultRef="76">
          <traceStyle color="#FFED1D2F" symbol="none" line-weight="1" line-style="Solid">lines</traceStyle>
        </trace>
        <trace resultRef="77">
          <traceStyle color="#FF662D91" symbol="circle" line-weight="1" line-style="None">lines</traceStyle>
        </trace>
      </traces>
      <graph-size width="381.06857242139" height="240" />
      <axes>
        <xAxis rank="1" legend-position="PlotBoundaryBottom" start="0" end="5">
          <axisLine position="origin" positionticmark="2" legendWidth="60.74" />
          <axisGrid>
            <gridFrequency>11</gridFrequency>
            <gridLabels display="true" />
            <gridLines />
            <tickMarks display="true" />
          </axisGrid>
          <axisLabel />
          <markers />
          <plotEquations>
            <plotEquation>
              <math>
                <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true" xmlns="http://schemas.mathsoft.com/math50">y</id>
              </math>
              <math>
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
          </plotEquations>
          <xyDomain scale-type="linear" auto-scale="true">
            <startValue>
              <placeholder xmlns="http://schemas.mathsoft.com/math50" />
            </startValue>
            <secondTickValue resultRef="78">
              <real xmlns="http://schemas.mathsoft.com/math50">0.5</real>
            </secondTickValue>
            <endValue resultRef="79">
              <real xmlns="http://schemas.mathsoft.com/math50">5</real>
            </endValue>
          </xyDomain>
        </xAxis>
        <yAxis rank="1" legend-position="PlotBoundaryLeft" start="-0.5" end="1">
          <axisLine position="origin" positionticmark="0" legendWidth="108.586666666667" />
          <axisGrid>
            <gridFrequency>7</gridFrequency>
            <gridLabels display="true" />
            <gridLines />
            <tickMarks display="true" />
          </axisGrid>
          <axisLabel />
          <markers />
          <plotEquations>
            <plotEquation>
              <math>
                <apply xmlns="http://schemas.mathsoft.com/math50">
                  <id labels="FUNCTION" xml:space="preserve">J0</id>
                  <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">y</id>
                </apply>
              </math>
              <math>
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
            <plotEquation>
              <math>
                <apply xmlns="http://schemas.mathsoft.com/math50">
                  <id labels="FUNCTION" xml:space="preserve">J0.sc</id>
                  <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">y</id>
                </apply>
              </math>
              <math>
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
          </plotEquations>
          <xyDomain scale-type="linear" auto-scale="true">
            <startValue resultRef="80">
              <real xmlns="http://schemas.mathsoft.com/math50">-0.5</real>
            </startValue>
            <secondTickValue resultRef="81">
              <real xmlns="http://schemas.mathsoft.com/math50">-0.25</real>
            </secondTickValue>
            <endValue>
              <placeholder xmlns="http://schemas.mathsoft.com/math50" />
            </endValue>
          </xyDomain>
        </yAxis>
      </axes>
    </xyPlot>
  </plot>
</region>

7. J1 の最初の 2 つのピークの座標を計算します。関数augmentおよびlocalmaxを使用して、指定した範囲内にあるピークを特定します。

関数 localmax には 2 列の入力行列が必要です。関数 augment を使用してこのような行列を作成します。

8. 関数matchを使用してピークの横方向の座標を求めます。最も正確な結果を得るためには、TOL の値を小さくします。

9. プロットにマーカーを追加して、最初の 2 つのピークをマーク付けします。

レンジ変数のステップは 0.1 です。これは、関数 match で識別される要素の添字の 1/10、つまり 1.8 と 8.5 でそれぞれピークが発生することを示します。

<region id="ID0E4PF6" actualWidth="590.46857242139" actualHeight="316.79999999999995" top="2246.4000000000005" left="0" height="316.79999999999995" width="590.46857242139" xmlns="http://schemas.mathsoft.com/worksheet50">
  <plot background-type="white" origin-positioning="true">
    <xyPlot>
      <title class="- topic/title " wwtype:type="Paragraph" xmlns:wwtype="urn:WebWorks-Type-Schema" />
      <legend />
      <traces>
        <trace resultRef="94">
          <traceStyle color="#FF662D91" symbol="none" line-weight="1" line-style="Solid">lines</traceStyle>
        </trace>
      </traces>
      <graph-size width="456.06857242139" height="259.2" />
      <axes>
        <xAxis rank="1" legend-position="PlotBoundaryBottom" start="-2" end="20">
          <axisLine position="origin" positionticmark="2" legendWidth="60.93" />
          <axisGrid>
            <gridFrequency>12</gridFrequency>
            <gridLabels display="true" />
            <gridLines />
            <tickMarks display="true" />
          </axisGrid>
          <axisLabel />
          <markers>
            <marker resultRef="95" display="true" color="#FFED1D2F" line-style="Dash" value="1.8">
              <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">X1</id>
            </marker>
            <marker resultRef="96" display="true" color="#FF068149" line-style="Dash" value="8.5">
              <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">X2</id>
            </marker>
          </markers>
          <plotEquations>
            <plotEquation>
              <math>
                <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true" xmlns="http://schemas.mathsoft.com/math50">x</id>
              </math>
              <math>
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
          </plotEquations>
          <xyDomain scale-type="linear" auto-scale="true">
            <startValue resultRef="97">
              <real xmlns="http://schemas.mathsoft.com/math50">-2</real>
            </startValue>
            <endValue>
              <placeholder xmlns="http://schemas.mathsoft.com/math50" />
            </endValue>
          </xyDomain>
        </xAxis>
        <yAxis rank="1" legend-position="PlotBoundaryLeft" start="-0.5" end="0.75">
          <axisLine position="origin" positionticmark="1" legendWidth="90.8833333333334" />
          <axisGrid>
            <gridFrequency>6</gridFrequency>
            <gridLabels display="true" />
            <gridLines />
            <tickMarks display="true" />
          </axisGrid>
          <axisLabel />
          <markers>
            <marker resultRef="98" display="true" color="#FFED1D2F" line-style="Dash" value="0.58151695173116513">
              <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">P1</id>
            </marker>
            <marker resultRef="99" display="true" color="#FF068149" line-style="Dash" value="0.27312196367405378">
              <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">P2</id>
            </marker>
          </markers>
          <plotEquations>
            <plotEquation>
              <math>
                <apply xmlns="http://schemas.mathsoft.com/math50">
                  <id labels="FUNCTION" xml:space="preserve">J1</id>
                  <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">x</id>
                </apply>
              </math>
              <math>
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
          </plotEquations>
          <xyDomain scale-type="linear" auto-scale="true">
            <startValue resultRef="100">
              <real xmlns="http://schemas.mathsoft.com/math50">-0.5</real>
            </startValue>
            <secondTickValue resultRef="101">
              <real xmlns="http://schemas.mathsoft.com/math50">-0.25</real>
            </secondTickValue>
            <endValue>
              <placeholder xmlns="http://schemas.mathsoft.com/math50" />
            </endValue>
          </xyDomain>
        </yAxis>
      </axes>
    </xyPlot>
  </plot>
</region>