studystat @ ウィキ

メニュー

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

**メニュー -[[トップページ]] -&treemenu2(title=MATLAB){[[MATLABとシェル]]|[[プロット0]]|[[グラフのプロット]]|[[グラフのプロット 2]]|[[グラフのプロット 3]]|[[数値計算]]|[[最適化（黄金分割法）]]|[[同じ列を持つ行列を作成]][[複数の X 軸、Y 軸の使用]]} -[[統計学関連]] -&treemenu2(title=Latex){[[Latex]]|[[図の出力]]} -&treemenu2(title=本の紹介){[[確率（本）]]|[[統計（本）]]|[[プログラミング（本）]]|[[その他（本）]]} -[[メニュー]] -[[@wikiで便利そうなコマンド]] -[[コマンドプロンプト]] ---- **リンク -[[MATLAB 公式ドキュメント>>http://www.mathworks.co.jp/help/ja_JP/techdoc/]] -[[初心者によるMATLABメモ>>http://sach1o.blog80.fc2.com/#]] -[[MATLABの極 >>http://beta-infinite.blogspot.com/]] -[[@wikiご利用ガイド>>http://atwiki.jp/guide/]] // リンクを張るには "[" 2つで文字列を括ります。 // ">" の左側に文字、右側にURLを記述するとリンクになります //**更新履歴 //#recent(20) &link_editmenu(text=ここを編集)

複数の 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)); }

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.

コマンドプロンプト

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(i)=d; miny(i)=f(d); end end plot(minx,miny,'r') %途中経過を赤ライン、最終結果を黒丸でプロット plot(minx(n),miny(n),'ko') fprintf('minimum f(%f) = %f\n',minx(n),miny(n)) %最終的な結果を表示 %関数のプロット x=-10:0.001:10; plot(x,f(x)) 結果 EDU>> golden_section f = @(x)3*(x-3).^2+2 [a b]=[-4.000000 8.000000] minimum f(3.026138) = 2.002050 黄金分割で使用した点を順に赤ラインでつないでいます。 #image(http://www47.atwiki.jp/studystat/?cmd=upload&act=open&page=%E6%9C%80%E9%81%A9%E5%8C%96%EF%BC%88%E9%BB%84%E9%87%91%E5%88%86%E5%89%B2%E6%B3%95%EF%BC%89&file=golden_section.jpg,width=400,height=400,http://www47.atwiki.jp/studystat/?cmd=upload&act=open&page=%E6%9C%80%E9%81%A9%E5%8C%96%EF%BC%88%E9%BB%84%E9%87%91%E5%88%86%E5%89%B2%E6%B3%95%EF%BC%89&file=golden_section.jpg,blank)

トップページ

2011-07-12T02:08:14+09:00

studystat@wiki（仮名）は MATLABを使っていて会得したtipsや、ハマった事柄をみんなで共有するwikiです。だいたい下記のようなスタイルで運営できたらと思います。（未確定） -何かを知ったら、とりあえず書いてみましょう。「こんなの作ってみた」とか「〜についてまとめてみた」とか何でもOKです。 -ドキュメントは完璧である必要は全くありません。ただ、参考文献（ウェブページの場合URL)は書いておきましょう。 -ソースにはコメントがあれば助かりますが、目的とソースと実行結果だけでも構いません。 MATLABについて MATLAB (マットラブ/マトラボ)は数値計算の商用ソフトウェアです。行列計算が得意で、さまざまな分野で利用されているそうです。類似のフリーのソフトとして octave や scilab があります。ですから、ここで紹介している内容はoctaveやscilabユーザーにも役に立つかもしれません。 .

数値計算

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) }