トップページ⇒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/
0210１３２人目の素数さん2006/03/16(木) 22:50:48
実験データをy=a/x+bの式にフィッティングしたいのですが、
データの中に{x.y}={0、100}を含んでいるので
無限式が出ると言われて上手くいきません。
どのようにすればいいのでしょうか。
0211１３２人目の素数さん2006/03/16(木) 23:20:23
>>210
0のでるデータを外す。
普通、実数は0で割れないでしょ。
02122102006/03/17(金) 20:11:16
>>211
ありがとうございます。やはりそうするしかないのでしょうか。
データ数が少ない、かつx=0付近が重要なのでそれは避けたいのですが。
0213１３２人目の素数さん2006/03/23(木) 09:38:11
>>210
自分で(フィッテングの)評価関数書いて、評価関数の最適化だけやらせるとか。
自分で書けば、「x=0付近が重要」と云うのも、それなりに盛り込めるでしょ？
0214１３２人目の素数さん2006/03/26(日) 15:03:14
0215１３２人目の素数さん2006/04/03(月) 21:43:37
そろそろ、新しいMathematicaでるのかな？
http://www.wolfram.co.jp/news/events/JMC2006/
0216１３２人目の素数さん2006/04/05(水) 17:20:18
ある数列が与えられていて、それがベータ関数の値であるときに、
元のベータ関数のパラメータはどうやって求めればいいんでしょうか。
V5.2を使っています。
0217１３２人目の素数さん2006/04/14(金) 21:03:53
wktk
0218１３２人目の素数さん2006/04/15(土) 12:01:53
age
0219１３２人目の素数さん2006/04/16(日) 01:39:42
0220１３２人目の素数さん2006/04/17(月) 01:42:18
Ver6ﾏﾀﾞｰ
0221１３２人目の素数さん2006/04/23(日) 22:59:53
ﾏｰﾀﾞﾀﾞﾖｰ
0222１３２人目の素数さん2006/04/29(土) 22:05:27
マスマスティかで不等式を2DPlotで描くにはどうすればいいの？
領域の色づけがヘルプを見ても分からない。
0223１３２人目の素数さん2006/05/01(月) 10:41:08
Mathematicaでアイパターン描こうとしているのですがなかなかうまくいきません
やり方としてはNRZのランダムパルスを発生させそれを離散フーリエ変換してナイキストフィルタのリストと乗算
それを逆フーリエ変換して同一平面上にプロットすれば良いと思うのですがほかにやり方あったら教えてください

>>222
filledplot
0224１３２人目の素数さん2006/05/10(水) 12:56:10
√(sin 2x＋2cos 5x)のグラフ書くために初めてmathematica使ったんですけど、
グラフがｘ軸に接しないのはなぜでしょうか？
0225１３２人目の素数さん2006/05/10(水) 22:20:26
>>224
接するグラフではないのだろう
0226１３２人目の素数さん2006/05/11(木) 16:36:46
FindRootの結果をグラフにしたいのですが、なんかスマートな方法ないですか？

普通にFindRootすると戻り値は
{hoge -> 51.6438}
みたいな文字列になってしまってるので困ってます。
0227１３２人目の素数さん2006/05/11(木) 21:16:10
>>226
ReplaceAllを使う
x /. FindRoot[x^2 == 2, {x, 1}]
みたいな感じで
0228１３２人目の素数さん2006/05/13(土) 12:15:14
>>224
Solve[(Sin[2x] + 2 Cos[5x])^0.5 == 0, x]
の結果、実数の範囲では
{x -> -1.5707963267948966`} {x -> 1.5707963267948966`}

Plotでは限界がある計算をさせてるみたい。
0229１３２人目の素数さん2006/05/13(土) 21:21:13
次のMathematicaの情報知ってる人いますか？
4.2から5のときのReduceはすごかったけどそれに匹敵するようなグレードアップ
があるといいなあ。
0230１３２人目の素数さん2006/05/14(日) 23:33:37
二つのベクトルから行列を作る方法を教えてください。
v . w^t のような感じで、v.Transpose[w] だとできません。
0231１３２人目の素数さん2006/05/15(月) 01:30:13
>>230
Outer[Times, w, v]
0232１３２人目の素数さん2006/05/18(木) 17:30:12
x^2(x-3)(2x+3)と２点で接する半径1.5の接点を求めたいのですがどうすればよろしいですか？
0233１３２人目の素数さん2006/05/20(土) 00:31:01
このソフトの最新バージョンナンバーはいくつですか？
0234１３２人目の素数さん2006/05/20(土) 00:34:29
5.2
0235１３２人目の素数さん2006/05/22(月) 16:22:35
ある曲線に、２点で接する円を書きたい時はどうすればいいんですか？
0236１３２人目の素数さん2006/05/24(水) 19:52:20
うちの教授が

「マテマチカ」

と呼ぶアプリはこれのことでしょうか？
0237１３２人目の素数さん2006/05/25(木) 16:19:16
ループを回して計算しているその途中経過をグラフィカルに表示するにはどうすればいいですか
計算が進むにつれ折れ線グラフがにゅーんと延びていくとかできればベストです。
0238１３２人目の素数さん2006/05/25(木) 19:41:40
>>236
その発音も間違いじゃない。
0239１３２人目の素数さん2006/05/26(金) 08:58:14
本屋さんの名前と混同してるんだろう
0240１３２人目の素数さん2006/05/27(土) 06:35:34
Mathematicaは線形連立方程式を数値的にではなく
代数的に解くことはできますか？
つまり
x=a+b
y=a-b
を解いて
a=(x+y)/2
b=(x-y)/2
のような答を、a,bに具体的な値を入れなくても
返してくれるのでしょうか？
0241１３２人目の素数さん2006/05/27(土) 06:40:20
できる
0242１３２人目の素数さん2006/05/27(土) 12:10:53
そろそろ、ver6.0が出るのかな？
0243１３２人目の素数さん2006/05/29(月) 23:55:30
n==0のとき、 f[x_,n_]:=1
それ以外のとき、f[x_,n_]:=x^n
となるようにf[x_,n_]を定義したいのですが、
どうしたらいいのでしょうか？
02442432006/05/30(火) 00:03:13
ちょっと調べたら見つかった
f[x_, n_] := If[n == 0, 1, x^n]
スレ汚しｽﾏｿ
0245１３２人目の素数さん2006/05/31(水) 14:40:52
出力がリストになるような関数をコンパイルするにはどうすればよいのでしょうか。
下記のようにうまくいきません


以下失敗例

In[1]:= test[X_, Y_] := {X + Y, X - Y}
In[2]:= comp = Compile[{X, Y}, test[X, Y]]
Out[2]= CompiledFunction[{X, Y}, ttest[X, Y], "-CompiledCode-"]
In[3]:= comp[1, 2]
CompiledFunction::cfse: コンパイルされた式(3.,-1.)はmachine-size real numberであるべきです.
CompiledFunction::cfex: "式2の外部評価においてエラーが発生しました．コンパイルなしで評価を続行します．"
Out[3]={3, -1}
0246１３２人目の素数さん2006/06/06(火) 21:04:38
>>242
もう機能満載で大きな変化はないんとちゃう？　
0247１３２人目の素数さん2006/06/07(水) 21:13:56
Mapleみたいに、外観が変わるんじゃない?
0248１３２人目の素数さん2006/06/10(土) 08:40:56
mapleはだんだん教育ツール化している気がする。
mathematicaはどうだろう？
0249１３２人目の素数さん2006/06/11(日) 15:07:04
計算ツールになってる。
0250１３２人目の素数さん2006/06/15(木) 01:51:12
チュートリアルの例文そのまま
result = NDSolve[{
y''[t] + 0.1 y'[t] == -Sin[y[t]], y[0] == 1, y'[0] == 0}, y, {t, 0, 50}]
をコピペしただけなのに
NDSolve::deqnの最初の引数にはFalseではなく，方程式か方程式のリストが必要です．＼"＼>"}]＼)
とかエラーが出やがるんですがどうしてでしょうか？
0251羽村2006/06/15(木) 16:20:27
>>250
Clear[y];
して再チャレンジ
0252UNI2006/06/16(金) 23:17:15
Mathematicaでアポロニウスの円を書きたいんですけど、どうすればよいでしょうか？
0253１３２人目の素数さん2006/06/17(土) 11:53:30
>>252

Print["アポロニウスの円"]
0254１３２人目の素数さん2006/06/17(土) 13:04:14
>>252
定義通りに式を書いて、それをPlotすればいいのでは？
0255１３２人目の素数さん2006/06/17(土) 20:57:27
噂に聞くMathematicaそのうち買おうと思ってますが、ver.6そろそろ出るそうでないですか。
メジャー ver upによってエンバグされたりすると悲しいので、今のうちに買おうかなと思ってるんですが、
何分はじめてなもんで、この判断が妥当かどうか微妙です。

皆様の過去のver up体験談(4 ==> 5で酷い目にあった or 全く問題なかった)などありましたら
聞かせて頂けないでしょうか
0256１３２人目の素数さん2006/06/17(土) 21:18:10
>>255
質問とはそれるけれど、
学生版以外は、Premier Serviceが無料で付くので、
ver.6に無料で変えれるよ。
0257１３２人目の素数さん2006/06/17(土) 21:52:50
>>255

ver4.0==>ver5.0(Win)で酷い目にあったこと
・立ち上げたときにピロロロンと鳴らなくなった
・境界条件に文字を含んだ簡単な微分方程式が解けなくなった
などなど
0258１３２人目の素数さん2006/06/17(土) 22:31:22
>>257
>立ち上げたときにピロロロンと鳴らなくなった
これは、環境設定のところで音を鳴らすことができる。

ver5.0はバグが多かったなぁ。
02592552006/06/17(土) 23:52:01
お返事ありがとうございました。

>>256
Academic(~20万)で買おうと思っとるんですが、どうもPremier Serviceは別売(~10万)のようです。
(今なら2万円程度おトクのようですが。。。)

>>257,258
むむ、やはり不具合などありましたか。

アドバイスありがとうございます。(比較的)枯れてる現バージョン購入の方向でいってみます
0260UNI2006/06/19(月) 22:46:12
お願いします！！Mathematicaを使って書けっていう問題が出ました・・・
問１　２点からの距離の比が一定な点の軌跡（アポロニウスの円）のグラフをかけ。
問２　１点と１直線から等距離にある点の軌跡（放物線）のグラフを書け。

コマンド教えてください！！
お願いします☆
0261GiantLeaves ◆6fN.Sojv5w 2006/06/20(火) 08:13:45
talk:>>260 軌跡がどうなるかあらかじめ計算して入力すればいいのではないのか？
0262ＵＮＩ2006/06/20(火) 15:40:06
Mathematica上で、計算していきたいんです!!
お願いします☆
0263１３２人目の素数さん2006/06/20(火) 20:08:33
\!\(Appolonius[m_, n_] := Module[{a, b, r}, \[IndentingNewLine]a =
1 + n\^2\/\(m\^2 - n\^2\); \[IndentingNewLine]b =
0; \[IndentingNewLine]r =
1 + n\^2\/\(m\^2 - n\^2\) - m\/\(m + n\); \[IndentingNewLine]Show[
Graphics[{Point[{0, 0}], Point[{1, 0}], Circle[{a, b}, r], Point[{a,
b}]}, AspectRatio -> Automatic]]]\)
0264UNI2006/06/25(日) 11:54:58
Mathematicaで、ｘ軸に垂直な直線を書きたいんですけど、どうすればいいですか？！
0265１３２人目の素数さん2006/06/25(日) 20:57:35
お前、それは中学校の数学で習う範囲だろ('A`)
小学生だってx=1とかx=0は知ってる人は知ってる
0266１３２人目の素数さん2006/06/25(日) 23:48:54
<<264

ｘ軸に垂直な直線、例えば、x=1 を書きたいとき、

<<Graphics`

ImplicitPlot[x==1,{x,-5,5},{y,-5,5},AspectRatio -> Automatic]
0267１３２人目の素数さん2006/06/26(月) 17:29:00
２０行１列のリストをＸ軸に、
もう一つの２０行１列のリストをＹ軸に、といったプロットって出来ますか？
簡単に出来そうに見えて出来ない・・・
0268１３２人目の素数さん2006/06/27(火) 02:18:25
>>267
x={{x[1]},...,{x[20]}};
y={{y[1]},...,{y[20]}};
ListPlot[Flatten[Transpose[x_List]],Flatten[Transpose[y_List]]]
でいいんじゃない？
02692672006/06/27(火) 12:48:02
>>268さん
レスありがとうございます。
やってみたのですが、うまくいきませんでした・・・。
前の書き込みでは説明不足だったのですが、
x={{x[1]},...,{x[20]}};
y={{y[1]},...,{y[20]}};
のところは、私がやろうとしていたのは、
Excelのデータから抜き出したデータの羅列、例えば
x={{1},{2},{3}};
y={{3},{4},{5}};
のような感じで、番号を振ることができないのでうまくいかないのかなと思います。
ヘルプから、これらのようなxとyを、
n={{1,3}{2,4}{3,5}};
といった行列に直せば良いことは分かっていたのですが、この変換がクセモノで。。

勉強不足と言われればそれまでですが、
妙案がございましたら、レスしていただけると幸いです。
02702682006/06/27(火) 13:00:27
>>269
そういうnを得たいんだったら、
n=Transpose[{Flatten[x],Flatten[y]}]
でどう？
02712672006/06/27(火) 13:28:27
>>268さん
素早いレスありがとうございます！
希望の通りの行列ができました！！
本当に、ありがとうございます！！

もっと精進していこうと思います。
0272１３２人目の素数さん2006/06/28(水) 03:08:19
アルキメデスのらせんを書く方法を教えて下さい。
もうさっぱりです。
0273１３２人目の素数さん2006/06/28(水) 08:23:47
極座標で
r=nθ
を表現すればいい。
0274１３２人目の素数さん2006/06/28(水) 21:09:38
すいません、極座標の表現の仕方がわかりません。

x:=Cos[θ]
y:=Sin[θ]

とおいて、置換で書いた見たんですが直線になってしまいます。
02752682006/06/29(木) 03:02:00
>>272
<<Graphics`Graphics`
PolarPlot[t,{t,0,2 Pi}];
0276１３２人目の素数さん2006/06/29(木) 19:31:46
>>275
書き方がわかりました。
ありがとうございました。
0277１３２人目の素数さん2006/06/30(金) 12:55:28
等高線に対しての流線の描き方がわかりません。
よろしくお願いします。
0278１３２人目の素数さん2006/07/04(火) 12:10:11
正の整数 N に対して，√N を根号の中ができるだけ小さな整数に
なるように変形して a√b の形に直すとき，この整数 a と b の
求め方を教えて下さい．よろしくお願いします．
0279１３２人目の素数さん2006/07/04(火) 23:50:17
>>277,278
かんがれ!
0280１３２人目の素数さん2006/07/06(木) 14:23:48
>>278
root[27]で実行すれば、3root[3]とかなるが、こういうことじゃないのか？
0281１３２人目の素数さん2006/07/12(水) 22:08:19
あるセルを評価したときに結果がとても膨大になってしまうことってありますよね。
そういうときヘルプの中のセルのように、見出しだけ表示しておいて、
横のところをダブルクリックすると、中身が出てくるようにしたいのですが、
どなたかやり方をご存知の方いらっしゃいますでしょうか。
どうかよろしくお願いいたしますm(_ _)m
0282１３２人目の素数さん2006/07/13(木) 23:07:49
>>280
その 3 と 3 とを取り出すにはどうしたらいいですか？
0283１３２人目の素数さん2006/07/13(木) 23:49:23
「Print」で出力した数値の羅列を、Excelに出力する方法が分からないよ〜え〜ん。。。
0284１３２人目の素数さん2006/07/14(金) 01:14:15
>>282
私も初心者なので分かりませんが。
f[n_]:=If[Head[root[n]===Integer,root[n],If[Head[root[n]===Power,1,root[n][[1]] ]]
でよいはず。

root[27]//FullForm は Times[3,Power[... なので root[27][[1]]=3.
0285１３２人目の素数さん2006/07/15(土) 10:59:08
>>283
tsvとかで別ファイルに書き出したらいいのでは？
0286１３２人目の素数さん2006/07/22(土) 01:53:35
Piのリストを作りたいのですがどういう関数使えばいいですか？
3.141592653を[3,1,4,1,5,9,2,6,5,3]みたいにしたいです
0287１３２人目の素数さん2006/07/22(土) 04:37:25
List
0288１３２人目の素数さん2006/07/22(土) 14:46:38
RealDigits
0289１３２人目の素数さん2006/07/22(土) 15:58:10
対数関数のグラフの表示方法を教えてください。
0290１３２人目の素数さん2006/07/22(土) 16:12:38
Plot
0291１３２人目の素数さん2006/07/22(土) 16:17:34
LogPlot[Exp[-x] + 4 Exp[-2x], {x, 0, 6}]
で表示できないのですが…
02922862006/07/22(土) 16:24:56
>>288
どもです！
0293１３２人目の素数さん2006/07/22(土) 16:45:05
対数関数のグラフの表示方法を教えてください。
LogPlot[Exp[-x] + 4 Exp[-2x], {x, 0, 6}]
で表示できないのですが…
0294１３２人目の素数さん2006/07/22(土) 16:58:21
>>293

<< Graphics`Graphics`
LogPlot[Exp[-x] + 4 Exp[-2x], {x, 0, 6}]
0295１３２人目の素数さん2006/07/22(土) 17:22:50
>>294
ありがとう御座います！

<< Graphics`Graphics`はパッケージを読み込むってことですよね？

LogLinearPlot[60-Log[x], {x, 10^0, 10^30}]
y=60-Log[x]
縦軸y（線形）と横軸x（対数）
で表示したいのですが、できませんでした。
どこが問題なのでしょうか？
0296１３２人目の素数さん2006/07/27(木) 23:52:08
複素数の計算とかはできないのでしょうか？
たとえば|z|^2をzで微分するとconj(z)になったりとかそういうやつです。
0297１３２人目の素数さん2006/07/28(金) 18:41:40
471
0298１３２人目の素数さん2006/07/31(月) 01:56:04
「ListPlot3Dはデフォルトで各データ点のx座標とy座標の値が、1から始まる連続する整数であると解釈する」
とあるのですが（ヘルプ）、これを変更する方法をご存知の方はいらっしゃいますか？
教えて頂けるとすごくうれしいです。

具体的に問題点を指摘しておきます。
Plot3D[x^2+y^2,{x,0,1},{y,0,1}]
と同じグラフィックスをListPlot3Dで描くということが課題です。
partition=24;ListPlot3D[Table[x^2+y^2,{x,0,1,1/partition},{y,0,1,1/partition}]]
を描くと、x軸・y軸ともに1から25までの整数値になります。これを0から1までの実数値にすることが課題です。
ちなみにTableのx^2+y^2を{x,y,x^2+y^2}にすればうまくいくかなと思ったのですがダメでした。
２次元（ListPlot）では{x,x^2}とすると思ったような描画ができるのに・・・拡張性が無いですね（改善してほしいです）。
0299１３２人目の素数さん2006/07/31(月) 10:03:02
i=1->25⇔i-1=0->24⇔x=(i-1)/24=0->1
0300１３２人目の素数さん2006/07/31(月) 10:38:17
>>299
>>298が指摘してる問題点はそういうことではないだろ
具体的にListPlot3Dで課題を実現してみようとすればわかる
03012992006/07/31(月) 11:54:24
>>300
機械は苦手なもので申し訳ありません。
ところで、定義域を{x,0,24}にするのと{x,1,25}にするのでは
目盛りは一緒なのですが、値は違うようです。
たとえば{x,-24,24}なんかにすると分かりやすいですが。

目盛りを変える方法まではわからないなぁ。
0302１３２人目の素数さん2006/07/31(月) 15:50:02
>>298
いろいろ試してみたけどむずかしいね
でも「デフォルト」って書いてあるくらいだから変更できるはず

>>301
日本語でおｋ
0303１３２人目の素数さん2006/07/31(月) 22:31:08
うーん、あとはTicks->ほにゃららで目盛りをいじる
ことぐらいしか思いつかん。
0304１３２人目の素数さん2006/07/31(月) 22:53:17
あー、わかった。

Help で MeshRange をみれ
0305１３２人目の素数さん2006/08/01(火) 00:38:00
楕円曲線における解（楕円曲線暗号）の求め方ってどうやるんでしょうか？
教えてください。
0306１３２人目の素数さん2006/08/01(火) 00:40:40
教員をやってます。バージョン４です。
初等幾何で表示する図形をマセマティカで表示使用と思っているのですが、
やはり、2次元の座標上に図形の方程式を使ってえがくしかないんでしょうか？
よろしくお願いいたします。
0307１３２人目の素数さん2006/08/01(火) 00:52:56
>>306
多分そうだと思う。
初等幾何の図を書くなら、ドローソフトを使った方が良いと思います。
03082982006/08/01(火) 04:16:47
>>304
ｷﾀ*･゜ﾟ･*:.｡..｡.:*･ﾟ(≧∀≦)ﾟ･*:.｡. .｡.:*･゜ﾟ･*!!!!
ｷﾀ*･゜ﾟ･*:.｡..｡.:*･ﾟ(≧∀≦)ﾟ･*:.｡. .｡.:*･゜ﾟ･*!!!!
ｷﾀ*･゜ﾟ･*:.｡..｡.:*･ﾟ(≧∀≦)ﾟ･*:.｡. .｡.:*･゜ﾟ･*!!!!
ｷﾀ*･゜ﾟ･*:.｡..｡.:*･ﾟ(≧∀≦)ﾟ･*:.｡. .｡.:*･゜ﾟ･*!!!!
ｷﾀ*･゜ﾟ･*:.｡..｡.:*･ﾟ(≧∀≦)ﾟ･*:.｡. .｡.:*･゜ﾟ･*!!!!

どうみても0-1区間です。
本当にありがとうございました。
0309１３２人目の素数さん2006/08/01(火) 23:03:59
>>306
そういう目的なら、シンデレラが良いかも。フリーじゃないけど、
書籍扱いで安い。
http://www.amazon.co.jp/gp/product/4431711201/
http://www.amazon.co.jp/gp/product/4431709665/
■ このスレッドは過去ログ倉庫に格納されています