方程式の解法

Size: px

Start display at page:

Download "方程式の解法"

ふみなひろなが
5 years ago
Views:

1 方程式の解法方程式 (f(x)) の解を求めるということはf(x)=0を意味しグラフ上では y=f(x) で X 軸との交点を求めることであるここではエクセルを用いて方程式の解を求める方法を3つ紹介するもちろん前回教わったニュートンラフソン法も行う 1. グラフを確認しながら求める方法エクセルの便利な点は数値を簡単にグラフ化できる点である中学校の数学でも教わったように方程式の解は X 軸と交差する点であるこれをグラフを確認しながら決定してみる次の方程式を解いてみるこの方程式を解くために y xlog e x 1のyの値の変化を見てみる xlog e x 1 0 図 1 A1-A18 に X 軸の値として 0.1 から 16 までの値を入れてみる B1 に =A1*LN(A1)-1 を入力することで y の値を計算する B1 を B2 から B18 にコピーする A1:B18 を選択してグラフ作成のアイコンを押し散布

2 図 2 図を選択し完了にするとグラフを作成することができるできた図 1 のグラフからわかるように f(x) の値は 2 付近で X 軸と交差していることがわかる実際の値を見てみると x=1 の時には負の値であるが x=2 では正の値になっているすなわち 1 と 2 の間に X 軸と交差していることがわかる X の値を 1 と 2 の間で変化させて同様に計算しグラフを作成してみる図 2 では X 軸との交点は丁度 1.7 と 1.8 の間であることがわかる E E-06 これを繰り返すと上の表のように 0 になる X を求めることができるこの場合には x= が解である 2. ニュートンラフソン法で求める ( 表計算 ) 次にニュートンラフソン法を用いてセル中で計算させて求めることを考えるその例が図 4 である

この例では A 列は計算回数である A2 には 0 が代入してあり A3 は =A2+1 として計算回数を表示するようにしてある B2 が初期値であるこの場合は 1 を初期値としている B3 はニュートンラフソン法の公式を代入してある x i 1 x i f ( xi ) xi ln xi 1 xi f '( x ) ln x 1 i 上式に基づいて B3 には

3 この例では A 列は計算回数である A2 には 0 が代入してあり A3 は =A2+1 として計算回数を表示するようにしてある B2 が初期値であるこの場合は 1 を初期値としている B3 はニュートンラフソン法の公式を代入してある x i 1 x i f ( xi ) xi ln xi 1 xi f '( x ) ln x 1 i 上式に基づいて B3 には =B2-(B2*LN(B2)-1)/(LN(B2)+1) としてある xi=b2 であるので B2 を用いて B3=xi+1 を計算しているこれを繰り返すために B4 以降にコピーすると繰り返し計算ができる C 列には xi+1-xi が表示されているこの結果を見ると 5 回目で誤差が 10-6 以下になっているそのときの解はグラフを見ながら求めた結果とよく一致している i

4 3. ニュートンラフソン法で求める (VBA) まず新しいワークシートへ移動しますそしてコマンドボタンを作成し名前をニュートンラフソン法に変更しましょうまた, シートの名前を方程式の開放に変更します. このボタンを押した状態でコマンドボタンを作成する作成したコマンドボタンの名前をニュートンラフソン法に変更するシートの名称を方程式の解法とするマクロの登録のウィンドウが開くので新規作成を押し, 起動したエディタで次の

5 プログラムを入力する Private Sub CommandButton1_Click() End Sub 非線形方程式.Show 次に VBA エディターで挿入 - ユーザーフォームで新しいユーザーフォームを作成しラベルテキストボックスコマンドボタンを利用して条件の入力画面を作成するラベルテキストボックスコマンドボタンさらにこのユーザーフォームのオブジェクト名を非線形方程式とするこのオブジェクト名の変更はエクセルの表中に作成したコマンドボタンで呼び出すプログラムに非線形方程式.show と記述したがここで作成したユーザーフォームを呼び出すためである次に計算実行のボタンをダブルクリックしてコード画面を表示しコードを入力するコード入力の際にテキストボックスの値を参照するこの場合テキストボックスの番号が重要であるこの例では以下のように対応している初期値 = : テキストボックス1 閾値 ( しきいち )= : テキストボックス2 最大反復回数 = : テキストボックス3 計算実行 : コマンドボタン1 終了 : コマンドボタン2

6 もし各自のテキストボックスの番号が異なる場合にはそれに対応してプログラムを読み直してほしいコード画面が現れたら次のプログラムを入力する

7 Private Sub CommandButton1_Click() Dim X1 As Double Dim eps As Double Dim ff As Double Dim df As Double Dim x As Double Dim max As Integer Dim No As Integer 1 Worksheets(" 方程式の解法 ").Activate Set WS = Worksheets(" 方程式の解法 ").Application WS.Cells(34, 2) = " 初期値 X1=" WS.Cells(34, 3) = TextBox1.Text X1 = CDbl(TextBox1.Text) WS.Cells(35, 2) = " 閾値 =" WS.Cells(35, 3) = TextBox2.Text eps = CDbl(TextBox2.Text) WS.Cells(36, 2) = " 最大反復回数 =" WS.Cells(36, 3) = TextBox3.Text max = CDbl(TextBox3.Text) No = 0 Do If No >= max Then Exit Do x = X1 No = No + 1 ff = x * Log(x) - 1 df = Log(x) + 1 X1 = x - (ff / df) WS.Cells(39 + No, 2) = "X1=" WS.Cells(39 + No, 3) = X1 WS.Cells(39 + No, 5) = No WS.Cells(39 + No, 7) = x - X1 Loop While (Abs(x - X1) > eps) End Sub WS.Cells(37, 2) = " 方程式の解 =" WS.Cells(37, 3) = X1 6

8 プログラムを解説する前にプログラムの流れを解説しておくニュートンラフソン法での解の計算は下記のような流れにより計算する変数定義 1 初期値設定 (x 0 閾値最大反復回数など ) 2,3,4 反復回数 (No) 最大反復回数 (max) Yes No xi+1 の計算 5 xi+1 の表示 No 閾値 (eps) xi+1-xi Yes 解の表示 6 終了

9 プログラムの解説 1 変数の定義プログラム中で用いる変数 ( 数値を入れておくメモリ ) は必ず定義されなければならないここでは2 種類の変数が定義されている Dim X1 As Double アイエオ変数の定義は上記のように記述されるア : 変数の定義を行いますという宣言イ : 定義する変数名変数名は自由につけることができるがアルファベットと数字および _( アンダーライン ) でその変数が何を意味するかを考えてつけるようにすると良いただし頭文字に数字やアンダーバーは使うことができない漢字やひらがなカタカナなども使うことができるが別のコンピュータに持って行った場合にエラーを起こす可能性があるので使うべきではないエ : イの変数をどの様な変数として定義するかという意味オ : 変数の種類変数の種類には次のようなものがあるデータの型データの特徴データの例 Integer 整数 ( 小数点のつかな 0,5, 100 い数値 ) Long 桁数の大きい整数 , (Integer も含む ) Single 実数 ( 小数点を含む数 -2.0, 12.8 値 ) Double 桁数がかなり大きい実数 (Single も含む ) , String 文字カレーうどん大学変数の型としては上記のもの以外にも Date, Byte 等があるここでは必要ないので割愛するプログラムで定義する変数は常にプログラムの一番最初で定義するそしてこの変数の有効範囲は定義した部分から End Sub まであるすなわち Private End Sub がその変数が有効な範囲であるプログラムではそれぞれのコマンドボタンやテキストボックスなどにそれぞれプログラムを記述することができるその単位をプロシージャと呼ぶが通常変数は各プロシージャ内でのみ有効であるということであるつまり別のプロシージャで同じ変数を定

10 義してもかまわないということである 2 利用するワークシートの定義 Worksheets(" 方程式の解法 ").Activate このコマンドは利用するワークシートを定義していますすなわち実際の操作で考えるとワークシートのタグをマウスでクリックすることです上記の例では方程式の解法と名付けたワークシートが選択されていますさらに選択したワークシートの名前を WS というな変数に代入していますこの変数は 1 では定義されていませんが実際にはオブジェクト変数というもので代入することで定義されていると考えてください 3 初期値の代入ここでは初期値閾値最大反復回数を各セルおよび変数に代入しています WS.Cells(34, 2) = " 初期値 X1=" 先頭の WS は用いるワークシートを示していますこの命令では 34 行目の 2 列目のセルに " 初期値 X1=" という文字列を代入しています前の授業では Range という命令を使いました Range の場合にはセルの位置の指定に A1 や B2 等の表現を用いましたしかしプログラムで行の位置や列の位置を変化させたい場合にはこれでは不便です Cells のセル指定では Cells( 行位置を示す数値列位置を示す数値 ) という指定の仕方をしますセル位置をすべて数値で表すことができるのでプログラム中でも扱いやすくなります具体的には 5 の中で出てくる 39+No のように演算形式で表現することもできます WS.Cells(34, 3) = TextBox1.Text この命令では 34 行目 3 列目のセルに TextBox1.Text の内容を代入していますつまり前の命令の右横のセルに TextBox1.Text の内容つまり初期値の値を代入しています X1 =CDbl(TextBox1.Text) さらに次の命令では X1 という変数に初期値を代入しています CDbl() は文字を数値 ( 倍精度 ) に変換する関数ですテキストボックスに入力された数字は文字として認識されます変数に文字を代入する場合には必ず数値に変換してください同様に 35,36 行目に閾値と最大反復回数を代入し eps と max にも代入しています

11 4 反復回数の初期化 No = 0 この命令では No という整数変数に0を代入していますこの変数は5で反復回数を数えるために使います計算が始まる前には必ず0でなければならないのでここで0にしています 5 ニュートンラフソン法の計算この部分が計算の本体です計算のフローチャートを見てくださいここに相当する部分は動作が複雑です反復回数が最大反復回数より大きいか小さいかの判断をする部分がありさらに計算値と前回の計算値の差が閾値より大きいかどうかの判断をする部分がありますさらに同じ計算を繰り返し行う動作をしなければなりません条件判断条件判断には次の命令を使います If 条件式 1 Then ElseIf 条件式 2 Then ElseIf 条件式 3 Then 実行文群 1 実行文群 2 実行文群 2 Else 実行文群 n End IF 判断文は条件式が満足されている場合 ( 真 ) にすぐ次の実行文を実行しそれ以外は実行せずに EndIf に飛び If 文を終了します

12 具体的には次のようになります If x<0 Then MsgBox x は負の数です ElseIf x<10 Then MsgBox x は 0 以上 10 未満の数です ElseIf x<20 Then Else End If MsgBox x は 10 以上 20 未満の数です MsgBox x は 20 以上の数ですこのプログラムで MsgBox というのが出ていますがこれはこの命令に続く文章を画面に表示する命令です x=-1 の場合この場合には最初の文の条件式 (x<0) が満足されますそのため画面に x は負の数ですと表示して終了します x=12 の場合この場合には x<0 は満足しませんからその次の実行文 (MsgBox x は負の数ですは実行しませんその次の条件式 (x<10) も満足しませんそして 3 番目の条件式 (x<20) は満足しますので画面に x は 10 以上 20 未満の数ですと表示して終了します x=35 の場合この場合には 3 つの条件式 (x<0, x<10, x<20) のどれも満足しませんそのため一番最後の Else の実行文 ( x は 20 以上の数ですを表示 ) を実行します

13 条件式に用いられる演算記号は次のようなものがあります記号意味 = 左辺と右辺が等しい > 左辺が右辺より大きい < 左辺が右辺未満 >= 左辺が右辺以上 <= 左辺が右辺以下 <> 左辺と右辺が等しくないまた条件式を複数つなぐ論理演算子 (And, Or, Not 等 ) も使うことができます繰り返し今回行うニュートンラフソン法では xi+1=xi-f(xi)/f (xi) を求めることにより数値解を求めますたとえば初期値 x0=1 から x1 を求めますさらに x1 を元に x2 を求めますさらに x2 x3 x4 x5 のように漸化式として一つ前の値から次の値を求めますこのような計算に対してプログラムする際にはどの様にするのでしょう x0 x1 に対刷るプログラムを書き x1 x2 に対するプログラムを書きさらにではどこまでプログラムを書けばよいのでしょう方程式によってはかなりの回数計算しなければならない場合があるはずですし簡単に解を得られる場合もありますこれまで考慮してプログラムをするには最初に自分で計算してみてプログラムを書く必要があり個々の方程式について別々のプログラムになってしまい効率が良くありませんはっきり言ってそんなことならプログラムを書く必要はないはずですあらゆるプログラム言語には繰り返し命令というのがあります VBA でも同じですコンピュータにとってあるルール ( プログラム ) に沿って繰り返しを行うことは得意なことです VBA では 3 種類の繰り返し命令があります Do.Loop ステートメント Do.Loop ステートメントは Do と Loop で囲まれた範囲の命令をある条件の間繰り返し実行するための命令ですその構文には 4 種類あります Do While 条件式 1 繰り返したい命令文 Loop Do Until 条件式繰り返したい命令文 Loop 3 Do 繰り返したい命令文 Loop While 条件式 2 Do 繰り返したい命令文 Loop Until 条件式 4

14 1~4の違いは2 種類ありますひとつは条件式の位置です条件式が Do の部分に記述されているもの (1と3) と Loop の部分に記述されているもの (2と4) ですこの2つの違いは条件判定をする位置です Do Loop ステートメントでは条件判定が行われその判定が True( 真 ) の場合には繰り返しを行い偽 (False) の場合には Do Loop ステートメントを終了しますそのため下のように Do の部分に条件判定がある場合には最初から判定が偽の場合には1 回もその繰り返しが実行されません x=-1 判定が偽のため 1 Do While x>0 Do Loop は実行されな 2 x=x-1 Loop x=-1 が表示される Debug.Print x これに対して Loop の部分に条件判定がある場合には最低でも1 回はその繰り返しが実行されます x=-1 最低でも1 回は実行さ 1 Do れる 2 x=x-1 Loop While x>0 x=-2 が表示される Debug.Print x もう一つの違いは条件判定の部分の While 1 2 と Until 3 4 ですこれは英語の意味と同じです While : 条件式が真の間繰り返しを実行する Untile : 条件式が真になったら繰り返しを終了する条件式以外でこの繰り返しを終了するには Exit Do という命令がありしばしば if 文と併せて用いられます For Next ステートメント For Nest ステートメントも Do Loop ステートメントと同様に繰り返しの命令です記述例を示します For カウンタ変数名 = 初期値 To 最終値 Step 増分繰り返すべき命令群 Next カウンタ変数名 For Next ステートメントでは繰り返し回数がわかっている場合に用いますカウンタ変数名というのは何回繰り返したかを判定するために用いる変数です

15 ( 必ずしも繰り返し回数とは一致しない ) 上の命令文の意味はカウンタ変数を初期値から初めて増分ずつ増加させ最終値になるまで繰り返せということを意味しています例を下に示します 1 2 x=0 For i=1 to 10 step 1 x=x+i next I MsgBox x i x(1) x(2) この例では i の値が 1 から1ずつ増加させ 10 になるまで繰り返しなさいということを示していますすなわち 10 回繰り返すということです x=0 1 For i=1 to 10 step 2 i x(1) x(2) 2 x=x+i next I MsgBox x この例は増分が 2 ですそのためカウンタ変数 i は 2 づつ増加しますもちろんカウンタ変数には明確な数値でなくても初期値と最終値に変数を指定することができます For Next ステートメントでも繰り返しの回数の途中であっても抜け出すことができますその命令は Exit For です While Wend ステートメント 3 つめの繰り返し命令は While Wend ステートメントです記述方法を示します While 条件式繰り返すべき命令群 Wend この命令では条件式が成り立っている間は繰り返し続けます逆に条件式が成立しなくなったら繰り返しの実行が終了します例で見てみましょう

16 intval=0 intadd=0 While intadd<=100 intval=intval+intadd intadd=intadd+1 Wend MsgBox intval このプログラムでは最初の 2 行で変数を 0 に初期化しています While Wend の条件式は intadd<=100 なので intadd という変数が 100 以下の時は繰り返します繰り返し命令では単に intadd の足し算を行っているだけですこの繰り返しの中ではカウンタ変数として intadd が用いられていますそして intadd=intadd+1 でカウンタを 1 つずつ増加させていますそして intadd が 100 まで繰り返しを行い 101 になったときに終了します While Wend ステートメントではこれを終了する命令がありませんそのため条件式および条件の変化をきちんと行わないと無限ループ ( いつまでたっても終わらないプログラム ) になってしまいがちです以上で基礎知識の説明は終わりましたので 5 の部分の説明を行います 1 Do 2 If No >= max Then Exit Do 3 x = X1 4 No = No ff = x * Log(x) df = Log(x) X1 = x - (ff / df) 8 WS.Cells(39 + No, 2) = "X1=" 9 WS.Cells(39 + No, 3) = X1 10 WS.Cells(39 + No, 5) = No 11 WS.Cells(39 + No, 7) = x - X1 12 Loop While (Abs(x - X1) > eps)

17 112 この間を繰り返す 2If No >= max Then Exit Do No:Do Loop の間を何回実行したかを示すカウンタ変数 Max: 入力した最大反復回数この判断文では繰り返しが最大反復回数を超えたら繰り返し操作を終了する (Exit Do) 3x=X1 x : 一つ前の方程式の数値解 X1 : x を用いて計算する新しい数値解この命令で新たに計算するにあたって x に前回計算した X1 の値を代入している 4No=No+1 カウンタ変数を1つ増やしている ( 何回計算したか確認のため ) 6 ff = x * Log(x) 1 y xlog e x 1 7 df = Log(x) + 1 y (y の導関数 ) f ( xi ) 8 X1 = x - (ff / df) xi 1 xi f '( x ) ( ニュートンラフソン法 ) 8 WS.Cells(39 + No, 2) = "X1=" X1= 回数前回との差 9 WS.Cells(39 + No, 3) = X1 の形式での結果の出力 10 WS.Cells(39 + No, 5) = No 11 WS.Cells(39 + No, 7) = x - X1 8~11では Cells の行の指定が 39+No となっています No は繰り返し変数ですつまり何行目に出力するかが繰り返し回数によって変わることを意味します 1 回目は40 行目 2 回目は41 行目ですこれにより計算過程を記録することができます 12 Loop While (Abs(x - X1) > eps) 条件式は Abs(x - X1) > eps です eps は閾値ですさらに Abs は絶対値を取る関数ですつまり前回と今回の計算結果の差が eps 以上であったら繰り返しを実行し続けることを意味していますプログラムの一番最後 (6 ページ目の 6) は最終結果の出力です i

Microsoft Word - VBA基礎(3).docx

Microsoft Word - VBA基礎(3).docx 上に中和滴定のフローチャートを示しましたこの中で溶液の色を判断する部分がありますこのような判断はプログラムではどのように行うのでしょうか判断に使う命令は IF 文を使います IF は英語でもし何々ならという意味になります条件判断条件判断には次の命令を使います If 条件式 1 Then ElseIf 条件式 2 Then ElseIf 条件式 3 Then 実行文群 1 実行文群 2 実行文群