[pear_error: message="Success" code=0 mode=return level=notice prefix="" info=""]
JUGEMテーマ:数学 たいていの積分はシンプソンの公式を用いることで、いい精度で計算することができます。従って、今回はシンプソンの公式を利用したプログラムを作成することで円周率の計算を行ってみようかと思います。理工系の皆様なら解ると思いますが、√(1-x^2)を0〜1の範囲で積分して、その値を4倍すれば円になるということを理解しているはずです。ですから、関数としてf(x)=4√(1-x^2)を積分すれば円周率を計算することができるのです。早速やってみましょう。プログラム例)Sub main() Dim N, DERU As LongDim...
Excel VBAによる数値解析(Excel 2007編) | 2010.01.03 Sun 05:49
JUGEMテーマ:数学 √2の近似値を求めるプログラムなら、ネット上でいっぱい紹介されているけれども、正直汎用性がありません。何故なら、あくまでx^2―2=0の解は解けても、方程式が変わると解けなくなるじゃないですか?もう少し一般的な方程式がとけるようなプログラムを作成することが大事だと思います。そこで、下記のようなプログラムを組んでみてはどうでしょうか?Sub main() Dim A As SingleDim ob1 As ObjectSet ob1 = Application.ThisWorkbook.Worksheets("Sheet1") A = InputBox("初期値を入力してください") Wi...
Excel VBAによる数値解析(Excel 2007編) | 2010.01.03 Sun 05:49
JUGEMテーマ:数学 Excelを使ってモンテカルロ法による円周率計算を行ってみました。アイディアは大変簡単なものです。下の図を見てください。図のように,まず、1×1の正方形の中に中心角が90°、半径1の扇形があるような図形を考えます。この図形の中に、乱数を利用して任意の座標(図の赤丸)を発生させます。この座標を統計的にものすごくたくさんとると、打ち出したすべての点と扇形の内部に存在する点の個数比が、確率的に正方形と扇形の面積比と一致するので、その事実を利用して円周率を定める手法がモンテカルロ法になります。円...
Excel VBAによる数値解析(Excel 2007編) | 2010.01.03 Sun 05:48
JUGEMテーマ:数学 正直、シンプソンの公式の方が精度がいいことから、あまり僕なんかは使わないんだけれども、せっかくだから紹介しておこうかと思います。Sub main() Dim N, DERU As LongDim A, B As SingleDim ob1 As ObjectSet ob1 = Application.ThisWorkbook.Worksheets("Sheet1") Do N = InputBox("分割数を入力してください。") DERU = 1 'Loopを抜けるためのパラメータ If N <= 0 Then MsgBox "分割数エラー" DERU = 0 End If Loop Until D...
Excel VBAによる数値解析(Excel 2007編) | 2010.01.03 Sun 05:47
JUGEMテーマ:数学 結局、ニュートン法のほうが収束自体ははやいのですけれども、二分法というやり方もあります。この解の探し方は、解の範囲を二つに絞ってから近似解を探すというものです。教科書に載っているから、二分法の宿題が出されている方のために、プログラムサンプル例を適当に作ったので紹介したいと思います。今回探す解はx^2−2=0としましょう。プログラムサンプル例)Sub main()Dim A, B As Single DoA = InputBox("解の範囲a〜bに対してaを入力してください")B = InputBox("解の範囲a〜bに対してbを入力してく...
Excel VBAによる数値解析(Excel 2007編) | 2010.01.03 Sun 05:47
JUGEMテーマ:数学 連立方程式を解く手法としてガウスジョルダン法というものがあげられます。例えば、三元の連立一次方程式を例に考えましょう。これらの方程式は一般的に、[a11,a12,a13][x] [b1][a21,a22,a23][y]=[b2][a31,a32,a33][z] [b3]といった行列形式で表すことができます。そこで、係数行列と定数のベクトルを組み合わせることで、新に下記のマトリックスを作りましょう。[a11,a12,a13,b1][a21,a22,a23,b2][a31,a32,a33,b3]この行列を作りましたら、次の操作に従って下記の形にすることを考えましょう。[1,0,0,c1][0,1...
Excel VBAによる数値解析(Excel 2007編) | 2010.01.03 Sun 05:46
JUGEMテーマ:数学 これまでに、数値的に積分する手段として、台形公式やシンプソンの公式という手法があること、そして数値的に方程式を解く手段として二分法やニュートン法等があると説明しました。しかしながら、これらはy=f(x)がどのような関数か具体的にわかっていないと使いようがない形式でプログラムを作成しています。そこで、これらの手法の適用範囲を拡大できるように、f(x)をラグランジュ補完法という手段で近似的に作る方法をプログラムで紹介したいと思います。(1)数値例http://www.st.chukyo-u.ac.jp/hatano/...
Excel VBAによる数値解析(Excel 2007編) | 2010.01.03 Sun 05:45
----------------------------------------------------------------------------------- リーマン予想に北野武さんもビックリ ここから ----------------------------------------------------------------------------------- NHK 魔性の難問〜リーマン予想・天才たちの闘い〜 ここから ----------------------------------------------------------------------------------- 『素数に憑かれた人たち ‐リーマン予想への挑戦‐』(ジョン・ダービーシャー,日経BP社) ここから ------------------...
奇跡の超魔術 | 2009.12.22 Tue 16:38
NHK特集は 魔性の難問 〜リーマン予想・天才たちの闘い 数学者が素数は美しいと表現していましたが、 リーマン予想は、素数の法則を解明する難問で、 150年間経ってもまだ解明されていないのだそうです。 リーマン予想は解決するのか? JUGEMテーマ:数学
時のわすれもの | 2009.11.15 Sun 23:28
JUGEMテーマ:数学数学と物理ってふと思った時に突然勉強したくなるあの数字が並んでるのみるとどうしてもゾクゾクしてしまうんだよね…
日向と海晴の交換日記と見せかけた気ままな日記 | 2009.09.03 Thu 10:06
全119件中 101 - 110 件表示 (11/12 ページ)