C言語
Q:コンパイラーとしてどんなものがありますか。
A:
コンパイラーには、
(1) Microsoft Visual C++
Development system Learning Edition Ver 6.0
アカデミックパッケージ(市価 11,000円程度) 学生、教員
(2) Borland C++ Builder
などがあります。
Q:C言語ではエスケープ文字としてバックスラッシュを使いますが、
日本語では円記号を使うのは何故ですか。
A:これは、バックスラッシュのASCIIコードが 16進数の5C(0x5c、5Ch)
であることと、日本語のJISコードにおいては円記号がASCIIコードの 5Ch
に割り当てられていることによります。
つまり、コードが一致しているからです。
Q:C言語では文字も数字も一緒くたに扱われているように思えます。
初心者には紛らわしいのですが。
A:
私たちが画面で見る文字も数字も、最小情報量であるビット(0 か 1)
の組み合わせで表現されます。
例えば、ASCIIコードでは、英字の A は41h(0x41)という値で表現されます。
最も機械に近いレベルで表現すると 01000001 です。
つまり、'A'も0x41も機械内部では 01000001 であり、画面上では 英字のAです。
したがって、C言語では、変数をXとして、X='A'とするのと、X=0x41
とするのは同じことになります。
紛らわしく感じるのは、画面上の数字や文字をそのまま受け取るからです。
次のように考えると良いでしょう。私たちは文字「A」を文字そのものとして
理解していますが、コンピュータが扱う最小情報量はビットで 0 か 1ですから、
「A」そのものは認識できません。そこで、「A」を 01000001 に当てはめ、機
械で「A」という文字を処理できるようにしているのです。ただし、0 と1の並
びは人間にわかりにくいので、通常は16進数で表現されます。Aは 01000001
ですから、16進数で41となります。
コンピュータでの表現 画面での英数字
------------------------------------------------------------
16進数の 48 65 6C 6C 6F 英字で H e l l o
16進数の 31 32 33 34 35 数字の 1 2 3 4 5
とにかく、私たちが使う文字や数字が機械の内部ではどのように扱われているか、
ということに興味を持つことです。
Q:よく「文字化けする」ということを耳にしますが、文字化けする原理を教えて下さい。
A:まず、C言語との関連で説明しましょう。例えば、エスケープ文字は円記号(¥)で、
そのASCIIコードは 5Chです。ここで、漢字の「表」を例に取ると、シフトJISでは16進数
の 955C です。もうお分かりですね。もしコンパイラなどが 955C を漢字と認識せず、95と5C
と別々に読むと、5Ch は '\' ですから文字として認識されません。したがって、955Ch は 95
だけとして読まれてしまい、「表」はどこかへ消えます。古いコンパイラ(Lattice C ver.2
など)ではよくあった話です。現在のCコンパイラは漢字を扱えますが、注意はしておく必要
があります。
次に、漢字コードの関係で文字化けが生じる原因を説明しましょう。漢字変換には、パソコン
ではJISコードやシフトJISコード、ワークステーションではEUCコードが使われます。
例えば、「大学」を例にとると、
大 学
シフトJISでは、91E5 8A77
JISでは、4267 3358
と16進数の値が異なります。このように、同一漢字でも16進数が異なることから、漢字を認識
する方法も異なるということになります。ですから、異なる機種の間で通信する場合はコード
を自動的に判別するソフトでないと文字化けすることがあります。自分のパソコンがどの漢字
コードを用いているか、通信の相手がどのコードを用いているかを調べておくことも重要です。
試しに、図書館を呼び出して、図書検索をEUCコードと異なるコートでやってみて下さい。
Q:通信ではファイルをバイナリで送ったり(アップロード)、受け取ったり(ダウンロード)
するのが安全といわれますが何故ですか。
A: Internet の Q&A に書いて有ります。そこをご覧下さい。
Q:マスクとはどのようなことですか。
A:例で示しましょう。
X & 0X00FF とは、数字X と 0X00FFとのビットごとのANDをとることです。
0X00FFをビットで示す(上段はビット、下段はビットフィールド)と、
0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1
---------------------------------------------------
15 14 13 12 11 10 9 8 7 6 5 4 3 2 1 0
となります。上位8ビットはゼロですから、Xと0X00FFとのANDをとった数も、
上位8ビットはゼロになります。つまり、上位8ビットはフィルターをかけ
て隠すことを意味しますので、マスクするというわけです。
その他の例
X & 0x7fff 15ビット目をマスク
X & 0x003f 6ビット目以上をマスク
X & 0x000f 4ビット目以上をマスク
(補足)
インターネットでは「サブネットマスク」が使われます。
サブネットマスクは十進数のオクテットで、
255.255.0.0
255.255.255.0
などと表されます。
10進数の255は、16進数の FF、2進数の 11111111 です。
10進数の 0 は、16進数の 00、2進数の 00000000 です。
これから、なぜ「マスク」というのか予想できますね。
「サブネットマスク」については、Internetのページに説明があ
ります。
Q:ビットシフトとはどのようなことですか。
A:例で示しましょう。
0XFF00 >> 8 の場合、右へ8ビットシフトすると、右側の8ビットは除か
れ、代わりに左側に8ビット分 0が詰められます。
図示しますと、
1111 1111 0000 0000 ・・・・ 0XFF00
----------------------
1111 1111 ・・・・ 右の8ビット( 0000 0000)が除かれる
0000 0000 ・・・・ 左に8ビット0が詰められる
----------------------
0000 0000 1111 1111 ・・・・ 結局、0X00FF となる。
Q:マスクやビットシフトはどのようなことに応用されますか。
A:例で示しましょう。
例題:構造体の中に「年月日」を示す16ビットのビットフィールド
が入っているものとします。
Borland C++ Ver3.0 の ffblk.ff_fdate を例にとっています。
ffblk.ff_fdate のビットフィールド・・・ 日:ビット 0-4
月:ビット 5-8
年:ビット 9-15
図示すると、
15 14 13 12 11 10 9 8 7 6 5 4 3 2 1 0
|------- 年 --------|---- 月 ---|---- 日 ------|
となります。
ここで、年(1980年以降の年数)をy、月(1-12)をm、日(1-31)をd、としましょう。
年を得るのは簡単ですね。ffblk.ff_fdate を9ビット右へシフトすればよい。
(「月日」の桁を削除し、残った「年」を右へ桁下げすると考えると分かり易い)
結局
y = ffblk.ff_fdate >> 9;
日を得るのはどうですか。ffblk.ff_fdate の 左側11ビットをマスクすればよい。
0000 0000 0001 1111 は 0X001F です。
結局、
d = ffblk.ff_fdate & 0X001F;
月を得るのは少々手強いです。月を得るのは、マスクとビットシフトの応用問題です。
まず、5ビット右へシフトさせて、日の5ビットを捨てましょう。
すると 00000年月 が イメージとしてのビットフィールドになります。
次に、月(下位4ビット)だけ残すために、上位12ビットをマスクして隠します。
マスクするのは上位12ビットですから、
0000 0000 0000 1111 との AND をとればよい。( 0X000F との AND をとる。)
結局、
m = ( ffblk.ff_fdate >> 5 ) & 0X000F
分かりましたか。初級者には、ビット操作は無理かもしれません。しかし、C言語の魅
力の一つはこういったビット単位の操作が可能なことにあるのです。
Q:最も簡単なプログラムを一つあげてCのプログラムを説明して下さい。
A:どの入門書でも、"Hello World" と表示するだけのものから始まります。
@最も簡単なプログラムは、以下のような5行のプログラムです。
=====================================================
#include "stdio.h"
void main(void)
{
printf("\n Hello World \n");
}
=====================================================
・Cプログラムは関数の集まりのような形をしています。
・ #include は、プリプロセッサにヘッダファイル(stdio.h)を取り込ませる命令です。
stdio.h は基本的な inputとoutputのきまりを指定したファイルです。
stdio.hの両端を右のようにくくります。<stdio.h> もしくは "stdio.h"
--------------------------------------------------------------------
注意: "<" と ">" は半角です!! ここでは全角になっていますが、こ
うしないと netscape navigator はタグ(tag)と受け取り、画面に表
示されなくなるからです。
--------------------------------------------------------------------
どのプログラムにもstdio.hは必要です。
プリプロセッサ制御命令にはセミコロン(;)はつけません。
main()が必要です。中身は{ と }の間に書きます。
printf()は画面に表示する命令です。表示の中身は "Hello World" です。
"\n" は改行を意味します。
ですから、"\n Hello World \n" は、まず改行し、"Hello World" を表示し、改行する
ことを意味します。文の最後にセミコロン(;)をつけます。
A次は、画面消去後 "Hello World" と表示させます。
=====================================================
#include "stdio.h"
void main(void)
{
printf("\x01b[2J");
printf("Hello World \n");
}
=====================================================
printf("\x01b[2J") は画面消去命令です。\x01b はエスケープコードです。
また、例えば画面を青くしたければ、
printf("\x01b[34m");
とします。これはDOSに対する命令の仕組みですので、こういう場合はこう命令
するのだと覚えるのがC言語上達の近道です。
"Hello World \n" と、"Hello World" の前に \n がないのは、画面を消去したの
で、改行する必要がないからです。
Q:ポインタが分からなくてC言語から落ちこぼれる人が多いと聞きますが。
A:そうらしいですね。
ポインタは変数や関数のアドレスをもつ変数です。
ポインタと配列との関係
pがポインタ変数で、その内容が "HELLO"であるとすると、以下のようになります。
p == "HELLO"
*p == p[0] == "H"
*(p+1) == p[1] == "E"
*(p+2) == p[2] == "L"
*(p+3) == p[3] == "L"
*(p+4) == p[4] == "O"
文字列へのポインタの例を示します。
void string(*st)
{
int i = 0;
int count = 0;
char *str;
str = st; /*文字列の代入(アドレスを一致させる)*/
while( str[i] != '\0' ) {
printf( "str[%d] = %c\n", i, str[i] );
/*1文字ずつプリントする*/
i++;
count++;
}
printf( "Number of Characters = %d \n", count );
/*文字数をプリントする*/
}
参考文献
Borland 『Borland C++ Turbo C++ ユーザーズガイド』
Microsoft 『Visual C++ ユーザーズガイド』
桜田・田口『Visual C++ 5.0 プログラミング入門』アスキー,1997.
平林雅英『新 ANSI C言語辞典』技術評論社,1997. \2,300+Tax.
Jan Walter, et al. 著/アクロバイト監訳『C++ パワフルテクニック大全集』
インプレス,1999. \3,980.
以下で、計量経済分析の Chow Testプログラムを示します。
以下のプログラムは、教育学部・人文学部1年生向けの講義「教養
原論」で構造変化を判断する手法を示すために書いたプログラムで
す。
含まれているのは、
統計の基礎: 平均、分散、共分散、決定係数
自由度、検定
行列の操作: 転置、積、逆行列
などです。
/******************************************************************
** CHOWG442.C (1995.6.5 Toshiyoshi FURUTA)
**
** OLS ESTIMATION
** y(i)=b1+b2*x2(i)+b3*x3(i)+・・・・+bk*xk(i)+e
** (i=1,2,・・・・,n)
** The Chow test
**
** Borland C++ ver3.0
*****************************************************************/*
/* data SSE df
** -------------------------------------
** COMBINED PERIODS n=n1+n2 SSE(1)=ee(1) n1+n2-k
** FORMER PERIODS n1 SSE(2)=ee(2) n1-k
** LATTER PERIODS n2 SSE(3)=ee(3) n2-k
**
** SSE(4)=SSE(2)+SSE(3) n1+n2-2k
** SSE(5)=SSE(1)-SSE(4) n-2k
** F_test
** F(k,n-2k)={SSE(5)/k)}/{SSE(4)/(n-2k)}
**
**
** References
** Jhonston,J., Econometric Methods,2nd ed.,McGraw-Hill, 1972.
** Theil,H., Principles of Econometrics, North-Holland, 1971.
**
** Test Data & Results:
** D.N.Gujarati,Basic Econometrics,2nd ed.,
** McGraw-Hill,pp.442-445.
**
******************************************************************/
#include <stdio.h>
#include <conio.h>
#include <math.h>
doubleee,ss,se,y2,dw,r2,adj_r2,ess,f_value;
double w[25][8], mean[8], var_cov[8][8], corr[8][8], ww[8][8], xx[8][8];
double b[8], e[25], t[8];
void datainput(void) ,estimation(int,int);
voidinverse(int), results(int, int);
voidpr_ww(int),pr_xxinv(int),pr_xy(int),pr_var_cov_yx(int);
voidpr_corr(int),pr_var_cov_b(int);
/*************************************************************************
** w[i][j] (i=0,1,・・・,n-1) ww[i][j] (i,j=0,1,・・・,k)
** (j=0,1,・・・,k)
**
** 0 1 2 ・・・ k 0 1 2 ・・・ k
** ------------------------ --------------------------
** 0| y(1) 1 x2(1) ・・・ xk(1) ** 0| y'y y'1 y'x1 ・・・ y'xk
** 1| y(2) 1 x2(2) ・・・ xk(2) ** 1| 1'y n 1'x1 ・・・ 1'xk
** | ・ ・ ** 2| x1'y x1'1 x1'x1 ・・・ ・
** | ・ w[i][j] ・ ** | ・ ・
** | ・ ** | ・ ww[i][j] ・
** n-1| y(n) 1 x2(n) ・・・ xk(n) ** n-1| y'xk xk'1 xk'x1 ・・・ xk'xk
**
***************************************************************************
** var_cov[i][j] (i,j=1,2,・・・,k)
**
** y x2 ・・・ xk
** -----------------
** y | 11 12 ・・・ 1k
** ・ | ・ ・ ・・・ ・
** ・ | var_cov[i][j]
** xk | k1 k2 ・・・ kk
***********************************************************/
void main(void)
{
int n,n1,n2,k;
int i,j,m;
/****************************************************
n: Data of the all periods
n1: Data of the former periods
n2=n-n1: Data of the latter periods
k: Number of Independent Variables
****************************************************/
n=18;
n1=9;
k=2;
clrscr();
datainput();/** w=[y x1 x2 ・・・・・・・・ xk] **/
printf("\n***** OLS : n=%d (1-%d) *****",n,n);
estimation(k,n);
results(k,n);
printf("\n***** OLS : n=%d (1-%d) *****",n1,n1);
estimation(k,n1);
results(k,n1);
n2=n-n1;
printf("\n***** OLS : n=%d (%d-%d) *****",n2,n1+1,n);
for(i=0; i<n2; i++)
for(j=0;j<=k; j++)
w[i][j]=w[i+n1][j];/** data transformation **/
estimation(k,n2);
results(k,n2);
}
/*------- ESTIMATION -----------*/
void estimation(int k,int n)
{
int i,j,m;
double sum;
/*------ w'w=[y x1・・・・xk]'[y x1・・・・xk] : ww[i][j] (i,j=0,1,・・・,k) -----*/
for(i=0; i<=k; i++) {
for(j=0; j<=k; j++) {
sum=0;
for( m=0; m<n ; m++ )
sum += w[m][i]*w[m][j];
ww[i][j]=sum;
}
}
/*------- x'x=w'w : xx[i][j] (i,j=1,2,・・・・・,k) --------*/
for(i=1; i<=k; i++)
for( j=1; j<=k; j++ )
xx[i][j]=ww[i][j];
pr_ww(k);
inverse(k);/* xx=ww^(-1) */
pr_xxinv(k);
pr_xy(k);
/**------- MEAN, DEVIATION : mean[]=[NUL y x2 x3 ・・・ xk]------*/
for(j=1; j<=k; j++)
mean[j]=0;
mean[1]=ww[1][0]/n;/** mean(y) = mean[1] **/
for(j=2; j<=k; j++)
mean[j]=ww[1][j]/n;/** mean(xi) = mean[i] (i=2,・・・,k) **/
/*----- VARIANCE-COVARIANCE MATRIX : var_cov[i][j] (i,j=1,2,・・・k)------*/
for(i=1; i<=k; i++)
for(j=1; j<=k; j++){
var_cov[i][j]=0;
corr[i][j]=0;
}
var_cov[1][1]= ww[0][0]/n-mean[1]*mean[1];
for(j=2; j<=k ; j++) {
var_cov[1][j]= ww[0][j]/n-mean[1]*mean[j];/** var_cov(y'x) **/
var_cov[j][1]=var_cov[1][j];
}
for(i=2; i<=k; i++) {
for(j=2; j<=k ; j++)
var_cov[i][j]= ww[i][j]/n-mean[i]*mean[j];/** var_cov(x'x) **/
}
pr_var_cov_yx(k);
/**------- CORRELATION MATRIX : corr[i][j] (i,j=1,2,・・・,k) -----*/
for(i=1; i<=k; i++){
for(j=1; j<=k ; j++)
corr[i][j]= var_cov[i][j]/sqrt(var_cov[i][i]*var_cov[j][j]);
}
pr_corr(k);
/*-------- b=[x'x]^(-1)x'y : b[i] (i=1,1,・・・,k) ------*/
for(i=1; i<=k; i++) {
sum=0;
for(j=1; j<=k; j++)
sum += xx[i][j]*ww[j][0];
b[i]=sum;
}
/*-------- e=y-xb : e[i] (i=0,1,・・・,n) --------*/
for(i=0; i<n; i++) {
sum=0;
for(j=1; j<=k; j++)
sum += w[i][j]*b[j];
e[i]=w[i][0]-sum;
}
ee=0;
for(i=0; i<n ; i++)
ee += e[i]*e[i];
ss=ee/(n-k);/*** unbiased estimator of sigma^2 */
se=sqrt(ss);/*** Standard Error of Estimation */
for(i=1; i<=k; i++)
t[i]=b[i]/sqrt(ss*xx[i][i]);/*** t-value */
pr_var_cov_b(k);/* var_cov(b) = ss*(xx)^(-1) */
sum=0;
for( i=0; i<n-1; i++ )
sum += (e[i+1]-e[i])*(e[i+1]-e[i]);
dw=sum/ee;/*** DW ratio **/
/*------- R^2 & adjusted R^2 -------*/
r2=1-ee/(n*var_cov[1][1]);
adj_r2=1.0-ss*(n-1)/(n*var_cov[1][1]);
/**********************************************************
** SSR(Regression Sum of Squares)=y'y-n*(mean_y)*^2 **
** SSE(Error Sum of Squares)=e'e **
** F={SSR/(k-1)}/SSE/(n-k)}={SSR/(k-1)}/(ee/(n-k)} **
** ={r2/(k-1)}/{(1-r2)/(n-k)} **
**********************************************************/
f_value=((n-k)*r2)/((k-1)*(1-r2));
}
/*--------------*/
void pr_ww(int k)
{
int i,j;
printf("\n -------- ww --------\n");
for( i=0; i<=k; i++ ) {
for( j=0; j<=k; j++ )
printf(" %14.5f ",ww[i][j]);
printf(" \n");
}
printf(" \n");
}
/*--------------*/
void pr_xxinv(int k)
{
int i,j;
printf("\n -------- [x'x]^(-1) --------\n");
for( i=1; i<=k; i++ ) {
for( j=1; j<=k; j++ )
printf(" %13.5f ",xx[i][j]);
printf(" \n");
}
printf(" \n");
}
/*--------------*/
void pr_xy(int k)
{
int i;
printf("\n -------- x'y --------\n");
for( i=1; i<=k; i++ ) {
printf(" %14.5f ",ww[i][0]);
}
printf(" \n");
}
/*--------------------*/
void pr_var_cov_yx(int k)
{
int i,j;
printf("\n --- variance covariance matrix : [y x2 x3 ・・・・・・ xk] ---\n");
for( i=1; i<=k; i++ ) {
for( j=1; j<=k; j++ )
printf(" %14.5f ",var_cov[i][j]);
printf(" \n");
}
printf(" \n");
}
/*--------------*/
void pr_corr(int k)
{
int i,j;
printf("\n --- correlation matrix : [y x2 x3 ・・・・ xk] ---\n");
for( i=1; i<=k; i++ ) {
for( j=1; j<=k; j++ )
printf(" %10.5f ",corr[i][j]);
printf(" \n");
}
printf(" \n");
}
/*------------------*/
void pr_var_cov_b(int k)
{
int i,j;
printf("\n --- variance covariance matrix : [b1 b2 ・・・・・・ bk] ---\n");
for( i=1; i<=k; i++ ) {
for( j=1; j<=k; j++ )
printf(" %10.5f ", ss*xx[i][j]); /* var_cov(b)=ss*(xx)^(-1) */
printf(" \n");
}
printf(" \n");
}
/*----------------- RESULTS -----------------*/
void results(int k,int n)
{
int i;
printf( "\n ----------- RESULTS -------------");
for( i=1; i<=k; i++)
printf( "\n b[%2d] = %11.5f",i,b[i]);
printf( "\n\n t-value ");
for( i=1; i<=k; i++)
printf( "\n b[%2d] = %11.5f",i,t[i]);
printf( "\n\n SSE(ee) =%11.5f",ee );
printf( "\n SE =%11.5f",se );
printf( "\n DW =%11.5f",dw );
printf( "\n R^2 =%11.5f",r2 );
printf( "\n ADJ.R^2 =%11.5f",adj_r2 );
printf( "\n F(%d,%d)=%11.5f",k-1,n-k,f_value);
printf("\n\n");
}
/*-------------- INVERSE MATRIX ------------------*/
void inverse( int order )
/** order: order of inverse matrix **/
{
/*--- omitted ---*/
}
/*-------------- DATA INPUT ------------------*/
void datainput(void)
{
/** 0 1 2 ・・・・・・・・ k w[i][j](i=0,1,2,・・・,n-1) **/
/** w=[y 1 x2 ・・・・・・・・ xk] (j=0,1,2,・・・,k) **/
/*** y[i] ***/
w[0][0]=0.36; w[1][0]=0.21; w[2][0]=0.08; w[3][0]=0.20; w[4][0]=0.10;
w[5][0]=0.12; w[6][0]=0.41; w[7][0]=0.50; w[8][0]=0.43; w[9][0]=0.59;
w[10][0]=0.90;w[11][0]=0.95;w[12][0]=0.82;w[13][0]=1.04;w[14][0]=1.53;
w[15][0]=1.94;w[16][0]=1.75;w[17][0]=1.99;
/*** x[i][j] ***/
w[0][1]=1.0; w[1][1]=1.0; w[2][1]=1.0; w[3][1]=1.0; w[4][1]=1.0;
w[5][1]=1.0; w[6][1]=1.0; w[7][1]=1.0; w[8][1]=1.0; w[9][1]=1.0;
w[10][1]=1.0;w[11][1]=1.0;w[12][1]=1.0;w[13][1]=1.0;w[14][1]=1.0;
w[15][1]=1.0;w[16][1]=1.0;w[17][1]=1.0;
w[0][2]= 8.8; w[1][2]= 9.4; w[2][2]=10.0; w[3][2]=10.6; w[4][2]=11.0;
w[5][2]=11.9; w[6][2]=12.7; w[7][2]=13.5; w[8][2]=14.3; w[9][2]=15.5;
w[10][2]=16.7;w[11][2]=17.7;w[12][2]=18.6;w[13][2]=19.7;w[14][2]=21.1;
w[15][2]=22.8;w[16][2]=23.9;w[17][2]=25.2;
}