トップページ⇒math
1001コメント292KB

〓 Mathematica 3 〓

■ このスレッドは過去ログ倉庫に格納されています
0001１３２人目の素数さん2005/07/15(金) 16:37:20
　　　　　　　　_....._{{ 〃
　　　　 　, - ' ,..､､.ヾ{{フ'⌒`ヽ､
　　　　／　 ,:', -‐‐` ´ '´⌒ヽ　ヾ:､
.　　　,'　　 ,'´　,ィ ,ｨ ,' , 　 ｀ヽ',　 ',-<
　 　 ,' 　　.i　 /|. /.| { i,　 i,　　}.　 }_,,）)
　　　!　|　 ! .,'-.{ !　!|; |`､.}ﾞ!.!　|.　 !　ヽ.
　　　',　',　|Vｧ=､ﾞ､ ｀ﾞ､!-_:ﾄ,ﾘ', l !　| 　 ﾞ',
　　　 ヽ､', l:!Kノ}. 　 　 f:_.)iﾞi: ﾘ !　l　ル　
　　　　 | l!iヽﾞ- '　, 　　.!__:ﾉ ﾞ ,ﾘ l ﾘ'´ 　 　　
.　　　　 ',|!!､ 　 　r‐┐ 　 ｀ ノ'. /,ｲ　　　　／￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣ 　
　　　　　 'i!ﾞ､ヽ､　ﾞｰ'　　_, ィ,:',:''´　　　 ＜　 Mathematica に関する話題はここに書いてね！
　　　　　　ﾞ:､ｨ､jヾｰ::: 'iﾍ　.ノ',ﾘ. 　 　 　　 ＼＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿
　　　,､- '´　ヽ､ﾞ､　　 {　｀＞"､ 　
　　/＼＼　　　　',　　 }　　 //｀ヽ

　★前スレ

〓Mathematica〓
http://science2.2ch.net/test/read.cgi/math/985023298/

〓 Mathematica 2 〓
http://science3.2ch.net/test/read.cgi/math/1078534285/
0034１３２人目の素数さん2005/07/29(金) 02:17:50
>>33
付属のMathematicaブックを順に読んでいくのが一番良いと思う。ただし、
プログラミングの素養がある程度あるという前提で。
0035１３２人目の素数さん2005/07/29(金) 11:58:26
{{0,0,0,0},{0,0,0,1},{0,0,1,0},{0,0,1,1},...,{1,1,1,1}}
のような配列を作りたいのですが、どうすればできますでしょうか？
よろしければヒントお願いします
0036１３２人目の素数さん2005/07/29(金) 15:15:36
BinDigits[b_]:=
Module[{s=ToString[BaseForm[b,2]]},
StringDrop[s,{Part[StringPosition[s,"\n"],1,1],-1}]]

BinDigitsN[b_,N_]:=Module[{l=Length[d=Characters[BinDigits[b]]]},
Join[Table[0,{N-l}],d]]

Table[BinDigitsN[n,4],{n,15}]
0037さくらﾀﾝ2005/07/29(金) 15:47:23
ちと改造

BinDigits[0]={0};
BinDigits[1]={1};
BinDigits[x_]:=Join[BinDigits[Quotient[x,2]],{Mod[x,2]}]

BinDigitsN[x_,n_]:=Module[{l=Length[d=BinDigits[x]]},Join[Table[0,{n-l}],d]]

Table[BinDigitsN[x,4],{x,15}]
0038さくらﾀﾝ2005/07/29(金) 15:54:48
ｽﾏｿ訂正
BinDigits[0]={};
BinDigits[x_]:=Join[BinDigits[Quotient[x,2]],{Mod[x,2]}]
でおｋだた
0039１３２人目の素数さん2005/07/29(金) 21:46:19
>>34
ありがとうございます。
mathematicaブックって洋書しか売っていないのでしょうか？
夏休みなので、学校のサーバに自宅からログインして使おうと
思っているので、ヘルプのマセマ本では使いにくいので
実際の本を手元に置いて使いたいのですが。
0040１３２人目の素数さん2005/07/30(土) 01:09:35
>>39
第4版は訳書がありますが、第5版は原書しかありません。
0041１３２人目の素数さん2005/07/30(土) 01:12:28
>>35
digits[n_] := Table[IntegerDigits[i, 2, n], {i, 0, 2^n - 1}];で
digits[4]とすれば求めるものが得られる。
0042１３２人目の素数さん2005/07/30(土) 22:23:53
>>40
ありがとうございました。
やっぱり、ないのですね。
第四版が古本屋で安く手に入らないか、さがしてみます。
0043１３２人目の素数さん2005/07/31(日) 01:42:52
http://documents.wolfram.com/v5/で読めるけどそれじゃ駄目なん?
0044１３２人目の素数さん2005/07/31(日) 02:08:55
>>43
>>39
0045１３２人目の素数さん2005/08/02(火) 19:48:52
グラフのアニメーションについて質問です。
1000x1000くらいの行列を記したファイルを連続で読み込んで、
そのグラフをアニメーションさせるにはどうしたらいいでしょうか。
配列にしてやってDoとShowで回せばいいと思うんですが、
ファイル"f1.dat", "f2.dat", "f3.dat"... とあった場合に
連続でImportさせるにはどうしたらよいのか分からないです・・・。
0046１３２人目の素数さん2005/08/02(火) 19:51:10
ファイルにチェインをつければいい
0047１３２人目の素数さん2005/08/03(水) 15:53:48
>>46
すみません、ファイルにチェインっていうのが
わからないんですがどういうことでしょうか？？
0048さくらﾀﾝ2005/08/03(水) 16:32:04
>ファイルにチェインをつければいい
俺も分からん、教えてんか。

とりあえず、
SetDirectory["データのあるディレクトリ"];
data = Map[Import, FileNames["*.dat"]];
で、できへん？
あ、でもファイル名の数字の部分をを"f01.dat"とかにせんと、順番が変なるわ。
0049さくらﾀﾝ2005/08/03(水) 16:34:51
数字の部分の桁数を揃えるって意味な。
0050１３２人目の素数さん2005/08/03(水) 17:35:33
>>48-49
> data = Map[Import, FileNames["*.dat"]];
Mapを使ってImportするとき、データは行列なんで"Table"で読み込みたいんですが
どうやって記述したらいいでしょうか。
諦めて手動でやろうと：
f01=Import["f01.dat","Table"];
f02=Import["f02.dat","Table"];
f03.....
として、ListDensityPlot[f01]、ListDensityPlot[f02]...で描いたグラフをアニメーションさせる、ってのは無理でしょうか。
何度もすみません。お願いします。
0051１３２人目の素数さん2005/08/03(水) 19:36:25
>>50
とりあえずTable[Import["f"<>ToString[i]<>".dat"],{i,3}]とすればとりあえず
{(f1.datのデータ),(f2.datのデータ),(f3.datのデータ),…}とはできる。
0052さくらﾀﾝ2005/08/03(水) 21:01:31
>>50
ファイルの拡張子がdatやったら、自動的にTable形式で読み込むとヘルプに書いとるから
特に指定する必要は無いと思うが、あえてやるんやったら
data = Map[Import[#, "Table"]&, FileNames["f*.dat"]];
て感じか。

>>51
そのやり方（・∀・）ｴｴﾅｧ!!
0053１３２人目の素数さん2005/08/07(日) 22:52:20
"Possible spelling error" って出力を抑える方法ってなかったでしょうか？
FAQ のような気がするんだけど、 みつからなかった。
0054352005/08/08(月) 14:08:45
>>36-38,41

遅くなりましたがありがとうございました。
大変勉強になりましたm(_ _)m
0055１３２人目の素数さん2005/08/09(火) 08:41:43
>>53
Off[General::spell];
0056１３２人目の素数さん2005/08/10(水) 23:08:50
age
0057１３２人目の素数さん2005/08/11(木) 00:34:02
この番組ピーターフランクル出てたっけ？
ときどき見てた
0058１３２人目の素数さん2005/08/11(木) 01:28:50
それはNHK教育の「マテマティカ」

これはPC用ソフトMathematica(マセマティカ)
0059１３２人目の素数さん2005/08/11(木) 17:48:08
>>57-58
ﾜﾗﾀw
0060１３２人目の素数さん2005/08/16(火) 17:42:49
複数のリストを対数軸で同じグラフにプロットしたいのですが・・・
MultipleListPlotでは対数軸は扱えないのでしょうか？？何か良い方法ないですかね。
0061１３２人目の素数さん2005/08/16(火) 17:45:18
　／⌒ヽ,　　,／⌒丶､　　　　　　　,-
　　　　　　　`,ヾ　　 ／　 　　,;;iiiiiiiiiii;､　　　＼　　　_ノｿ´
　　　　　　　　iｶ ／　　 　,;;´ 　;lllllllllllllii､　　　 ＼　ｉｶ
　　　　　　　　iｻ'　 　 　,;´　 ,;;llllllllllllllllllllii､　　　　ｆｻ
　　　　　　　　 !ｶ､._　 ,=ゞiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiii!!　＿_ｆｶﾍ.
　　　　　　　/ 　`ヾｻ;三ミミミミﾐご彡彡彡ミヾｻ`´　'i､
　　　　　　 i' 　　,.＿Ξミミミミミミき彡/////ii＿　　　|
　　　　　　 |　　;ｶ≡|ヾヾヾミミミミミぶ､//巛iﾘ≡ｶi　　|
　　　 　 　 |　　iｻ　 |l lヾヾｼヾミミﾐﾐり|iｉ//三iﾘ　`ｻi　 |
　　　　 　 |　　,ｶ ,ｶll|l l lヾリﾘﾘﾘﾘ川川|爪ミミiﾘllｶ､ｶi　 |
　 　 　 　 |　　;iｻ,ｻ |l l l リﾘ川川川川|爪ミﾐiiﾘ ｻi ｻi　 |
　 　 　 　 | 　 iｶ ;ｶ, |l l ﾘﾘﾘﾘ川川川川l爪ﾐﾐilﾘ ,ｶi ｶi　　|
　　　　 　 | 　iｻ ;ｻ, |ﾘ ﾘﾘ川川川川川l爪ミﾐiﾘ ,ｻi ｻi　　|
　　　　 　 |　 iｻ ;iｶ, | ﾘ彡彡川川川川|爪ﾐﾐｉﾘ　,ｶi :ｻ、　|
　　 　　　　,i厂　ｉｻ, |彡彡彡彡ﾉ|川川|爪ﾐミﾘ　,ｻi　`ヘ､
　　　　　 ,√　　,:ｶ,　|彡彡彡彡ﾉ川川|ゞﾐﾐﾐﾘ 　,ｶi　　 `ヾ
　　　　　´　　 　;ｻ, 　|彡彡彡彡川川ﾘゞﾐミﾘ　　,ｻi
　　　　　　　　　;ｶ,　　|彡彡彡彡リﾘﾘミミミｼ　 　,ｶi
　　　　　　　 　,;ｻ,　 　|彡彡ノリリﾘﾘミミミｼ　　　 ,ｻi
　　　　　　　 ;ﾒ'´　　　　i彡ノリﾘﾘﾘリゞﾐミｼ　　　　　`ﾍ､
　　　　　　　;ﾒ　　　　　　ヾリリﾘﾘノ巛ゞシ　　　　　　 `ﾍ､
　　　　　　;ﾒ　　　　　　　　｀`十≡=十´　　　　　　　　　`ﾍ､
0062１３２人目の素数さん2005/08/16(火) 21:13:28
nyでみつけた
感謝
0063GiantLeaves ◆6fN.Sojv5w 2005/08/16(火) 22:08:45
talk:>>62 とりあえずPlot関数とIntegrate関数とFactor関数を自作せよの刑。
0064１３２人目の素数さん2005/08/17(水) 02:02:08
Mathematica高すぎだよね
0065１３２人目の素数さん2005/08/17(水) 07:55:51
買う人が限られるからな。
0066GiantLeaves ◆6fN.Sojv5w 2005/08/17(水) 08:25:25
入学さえすれば無料でMathematicaを使えるケースもある。
0067１３２人目の素数さん2005/08/18(木) 17:33:42
Position[{a, b, c, d, e}, a]

の出力は{{1}}となりますが，1と出力させることはできないのでしょうか？

{{Pos}}=Position[{a, b, c, d, e}, a];
Pos

とすれば一応できますがなんか見苦しいもので。。。
0068１３２人目の素数さん2005/08/18(木) 18:15:19
Position[{a, b, c, d, e}, a][[1, 1]]
0069１３２人目の素数さん2005/08/21(日) 21:49:19
複数の式を1つの関数にまとめたいのですがどうすればよいでしょうか？
f[x_,y_]:=[expr1;expr2;]としてみたのですができませんでした
0070１３２人目の素数さん2005/08/21(日) 22:49:18
>>69
「複数の式を1つの関数にまとめる」の意味がわからないが、
複数の手続きを連続して実行したいなら、

f[x_,y_] := (expr1; expr2;)
0071１３２人目の素数さん2005/08/22(月) 20:47:11
{{1},{2},{3},...,{n}}
のようなリストを
{1,2,3,...,n}
としたいのですがどうすればよいですか？
Joinを使おうとしたのですがnが大きいとすべて列挙するのが難しいので・・

よろしくお願いします
0072GiantLeaves ◆6fN.Sojv5w 2005/08/22(月) 21:41:33
talk:>>71 Flatとかいう関数が無かった？あるいはPart関数とリストに要素を追加する関数を組み合わせるか。
0073GiantLeaves ◆6fN.Sojv5w 2005/08/22(月) 21:44:02
talk:>>71 http://documents.wolfram.co.jp/v5/Tour/MathematicaAsAProgrammingLanguage.ja.html
0074１３２人目の素数さん2005/08/23(火) 10:13:35
ありがとうございました
0075１３２人目の素数さん2005/08/25(木) 13:57:06
基本的なことかもしれませんが・・・
x=n(n=const)
のグラフはどのように描けばよいのでしょうか？
Plot[x = 3, {x, 0, 5}]
としてみたのですがエラーが出てy=3が描画されてしまいます
0076１３２人目の素数さん2005/08/26(金) 21:54:12
ParametricPlot[{x=3,y=t},{t,0,5}]
こんな感じで、いかがでしょう?
0077１３２人目の素数さん2005/08/30(火) 19:20:33
うまいね
0078１３２人目の素数さん2005/08/31(水) 00:29:31
age
0079１３２人目の素数さん2005/09/08(木) 21:52:47
sage
0080１３２人目の素数さん2005/09/09(金) 00:25:59
age
0081１３２人目の素数さん2005/09/09(金) 02:34:09
ﾈﾀくれ
0082１３２人目の素数さん2005/09/09(金) 12:12:43
3次元で r = θ の極座標を描きたい。
0083１３２人目の素数さん2005/09/10(土) 01:23:23
>>82
SphericalPlot3D?
0084１３２人目の素数さん2005/09/16(金) 12:53:36
分散共分散行列を1発で求める関数はありますでしょうか？
0085１３２人目の素数さん2005/09/16(金) 16:44:11
CovarianceMatrix or CovarianceMatrixMLE
0086１３２人目の素数さん2005/09/18(日) 16:31:55
質問させてください
1×n行列を転置してn×1行列にしたいのですが,
Transpose[matrix]
を使ってもエラーが出てできませんでした
どうすればよいのでしょうか？
0087１３２人目の素数さん2005/09/18(日) 18:36:47
In[1]:=Transpose[{{1,2,3}}]
Out[1]={{1},{2},{3}}
0088１３２人目の素数さん2005/09/18(日) 22:15:06
>>87
ありがとうございました
参考になりました

もう1つ質問させてください
リストの正数や負数の数をカウントしたいのですがどうすればよいのでしょうか？
Count[list , _>0]やCount[list , # > 0 &]
としてみたのですがうまくいきませんでした
0089１３２人目の素数さん2005/09/19(月) 08:06:03
Count[{0, 1, -2, 3, 4, -5}, _?Positive]
Count[{0, 1, -2, 3, 4, -5}, _?(# > 0 &)]
0090１３２人目の素数さん2005/09/19(月) 19:48:04
ある計算とその結果を出力させる関数1を別の関数2から実行するとき,
一時的に関数1の出力をなくすようなことはできますか？
0091１３２人目の素数さん2005/09/20(火) 09:04:00
age
0092１３２人目の素数さん2005/09/30(金) 02:57:30
mage
0093１３２人目の素数さん2005/09/30(金) 13:25:36
agenai
0094１３２人目の素数さん2005/10/02(日) 14:20:05
今4.1を使っているのですが、5.2にバージョンアップするだけの
価値はありますか？
0095１３２人目の素数さん2005/10/02(日) 15:29:26
公式ウェブページ見て、便利そうと思わなければ特に変えなくてもいいと思う。

地味なところでは計算精度・速度の向上もあるが、100*100の行列の演算くらいの大規模計算で無い限り差が分からないと思う。
0096942005/10/03(月) 17:56:08
>>95
ありがとうございます。
あまり利点がなさそうなのでうｐは当分見送ることにします。
0097１３２人目の素数さん2005/10/04(火) 00:42:27
>>96
取り敢ず体験版使ってみれば？
日本語化されてるし。
0098942005/10/05(水) 04:17:31
>>97
使ってみました。
特に面白そうなのはなかったんですがReduceだけは便利そうですね。
Sin(x)==Cos(x)みたいなのが解けるっていうのにはちょっと惹かれました。
0099１３２人目の素数さん2005/10/07(金) 12:38:29
リストをランダムに並び替えるにはどのようにしたら良いでしょうか？
アドバイスよろしくお願いします
0100１３２人目の素数さん2005/10/07(金) 16:18:44
NormalPValueは平均0分散1のときのp値を求めることができるようですが，
任意の平均，分散の正規分布のp値を求めたいときはどうすればいいですか？
0101１３２人目の素数さん2005/10/15(土) 06:41:18
　
0102１３２人目の素数さん2005/10/18(火) 17:11:12
AxesLabelで軸のラベルを書くとy軸のラベルはy軸と直角に表示されますが，
これを回転してy軸と水平に表示させる方法はありますか？
0103１３２人目の素数さん2005/10/19(水) 18:54:43
ListPlotで、離散的な関数を図示しているのですが、
Frame-> Trueの際にTicksを有効にしたいのですが、
なにか名案があったらお教え下さい。
0104１３２人目の素数さん2005/10/20(木) 20:00:43
age
0105１３２人目の素数さん2005/10/23(日) 23:52:15
N-joy the noiZ


sR=441*4;siZe=sR*75;timeLine=Range[siZe];t=n/siZe;
a=0.41+0.59*Exp[-3*(2t-1)^2]*(Sin[Pi*((6*Abs[t-0.5])!)])^2; mL=5;mR=11;x0=0.98;
yL=Cos[mL*ArcCos[x]];fL=Compile[{x,n},Evaluate[x+a*(yL-x)]];
dL=FoldList[fL,x0,timeLine];
ListPlot[dL,PlotStyle->{RGBColor[0.05,0.2,0.85]}];
yR=Cos[mR*ArcCos[x]];fR=Compile[{x,n},Evaluate[x+a*(yR-x)]];
dR=FoldList[fR,x0,timeLine];
ListPlot[dR,PlotStyle->{RGBColor[0.85,0.2,0.1]}];
ListPlay[{dL,0.5*dR},SampleRate->sR,SampleDepth->16]
0106１３２人目の素数さん2005/10/25(火) 02:26:22
ﾜﾚﾏﾃｨｶはどこで拾えますか？（MathID とか聞いてこないやつ）
0107１３２人目の素数さん2005/10/25(火) 18:34:26
氏ね
0108１３２人目の素数さん2005/10/25(火) 19:42:44
反応するなよ
0109１３２人目の素数さん2005/10/25(火) 19:57:35
みんな市んじまえ
0110１３２人目の素数さん2005/10/26(水) 10:35:52
>>109
どうした。悩みがあるなら聞くぞ？
0111１３２人目の素数さん2005/10/26(水) 20:57:03
>>99>>100>>102>>103
に答えてくれないかな、いや、ください。
0112１３２人目の素数さん2005/10/29(土) 16:59:19
x^2 + y^2 = r^2
を
x = rcosθ
y = rsinθ
みたいに媒介変数表示に変換したいです。
上記のような単純な円だといいのですが、
任意の楕円
ax^2 + bxy + cy^2 + dx + ey = 1
とかから媒介変数に変換するにはどうしたらいいんでしょうか？
そもそもmathematicaでできますか？
0113１３２人目の素数さん2005/10/30(日) 01:16:21
>>111
　>>99
　　RandomRelist[x_List] := Block[{n = x, p}, Do[p = Random[Integer, {1, i}]; n[[{p, i}]] = n[[{i, p}]], {i, Length[x]}]; n]
　>>100
　　HypothesisTests.mを読む
　>>102
　　RotateLabel
0114１３２人目の素数さん2005/10/31(月) 05:23:47
age
0115１３２人目の素数さん2005/11/01(火) 12:00:22
ver. 3 くそ
ver. 4 本来３だったもの
ver. 4.1 バグ取り
ver. 4.2 新たなバグ付き
ver. 5 どうでもいい
ver. 5.1 バグ取り
ver. 5.2 新たなバグ付き
ver. 6.0 さらなるバグ付き
0116１３２人目の素数さん2005/11/01(火) 20:49:47
でﾜﾚﾏﾃｨｶは？
0117１３２人目の素数さん2005/11/04(金) 22:35:23
放流するか
0118１３２人目の素数さん2005/11/04(金) 23:08:31
ny使い方わからん
0119１３２人目の素数さん2005/11/05(土) 18:30:42
>>117
Student版買ったのに・・・
登録とかしなくて使えんの？
0120１３２人目の素数さん2005/11/06(日) 01:53:51
nyはいかん．矢風億がいい．
0121１３２人目の素数さん2005/11/06(日) 06:20:31
>>120
億漁ったけどイマイチだな・・・
以前なら自サイトへのリンクとかあって色んなソフト販売してる奴らが多かったから
直接聞くことができたが、今は取り締まりが厳しいのかな。
ライムやっても罠入りばっかだしな・・・。
0122名無しさん＠Mathematica2005/11/06(日) 13:04:32
>>117
お使いのはVersion4や4.1とか4.2でしょうか
もしそうでしたら[email protected]へご連絡いただければ
便利なつーるをおわたしします
0123１３２人目の素数さん2005/11/07(月) 01:29:29
変数の値が、0.000000000000001
だった場合に、これを 0.000000000000001　
と出力させる方法はないですか？？
それか、1*10^(-14) と出力する方法はないでしょうか。

mathematica　の出力を 他のプログラムに食わせることを考えてます。
CForm が期待してるものに近いんだけど、そのまま使えないです。いまは、、CForm の出力に、
ラッパーかませてます。おしえてくださいませ。
0124１３２人目の素数さん2005/11/07(月) 07:26:44
5.1の日本語が使える鍵genはまだぁ〜
0125１３２人目の素数さん2005/11/07(月) 13:55:24
>>123
0.000000000000001 だったら

fooForm[n_] := NumberForm[n, ExponentFunction -> (Null &)]

1*10^(-14) なら

barForm[n_] := NumberForm[n,
NumberFormat -> (SequenceForm[#1, "*", #2, "^", #3] &),
ExponentFunction -> Identity
]
0126１３２人目の素数さん2005/11/08(火) 01:20:26
>>124
汗水垂らして働いて買えよ。
0127１３２人目の素数さん2005/11/09(水) 18:52:43
鼻水垂らして割れ物あさってまつ
0128１３２人目の素数さん2005/11/10(木) 19:35:06
NDsolveで得られた補完関数ってどうやって微分するのでしょうか？
0129１３２人目の素数さん2005/11/14(月) 19:11:11
D
0130１３２人目の素数さん2005/11/16(水) 22:11:26
http://www.pcsoftoff.com/2bu.htm
Mathematica5.1大安売り
0131１３２人目の素数さん2005/11/17(木) 00:08:24
>>130
IPアドレス 219.101.229.147
ホスト名 www.pcsoftoff.com
IPアドレス
　割当国　※ 日本 (JP)

国内IPか。へぇ。
0132１３２人目の素数さん2005/11/17(木) 02:34:42
>>130
ACCSに通報しますた。
0133１３２人目の素数さん2005/11/17(木) 13:59:11
日本語も使えるのかな？
多分使えないモノだなwww
■ このスレッドは過去ログ倉庫に格納されています