数値解析法及び演習 第三回

DOループの応用
  1. 今日の目標
  2. 異なる型の混合演算について
  3. gnuplotの基本
  4. 配列
  5. 最小二乗法
  6. 今日の宿題

[1] 今日の目標

  1. 数の型の理解を深める
  2. 配列の概念を学ぶ

[2] 異なる型の混合演算について

FORTRANでは実数と整数を区別して変数に格納します. 正しい型の変数を用意しておかないと計算結果がおかしくなります. 変数の型には注意をはらうこと!

  • 実数は必ず小数点がつく! 3.D0は実数，3は整数． 整数型と実数型の変数を混ぜて使うときには注意が必要である．以下に考えられるパ ターンを挙げる． 
          implicit none
      integer m,n,l
      real*8 a,b,c,d,e

      l=5/2
      a=5/2
      m=5.D0/2.D0
      b=5.D0/2.D0
      n=5/2.D0
      c=5/2.D0
      e=10

      write(*,'(3(I2,x),2(D24.16,x))') l,m,n,a,b
      write(*,'(3(D24.16,x))') c,d,e
      end
    結果はどうなったであろうか? dの値は？

    • [3] gnuplotの基本

    gnuplot(ニュープロットとよむ)は，UNIXシステムにおける簡便なグラフ描画 ソフトである． 起動は 
    gnuplot
    で行われる．
    最も基本的なデータは 
    -1.0   -3.2
 0.0   -0.9
 1.0    1.2
 2.0    2.8
 3.0    5.0
    といった形で与えられる．左側のデータをx座標，右側のデータ(ファイル名を plot.datとする)をy座標とおもってxy平面上にプロットするときは，gnuplotを起動した画面で 
    plot 'plot.dat'
    とするとグラフが描かれる．さらに 
    plot 'plot.dat' w l
    とするとデータ点を線で結ぶ．また，gnuplotにはいろいろな関数が用意されていて， 
    plot sin(x)
    とするとサインカーブを描いてくれる．他にもワザはたくさんあり
    GNUPLOT入門（kawanoさんのページ）
    等に出ているので参考にしてください．

    [4] 配列

    通常の変数（整数型，実数型ともに）は数値を一つしか格納できない．これでは多くのデータを扱うときに困る．そこで「配列」とよばれるものが用意されている．配列は複数の要素からなり，要素の番号を指定することで複数の数値を格納することができる．

    配列の宣言:
    宣言文 配列名（寸法の下限：寸法の上限）

    例1 
    real*8 a(10)
    a(1)からa(10)まで要素を10個持つ倍精度実数型配列を宣言

    例2 
    integer i(-1:8)
    i(-1)からi(8)まで要素を10個もつ実数型配列を宣言

    例： 
           implicit none 
       integer i
       real*8 a(10)

       a(1)=1.D0
       a(2)=2.D0
       a(3)=3.D0


       do i=1,10
          write(*,*) i,a(i)
       enddo

       end
    a(4)以降の要素はどう表示されたでしょうか？

    練習問題1

    次のプログラムを改造して，

    それぞれの総和を求めなさい． 

           implicit none 
       integer i
       real*8 a(10),b(10),sum(3)

       do i=1,3
          sum(i)=0.D0
       enddo
    
       do i=1,10
          a(i)=1.D0*i
          b(i)=2.D0*i
       enddo

       **ここに必要なプログラムを書き足して下さい**

       end

    練習問題2

    上のプログラムを改造して，

    を求めなさい．

    [5] 最小二乗法

    二つの量を測定したときに，そのあいだに線形の関係が期待されることがよくある． もっともありそうな傾き(B)と切片(A)を求めるのが最小二乗法である．
    最小二乗法の解説

    最小二乗法では, 測定されたデータの組(x,y)にフィットする直線を求める. データxは正確に求められているとし, またデータyのバラツキは正規分布をしているものとする. この場合, フィットする直線をy=A+Bxとすると, A,Bはそれぞれ

    $\displaystyle A$ $\textstyle =$ $\displaystyle {\sum x_i^2\sum y_i-\sum x_i\sum x_i y_i\over N\sum x_i^2-(\sum x_i)^2}$  
    $\displaystyle B$ $\textstyle =$ $\displaystyle {N\sum x_iy_i-\sum x_i\sum y_i\over N\sum x_i^2-(\sum x_i)^2}$  

    となる. したがって切片Aと傾きBを算出するためには次の和をもとめればよい．

    練習問題3

    次のx，yのデータをファイルから読み込ませ，最小二乗法によりフィットした 直線の傾きと切片を求めるプログラムを作成せよ. AとBを求めたら, その直線 とデータをgnuplotで描き, フィッテイングが正しいことを確かめよ. 
    
-1.0   -3.2
 0.0   -0.9
 1.0    1.2
 2.0    2.8
 3.0    5.0

    [6] 今日の宿題

    練習問題3のプログラム, 答えを送ること. 〆切は10/22.


    日程表へ戻る <<