なかのにっき

最新版はこちら。 突っ込みは各日付の BBS エントリのほか、 メール (nakano@st.seikei.ac.jp) や フォーム からどうぞ。 なおスパム除けのため、BBS 機能には 緩い認証を入れて います。 検索エンジンから来た方は、エンジンの方のキャッシュを見るか、 下の簡易検索を試してみてください。


hns - 日記自動生成システム - Version 2.19.5

先月 2002年12月 来月
1 2 3 4 5 6 7
8 9 10 11 12 13 14
15 16 17 18 19 20 21
22 23 24 25 26 27 28
29 30 31
Namazu for hns による簡易全文検索
詳しくは 詳細指定/ヘルプを参照して下さい
検索式:

2002年12月21日() [n年日記]

#1 [debian] logcheck

たまたま /var/tmp/logcheck を消しちゃったら、動作しなくなった。 うーん、これパッケージ supplied なディレクトリなのか。 なんで /var/lib 以下でやっちまわないんやろか。 せめてなければ作ってほしいような気が。

#2 [labo] MgO フィッティング用プログラム (2)

基板屈折率の扱いに悩み中。 どうせ 400〜800nm なので k=0 と置いて透過データから作ってしまえばいいのだが、 生データに対するフィルタと補間をどうするのが一番 手抜きできる スマートか。 本をパラパラ眺めたりしているだけで一向に進まないので、 後回しにして他のとこ実装しちゃうか。

cvs import:

thory.f と amoeba の簡単な interface を書いて、 ひとまず MgOfit というモジュールで import した。

ガラス透過率測定:

プログラミングと並行して。 データは ~/project/MgOfit/tglass 以下。 0F, 1R, 2F, 3R と向きを反転させつつ 4 回測った。 一応スタンダードに電源 on から 1h 待って測定。 以下気になったこと。 ちなみに gnuplot で img2dat の結果を見るときは
plot 'ふがほげ.dat' every ::1:1
と every を使えば変換せずにすむ。

使った Maxell のフロッピーでリードエラーが多発。 分光計のドライブとは相性が悪いようだ。 データを移動後あらゆる廃棄。

ガラス透過率の平滑化:

Excel で作業。 ファイルは ~/project/MgOfit/tglass/tglass.xls。
Function movave(cc As range, width As Variant) As Variant
Dim ccrow As Integer
Dim cccol As Integer
Dim sum As Variant

    ccrow = cc.Row
    cccol = cc.Column
    
    sum = 0
    For irow = ccrow - width To ccrow + width
        sum = sum + Cells(irow, cccol)
    Next
    movave = sum / (2 * width + 1)
End Function
とかいう関数を作って眺めてみる。 ざっと見た感じ、15〜25 nm 以上の移動平均を取らないと、 ノイズが残る感じか。 もうちょい取得時のバンドパスを広くして、 積分時間も大きくするほうがいいのかもしれないな。

Excel でテキスト保存:

うへ、データをテキスト保存しようとすると、 さっきの関数で計算したセルが 0 クリアされる。 しょうがないので別シートに手動で「値貼り付け」してから保存。 せめてマクロ化しようと思ったけど、 マクロ実行時にやっぱり 0 クリアされる。だめだこりゃ。死ネ。

スプライン関数:

4 本のデータを波長ごとに平均し、71 点 (=28.4 nm) の移動平均をとり、 400〜800 の区間を 20nm おきにサンプリング、 この生データを節点としてスプライン関数をつくることにした。 D 論の時に作ったルーチンを引っ張り出すために以下を実施。
cvshome -d sparse HISPUT/contrib/sparse
cvshome -d spline HISPUT/lib/nakano
cd spline
cvs tag -b for_MgOfit_20021221
cvs up -r for_MgOfit_20021221
うりうり書いて nglass.F 完成。 リポジトリの tglass/ が整理したガラスの透過率データと DATA 文を作るためのスクリプト・Makefile で、 nglass/ 以下がテストルーチン。

spline1d.F に与える x のデータが昇順でないといけないことに気づかずにしばらく嵌まる。 ソースのコメントには書いてあったけど、エラーにしないとな。 x が等間隔じゃなくても OK にするのも含め、いつか直さんと。

完成:

amoeba へのインターフェース (fom.F) に、 初期化と結果出力の各ルーチンを書く。 ひとまず動くバージョンができた。 生データから こないだ作った前処理 とかを含めて一括処理するスクリプト書きはまた今度にしよう。 fortran の書式付き read と実際のデータとの擦り合わせがちと面倒くさいな。

#3 [linux] make の日本語カタログ

Makefile:19: *** 分離記号を欠いています (8 個の空白でしたが、TAB のつもりですか?).  中止。
すげえむかついた(わら > 「つもりですか?」

#4 電話

坂上から。 「牡蛎をもらったので明日食いに来い」 というお誘い。
コメント [全部読む/投稿する]

以上、1 日分です。
タイトル一覧
カテゴリ分類
book
dept
issp
labo
paper
snap
stock
vsj
Powered by hns-2.19.5, HyperNikkiSystem Project

中野武雄 (NAKANO, Takeo) <nakano@st.seikei.ac.jp> Since 1999-10-07
RSS feed, 更新時刻, LIRS エントリ, アクセス制御 (解説)

中野のホームページへ