Ｃ言語第８回

PDFﾌｧｲﾙの印刷

【学内】PDFファイルリーダーで印刷可能です。PDFファイルでの配布のｐｄｆのリンクをクリックしてください。現われる講義録をで印刷してください。

	PDFファイル閲覧・印刷にはPDFファイルリーダーが必要です。
	無償で使用できる高性能なPDF作成ソフト

第８回：数値シミュレーション２（PDFファイルでの配布：c8.pdf）（Web版は更新されていません、pdfをご覧下さい）

【１】２階の微分方程式と差分方程式

微分方程式を、

とする。これを２つの１階の微分方程式に変更する。

それぞれの１階の微分方程式にRunge-Kutta法を適用する。

【２】ＲＬＣ回路のシミュレーション

図のようなＲＣ回路にコイル（Ｌ）が入っている。スイッチＳで端子を電池（Ｅ）側から切り替えると、コンデンサに生じた電位が徐々に消滅していく。そのときのコンデンサの電圧変化は、時間とともに変化する。電流の向きを、ＲからＬの方向にとると、

となる。これを２つの１階の方程式にするには、

と置くと

になる。初期条件として

にする。

【３】Ｃ言語プログラム

数値として

を与える。

になる。

【４】表計算ソフトの自動立ち上げ（Windows用のみ）

計算結果を画面表示と同時にファイルを作成し、作成したファイルを自動的に表計算ソフトで表にして表示するようにする。

このプログラムが実行されるドライブ名の取得：GetCurrentDirectory関数
- GetCurrentDirectory(整数, 番地)
  - GetCurrentDirectory(文字配列の個数, 取得したファルダ名の格納用番地)
    - char name[100]:
    - GetCurrentDirectory(100, name)
    - name[0]がドライブ名の文字（コード）になっている
ｃsvファイルを作成：fopen関数
- fopen(番地, 決まった文字列);
- fopen(csvファイル名の格納番地, "wt");
  - FILE * file_open;
  - file_open= fopen(filename, "wt");
  - file_openの数値が零（NULL）だと作成に失敗
ｃsvファイルへの書き込み：fprintf関数
- fprintf(FILE * 型番地, ...);
  - printf("%5.2f %10.3e\n", t, x1);
  - fprintf(file_open, "%5.2f, %10.3e\n", t, x1);
ｃsvファイルを閉じる：fclose関数
- fclose(FILE * 型番地)
  - fclose(file_open)
表計算ソフトでｃsvファイルを開く：ShellExecute関数
- ShellExecute((HWND)0, "open", 番地, NULL, NULL, SW_SHOWNORMAL)
  - ShellExecute((HWND)0, "open", csvファイル名の格納番地, NULL, NULL, SW_SHOWNORMAL)
  - 戻り値は特別でHWND型になっている。戻り値が３２以下だとcsvファイル用表計算ソフトの起動に失敗