情報工学Ⅱ （第９回） 月曜４限 担当：北川 晃.

Presentation on theme: "情報工学Ⅱ （第９回） 月曜４限 担当：北川 晃."— Presentation transcript:

1 情報工学Ⅱ （第９回） 月曜４限 担当：北川 晃

2 配列とは？ 表を一つの変数のように扱う 一次元の配列： → a(i) ベクトル 二次元の配列： → b(i,j) 行列

3 配列の引数には，整数型の式，関数なども許される
配列の表記法 配列内の個々のデータ（配列要素）は， 配列名に添え字を付けて表す． alpha(16) n(5,1) p(i,j) b(i+1) x(2*n-1) k(l(m)) 配列の引数には，整数型の式，関数なども許される 配列の宣言 一次元， 単精度実数型 の配列の宣言 a(0)～a(10) Dim a(10) As Single Dim b(8,8) As Integer Dim c(3,5) As String 二次元，整数型の配列の宣言 b(0,0)～b(0,8) b(1,0)～b(1,8) … b(8,0)～b(8,8) 配列の番号は0から始まる 使わない要素は0と見なされるだけなので， 少し余裕を見て大きめに宣言するとよい

4 配列の操作 配列Aの内容をBにコピーする（下限1，上限N）． For i As Integer = 1 To n b(i)=a(i)
Next 行列A，Bの和を作ってCに代入する（M行N列とする）． For i as Integer = 1 To m For j As Integer = 1 To n c(i,j) = a(i,j) + b(i,j) Next j Next i

5 例題:九九の表 九九の表を書き出すプログラムを作れ． 1 2 3 4 5 6 7 8 9

6 九九の表：プログラム例 行ごとに改行する

7 行列の積：プログラム例 Console.Title = "九九の表" Dim i, j, kuku(9, 9) As Integer
For j = 1 To 9 For i = 1 To 9 kuku(i, j) = i * j If i <> 9 Then Console.Write("{0,4}", kuku(i, j)) Else Console.WriteLine("{0,4}", kuku(i, j)) End If Next 表の右端で改行する

8 ベクトルの内積 ベクトルAとベクトルBの内積を計算する（下限1，上限N）． s = 0.0
For i As Integer = 1 To n s = s + a(i) * b(i) Next

9 行列の積の成分計算 行列Aと行列Bの積を計算する（下限1，上限N）． 第(i,j)成分は・・・ s = 0.0
For k As Integer = 1 To N s = s+a(i,k)*b(k,j) Next k

10 プログラミング演習 3×3の正方行列A，Bを外部ファイルからデータとして 読み込み，その積C=ABを計算するプログラムを作れ．
data_a.txt -2 1 4 3 2 -1 1 2 -3 -1 5 data_b.txt

11 行列の積：出力例 行ごとに改行する

12 Streamクラスを用いたファイルの接続
プログラム中における ファイルの名前 … Dim インスタンス名 As New IO.StreamReader( _ “ファイルの絶対パス”, _ System.Text.Encoding.Default) s = インスタンス名. インスタンス名.Close() ReadToEnd：ファイルの内容すべてを取得 ReadLine：ファイルから１行読み込む Read ：ファイルから１文字読み込む ファイルを閉じる

13 外部ファイルからのデータの読み込み：例文
Dim a, b, c As Single Dim ReadText As　New IO.StreamReader( _ “D:\read_test\data.txt”, _ System.Text.Encoding.Default) a = ReadText.ReadLine() b = ReadText.ReadLine() c = ReadText.ReadLine() Console.WriteLine(“a={0}”,a) Console.WriteLine(“b={0}”,b) Console.WriteLine(“c={0}”,c) ReadText.Close() １行ずつ 読み込み ファイル名”data.txt”として “D:\read_test\”に保存 1 ファイル を閉じる

14 外部ファイルからのデータの読み込み：出力例

15 行列の積：プログラム例 Dim ReadText1 As New IO.StreamReader( _
"D:\...\data_a.txt", System.Text.Encoding.Default) Dim ReadText2 As New IO.StreamReader( _ "D:\...\data_b.txt", System.Text.Encoding.Default) Dim n As Integer = 3 Dim a(10, 10), b(10, 10), c(10, 10) As Single For i As Integer = 1 To n For j As Integer = 1 To n a(i, j) = ReadText1.ReadLine() b(i, j) = ReadText2.ReadLine() Next ReadText1.Close() ReadText2.Close() 二つの外部 ファイルを接続 配列の宣言 （サイズは大きめに） 行列データの 読み込み ファイルを閉じる

16 行列の積：プログラム例（続き） For i As Integer = 1 To n For j As Integer = 1 To n
c(i, j) = 0.0 For k As Integer = 1 To n c(i, j) = c(i, j) + a(i, k) * b(k, j) Next If j <> n Then Console.Write("{0,10}", c(i, j)) Else Console.WriteLine("{0,10}", c(i, j)) End If 積の(i,j)成分の計算 行の右端で改行する

17 データの統計処理に関する諸問題 データの分布の様子 データの期待値（平均値） データの分散，標準偏差（ばらつきの度合い）
回帰直線へのフィッティング（最小二乗法） 一様乱数の発生 正規乱数への変換 モンテカルロシミュレーション

18 例題:度数分布 n個のデータx1, x2, …, xnが与えられている． ‘normal_data.txt’よりデータを逐次読み込んで，
次の区分に従って度数分布を調べよ． x＜-2.25 -2.25≦x＜-1.75 -1.75≦x＜-1.25 -1.25≦x＜-0.75 -0.75≦x＜-0.25 -0.25≦x＜ 0.25 0.25≦x＜ 0.75 0.75≦x＜ 1.25 1.25≦x＜ 1.75 1.75≦x＜ 2.25 2.25≦x

19 度数分布：考え方 それぞれ領域を以下のように定める． x＜-2.25：領域0 -2.25≦x＜-1.75：領域1
データを読み込んで，領域iに当てはまるとき， c(i)の値を1増やす． p(i)=c(i)/nの値を計算し，書き出す． x＜-2.25：領域0 -2.25≦x＜-1.75：領域1 -1.75≦x＜-1.25：領域2 … ：… 1.75≦x＜ 2.25：領域9 2.25≦x ：領域10

20 出力例：度数分布

21 度数分布：プログラム例 読み込んだ値がどの領域に属するか判定
Dim n, c(100) As Integer, x, p(100) As Single Dim ReadText As New IO.StreamReader( _ "D:\…\normal_data.txt", System.Text.Encoding.Default) For i As Integer = 0 To 10 c(i) = 0 Next n = ReadText.ReadLine() For i As Integer = 1 To n x = ReadText.ReadLine() If x < Then c(0) = c(0) + 1 ElseIf <= x And x < Then c(1) = c(1) + 1 ElseIf <= x And x < Then c(2) = c(2) + 1 読み込んだ値がどの領域に属するか判定

22 度数分布：プログラム例（続き） 読み込んだ値がどの領域に属するか判定
ElseIf <= x And x < Then c(3) = c(3) + 1 ElseIf <= x And x < Then c(4) = c(4) + 1 ElseIf <= x And x < 0.25 Then c(5) = c(5) + 1 ElseIf 0.25 <= x And x < 0.75 Then c(6) = c(6) + 1 ElseIf 0.75 <= x And x < 1.25 Then c(7) = c(7) + 1 ElseIf 1.25 <= x And x < 1.75 Then c(8) = c(8) + 1 ElseIf 1.75 <= x And x < 2.25 Then c(9) = c(9) + 1 Else c(10) = c(10) + 1 End If Next 読み込んだ値がどの領域に属するか判定

23 度数分布：プログラム例（続き） 各領域の規格化度数を計算 For i As Integer = 0 To 10
p(i) = c(i) / n Console.WriteLine("p({0})={1}", i, p(i)) Next

24 出力例：度数分布（エクセルによる表示） 釣り鐘型（正規分布）