〓 Mathematica 3 〓
■ このスレッドは過去ログ倉庫に格納されています
0001132人目の素数さん
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/
0022ニコちゃん ◆9COufv/UGU
2005/07/25(月) 20:46:49処理スピード、機能は同等。
ハードコピーのMathematica Bookがついてこない。
プリントアウトすると、下に"Printed by Mathematica for Students"がつく。
学生の証明書が必要。
か。 しかし、ワレマティカって言われるくらいの代物だし(以下自粛。。
0023132人目の素数さん
2005/07/25(月) 21:13:34> しかし、ワレマティカって言われるくらいの代物だし(以下自粛。。
詳しく。
何がいけないんですか?
0024132人目の素数さん
2005/07/25(月) 23:41:49一番の違いはインストール時にMathIDが要求されること。そのため異なるPCにインストールするには新たにMathIDをWolfram Researchに請求する必要がある。
Pro版はMathIDは必要なく原理的には何台でもインストールできる。(もちろんライセンス上問題はあるが…)
複数台にインストールできるからといってMathematicaをMath linkを使って並列計算させようしてもPro版だと確か4台までしかLinkできなかったはず。それ以上を望むならさらに上位バージョンが必要。
0025132人目の素数さん
2005/07/26(火) 03:00:450026132人目の素数さん
2005/07/26(火) 11:54:21なんにも自分で調べないのに5万や10万を払ってソフトを買うわけがない。
0027132人目の素数さん
2005/07/27(水) 07:04:49「3x+1問題
まず,(1以上の)自然数 n に対して次のような関数 f(n)
f(n) = n/2 : n が偶数の場合
f(n) = 3n+1 : n が奇数の場合
を考える.
これに対し,任意の自然数 n に対して,関数 f を繰り返し適用するといつか必ず 1 になるという「予想」がある.
この予想は真か偽か… というのがこの問題である.
さて,この予想の真偽を体感すべく,上の f(n) を Mathematica で 定義し,
n=1 から n=100 程度まで予想が成り立つかどうか試してみよ.
また、fm(n) = 1 なる最小の m を K(n) と書いたとき,
K(n) がいくつになるか,n=1 から n=100 程度まで調べ,グラフにせよ.」
f(n)の定義までは出来たのですが、K(n)の定義のしかたがわかりません。
どなたかご教授をお願いします。
0028132人目の素数さん
2005/07/27(水) 08:55:590029132人目の素数さん
2005/07/27(水) 09:32:58f[n_]:=n/2/;EvenQ[n]
f[n_]:=3n+1/;OddQ[n]
K[n_]:=Length[NestWhileList[f,n,#!=1&]]-1
ListPlot[Table[K[n],{n,100}]];
0030132人目の素数さん
2005/07/27(水) 10:40:51ageなくても見てる奴は見てる。2ちゃんねるブラウザがある奴なら
0031132人目の素数さん
2005/07/28(木) 16:28:20それ何?
0032132人目の素数さん
2005/07/28(木) 17:33:582ちゃんねるブラウザ の検索結果 約 35,200 件
0033132人目の素数さん
2005/07/29(金) 00:04:490034132人目の素数さん
2005/07/29(金) 02:17:50付属のMathematicaブックを順に読んでいくのが一番良いと思う。ただし、
プログラミングの素養がある程度あるという前提で。
0035132人目の素数さん
2005/07/29(金) 11:58:26のような配列を作りたいのですが、どうすればできますでしょうか?
よろしければヒントお願いします
0036132人目の素数さん
2005/07/29(金) 15:15:36Module[{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:23BinDigits[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:48BinDigits[0]={};
BinDigits[x_]:=Join[BinDigits[Quotient[x,2]],{Mod[x,2]}]
でおkだた
0039132人目の素数さん
2005/07/29(金) 21:46:19ありがとうございます。
mathematicaブックって洋書しか売っていないのでしょうか?
夏休みなので、学校のサーバに自宅からログインして使おうと
思っているので、ヘルプのマセマ本では使いにくいので
実際の本を手元に置いて使いたいのですが。
0040132人目の素数さん
2005/07/30(土) 01:09:35第4版は訳書がありますが、第5版は原書しかありません。
0041132人目の素数さん
2005/07/30(土) 01:12:28digits[n_] := Table[IntegerDigits[i, 2, n], {i, 0, 2^n - 1}];で
digits[4]とすれば求めるものが得られる。
0042132人目の素数さん
2005/07/30(土) 22:23:53ありがとうございました。
やっぱり、ないのですね。
第四版が古本屋で安く手に入らないか、さがしてみます。
0043132人目の素数さん
2005/07/31(日) 01:42:520044132人目の素数さん
2005/07/31(日) 02:08:55>>39
0045132人目の素数さん
2005/08/02(火) 19:48:521000x1000くらいの行列を記したファイルを連続で読み込んで、
そのグラフをアニメーションさせるにはどうしたらいいでしょうか。
配列にしてやってDoとShowで回せばいいと思うんですが、
ファイル"f1.dat", "f2.dat", "f3.dat"... とあった場合に
連続でImportさせるにはどうしたらよいのか分からないです・・・。
0046132人目の素数さん
2005/08/02(火) 19:51:100047132人目の素数さん
2005/08/03(水) 15:53:48すみません、ファイルにチェインっていうのが
わからないんですがどういうことでしょうか??
0048さくらタン
2005/08/03(水) 16:32:04俺も分からん、教えてんか。
とりあえず、
SetDirectory["データのあるディレクトリ"];
data = Map[Import, FileNames["*.dat"]];
で、できへん?
あ、でもファイル名の数字の部分をを"f01.dat"とかにせんと、順番が変なるわ。
0049さくらタン
2005/08/03(水) 16:34:510050132人目の素数さん
2005/08/03(水) 17:35:33> data = Map[Import, FileNames["*.dat"]];
Mapを使ってImportするとき、データは行列なんで"Table"で読み込みたいんですが
どうやって記述したらいいでしょうか。
諦めて手動でやろうと:
f01=Import["f01.dat","Table"];
f02=Import["f02.dat","Table"];
f03.....
として、ListDensityPlot[f01]、ListDensityPlot[f02]...で描いたグラフをアニメーションさせる、ってのは無理でしょうか。
何度もすみません。お願いします。
0051132人目の素数さん
2005/08/03(水) 19:36:25とりあえずTable[Import["f"<>ToString[i]<>".dat"],{i,3}]とすればとりあえず
{(f1.datのデータ),(f2.datのデータ),(f3.datのデータ),…}とはできる。
0052さくらタン
2005/08/03(水) 21:01:31ファイルの拡張子がdatやったら、自動的にTable形式で読み込むとヘルプに書いとるから
特に指定する必要は無いと思うが、あえてやるんやったら
data = Map[Import[#, "Table"]&, FileNames["f*.dat"]];
て感じか。
>>51
そのやり方(・∀・)エエナァ!!
0053132人目の素数さん
2005/08/07(日) 22:52:20FAQ のような気がするんだけど、 みつからなかった。
0055132人目の素数さん
2005/08/09(火) 08:41:43Off[General::spell];
0056132人目の素数さん
2005/08/10(水) 23:08:500057132人目の素数さん
2005/08/11(木) 00:34:02ときどき見てた
0058132人目の素数さん
2005/08/11(木) 01:28:50これはPC用ソフトMathematica(マセマティカ)
0059132人目の素数さん
2005/08/11(木) 17:48:08ワラタw
0060132人目の素数さん
2005/08/16(火) 17:42:49MultipleListPlotでは対数軸は扱えないのでしょうか??何か良い方法ないですかね。
0061132人目の素数さん
2005/08/16(火) 17:45:18`,ヾ / ,;;iiiiiiiiiii;、 \ _ノソ´
iカ / ,;;´ ;lllllllllllllii、 \ iカ
iサ' ,;´ ,;;llllllllllllllllllllii、 fサ
!カ、._ ,=ゞiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiii!! __fカヘ.
/ `ヾサ;三ミミミミミご彡彡彡ミヾサ`´ 'i、
i' ,._Ξミミミミミミき彡/////ii_ |
| ;カ≡|ヾヾヾミミミミミぶ、//巛iリ≡カi |
| iサ |l lヾヾシヾミミミミり|ii//三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
,;サ, |彡彡ノリリリリミミミシ ,サi
;メ'´ i彡ノリリリリリゞミミシ `ヘ、
;メ ヾリリリリノ巛ゞシ `ヘ、
;メ ``十≡=十´ `ヘ、
0062132人目の素数さん
2005/08/16(火) 21:13:28感謝
0063GiantLeaves ◆6fN.Sojv5w
2005/08/16(火) 22:08:450064132人目の素数さん
2005/08/17(水) 02:02:080065132人目の素数さん
2005/08/17(水) 07:55:510066GiantLeaves ◆6fN.Sojv5w
2005/08/17(水) 08:25:250067132人目の素数さん
2005/08/18(木) 17:33:42の出力は{{1}}となりますが,1と出力させることはできないのでしょうか?
{{Pos}}=Position[{a, b, c, d, e}, a];
Pos
とすれば一応できますがなんか見苦しいもので。。。
0068132人目の素数さん
2005/08/18(木) 18:15:190069132人目の素数さん
2005/08/21(日) 21:49:19f[x_,y_]:=[expr1;expr2;]としてみたのですができませんでした
0070132人目の素数さん
2005/08/21(日) 22:49:18「複数の式を1つの関数にまとめる」の意味がわからないが、
複数の手続きを連続して実行したいなら、
f[x_,y_] := (expr1; expr2;)
0071132人目の素数さん
2005/08/22(月) 20:47:11のようなリストを
{1,2,3,...,n}
としたいのですがどうすればよいですか?
Joinを使おうとしたのですがnが大きいとすべて列挙するのが難しいので・・
よろしくお願いします
0072GiantLeaves ◆6fN.Sojv5w
2005/08/22(月) 21:41:330073GiantLeaves ◆6fN.Sojv5w
2005/08/22(月) 21:44:020074132人目の素数さん
2005/08/23(火) 10:13:350075132人目の素数さん
2005/08/25(木) 13:57:06x=n(n=const)
のグラフはどのように描けばよいのでしょうか?
Plot[x = 3, {x, 0, 5}]
としてみたのですがエラーが出てy=3が描画されてしまいます
0076132人目の素数さん
2005/08/26(金) 21:54:12こんな感じで、いかがでしょう?
0077132人目の素数さん
2005/08/30(火) 19:20:330078132人目の素数さん
2005/08/31(水) 00:29:310079132人目の素数さん
2005/09/08(木) 21:52:470080132人目の素数さん
2005/09/09(金) 00:25:590081132人目の素数さん
2005/09/09(金) 02:34:090082132人目の素数さん
2005/09/09(金) 12:12:430083132人目の素数さん
2005/09/10(土) 01:23:23SphericalPlot3D?
0084132人目の素数さん
2005/09/16(金) 12:53:360085132人目の素数さん
2005/09/16(金) 16:44:110086132人目の素数さん
2005/09/18(日) 16:31:551×n行列を転置してn×1行列にしたいのですが,
Transpose[matrix]
を使ってもエラーが出てできませんでした
どうすればよいのでしょうか?
0087132人目の素数さん
2005/09/18(日) 18:36:47Out[1]={{1},{2},{3}}
0088132人目の素数さん
2005/09/18(日) 22:15:06ありがとうございました
参考になりました
もう1つ質問させてください
リストの正数や負数の数をカウントしたいのですがどうすればよいのでしょうか?
Count[list , _>0]やCount[list , # > 0 &]
としてみたのですがうまくいきませんでした
0089132人目の素数さん
2005/09/19(月) 08:06:03Count[{0, 1, -2, 3, 4, -5}, _?(# > 0 &)]
0090132人目の素数さん
2005/09/19(月) 19:48:04一時的に関数1の出力をなくすようなことはできますか?
0091132人目の素数さん
2005/09/20(火) 09:04:000092132人目の素数さん
2005/09/30(金) 02:57:300093132人目の素数さん
2005/09/30(金) 13:25:360094132人目の素数さん
2005/10/02(日) 14:20:05価値はありますか?
0095132人目の素数さん
2005/10/02(日) 15:29:26地味なところでは計算精度・速度の向上もあるが、100*100の行列の演算くらいの大規模計算で無い限り差が分からないと思う。
0097132人目の素数さん
2005/10/04(火) 00:42:27取り敢ず体験版使ってみれば?
日本語化されてるし。
009894
2005/10/05(水) 04:17:31使ってみました。
特に面白そうなのはなかったんですがReduceだけは便利そうですね。
Sin(x)==Cos(x)みたいなのが解けるっていうのにはちょっと惹かれました。
0099132人目の素数さん
2005/10/07(金) 12:38:29アドバイスよろしくお願いします
0100132人目の素数さん
2005/10/07(金) 16:18:44任意の平均,分散の正規分布のp値を求めたいときはどうすればいいですか?
0101132人目の素数さん
2005/10/15(土) 06:41:180102132人目の素数さん
2005/10/18(火) 17:11:12これを回転してy軸と水平に表示させる方法はありますか?
0103132人目の素数さん
2005/10/19(水) 18:54:43Frame-> Trueの際にTicksを有効にしたいのですが、
なにか名案があったらお教え下さい。
0104132人目の素数さん
2005/10/20(木) 20:00:430105132人目の素数さん
2005/10/23(日) 23:52:15sR=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]
0106132人目の素数さん
2005/10/25(火) 02:26:220107132人目の素数さん
2005/10/25(火) 18:34:260108132人目の素数さん
2005/10/25(火) 19:42:440109132人目の素数さん
2005/10/25(火) 19:57:350110132人目の素数さん
2005/10/26(水) 10:35:52どうした。悩みがあるなら聞くぞ?
0111132人目の素数さん
2005/10/26(水) 20:57:03に答えてくれないかな、いや、ください。
0112132人目の素数さん
2005/10/29(土) 16:59:19を
x = rcosθ
y = rsinθ
みたいに媒介変数表示に変換したいです。
上記のような単純な円だといいのですが、
任意の楕円
ax^2 + bxy + cy^2 + dx + ey = 1
とかから媒介変数に変換するにはどうしたらいいんでしょうか?
そもそもmathematicaでできますか?
0113132人目の素数さん
2005/10/30(日) 01:16:21>>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
0114132人目の素数さん
2005/10/31(月) 05:23:470115132人目の素数さん
2005/11/01(火) 12:00:22ver. 4 本来3だったもの
ver. 4.1 バグ取り
ver. 4.2 新たなバグ付き
ver. 5 どうでもいい
ver. 5.1 バグ取り
ver. 5.2 新たなバグ付き
ver. 6.0 さらなるバグ付き
0116132人目の素数さん
2005/11/01(火) 20:49:470117132人目の素数さん
2005/11/04(金) 22:35:230118132人目の素数さん
2005/11/04(金) 23:08:310119132人目の素数さん
2005/11/05(土) 18:30:42Student版買ったのに・・・
登録とかしなくて使えんの?
0120132人目の素数さん
2005/11/06(日) 01:53:510121132人目の素数さん
2005/11/06(日) 06:20:31億漁ったけどイマイチだな・・・
以前なら自サイトへのリンクとかあって色んなソフト販売してる奴らが多かったから
直接聞くことができたが、今は取り締まりが厳しいのかな。
ライムやっても罠入りばっかだしな・・・。
■ このスレッドは過去ログ倉庫に格納されています