「微分」 現場数学(その1)

更新日

投稿日

♦ 数学の教科書に書いてある通りにプログラムしてもだめ!では計算できない!

 サイン・コサイン何になる…とか嘯(うそぶ)いていた学生時代。しかし「現場に出て数学の重要性をひしひしと感じているのに、確実な理解が出来ていない…」と不安を抱いている技術者の皆さん。そこで、今回から18回にわたり「本当に意味のある“現場数学”とは何か」について私ともう一度考えてみませんか?

 

1.“現場数学”では、hは有限でなければならない

 高校で微分とか積分を習った時、先生が今まで見たことのない如何(いか)にも新しい式であると黒板に書く様に感激した経験を思い出します。しかし、これでは本当は計算できないのです。「現場数学」では、hは有限でなければなりません。となると…小さい、と思われる数、例えば0.00001とかにすれば良いということでしょうか?いえいえ、全くそうではないのです。この辺りから疑問を感じていた人は多いのです。でも、高校の先生に聞くと、ムニャムニャとした返事。それで数学が嫌いになった人がいるほどです。どうなっているのでしょうか?

 実はhは小さいどころか無限大でも良いのです!これは極端な話です。でも、高校数学で重要なことを教えていないことは確かなのです。また、大学に入ると微積分は途端にε-δの話になり、余計に誤魔化されてしまいます。実は、最初に考えなければならないのは対称性の問題なのです。通常の微分係数の記述は、

    

 です。しかし、これでは、xの位置から正方向のみの変化を見ることになり、その位置で増加している関数なら当然、x点での値より大きくなります。それにも関らず なのだから、正しい値に近づくはずと習うのです。これは一種の詭弁(きべん)です。二次関数y=x2にこの微分係数の式を当てはめてみましょう。

    

 となります。hだけの誤差があり、 で正しい2xになるという論法です。実は、この公式は前進公式と呼ばれます。xから正方向だけをみているからです。x点での微分係数を求めたいのですから、-h方向もみないといけません。それは後退公式と呼ばれる次式です。

    

 この式の誤差は-hです。ということは(2)式と(3)式の平均をとれば、hの大きさに依(よ)らず微分係数は2xと厳密に求まります。

    

 つまり微分係数はhに依らず、極端にいえばhは無限大でも良い式が作れたのです。この例では、2次関数の微分は一次関数ですから、hに依らず誤差が消えてしまいました。3次関数でも、何次関数でも同様に誤差を消してしまう公式が作れます。ぜひ、自分の手を動かして試して下さい。

2. より本質理解し、数式を上手に扱う

 この例題は、教科書に書いてある通りの式をそのままプログラムして、数値計算してはいけないことを如実に示しています。「現場数学」では、数値をなるべく正確にしかも効率良く計算しなければなりません。それを身につけるためには、このような高校の教科書に載っていなかった方策を勉強する必要があります。これは、頭の良...

♦ 数学の教科書に書いてある通りにプログラムしてもだめ!では計算できない!

 サイン・コサイン何になる…とか嘯(うそぶ)いていた学生時代。しかし「現場に出て数学の重要性をひしひしと感じているのに、確実な理解が出来ていない…」と不安を抱いている技術者の皆さん。そこで、今回から18回にわたり「本当に意味のある“現場数学”とは何か」について私ともう一度考えてみませんか?

 

1.“現場数学”では、hは有限でなければならない

 高校で微分とか積分を習った時、先生が今まで見たことのない如何(いか)にも新しい式であると黒板に書く様に感激した経験を思い出します。しかし、これでは本当は計算できないのです。「現場数学」では、hは有限でなければなりません。となると…小さい、と思われる数、例えば0.00001とかにすれば良いということでしょうか?いえいえ、全くそうではないのです。この辺りから疑問を感じていた人は多いのです。でも、高校の先生に聞くと、ムニャムニャとした返事。それで数学が嫌いになった人がいるほどです。どうなっているのでしょうか?

 実はhは小さいどころか無限大でも良いのです!これは極端な話です。でも、高校数学で重要なことを教えていないことは確かなのです。また、大学に入ると微積分は途端にε-δの話になり、余計に誤魔化されてしまいます。実は、最初に考えなければならないのは対称性の問題なのです。通常の微分係数の記述は、

    

 です。しかし、これでは、xの位置から正方向のみの変化を見ることになり、その位置で増加している関数なら当然、x点での値より大きくなります。それにも関らず なのだから、正しい値に近づくはずと習うのです。これは一種の詭弁(きべん)です。二次関数y=x2にこの微分係数の式を当てはめてみましょう。

    

 となります。hだけの誤差があり、 で正しい2xになるという論法です。実は、この公式は前進公式と呼ばれます。xから正方向だけをみているからです。x点での微分係数を求めたいのですから、-h方向もみないといけません。それは後退公式と呼ばれる次式です。

    

 この式の誤差は-hです。ということは(2)式と(3)式の平均をとれば、hの大きさに依(よ)らず微分係数は2xと厳密に求まります。

    

 つまり微分係数はhに依らず、極端にいえばhは無限大でも良い式が作れたのです。この例では、2次関数の微分は一次関数ですから、hに依らず誤差が消えてしまいました。3次関数でも、何次関数でも同様に誤差を消してしまう公式が作れます。ぜひ、自分の手を動かして試して下さい。

2. より本質理解し、数式を上手に扱う

 この例題は、教科書に書いてある通りの式をそのままプログラムして、数値計算してはいけないことを如実に示しています。「現場数学」では、数値をなるべく正確にしかも効率良く計算しなければなりません。それを身につけるためには、このような高校の教科書に載っていなかった方策を勉強する必要があります。これは、頭の良くなる本などで、特別な例を覚えておいて計算を早くする、という類のものではありません。より本質を理解して、数式を上手に扱う技法なのです。計算をするのに、そろばんを使おうとスーパーコンピューターを使おうと、どれも桁数は有限です。 は何倍精度を持ってしても厳密には成り立たないのです。そうではなく「現場数学」では、そもそも誤差を発生させなくする定式化、なるべく小さくする定式化をすることが肝要なのです。

 

 次回に続きます。

   続きを読むには・・・


この記事の著者

川添 良幸

市販の材料設計シミュレーションプログラムでは満足できない御社技術者に、本当に意味のある物理の基本に基づいた設計法を伝授します。我々の計算と実験結果が合わない場合は実験の方に問題があると言えるレベルを達成しています。

市販の材料設計シミュレーションプログラムでは満足できない御社技術者に、本当に意味のある物理の基本に基づいた設計法を伝授します。我々の計算と実験結果が合わな...


「SQC一般」の他のキーワード解説記事

もっと見る
~ 「定幅図形」マンホールのフタはなぜ丸い? 現場数学(その17)

  ♦ ルーローの三角形と定幅図形 1.自動車のタイヤは丸くないといけない?  ピラミッドを造った時、大きな石をどうやって...

  ♦ ルーローの三角形と定幅図形 1.自動車のタイヤは丸くないといけない?  ピラミッドを造った時、大きな石をどうやって...


~級数展開の不思議-桁落ち 現場数学(その2)

♦ 発想の転換で“桁落ち”回避  有限桁数の計算機を使って計算することによって起こる深刻な問題の例を挙げましょう...

♦ 発想の転換で“桁落ち”回避  有限桁数の計算機を使って計算することによって起こる深刻な問題の例を挙げましょう...


DPMO(Defects Per Million Opportunity) 工程能力の理解 (その3)

  【工程能力の理解 連載目次】  工程能力とは  工程能力指数-Cp/Cpk  ...

  【工程能力の理解 連載目次】  工程能力とは  工程能力指数-Cp/Cpk  ...


「SQC一般」の活用事例

もっと見る
統計教育が必須な検査担当者

 測定や検査結果の値で合格判定を行い、製品性能の品質保証をする事は製造業に取っては当たり前と言える日常的作業です。一方で検査や測定の正確さや信頼性を担保す...

 測定や検査結果の値で合格判定を行い、製品性能の品質保証をする事は製造業に取っては当たり前と言える日常的作業です。一方で検査や測定の正確さや信頼性を担保す...


信頼できる結果を得るために統計的手法の活用を

昔の事になりますが、配属された職場の先輩の下である実験を手伝いました。先輩は図1のような計画をたてました。2つの量的因子A、Bの効果を確認する実験でした。...

昔の事になりますが、配属された職場の先輩の下である実験を手伝いました。先輩は図1のような計画をたてました。2つの量的因子A、Bの効果を確認する実験でした。...


飛行機事故とセンサー・フィードバックの分布

   ボーイング 737 Max 8機が短い期間に2回の墜落死亡事故を起こし、多くの航空会社が同機種の運航を停止する事態となっています。多くの...

   ボーイング 737 Max 8機が短い期間に2回の墜落死亡事故を起こし、多くの航空会社が同機種の運航を停止する事態となっています。多くの...