- 292 名前：デフォルトの名無しさん mailto:sage [2008/05/31(土) 22:34:18 ]
- [1] ゼミ実験
[2]
○実験１、マルコフ情報源の生成プログラムの作成
プログラムは
double markov(int K ,double alfa, double *a_coeff, double *mem);
のように、次数Ｋとαの値および長さＫの係数配列a_coeffとメモリ配列memのポインタを引数としマルコフ情報源出力Ｘｎを返す関数文で書く。
main文では、Kとαの値の設計および配列a_coeffとmemの確保と配列の値の設定を行う。
なお、memは全零に初期設定する。
main文で関数marukovを呼出す度にマルコフ情報源の系列の値が順次得られ、N個のデータをファイルに記憶する。
なお、一様乱数からガウス関数への変換は(1)または(2)の方法を用いる。
○実験２、マルコフ情報源の生成実験
上で作成したプログラムを用いて、K=１の１次マルコフ情報源についてのみ系列の生成実験を行う。
ここで、係数a1は安定条件より|a1|＜1の範囲に限定する必要があり、実験では、a1=±0.1、±0.9の４種類についてマルコフ情報源出力系列の1024個のデータ（最初の１番目のデータは捨てる）をファイルに出力され、それらのグラフを描く。
なお、係数αは出力系列Xnの分散σx＾2=1になるように決定する。
問題【１】実験１のマルコフ情報源の生成プログラムを示しなさい。
問題【２】実験２のマルコフ情報源出力系列を図示し、AR係数a1=±0.1、±0.9の場合について、系列の違いを定性的に説明しなさい。
また、その違いの出る理由をK=1のときの式(2)のファイルの周波数応答F（e＾jω)=α/(1+a1･e＾-jω)より説明しなさい。
Uｎ=(1./(RAND_MAX+1.))*rand(); ････(1)式
F(z)=α/(1+a1･z＾-1+a2･z＾-2+･･････+ak･z＾-K) ･････(2)式
[3] 環境
[3.1] OS：ＸＰ
[3.2] コンパイラ名とバージョン： Vidual Studio2005
[3.3] 言語：C++
[4] 期限： 2008年6月1日24:00
[5] その他の制限：
急ぎですいません。問題【２】実験２は図示があるので問題【１】だけでもいいのでお願いします。
|
 |
