〓 Mathematica 七 〓 [転載禁止]©2ch.net
レス数が1000を超えています。これ以上書き込みはできません。
b^2/a^4 を (b/a^2)^2 に変形するみたいなのってどういうコマンド使えば良いんだろ? 数値最小化で
NMinimize[目的関数、制約条件]、変数]
というとき、たとえば制約条件がλ_1>0、λ_2>0、・・・、λ_100>0のとき、
λ_k>0,{k,1,100}と表したいのですが、うまくいきません。
制約条件の数が少なければλ_1>0、λ_2>0、・・・でよいですが、
1個ずつ入れるには多すぎるのです。
ご教示の程よろしくお願いいたします。 自己解決したでござるよ
Tableを使って
Table[λ_k, k] >= 0, {k, 1, 100}]}
でござるな {(a : Blank[1]) > 0, (a : Blank[2]) > 0, (a : Blank[3]) >
0, (a : Blank[4]) > 0, (a : Blank[5]) > 0, (a : Blank[6]) >
0, (a : Blank[7]) > 0, (a : Blank[8]) > 0, (a : Blank[9]) >
0, (a : Blank[10]) > 0} Sum[λ_k,{k,1,10}]=15
λ_k は自然数とする。
このとき解の個数を求めよ >>727
Flatten[Permutations /@ IntegerPartitions[15, {10}], 1] // Length
-> 2002
かな。 >>728
素晴らしいですね。
わたしは頭を使うのが苦手で
n100 = 10; A = Array[a, {n100}];
Solve[({(A // Total) == 15, Table[a[ll] > 0, {ll, 1, n100}]}) //
Flatten, A, Integers] // Length
でやりました。
貴君のほうが早く答えが出ますね。
ありがとうございました。 Enterキーのみで評価開始しちゃうときがあるんだけど、
これを止める方法誰かおしえてー plotArealElement[r_, phi_, dr_, dphi_] := Module[
{g1, g2, x, y, u, v, plr},
g1 = ParametricPlot[s {Cos[phi], Sin[phi]}, {s, r, r + dr}, {psi, phi, phi + dphi}];
u = r {Cos[phi], Sin[phi]};
v = r {Cos[phi + dphi], Sin[phi + dphi]};
x = v - u;
v = (r + dr) {Cos[phi], Sin[phi]};
y = v - u;
plr = Options[g1, PlotRange];
g2 = ParametricPlot[s x + t y, {s, 0, 1}, {t, 0, 1}, plr];
Show[g1, g2]
] >>744
g1とg2を同じPlotRangeで描きたいと思い、↑のコードを書きましたが、
エラーが出てしまいます。
どうすればいいでしょうか? >>745
g2 = ParametricPlot[s x + t y, {s, 0, 1}, {t, 0, 1}, plr];
↓
g2 = ParametricPlot[s x + t y, {s, 0, 1}, {t, 0, 1}];
で動作確認したら 例えば、
Subscript[x, 1]
を変数として使いたいのですが、
別のところで、
x = 3
などとやるとエラーが出ます。
どうすればいいでしょうか? 私の環境では問題なく使えますが、、。
例えば
Subscript[x, 1]=123
とすれば
Subscript[x, 1]
は123を返してきます。
しかし、ただのxで参照したら別の変数でしょう? あ、
f[a]=321
とすれば
f[a]
は321を返すんだから説明になってなかったですね。
でも
x_1(ここでは表現が難しい、、)は123を返してきます。 なぜ以下の関数は引数のリストの要素の和を計算するのでしょうか?
sumList[lis_]:=First[lis//.{x_, y___} -> x + {y}] >>750
sumListの右辺で最初にx+{y}が計算されて{x+y}となるみたいだ。
yはSequenceが入る。
例えば
x+Sequence[a,b,c]は
x+a+b+c
となるようだからlistの要素の総和になる。 3 + {Sequence[1, 2]}
を計算すると、
{4, 5}
になりました。 >>751
>sumListの右辺で最初にx+{y}が計算されて{x+y}となるみたいだ。
これはなぜなのでしょうか? {x,y,z,w}//.{x_, y___} -> x + {y}]
.....{x+y+{z,w}}---> x+y+z+w >>753
https://reference.wolfram.com/language/ref/Rule.html
によればRule“->”の右辺は即事に評価すると書いてあるから、
x+{y}はx,yに代入される前に評価されて{x+y}になり
そのあと代入されて、計算が続くということだろう。 Poor man's Mathematica used on Raspberry Pi FindEquationalProof
いつの間にか定理の自動証明ができるレベルになってたんですね
群の証明なんかが自動的にできてて驚いた 金持ちの数学者>>>>>>>>>。。。。。>>>貧乏人の数学者>>>>>2chの数学者 なんか久しぶりに数学勉強したくなってきた。
高校1年生の2次関数勉強中。サポートよろ! http://tvcap.info/2019/2/5/n2i9190205-2243560390.jpg
http://tvcap.info/2019/2/5/n2i9190205-2244010380.jpg
http://tvcap.info/2019/2/5/n2i9190205-2244040827.jpg
http://tvcap.info/2019/2/5/n2i9190205-2244240232.jpg
http://tvcap.info/2019/2/5/n2i9190205-2244110703.jpg
http://tvcap.info/2019/2/5/n2i9190205-2244390628.jpg
http://tvcap.info/2019/2/5/n2i9190205-2244540005.jpg
http://tvcap.info/2019/2/5/n2i9190205-2245000222.jpg
http://tvcap.info/2019/2/5/n2i9190205-2245040254.jpg
http://tvcap.info/2019/2/5/n2i9190205-2245180265.jpg
http://tvcap.info/2019/2/5/n2i9190205-2245250755.jpg
http://tvcap.info/2019/2/5/n2i9190205-2245300342.jpg 【超悪質!盗聴盗撮・つきまとい嫌がらせ犯罪者の実名と住所を公開】
@井口・千明(東京都葛飾区青戸6−23−16)
※盗聴盗撮・嫌がらせつきまとい犯罪者のリーダー的存在/犯罪組織の一員で様々な犯罪行為に手を染めている
低学歴で醜いほどの学歴コンプレックスの塊/超変態で食糞愛好家である/醜悪で不気味な顔つきが特徴的である
A宇野壽倫(東京都葛飾区青戸6−23−21ハイツニュー青戸202)
※色黒で醜く太っている醜悪黒豚宇野壽倫/低学歴で人間性が醜いだけでなく今後の人生でもう二度と女とセックスをすることができないほど容姿が醜悪である
B色川高志(東京都葛飾区青戸6−23−21ハイツニュー青戸103)
※色川高志はyoutubeの視聴回数を勝手に短時間に何百何千時には何万回と増やしたり高評価・低評価の数字を一人でいくつも増やしたり減らしたりなどの
youtubeの正常な運営を脅かし信頼性を損なわせるような犯罪的業務妨害行為を行っています
※色川高志は現在、生活保護を不正に受給している犯罪者です/どんどん警察や役所に通報・密告してやってください
【通報先】
◎葛飾区福祉事務所(西生活課)
〒124−8555
東京都葛飾区立石5−13−1
рO3−3695−1111
C清水(東京都葛飾区青戸6−23−19)
※低学歴脱糞老女:清水婆婆 ☆☆低学歴脱糞老女・清水婆婆は高学歴家系を一方的に憎悪している☆☆
清水婆婆はコンプレックスの塊でとにかく底意地が悪い/醜悪な形相で嫌がらせを楽しんでいるまさに悪魔のような老婆である
D高添・沼田(東京都葛飾区青戸6−26−6)
※犯罪首謀者井口・千明の子分/いつも逆らえずに言いなりになっている金魚のフン/親子孫一族そろって低能
E高橋(東京都葛飾区青戸6−23−23)
※高橋母は夫婦の夜の営み亀甲縛り食い込み緊縛プレイの最中に高橋親父にどさくさに紛れて首を絞められて殺されそうになったことがある
F長木義明(東京都葛飾区青戸6−23−20) ※日曜日になると風俗店に行っている Mathematicaは手続き型ではなく関数型でコードを書いたほうが高速化するらしいけど、並列計算させたいとき
関数型で書けばうまくいきますかね? そうだよ。関数型で書いて最後に //Parallelize を付けるだけで分散処理してくれる。 最後に //Parallelize を付けるだけ
ってのはMathematicaじゃないような気もするけど gridのこと?
自マシンのマルチコアプロセッサー利用なら何もいらない
正規ライセンスあればLightweight Gridも利用できる これはマセマティカがタダになったという理解でいいの? >>777
入力はキャラクタベースのTexでやるってこと? Potential[x_] = 0.5*k*x^2;
Threshold[n_] = 10*n;
If[Potential[x] >= threshold[n],
Force[x_, v_] := -k*x-G*v;];
のように、条件にあった時だけIf文内を実行したいのですが、If文の条件判定がうまく機能していないようです。
条件分に変数のある関数が入っていたらまずいのでしょうか。
下のように、条件文の前で関数の値を置き直せば問題ないのでしょうか。
こうすると動くようになったと思いますが、実際のものは煩雑でして正しく動作しているのか分かりません…。
よろしくお願い致します。
Potential[x_] = 0.5*k*x^2;
Threshold[n_] = 10*n;
Potential_value = Potential[x_];
Threshold_value = Threshold[n_];
If[Potential_value >= Threshold_value,
Force[x_, v_] := -k*x-G*v;]; >>782
>If[Potential[x] >= threshold[n],
>Force[x_, v_] := -k*x-G*v;];
よくわからないけど、If文の条件が実際に評価されるときに(x,nに)値が入ってるんですか?そうでないと関数Fは未定義のままですよね?
それとThresholdはMathematica自身の関数として定義済みでは?
それと、Potential_value =と_valueをつけると変数にセットできないんじゃありませんか?私の使ってるバージョンではできません。 x、nに値は入っています。
初心者でトンチンカンな事を聞いてるみたいですみません。
>それとThresholdはMathematica自身の関数として定義済みでは?
これ、そもそも○○[x_]などと定義するのが間違ってるって事でしょうか?
xに値が入るなら○○[x_]と定義する必要がない、と。
>それと、Potential_value =と_valueをつけると変数にセットできないんじゃありませんか?私の使ってるバージョンではできません
たしかにできないみたいです。
すみません。
x,nをループで回して、
Potential = 0.5*k*x^2;
Threshold = 10*n;
If[Potential >= threshold,
Force[v_] := -k*x-G*v;];
というようにすればいいのですかね。 ループの中で
Force[v_] := -k*x-G*v;
と何度も定義し直すのは何か意図があるのですか?
素直にループの外で
Force[v_]:=If[ 0.5 k x^2>=10 n, -k x-G v ,必要なら何か定義...]
と定義しておくのはダメなのですか? >>784
Mathematicaの入門書を読むことをお勧めします。 条件変えるループの中で定義したいんです
皆さんどうされてるんでしょう
本読み直します ヒントになるかわかりませんが、
関数の引数に関数を与えることもできます。 どこのがっこんだよ
くだらんものをだすな 資源の浪費だ >>792
あ、スマホでは出ないんだな。PCで見てくれ Clear[p,q];
Sum[If[q == p, q, 0], {q, 1, 3}] /. p -> 2
の結果として 2 になること期待したんだけれど,
qが返ってきてしまった…
なして? If[1 == p, q, 0] + If[2 == p, q, 0] + If[3 == p, q, 0] /. p -> 2
だから q だろうね。
Sum[If[q == p, q, 0] /. p -> 2, {q, 1, 3}]
だと2だよ。 If の属性に HoldRest が設定されているからだな。
Sum は Block を使って変数を置き換えていく。Block は変数を置き換えながら内側の式を評価するが、HoldRest のせいで第2引数の q は評価されずにそのまま残る。
第1引数は評価され置き換わるが、If は第1引数が True または False になるまで開かれないから、結果 >>798 の1行目のように展開され、答えは q となる。 0130
ふうL@Fu_L12345654321
学コン1傑いただきました!
とても嬉しいです!
https://pbs.twimg.com/media/D-IuUuqVUAALnAB.jpg
https://twitter.com/Fu_L12345654321/status/1144528199654633477
https://twitter.com/5chan_nel (5ch newer account) >>797
Sum[If[q == p, Evaluate@q, 0], {q, 1, 3}] /. p -> 2
でも2が返るよ。
昔lispやってたので馴染めないんだが、、、。 松坂和夫著『解析入門中』を読んでいます。
周囲の長さが一定 2*s である三角形の面積の最大値を求めよ。
ヘロンの公式より、
S = sqrt(s * (s - x) * (s - y) * (s - z))
φ(x, y, z) = x + y + z - 2*s
(grad φ)(x, y, z) = (1, 1, 1) ≠ (0, 0, 0)
0 < x < y + z = 2*s - x
0 < y < z + x = 2*s - y
0 < z < x + y = 2*s - z
だから、
0 < x < s
0 < y < s
0 < z < s
でなければならない。
φ(x, y, z) = 0
0 < x < s
0 < y < s
0 < z < s
という条件下で、
f(x, y, z) = (s - x) * (s - y) * (s - z)
を最大化する (x, y, z) を求めて、面積 S を計算すればよい。
↓は制約条件を満たす点の集合、いわゆる実行可能領域です。
https://imgur.com/XzhCaP7.jpg https://imgur.com/Fwn3FfO.jpg
↑実行可能領域の正三角形の重心で面積が最大になることが分かりますね。 f が具体的にどんな関数なのかは教えずに、 f が t の関数であることをMathematicaに教えるにはどうすればいいのでしょうか? iとjが添え字で、x[i,j]を総和記号で足し合わせるとき、 j≠iのときは足さないのはどうやって入力するの?
n m
Σ Σ x[i,j]
i=1 j=1
j≠i
みたいなやつ。意味わかります? >>810 レスを有り難うございます。
ただ、それはわかってますが、もうちょっとスマートに書きたいということです。 ラテン語またはギリシア語の「マテーマタ」(古希:Μαθηματα) は「学ばれるべき事ども(Lesson)」という意味である。 対角成分が発散してないなら引くかな
その方が早そうだし
発散するときは
x[i_, j_] := 0 /; i == j;
x[i_, j_] := X[i, j];
Sum[x[i, j], {i, 1, 3}, {j, 1, 3}]
Out=X[1, 2] + X[1, 3] + X[2, 1] + X[2, 3] + X[3, 1] + X[3, 2]
とか? >>813 お返事が遅くなってすみません。有り難うございます。それでうまくいきました。 Wolfram Languageの実践的な解説をしているサイトを教えてください。
Mathematicaの本を見ても基本的なことしか書いていないですよね。 合成関数について質問です。
f[x_, y_] := x^2 + y^3
g[t_] := {t^2, 3*t+1}
とします。
f[g[t]] とはできないですよね?
f[v_] := v[[1]]^2 + v[[2]]^3
g[t_] := {t^2, 3*t+1}
としないといけないんですか? 引数がx,yの2個の場合しか定義してないのに
1個しか与えてないのだから未定義扱いになるのは当然では? Mathematicaの年間ライセンスが三月十四日円周率キャンペーンで安くなってるけど買おうかな? 斎藤正彦著『斎藤正彦線型代数学』を読んでいます。
正規行列はユニタリ行列によって対角化できるという定理が書いてあります。
Wolfram Language 12で、 A = U*D*U^* となる U, D を計算する関数を作りました。
orthogonalizationU[A_] := Module[
{eigenvs, U, DIA},
eigenvs = Eigenvalues[A];
DIA = DiagonalMatrix[eigenvs];
eigenvs = DeleteDuplicates[Eigenvalues[A]];
U = Transpose @ ((Flatten[#, 1])& @ (Orthogonalize /@ (NullSpace /@ ((A - #*IdentityMatrix[Length[A]])& /@ eigenvs))));
{U, DIA}
] 斎藤正彦著『斎藤正彦線型代数学』を読んでいます。
以下の定理が書いてありますが、 A = U*T (T が上三角行列)の場合しか証明されていません。
2.5.17 【定理】
任意の n 次正則行列 A はユニタリ行列 U と、対角成分が正の実数であるような上三角行列(下三角行列でもいい) T の積 U*T (T*U でもいい)として
一意的に表わされる。 他のケースについて書いておきます:
(2)
A^T = U * T (T は上三角行列)
と分解し、両辺の転置を取ると、
A = T^T * U^T (T^T は下三角行列)
(3)
A の列たちをリバースした行列を A' とする。
A' = U * T (T は上三角行列)
と分解する。
T の列たちをリバースした行列を T' とする。
A = U * T' (T' は左上三角行列)
である。
U の列をリバースした行列を U' とする。
T' の行をリバースした行列を T'' とする。
A = U' * T'' (T'' は下三角行列)
である。
(4)
(3)の手順にしたがって、
A^T = U * T (T は下三角行列)
と分解する。
A = T^T * U^T (T^T は上三角行列)
(2)〜(4)のケースの一意性についても
A = U*T (T が上三角行列)のタイプの分解の一意性から自明ですよね。 Wolfram Language 12で、
4つのタイプの分解を行う関数を作りました。
gs[U_, a_] := Append[U, Simplify[normalize[a - Sum[(a.u)*u, {u, U}]]]];
gramschmidt[A_] := Module[{U = {}}, Do[U = gs[U, A[[i]]], {i, 1, Length[A]}]; U];
qrdecom[A_] := Module[{U, T},
U = gramschmidt[A];
T = Table[If[j > k, 0, A[[k]].U[[j]]], {j, 1, Length[A]}, {k, 1, Length[A]}];
{Transpose[U], T}];
qrdecom2[A_] := Module[
{B, U, T},
B = Reverse[A];
U = gramschmidt[B];
T = Table[If[j > k, 0, B[[k]].U[[j]]], {j, 1, Length[B]}, {k, 1, Length[B]}];
{Transpose[Reverse[U]], Reverse[Transpose[Reverse[Transpose[T]]]]}];
qrdecom3[A_] := Module[{B, U, T},
B = Transpose[A];
U = gramschmidt[B];
T = Table[If[j > k, 0, B[[k]].U[[j]]], {j, 1, Length[B]}, {k, 1, Length[B]}];
{Transpose[T], U}];
qrdecom4[A_] := Module[
{B, C, U, T},
B = Transpose[A];
C = Reverse[B];
U = gramschmidt[C];
T = Table[If[j > k, 0, C[[k]].U[[j]]], {j, 1, Length[C]}, {k, 1, Length[C]}];
{Transpose[Reverse[Transpose[Reverse[Transpose[T]]]]], Reverse[U]}]; >>820,821,822,823
こいつは数学の本スレで連投を繰り返している荒らしです
レスしないでください Wolfram Engine+Jupyter notebookしか使ったことないから、本物のMathematicaがどんな感じなのかわからん 正規MathematicaからWolfram Engine呼び出して使えば・・・・・ フロントエンドカーネルと並列処理のローカルカーネルとしては使える
並列処理のリモートカーネルとして使えないものか RpiのMathematicaフロントエンドを高性能のWindowsマシンにインストールしたEngineにリモート接続してフル機能使えるやん Jyupiter notebokk+Wolfram Engineより、Mathematica FE on RPi+Wolfram Engineのほうが遙かに幸せになれる ラズパイのMathematicaって仮想マシンでも動くの? Raspberry Pi Desktop X86を仮想化作業中
ここに入っていればRaspberry Piも要らないね
Mathematica 6のフロントエンドで動けば助かったけど、6.0 FE+12.1 Engineの組み合わせはちょっと無理っぽかった
Macで動いたって言うブログ見たけれど Engineをフロントエンドと並列計算用サブカーネルに使うことは簡単にできる
並列計算用サブカーネルにのみ使って、正規の8カーネル+Engine 8カーネルで16カーネル並列計算できないか
なんとなくやり方はわかったけど、これは後でやってみる やっぱx86版にはMathematicaはないみたいね 5つの8カーネル並列計算プロセスを1台にインストールされたMathematicaフロントエンドで制御して計算する フロントエンド落ちたら全部おじゃんになる
やりすぎたか これはMathematicaタダという理解でいいの? 単一のマシンでなくてもいいなら
無償のラズパイMathematicaフロントエンドか別に使えるバージョン新しめのMathematicaフロントエンドを持っているなら Jyupiter Notebookでもタダで使えてた
でもMathematicaのフロントエンドを使えば高額な有償Mathematicaと同じになっちゃう夢を見た フロントエンドにリモートMathKernelとの通信機能が実装されていれば
Wolfram EngineのMathKernelはMathematicaパッケージのMathKernelと全く同じで、特に機能制限されてない ラズパイってアームだよな
スマホでMathematicaとか未来だ ARM Raspbianをエミュレータで動かすことも可能なのか Raspberianは動いても、エミュ上ではMathematicaは動かないように細工されているらしい Jyupiter Notebookをフロントエンド、Wolfram Engineをバックエンドとして開発に使ってもよいとWolframが
許諾しているのなら、Jyupiter NotebookをMathematica Notebookに置き換えても構わないよな?
それを許諾しないとは書いてない 結局、ラズパイ買ってきてRasbianの中のMathematicaを色々設定して、高スペックサーバーにインストールしたWolfram Engineを呼び出すのが良いってこと? 良いか悪いかは言わない
ネットワーク経由のリモート計算でならロハでもできるというだけ
それが無理ならおとなしくローカルでJyupiter Notebook Jyupiter Notebookで機能限定がなくフルスペックの計算ができるのは確認できていた
ただJyupiter NotebookとMathematica Notebookの使用感が全く異なっており、ファイル互換性もない
なんとかMathematica NotebookをフロントエンドにしてNotebookを統一したいというのが夢のきっかけだった 12.0から12.1にしてみたが違いがわからん
何か変わったか? Raspberry Pi 4B 4GでMathematicaを使ってるけど普通に使える。ベンチマークだと、15年前のPowerMacと同じくらいの性能。少し不満なのはできた3Dグラフィックをマウスで動かすのが遅くてうまくいかないことと、入力補完機能がないことくらい。 それでうまくいくシステムもあろうが、Windowsのリモートカーネルは別の方法で簡単に動く Raspiのmathematicaで1秒を超えるような計算をしたことがほとんどない
3Dグラフィックが遅いのも入力支援機能がないのもフロントエンドの問題だから解決しない rshの代わりにWnRM使ってリモートカーネルを並列計算用のWorkerカーネルとして動かそうとしたけど今一歩うまくいかない 今のところ出来ているのは一つのNotebookフロントエンドから多数のマシンにそれそれ独立したマスターカーネルを立ち上げて、
同じマシンで8 x Workerカーネル動かして計算すること
マスターカーネルを一つにして、N個のマシンにWorkerカーネルを分散して、8N Workerカーネルで並列計算してパフォーマンス
を見てみたいのだが、うまくいってない つまり、N個のCPUがあるとき
N x 8 Wokerカーネルの並列計算は出来ているが
1 x 8N workerカーネルの並列計算が出来てない
今のところSurface Pro7も使ってN=6
Xeon E5-2643v3 (実6コア/仮想12スレッド) x2
Core i7-8700T(実6コア/仮想12スレッド)x2
Core i7-4790K(実4コア/仮想8スレッド)x1
Core i5-1035G4(実4コア/仮想8スレッド)x1
それぞれのコアで8 Workerカーネル動かして並列計算してる >>853
てことは初代iMacよりも性能上なのか
学生はStudent版買うよりこっちのほうがいいな なんかコールバックしてくるポートがおかしいんだよな
バグなのかな? Lightweight Grid利用だけど、3台の実マシン(12C/24T+6C/12T+6C/12T)で24カーネル並列計算
もっとコアが欲しい ローカルにあるWolfram Engineのカーネルをリモートカーネルのように起動して起動カーネル数増やすことも出来る Mathematicaの入っているPCでremoteカーネルにローカルのWolfram engineのカーネルを起動するようにハックすれば、
1台のPCで8カーネルの制約を超えた並列計算が可能 32カーネル並列計算
並列カーネルの状態表示が壮観 ローカルカーネルをリモートカーネル的に使うことはできるが、Windowsでlightweight gridなしにリモートカーネルを使うのは
難しい
引数がどうにもうまく渡せない 数学掲示板群 ttp://x0000.net/forum.aspx?id=1
学術の巨大掲示板群 - アルファ・ラボ ttp://x0000.net
数学 物理学 化学 生物学 天文学 地理地学
IT 電子 工学 言語学 国語 方言 など
PS 連続と離散を統一した!
ttp://x0000.net/topic.aspx?id=3709-0
微分幾何学入門
ttp://x0000.net/topic.aspx?id=3694-0 ようやく成功した
これでlightweight gridに頼らなくても済みそうだ 環境設定→並列処理→リモートカーネル→カスタム起動コマンドを使用
・ローカルカーネル x8
・ローカルにあるWolfram engineのカーネル x8
・2つの別マシンのカーネル x8x2
を起動して32カーネル クライアント・サーバーともWindows10
Windows10オプションのOpenSSH(クラ・サバ)を使用 要するにもうMathematicaはタダということなのかね? スキルあるIT技術者には弄るオモチャとして開放するよ
(フロントエンドや数学的な)サポートはしないけど
ってとこじゃない?。 >>875
完全タダは無理じゃない?
自分は正規ライセンスのMathematicaにWolfram Engineを組み合わせてできる範囲の並列分散処理をやってる
ラズパイ買えば、含まれる無償MathematicaとWolfram Engineの組み合わせでできることは拡大するが、ラズパイ
のコストはかかる
また、この方法はバックエンドの処理が重い場合には有効だけれども、複雑なグラフの描画とかフロントエンドの処理
が重い場合にはあまり意味がない フロントをラズパイで
Wolfram EngineをPCで
ってパターンならラズパイ分の値段だけでいける? >>878
いける
それにメリットを感じるかどうかは、やりたいことによる フロントエンド ― マスターカーネル ― スレーブカーネル
で、(スレーブカーネルは並列処理のときのみ)
フロントエンド をラズパイ上のMathematicaで、 マスターカーネル ― スレーブカーネル をPC上のWolfram Engineで
分散処理することができる
また、フロントエンド ― マスターカーネル をラズパイ上のMathematicaで、スレーブカーネルをPC上のWolfram Engineで
分散処理させることもできる
並列処理計算なんてしない人が多いんだろうけど
その場合はマスターカーネルをPCに任せるという話になる。マスターカーネルでの処理が重ければ分散処理の意義が出てくる
だろう >>873
これ、ParallelTableとかでシームレスに32並列行けるの? >>881
フロントエンド をラズパイ上のMathematicaで、 マスターカーネル―スレーブカーネル をPC上のWolfram Engineで分散処理は、
https://community.wolfram.com/groups/-/m/t/94917
https://support.wolfram.com/ja/12495?src=mathematica
http://www3.econ.fukuoka-u.ac.jp/usage/doku.php/windows
を基に設定した。これらを理解したうえで設定すればできた。
フロントエンド―マスターカーネル をラズパイ上のMathematicaで、スレーブカーネルをPC上のWolfram Engineで分散処理は、
Windows10のオプション機能であるOpenSSHを使えば実現できる >>883
ありがとうございます。できるかどうか分かりませんが、とりあえず読んでみようと思います。 そのリンクには、カーネルをリスンモードにしてフロントエンドからコネクトする方法と、フロントエンドをリスンモードにしてカーネルから
コネクトする方法が書かれているが、自分は前者で実現した。
PC側でカーネルをリスンモードで起動するバッチファイルを作成し、先にカーネルを起動する。フロントエンドのカーネル環境設定
でフロントエンドからPCに接続するようにMLOpenの引数を設定したうえで、マスターカーネルを起動する。
LinkNameは自分で適当にポートを2つ決め打ちして両方で合わせればいい。
やろうと思えば、PC上のカーネルをOpenSSHなどでリモート起動することもできるはず。 ラズパイを仮想マシンで動かせばPCだけで完結できる? ARMのエミュではMathematicaは起動しないように制限かけているらしい
残念 このコードだけ見て何をしているか分かりますか?
解説お願いします
seedList =
With[{vertices = .1 N[PolyhedronData["Icosahedron"][[1, 1]]]},
Join[Map[{#, 2} &, vertices],
Map[{# + {1, 1, 1}, -2} &, vertices]]]; Wolfram Engineの最新版についての案内のメールが来たけど、アップデートしたほうがいいの? Wolfram Engineを最新版にアップデートしたら、Jupyterから使えなくなってしまいました。kernel errorが出ます。どうすればいいですか? anacondaをアンインストールして再インストールしてもだめでした。
jupyter notebookを起動してNEWを押すとWolframがなぜか選択できるようになっています。
どうすればいいですか? jupyter notebookを起動してNEWを押したときにWolfram Engine 12.0とWolfram Engine 12.1の両方が現れます。古い方のWolfram Engine 12.0の表示を消すにはどうすればいいですか? aを正定数として
x^2 + y^2 <= a^2の範囲にある円筒面x^2 + z^2 =a^2を表示したい
Plot3Dで試してみたが(↓)違う図が表示される
なんて書けばいい?
a=1
Plot3D[x^2 + z^2 == a^2, {x, -a, a}, {z, -a, a},
RegionFunction -> Function[{x, y, z}, x^2 + y^2 <= a^2],
BoxRatios -> Automatic] 参考
a=1
RegionPlot3D[x^2 + y^2 < a^2, {x, -2, 2}, {y, -2, 2}, {z, -2, 2}] こんにちは
例えば
FinancialData["MITSUBISHI"]
とやると,三菱商事の株価(現在値)が出ますが
企業名と,上のようなMathematica上のシンボルの紐付け一覧はどこを見ればわかるでしょうか? Manipulate関数について
↓こうするとプロットされるのに
Manipulate[ Plot[ Sin[x-a], {x,0, 2π} ], {a,0,2π} ]
↓こうすると何もプロットされない理由を教えてください。
foo[x_, a_] := Sin[x-a];
Manipulate[ Plot[ foo[x,a], {x,0, 2π} ], {a,0,2π} ]
Manipulate の内側に全てを書くのは厳しい時があるので、解決策もお願いします。 foo[x_,a_] := Sin[x-a];
n = 10;
With[ {tbl = Table[ Plot[ foo[x, 2π(k-1)/n], {x,0, 2π} ], {k,n} ]},
Manipulate[ tbl[[k]], {k,1,n,1} ]]
回りくどいが Withで見かけ上の局所変数を導入するとプロットできるようになった。
これが ModuleやBlockだと上手くいかない。どういう理屈なのかは分からない。 >>904
12.0でやってみたが、どっちも上手く行くけどなあ。 >>906 試してくれてありがとう。
こちらの環境は
Wolfram Language 12.1.1 Engine for Mac OS X x86 (64-bit)
もしかしたら「正しい挙動」にアップデート修正されてしまったのかも知れない。
これまでのコードが動かなくて困る人出てきそうだけどいいのかな。 Wolfram Language 12.1.1 Engine
フリーのやつか。製品版だと問題ないとか? そうフリーのやつです。 でもコアの部分は製品版と共通のはずです。
わざわざそんなコード解釈部分をいじって差をつけたりしないんじゃないかと。 フロントエンドの問題じゃないの?
JupyterだとしたらManipulateでのトラブルはいくつか見た > JupyterだとしたらManipulateでのトラブルはいくつか見た
これは出力を 〜.nb ファイルにExportして Wolfram Playerで閲覧しないとマウスでグリグリ動かせない (Plot3Dとかも同様)
それって面倒じゃん?そういう類の話ではないかな。 それについての不満は無いです。
そして wolframscript のREPLで試しても同じなので、Jupyterの問題ではなさそう。
複合要因の「おま環」(問題が起きてるのは、おまえの環境だけだ) の可能性もあるので、深入りしても時間の無駄かもしれません。 JupyterでWolframEngineを使っていますが,グラフの色を例えば,Redと指定しても真っ赤な色にはならず,ピンクのような色になってしまいます.どうすれば正しく色を表示できますか? Jupyterは知らないが
ふつうはRGBとかHueとかで指定する >>916
ありがとうございます.
RGBでRを1,Gを0,Bを0と指定しても,ピンクのような色になってしまいます.
Jupyter Notebookがおかしいのではないかと思っております. Raspi4のUbuntu64でadd architecture armhfしてラズパイ用のmathematicaが実行できた。
ただし、keygen必要。
Out[2]//InputForm=
{"MachineName" -> "ubuntu", "System" -> "Linux ARM (32-bit)",
"BenchmarkName" -> "WolframMark", "FullVersionNumber" -> "12.1.1",
"Date" -> "December 24, 2020", "BenchmarkResult" -> 0.206,
"TotalTime" -> 67.311, "Results" -> {{"Data Fitting", 2.757},
{"Digits of Pi", 1.616}, {"Discrete Fourier Transform", 10.401},
{"Eigenvalues of a Matrix", 3.524}, {"Elementary Functions", 5.119},
{"Gamma Function", 2.196}, {"Large Integer Multiplication", 2.28},
{"Matrix Arithmetic", 1.637}, {"Matrix Multiplication", 10.93},
{"Matrix Transpose", 4.91}, {"Numerical Integration", 2.502},
{"Polynomial Expansion", 0.339}, {"Random Number Sort", 1.19},
{"Singular Value Decomposition", 8.676}, {"Solving a Linear System", 9.234}}} と煽ってみたが、qemu使えばM1 MacでUbuntu Arm 64bit 動くのか。
https://qiita.com/bellx2/items/c06e1c00ad9de77a4fb9
M1 MacでMathematica使えるじゃまいかw Plot3DでMathematica5あたりのカラーリングで表示させるにはどうしたらいいの?Mathematica9とかグラフが黄色みたくてキモイ昔は赤青色味でキレイだったのに >>921
無理だろ。>>918もラズパイ上でやってるし X86のMathematicaをM1 Macで走らせたベンチは出てるな。
https://community.wolfram.com/groups/-/m/t/2112731#
2.59だから上のラズパイ4より10倍速いw コレなんかも見るとMacとラズパイでUbuntu Arm動くならアプリの互換性ありそうなんだがなぁ
https://qiita.com/kose3/items/af9edc9c40c9ae8fc5c3 MacのM1で動かしてる人いる?
速度どのくらい? 簡単な質点の運動の計算をしてるんだけど、
加わる力の定義をSin[]やExp[]などの入った式でして、
それを1000点ごとに
list0 = Table[0, {t, 0, last, sample*dt}];
というような形でリストにして
Export["***.dat", list0, "Table"]でデータとして書き出してるんだけど、
数式のまま吐き出されることがあります。
数値で吐き出させるためにはどうするんでしょうか?
いろいろ試すと、力の定義の中の数値の桁数が多いと、数式で吐き出される気がします。 Export["***.dat", N[list0], "Table"]
ってやってもそうなるの?
とりあえずExportの直前でlist0を表示させて,
どうなってるのかチェックしたらええんじゃないだろうか.
有理数になってる場合なんかは分数のまま出力されたりする. ありがとうございます。
試してみます。
煩雑な分数のまま書き出される場合があります。
で、グラフ描画ソフトで読み込めなくて困っていました。 >>929
N[list0]で解決しました。ありがとうございました。
基礎すっ飛ばして(昔やったかもしれませんが忘れた)実際に使われてるプログラムを変形して使ってると困ったことになりますね。 創造の館のおかげでスピーカーの能率は高いほうがいいというのがわかってよかった。 某スレの問題を解こうとして書いた式だけど、これってForとか制御構造を使わないと書けないのかな?
For[n=2,n<=30240,n++,d=Divisors[n];l=Length[d];For[i=1,i<=l-12,i++,e=Extract[d,i]; s=Select[e<=# && #<=2e &][d]; If[Mod[n,2e]==0 && Length[s]>=13,Print[n,s],0]]]
期待する処理:2以上30240以下のnについて、nの約数eと2eの間に(eと2eを含めて)nの約数が13個以上あるものを見つけてプリントする バレるか否かはともかくとして、Mathematicaでの計算結果を論文とかに発表したり、計算結果をプログラムとかに組み込んだりするのに必要なライセンスで、一番安いやつってどれになるんでしょうか? 大学のライセンスがせっかく使えるので何かしたいのですが検算ぐらいにしか使えてません ffmpegインストールし、システム環境変数Pathに実行ファイルのパスを登録しているのに
VideoGeneratorやVideoTrimなどを使うと
制限付きのFFmpegを使用しています.より完全なコーデックサポートのためにはFFmpegをインストールしてください.
になるのですが、なぜだと思われますか???
PowerShellでffmpegで実行されるのでパスは通っています
OSはWin7です iphoneでwolfram player動くみたいだけどこれ編集ができないみたい
iPhoneでノートブック編集する方法ってある? ここで書くのもあれなんですがMapleとどっちがいいんですか?
MathematicaにできてMapleにできないあるいはその逆とか Mapleのほうが普通のプログラミング言語っぽいですよね。 Mapleは手続き型
Mathematicaは関数型
ですよね。
Mapleのほうが使いやすいと思います。 RemoteObjectという分散処理メソッドが取り入れられたようだが、使い方がわからん Plot3D[x^2 y^3 + (x - 1)^2 y, {x, -2, 2}, {y, -2, 2},
PlotRange -> 1.5, BoxRatios -> 1]
といていたら、講師【女性)ににらみつけられ、助教【男性)に教室から出されました。
理由がよくわかりません。
教えてください。
悩める童貞大学生 PlayerとEngineの違いって何かわかる人いる?
WolframScriptの有無はわかるんだけど、Jupyter-notebookとか自前のフロントエンドさえ用意すれば
同じに使えるのかな?
playerにはWolframScriptがないため定期的なライセンス認証アップデートも不要なのかな?
ご存じの方いれば教えてください どうやら>>948の推測は正しいらしい
1台のPCでMathematica/Wolfram Ensine/Wolfram Playerの8x3=24カーネル起動もできたわ >>948
はい、可能です。その方法を何故かWolframの公式ツイアカが紹介していたこともあったw >>950
なんだかPlayerとEngineの切り分けが難しいね
カーネルにMathematica Notebookを追加したもの=Matthematica
カーネルにPlayer(ライセンス認証不要)とWoframScript(ライセンス認証必要)を追加したもの=Wolfram Engine
カーネルにPlayer(ライセンス認証不要)を追加したもの=Wolfram Player
という感じなんだろうか?
カーネルとPlayerは商用以外無償でライセンス認証不要
開発環境であるMathematica Notebook/Wolfram Scriptはライセンス認証必要で
Mathematica Notebookは有償Wolfram Scriptは無償でライセンス認証不要
って感じかな?
開発環境を自前で用意すれば、ライセンス認証も不要なPlayerがお勧めなんだろうか
ファイルサイズを比較すると、PlayerはフルスペックのNotebookに機能制限かけたものっぽいね 訂正
カーネルとPlayerは商用以外無償でライセンス認証不要
開発環境であるMathematica Notebook/Wolfram Scriptはライセンス認証必要で
Mathematica Notebookは有償でWolfram Scriptは商用以外無償 価値のほとんどは本来カーネルにあるのだろうけど、開発環境側に転嫁してコスト回収しているというのが最初は驚きだった Playerはiphoneで動くということはiphoneでmathematica動くようにできるのか? >>953
開発環境とか統一性とか、Jypterと比べるとよくできていると思う。 ハックして76スレッド並列計算しているよ
6PC+2EWS Playerは定期的な再アクチが必要ないから便利に使える VMWareのはそれ単体で一応使えるけど、Wolfram Playerだけで実用的な計算は難しくね?裏技はあるにしても。 「同じ」という言葉から受ける印象のあいまいさの範囲の問題
カーネル部分は無償、インターフェースを有償にしている部分を同じと表現した PlayerハックすればMathematicaの機能全部使えそうだな 初心者向けのガイドブックを
読んでるのですが、ほとんどの例が
notebookで1行ごと実行する形式になってます。
例えば
ある数値データのファイルを読み取って
計算した結果をファイルに出力する
例えば結果をレポート形式で出力するような
プログラムを作りたい場合、
どう言うに流れになるのでしょうか?
毎回ノートブックを開いて行ごと順番に
実行するしかないのでしょうか? >>965
なんでよ
一つのセルにセミコロン;で区切ってコマンド書けば、そのセルを実行すれば順次実行になるでしょ 『線形代数とMathematica』という本に以下の関数があります。
squareMatrixQ[a_?MatrixQ] := Equal @@ Dimensions[a]
この関数の引数をa_?MatrixQと制限しています。
squareMatrixQというのは引数が正方行列かどうかをチェックするための関数です。
ですから、数の1を入力すればFalseを返さないといけないと思います。
ところが、この関数は引数を制限しているため、数の1を入れて評価することが出来ません。
***Qという関数をこのように書くのは良くないと言えるでしょうか? squareMatrixQはMatrixQがFalseの時の値を持っていないから、
Falseを返したいなら、書き方が少し違いますが
squareMatrixQ[a_] := False /; !MatrixQ[a]
squareMatrixQ[a_] := Equal @@ Dimensions[a] /;MatrixQ[a] >>970
やはり、 **Qという関数だから、どんな引数が来てもTrueかFalseか判定できないと駄目ですよね。 Wolfram Programming Labが廃止ってマ? SquareMatrixQ[1] False
SquareMatrixQ[{{1}}] True
でなにかもんだいあり? >>94
Plot3D[x^2 y^(9) + (x - 1/2)^2 y, {x, -2, 2}, {y, -2, 2},
PlotRange -> 1.5, BoxRatios -> 1]
をみせなさい。
上品にしかしちゃんとね SquareMatrixQ[m]
gives True if m is a square matrix, and False otherwise. G を群とする。
H1 と H2 は G の正規部分群で、互いに同型であるとする。
G/H1 と G/H2 は常に同型であるか? G を群とする。
H1 と H2 は G の正規部分群で、 G/H1 と G/H2 は互いに同型であるとする。
H1 と H2 は常に同型であるか? .ィ:.l ヽ: : : : : : : : l ヽ: : : : : : : : \
/: : :l ヽ.: : : : : :l ヽ: : : : : : : : : ヽ
,/: : : ヘ l: : : : : l ヽ: : : : : : : : : :ヽ
//: : : :./ヽ l: : : : : lヽ , l : : : : : : : : : : ヽ
/:/: : : : :' l: : : : : :l li: : : : : : :i: : : : :ヘ
/:/: : : : :,' 、: : : : :.l li: : : : : : :l : : : : :.l
i: i: : : : : l 、: : : : :l li: : : : : : :l: : : : : :l
l: l: : : : :.l ヽ - 、 ヽ: : : :ヘ __',ヘ: : : .l: : l : : : : :.l
l: l : : : : l二‐- ` ヽ: : : ヘ´__ヘ l: : : l: : l : : : : :.l
∨: : : : :l┬┰┬ ヽ: : .ヘ ┬┰┬l: : : l: : l : : : : :.l
∧ : : :.Cゝ__.ノ ヽ: : :ヽゝ__.ノl: : : l: : l : : : : :.l
/: : l: : : ヘ ヽ: : \_ ィ: : : l: : l: : : : : :l
./: : : l: : : : :、 ' ` ー―l: : : l: : l: : : : : :l
/: : : :ハ: : : : iヽ l: : : l: : l: : : : : :lKGB はあれすね
/: : : : /:ヘ: : : : ヘ ヽ ー - /l: : : l: : lヽ: : : : l
,/: : : : /, ´ヽ: : : :.ヘ |` 、 ィ彡'イl: : : l: : l iヽ: l
/: : : : // 、: : : :∨i 、 ヽー≦三三ヲ: : : ハ: : l .li `<
': : : : :/: : :l 、: : : i :l ヽ、 `ー‐''"´ /: :/l:ヘ: : l li `ヽ
: : : /: : : : : :l ;: : : :il 、ヽ、 /: :/ l ヘ: :.l li ヽ 新しいバージョンにすると、
(OSを新しいのにしたせいもあるかもしれません)
πがΠって表示されるようになったんだけど、
数式を表現できて、かつ、πってちゃんと表示できるフォントってありますか? グラフィック関連を系統的に覚えるのにいい本やサイトを教えてください. 例題で学ぶMathematica(グラフィックス編)
白石修二 高え~
matlabでさえ個人用は8000円にご時世に。
しかも文法がめっちゃ独特、prologかよ。
グラフを思い通りに書くのにも一苦労 >>986
サンクス
結構、洋書のmathematica本何冊も買って勉強したけど、どう見ても同じ事するならMapleの方が使い易そうだわ。
Mapleの方が見た目がダサいけど、安いしいいね。存在を忘れてた。 >>1
マセマティカっていつからコピーできなくなってたのかな
シリアルだけじゃ動かない >>977
駄目に決まってろ
Z2xZ4でZ2への全射2種類考えるだけていうか
すぐにでもそれを思わないのは何故? Mathematicaは文法が独自だから、
そんなものによりかかって、仕事をしてしまうとベンダーロックインされるだけだぞ。
実際、目玉が飛び出る程高い値段だろ。
なるべくソースコードも公開されているフリーなソフトを公共財として使い、
開発者をなんらかの形でサポート(たとえば金銭の寄贈)をして開発を継続
してもらい、しだいにソフトを向上してもらう方が長期的には良いんだよ。 ラズパイ買えラズパイ。
まあ、ラズパイも価格が暴騰して、ある程度Mathematicaが動く8gb版は2万以上するけど。 Mathematica8のstudent版買った時は既にオンライン認証方式だった.
少なくとも5.2はオフラインでいけた希ガス 値上がりしてたんだな
Mathematica13 ホームエディション 51,040 円 >Mathematica13 ホームエディション 51,040 円
それって大学の任期付きの助手でもその値段で買える? 教職員だと買い取り価格が30万円ぐらいしなかったか?
趣味用だと、仕事(たとえば給料を貰ってする研究や授業)には使えないはずだよ。 このスレッドは1000を超えました。
新しいスレッドを立ててください。
life time: 2772日 19時間 56分 3秒 5ちゃんねるの運営はプレミアム会員の皆さまに支えられています。
運営にご協力お願いいたします。
───────────────────
《プレミアム会員の主な特典》
★ 5ちゃんねる専用ブラウザからの広告除去
★ 5ちゃんねるの過去ログを取得
★ 書き込み規制の緩和
───────────────────
会員登録には個人情報は一切必要ありません。
月300円から匿名でご購入いただけます。
▼ プレミアム会員登録はこちら ▼
https://premium.5ch.net/
▼ 浪人ログインはこちら ▼
https://login.5ch.net/login.php レス数が1000を超えています。これ以上書き込みはできません。