studystat @ ウィキ

複数の X 軸、Y 軸の使用

2013-01-08T08:45:46+09:00

複数の X 軸、Y 軸の使用

図の出力

2013-01-04T19:39:08+09:00

図が思い通りの位置に出ない場合の対処方法.

同じ列を持つ行列を作成

2012-11-12T21:23:18+09:00

同じ行を持つ行列を作成同じ列Aを持つ行列Bを作成 #blockquote{ A=[ 1 2 3 ]; B=A(ones(1,5),:); } 同じ列Aを持つ行列Bを作成 #blockquote{ A=[ 1 2 3 ]'; B=A(:,ones(1,5)); }

コマンドプロンプト

2012-01-09T05:15:41+09:00

*windowsによるフォルダ内のファイル名の取得方法 cd 取得したいフォルダのパス dir /b > 適当な名前.txt cd カレントディレクトリに移動 dirファイル名を一覧表示 dir /b ファイル名を一覧表示(拡張子なし)

グラフのプロット 1

2011-10-03T20:26:13+09:00

[[グラフのプロット]] aからbまでをプロット x軸 xlim([a b]); y軸 ylim([a b]); z軸 zlim([a b]);

最適化（黄金分割法）

2011-07-25T21:37:53+09:00

**黄金分割法黄金分割法を使って関数の極小値を求めます。（後でもう少し説明を加えます） close all axis([1,5,1,10]) hold on f=@(x) 3*(x-3).^2+2 %対象の関数 n=10; %ループ回数 %極小値を含む区間を決める %(この方法は囲い込み法と呼ばれるらしい) x(1)=10; %初期値 h=1; %ステップサイズ x(2)=x(1)+h; if f(x(2))>f(x(1)) %f(x(2))>f(x(1))ならx(1)とx(2)を入れ替え、ステップhの方向を変える change=x(2); x(2)=x(1); x(1)=change; h=-h; end for i=1:n h=2*h; x(3)=x(2)+h; if f(x(3))<=f(x(2)) x(1)=x(2); x(2)=x(3); else break; end end if h>0 a=x(1); b=x(3); else a=x(3); b=x(1); end fprintf('[a b]=[%f %f]\n',a,b) %決定された区間を出力 %黄金分割法 p=(5^(1/2)-1)/2; c=a+(b-a).*p.^2; d=a+(b-a).*p; for i=1:n if f(c)<=f(d) b=d; d=c; c=a+(b-a).*p.^2; minx(i)=c; %分割の様子をプロットしたいので結果を配列に格納しておきます miny(i)=f(c); else a=c; c=d; d=a+(b-a).*p; minx

数値計算

2011-07-12T00:15:54+09:00

*数値積分積分の方法 MATABで積分計算をする方法 quad という数値積分のコマンドがありますが、今回は台形を使った数値積分を作ってみました. #blockquote{ %台形則を使った数値積分のコマンド % S=int_a^b f(x)dx n=10000; a=-10; b=10; h=(b-a)/n; F = @(x)x.^2; S=0; for i=1:n S=S+( F(a+(h*(i-1)) ) + F(a+(h*i) ) )*(h/2) ; end S S2=quad(F,a,b) }

Latex

2012-01-09T05:48:41+09:00

*texで作成したPDFにフォントを埋め込む方法 \tex\share\texmf\fonts\map\dvipdfmx\base にある dvips dvipdfm % Ryumin and GothicBBB found in PostScript printers: 以下にある次の4行をコメントアウトし, %rml H Ryumin-Light %gbm H GothicBBB-Medium %rmlv V Ryumin-Light %gbmv V GothicBBB-Medium 以下の4行を加える. rml H :0:msmincho gbm H :0:msgothic rmlv V :0:msmincho gbmv V :0:msgothic *新しいstyファイルを追加する方法 windowsの場合 [[コマンドプロンプト]]にて (1) 以下のコマンドを実行して、パスを調べる kpsewhich platex.ltx 例: c:/usr/local/share/texmf/tex/platex/mytexstyle/cite.sty (2) styファイルを置いて (3) 以下のコマンドを実行 mktexlsr *脚注を最後にする方法以下のコマンドをプリアンブルに追加 \usepackage{endnotes} ¥let¥footnote=¥endnote 脚注を出したいところに \theendnotes と書けばOK.

グラフのプロット 3

2011-07-07T04:28:21+09:00

3つめのグラフ描画例は, 金融データを扱います. 今回のcsvデータの日付はシリアルコードではなく数値で記録されています。したがって, 数値のデータから文字列を生成し, グラフのx軸へ反映させます. 詳しい説明は後でするつもりです. 例3 株価のグラフ作成 (データは　http://finance.yahoo.co.jp/　より取得) #blockquote{ % A example of Plot 3 (Stock Price) % 2011/07/07 clear all; %--data-import--% % X= [Year M_D price volume] all data are recorded as number X= csvread('JS7203d.csv'); year=X(:,1); M_D=X(:,2); price=X(:,3); volume=X(:,4); L=length(price); %--data label--- % separate manth and day MM=floor(0.01*M_D); DD=M_D-MM*100; F=[strcat(int2str(MM(1)),'/',int2str(DD(1)))]; for i=2:L S=strcat(int2str(MM(i)),'/',int2str(DD(i))); F = char(F,S); end date=cellstr(F); %--set-number of date label--% n=5; %--plot-pre-setting-% h=floor ( ( L-1 )/( n-1 ) ); M=h*n; lab=date(1:h:M,1); %---plot---% plot(price,'b'); %--plot-setting-% legend('price (yen)'); set(gca,'xtick',1:h:L); set(gca,'xticklabel',lab); ylabel('price', 'FontSi

統計（本）

2011-07-07T16:47:06+09:00

このページでは統計学を勉強するうえで参考になると思われる書籍を紹介していきます。いまのところ、とりあえずタイトルだけ。 *統計学 **入門統計学入門 (基礎統計学) 東京大学教養学部統計学教室 (1991) 東京大学出版会 #amazon(4130420658) コア・テキスト統計学　大屋幸輔 (2003)　新世社 #amazon(4883840506) **数理統計学入門数理統計学 P.G.ホーエル (著), 浅井晃 (訳), 村上正康 (訳) (1978) 培風館 #amazon(4563008281) Linear Statistical Inference and Its Applications, C.R. Rao (1973) John Wiley & Sons #amazon(0471218758) **学術誌 Annals of Statistics Biometrika Journal of the American Statistical Association Journal of the Royal Statistical Society Series B Journal of Statistical Planning and Inference Journal of Statistical Computation and Simulation Computational Statistics & Data Analysis Communications in Statistics - Theory and Methods Communications in Statistics - Simulation and Computation