- 1 名前：１３２人目の素数さん mailto:sage [2015/02/27(金) 10:38:04.85 ID:EbhLAWFK.net]
- _....._{{ 〃
, - ' ,..､､.ヾ{{フ'⌒`ヽ､
／ ,:', -‐‐` ´ '´⌒ヽ ヾ:､
. ,' ,'´ ,ィ ,ｨ ,' , ｀ヽ', ',-<
,' .i /|. /.| { i, i, }. }_,,）)
! | ! .,'-.{ ! !|; |`､.}ﾞ!.! |. ! ヽ.
', ', |Vｧ=､ﾞ､ ｀ﾞ､!-_:ﾄ,ﾘ', l ! | ﾞ',
ヽ､', l:!Kノ}. f:_.)iﾞi: ﾘ ! l ル
| l!iヽﾞ- ' , .!__:ﾉ ﾞ ,ﾘ l ﾘ'´
. ',|!!､ r‐┐ ｀ ノ'. /,ｲ ／￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣
'i!ﾞ､ヽ､ ﾞｰ' _, ィ,:',:''´ ＜ Mathematica に関する話題はここに書いてね！
ﾞ:､ｨ､jヾｰ::: 'iﾍ .ノ',ﾘ. ＼＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿
,､- '´ ヽ､ﾞ､ { ｀＞"､
/＼＼ ', } //｀ヽ
過去スレ
〓Mathematica〓
science2.2ch.net/test/read.cgi/math/985023298/
〓 Mathematica 2 〓
science3.2ch.net/test/read.cgi/math/1078534285/
〓 Mathematica 3 〓
science6.2ch.net/test/read.cgi/math/1121413040/
〓 Mathematica 四 〓
kamome.2ch.net/test/read.cgi/math/1197270001/ (鯖飛びでログ消滅)
mimizun.com/log/2ch/math/1197270001/
〓 Mathematica 5 〓
kamome.2ch.net/test/read.cgi/math/1285859504/
〓 Mathematica 伍 〓 (実質6)
wc2014.2ch.net/test/read.cgi/math/1320969748/
- 754 名前：￥ mailto:sage [2018/04/06(金) 01:05:11.63 ID:I+Mybrk/.net]
- ￥
- 755 名前：￥ mailto:sage [2018/04/06(金) 01:05:45.11 ID:I+Mybrk/.net]
- ￥
- 756 名前：￥ mailto:sage [2018/04/06(金) 01:06:02.43 ID:I+Mybrk/.net]
- ￥
- 757 名前：￥ mailto:sage [2018/04/06(金) 01:06:32.25 ID:I+Mybrk/.net]
- ￥
- 758 名前：￥ mailto:sage [2018/04/06(金) 01:06:48.35 ID:I+Mybrk/.net]
- ￥
- 759 名前：￥ mailto:sage [2018/04/06(金) 01:07:24.79 ID:I+Mybrk/.net]
- ￥
- 760 名前：１３２人目の素数さん [2018/06/08(金) 08:18:02.56 ID:oMit71Pl.net]
- 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]
]
- 761 名前：１３２人目の素数さん [2018/06/08(金) 08:18:58.85 ID:oMit71Pl.net]
- >>748
g1とg2を同じPlotRangeで描きたいと思い、↑のコードを書きましたが、
エラーが出てしまいます。
どうすればいいでしょうか？
- 762 名前：１３２人目の素数さん mailto:sage [2018/06/09(土) 03:34:04.34 ID:JHZIvIjE.net]
- >>749
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}];
で動作確認したら
- 763 名前：１３２人目の素数さん [2018/09/21(金) 17:37:02.38 ID:0uIdegM1.net]
- 例えば、
Subscript[x, 1]
を変数として使いたいのですが、
別のところで、
x = 3
などとやるとエラーが出ます。
どうすればいいでしょうか？
- 764 名前：１３２人目の素数さん mailto:sage [2018/09/21(金) 19:33:31.11 ID:3x53asMk.net]
- 私の環境では問題なく使えますが、、。
例えば
Subscript[x, 1]=123
とすれば
Subscript[x, 1]
は123を返してきます。
しかし、ただのxで参照したら別の変数でしょう？
- 765 名前：１３２人目の素数さん mailto:sage [2018/09/21(金) 19:38:24.31 ID:3x53asMk.net]
- あ、
f[a]=321
とすれば
f[a]
は321を返すんだから説明になってなかったですね。
でも
x_1(ここでは表現が難しい、、）は123を返してきます。
- 766 名前：１３２人目の素数さん [2018/10/15(月) 17:26:57.09 ID:rkpLwJKB.net]
- なぜ以下の関数は引数のリストの要素の和を計算するのでしょうか？
sumList[lis_]:=First[lis//.{x_, y___} -> x + {y}]
- 767 名前：１３２人目の素数さん mailto:sage [2018/10/17(水) 09:40:44.17 ID:q7AoWChQ.net]
- >>754
sumListの右辺で最初にx+{y}が計算されて{x+y}となるみたいだ。
yはSequenceが入る。
例えば
x+Sequence[a,b,c]は
x+a+b+c
となるようだからlistの要素の総和になる。
- 768 名前：１３２人目の素数さん [2018/10/18(木) 19:08:49.50 ID:rmFVsQv+.net]
- 3 + {Sequence[1, 2]}
を計算すると、
{4, 5}
になりました。
- 769 名前：１３２人目の素数さん [2018/10/18(木) 19:09:39.48 ID:rmFVsQv+.net]
- >>755
>sumListの右辺で最初にx+{y}が計算されて{x+y}となるみたいだ。
これはなぜなのでしょうか？
- 770 名前：１３２人目の素数さん mailto:禁止 [2018/10/18(木) 22:12:06.14 ID:98PMbrHl.net]
- {x,y,z,w}//.{x_, y___} -> x + {y}]
.....{x+y+{z,w}}---> x+y+z+w
- 771 名前：１３２人目の素数さん mailto:sage [2018/10/19(金) 08:00:36.08 ID:u3e62vYe.net]
- >>757
https://reference.wolfram.com/language/ref/Rule.html
によればRule“->”の右辺は即事に評価すると書いてあるから、
x+{y}はx,yに代入される前に評価されて{x+y}になり
そのあと代入されて、計算が続くということだろう。
- 772 名前：１３２人目の素数さん [2018/10/23(火) 03:53:40.31 ID:c7NxGayY.net]
- Poor man's Mathematica used on Raspberry Pi
- 773 名前：１３２人目の素数さん mailto:sage [2018/10/28(日) 00:03:19.87 ID:XtA7LFXv.net]
- FindEquationalProof
いつの間にか定理の自動証明ができるレベルになってたんですね
群の証明なんかが自動的にできてて驚いた
- 774 名前：１３２人目の素数さん [2018/11/02(金) 22:27:32.97 ID:WKAh5t24.net]
- 金持ちの数学者＞＞＞＞＞＞＞＞＞。。。。。＞＞＞貧乏人の数学者＞＞＞＞＞２ｃｈの数学者
- 775 名前：１３２人目の素数さん mailto:sage [2018/12/03(月) 00:58:26.74 ID:v688weAA.net]
- なんか久しぶりに数学勉強したくなってきた。
高校1年生の2次関数勉強中。サポートよろ！
- 776 名前：１３２人目の素数さん mailto:sage [2018/12/03(月) 18:27:35.82 ID:1vfhx/Vq.net]
- 何故それをここに書く
- 777 名前：１３２人目の素数さん [2019/02/05(火) 22:45:31.87 ID:zhX1PF9S.net]
- ｈttp://tvcap.info/2019/2/5/n2i9190205-2243560390.jpg
ｈttp://tvcap.info/2019/2/5/n2i9190205-2244010380.jpg
ｈttp://tvcap.info/2019/2/5/n2i9190205-2244040827.jpg
ｈttp://tvcap.info/2019/2/5/n2i9190205-2244240232.jpg
ｈttp://tvcap.info/2019/2/5/n2i9190205-2244110703.jpg
ｈttp://tvcap.info/2019/2/5/n2i9190205-2244390628.jpg
ｈttp://tvcap.info/2019/2/5/n2i9190205-2244540005.jpg
ｈttp://tvcap.info/2019/2/5/n2i9190205-2245000222.jpg
ｈttp://tvcap.info/2019/2/5/n2i9190205-2245040254.jpg
ｈttp://tvcap.info/2019/2/5/n2i9190205-2245180265.jpg
ｈttp://tvcap.info/2019/2/5/n2i9190205-2245250755.jpg
ｈttp://tvcap.info/2019/2/5/n2i9190205-2245300342.jpg
- 778 名前：低学歴脱糞老女・清水婆婆の連絡先：葛飾区青戸６－２３－１９ [2019/03/03(日) 11:46:45.53 ID:KV/cokeJ.net]
- 【超悪質！盗聴盗撮・つきまとい嫌がらせ犯罪者の実名と住所を公開】
①井口・千明（東京都葛飾区青戸６－２３－１６）
※盗聴盗撮・嫌がらせつきまとい犯罪者のリーダー的存在／犯罪組織の一員で様々な犯罪行為に手を染めている
低学歴で醜いほどの学歴コンプレックスの塊／超変態で食糞愛好家である／醜悪で不気味な顔つきが特徴的である
②宇野壽倫（東京都葛飾区青戸６－２３－２１ハイツニュー青戸２０２）
※色黒で醜く太っている醜悪黒豚宇野壽倫／低学歴で人間性が醜いだけでなく今後の人生でもう二度と女とセックスをすることができないほど容姿が醜悪である
③色川高志（東京都葛飾区青戸６－２３－２１ハイツニュー青戸１０３）
※色川高志はyoutubeの視聴回数を勝手に短時間に何百何千時には何万回と増やしたり高評価・低評価の数字を一人でいくつも増やしたり減らしたりなどの
youtubeの正常な運営を脅かし信頼性を損なわせるような犯罪的業務妨害行為を行っています
※色川高志は現在、生活保護を不正に受給している犯罪者です／どんどん警察や役所に通報・密告してやってください
【通報先】
◎葛飾区福祉事務所（西生活課）
〒１２４－８５５５
東京都葛飾区立石５－１３－１
℡０３－３６９５－１１１１
④清水（東京都葛飾区青戸６－２３－１９）
※低学歴脱糞老女：清水婆婆 ☆☆低学歴脱糞老女・清水婆婆は高学歴家系を一方的に憎悪している☆☆
清水婆婆はコンプレックスの塊でとにかく底意地が悪い／醜悪な形相で嫌がらせを楽しんでいるまさに悪魔のような老婆である
⑤高添・沼田（東京都葛飾区青戸６－２６－６）
※犯罪首謀者井口・千明の子分／いつも逆らえずに言いなりになっている金魚のフン／親子孫一族そろって低能
⑥高橋（東京都葛飾区青戸６－２３－２３）
※高橋母は夫婦の夜の営み亀甲縛り食い込み緊縛プレイの最中に高橋親父にどさくさに紛れて首を絞められて殺されそうになったことがある
⑦長木義明（東京都葛飾区青戸６－２３－２０） ※日曜日になると風俗店に行っている
- 779 名前：１３２人目の素数さん mailto:sage [2019/03/29(金) 23:09:53.69 ID:fHKb9Xcp.net]
- ver12マダー
- 780 名前：１３２人目の素数さん mailto:sage [2019/04/17(水) 06:08:16.62 ID:o2zp+n+z.net]
- 出たね。英語版は。
- 781 名前：１３２人目の素数さん mailto:sage [2019/04/20(土) 23:47:30.76 ID:ucNxwlZc.net]
- ユークリッド幾何の自動定理証明っぽい機能があるな
- 782 名前：１３２人目の素数さん [2019/04/23(火) 09:51:04.81 ID:6y6+1NWL.net]
- Mathematicaは手続き型ではなく関数型でコードを書いたほうが高速化するらしいけど、並列計算させたいとき
関数型で書けばうまくいきますかね？
- 783 名前：１３２人目の素数さん mailto:sage [2019/04/24(水) 13:32:48.98 ID:b+zqCTiq.net]
- そうだよ。関数型で書いて最後に //Parallelize を付けるだけで分散処理してくれる。
- 784 名前：１３２人目の素数さん [2019/04/25(木) 11:28:57.48 ID:LeBllpmo.net]
- >>771
サンキュー
頑張ってみる
- 785 名前：１３２人目の素数さん [2019/04/25(木) 13:33:39.52 ID:LeBllpmo.net]
- 最後に //Parallelize を付けるだけ
ってのはMathematicaじゃないような気もするけど
- 786 名前：１３２人目の素数さん [2019/04/25(木) 14:35:39.01 ID:9SdOWbtw.net]
- Gurdがいるんじゃないの？
- 787 名前：１３２人目の素数さん [2019/04/25(木) 17:17:29.74 ID:LeBllpmo.net]
- gridのこと？
自マシンのマルチコアプロセッサー利用なら何もいらない
正規ライセンスあればLightweight Gridも利用できる
- 788 名前：１３２人目の素数さん mailto:sage [2019/05/25(土) 01:07:49.66 ID:pma/hkB5.net]
- 無料化キター！
https://github.com/WolframResearch/WolframLanguageForJupyter
でnotebook代替もできるぞ！
- 789 名前：１３２人目の素数さん mailto:sage [2019/05/25(土) 01:52:51.18 ID:7q7gZ5Hj.net]
- https://twitter.com/wolframjapan/status/1131648749225369600?s=21
(deleted an unsolicited ad)
- 790 名前：１３２人目の素数さん mailto:sage [2019/05/26(日) 16:28:52.98 ID:TO+h+abu.net]
- https://nbviewer.jupyter.org/github/genkuroki/msfd28/blob/master/Free%20Wolfram%20Engine.ipynb
- 791 名前：１３２人目の素数さん mailto:sage [2019/06/07(金) 00:14:43.89 ID:KZdus4b3.net]
- まじかよ
- 792 名前：１３２人目の素数さん [2019/06/20(木) 15:23:01.20 ID:4cvo3/Hk.net]
- 数式もwysWygできる？
- 793 名前：１３２人目の素数さん mailto:sage [2019/06/20(木) 17:31:22.83 ID:BTA5nRWW.net]
- できるよ
TeXアウトプットもできる
- 794 名前：１３２人目の素数さん mailto:sage [2019/06/20(木) 22:22:20.22 ID:2OU9iC7f.net]
- これはマセマティカがタダになったという理解でいいの？
- 795 名前：１３２人目の素数さん [2019/06/21(金) 16:56:40.28 ID:p8Ey+i/x.net]
- >>781
入力はキャラクタベースのTexでやるってこと？
- 796 名前：１３２人目の素数さん mailto:sage [2019/06/25(火) 10:50:36.14 ID:Pn3mqXKf.net]
- >>442
出尽くしってことだろうね
- 797 名前：１３２人目の素数さん [2019/06/29(土) 16:32:34.25 ID:DHiuKlHq.net]
- 〓 Mathematica 七 〓
ふうL@Fu_L12345654
- 798 名前：321
学コン1傑いただきました！
とても嬉しいです！
https://pbs.twimg.com/media/D-IuUuqVUAALnAB.jpg
https://twitter.com/Fu_L12345654321/status/1144528199654633477
(deleted an unsolicited ad) [] - [ここ壊れてます]
- 799 名前：１３２人目の素数さん mailto:sage [2019/07/01(月) 16:09:26.13 ID:zLqAoY+7.net]
- 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;];
- 800 名前：１３２人目の素数さん mailto:sage [2019/07/01(月) 19:11:54.24 ID:5cPmcjg8.net]
- >>786
>If[Potential[x] >= threshold[n],
>Force[x_, v_] := -k*x-G*v;];
よくわからないけど、If文の条件が実際に評価されるときに(x,nに）値が入ってるんですか？そうでないと関数Fは未定義のままですよね？
それとThresholdはMathematica自身の関数として定義済みでは？
それと、Potential_value =と_valueをつけると変数にセットできないんじゃありませんか？私の使ってるバージョンではできません。
- 801 名前：１３２人目の素数さん mailto:sage [2019/07/01(月) 22:59:03.83 ID:wgvgoqOx.net]
- 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;];
というようにすればいいのですかね。
- 802 名前：１３２人目の素数さん mailto:sage [2019/07/02(火) 08:20:06.64 ID:AbMvpQaB.net]
- ループの中で
Force[v_] := -k*x-G*v;
と何度も定義し直すのは何か意図があるのですか？
素直にループの外で
Force[v_]:=If[ 0.5 k x^2>=10 n, -k x-G v ,必要なら何か定義...]
と定義しておくのはダメなのですか？
- 803 名前：１３２人目の素数さん mailto:sage [2019/07/02(火) 08:21:12.70 ID:AbMvpQaB.net]
- >>788
Mathematicaの入門書を読むことをお勧めします。
- 804 名前：１３２人目の素数さん mailto:sage [2019/07/02(火) 09:12:10.45 ID:ZLe7Frjb.net]
- 条件変えるループの中で定義したいんです
皆さんどうされてるんでしょう
本読み直します
- 805 名前：１３２人目の素数さん mailto:sage [2019/07/02(火) 10:22:56.70 ID:AbMvpQaB.net]
- ヒントになるかわかりませんが、
関数の引数に関数を与えることもできます。
- 806 名前：１３２人目の素数さん [2019/07/04(木) 00:31:56.85 ID:WjmhsYjy.net]
- 3200
ふうL@Fu_L12345654321
学コン1傑いただきました！
とても嬉しいです！
https://pbs.twimg.com/media/D-IuUuqVUAALnAB.jpg
https://twitter.com/Fu_L12345654321/status/1144528199654633477
(deleted an unsolicited ad)
- 807 名前：１３２人目の素数さん [2019/07/16(火) 17:48:13.35 ID:Z38S986h.net]
- どこのがっこんだよ
くだらんものをだすな 資源の浪費だ
- 808 名前：１３２人目の素数さん mailto:sage [2019/07/16(火) 18:49:37.97 ID:lFcROgSU.net]
- m.wolframalpha.com/input/?i=pikachu++curve
↑のθ(～)は何の演算子ですか？
- 809 名前：１３２人目の素数さん mailto:sage [2019/07/17(水) 07:29:54.91 ID:0FwmSYn+.net]
- >>795
右下見れ
- 810 名前：１３２人目の素数さん mailto:sage [2019/07/17(水) 08:59:41.06 ID:hcq7HK7z.net]
- どこの右下ですか？
- 811 名前：１３２人目の素数さん mailto:sage [2019/07/17(水) 09:00:40.70 ID:iozGF0zU.net]
- >>796
あ、スマホでは出ないんだな。PCで見てくれ
- 812 名前：１３２人目の素数さん mailto:sage [2019/07/17(水) 09:03:26.83 ID:iozGF0zU.net]
- https://i.imgur.com/9VPTZr7.jpg
- 813 名前：１３２人目の素数さん mailto:sage [2019/07/17(水) 10:01:55.79 ID:hcq7HK7z.net]
- ありがとうございます
- 814 名前：１３２人目の素数さん mailto:sage [2019/07/17(水) 16:05:53.82 ID:A5lmFvCD0.net]
- Clear[p,q];
Sum[If[q == p, q, 0], {q, 1, 3}] /. p -> 2
の結果として 2 になること期待したんだけれど，
qが返ってきてしまった…
なして？
- 815 名前：１３２人目の素数さん mailto:sage [2019/07/17(水) 18:52:19.50 ID:BWCFMYci.net]
- 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}]
だと２だよ。
- 816 名前：１３２人目の素数さん mailto:sage [2019/07/19(金) 07:51:55.90 ID:h7ko1Drl.net]
- If の属性に HoldRest が設定されているからだな。
Sum は Block を使って変数を置き換えていく。Block は変数を置き換えながら内側の式を評価するが、HoldRest のせいで第2引数の q は評価されずにそのまま残る。
第1引数は評価され置き換わるが、If は第1引数が True または False になるまで開かれないから、結果 >>802 の1行目のように展開され、答えは q となる。
- 817 名前：１３２人目の素数さん [2019/07/20(土) 11:01:21.03 ID:bSAoQnjE.net]
- 0130
ふうL@Fu_L12345654321
学コン1傑いただきました！
とても嬉しいです！
ｈｔｔｐｓ://pbs.twimg.com/media/D-IuUuqVUAALnAB.jpg
ｈｔｔｐｓ://twitter.com/Fu_L12345654321/status/1144528199654633477
(deleted an unsolicited ad)
- 818 名前：１３２人目の素数さん mailto:sage [2019/07/24(水) 13:02:02.44 ID:Hi3z86NN.net]
- >>801
Sum[If[q == p, Evaluate@q, 0], {q, 1, 3}] /. p -> 2
でも2が返るよ。
昔lispやってたので馴染めないんだが、、、。
- 819 名前：１３２人目の素数さん [2019/11/14(木) 23:23:19.39 ID:UAZBxqyu.net]
- 松坂和夫著『解析入門中』を読んでいます。
周囲の長さが一定 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
- 820 名前：１３２人目の素数さん [2019/11/14(木) 23:24:35.74 ID:UAZBxqyu.net]
- ↑の画像では、
s = 1 としました。
- 821 名前：１３２人目の素数さん [2019/11/14(木) 23:28:37.21 ID:UAZBxqyu.net]
- https://imgur.com/Fwn3FfO.jpg
↑実行可能領域の正三角形の重心で面積が最大になることが分かりますね。
- 822 名前：１３２人目の素数さん mailto:sage [2019/11/18(月) 22:50:48 ID:DwjlzO99.net]
- 数学でMATLABとかは使わないんですか？
- 823 名前：１３２人目の素数さん mailto:sage [2019/11/22(金) 23:40:53.18 ID:yWUYHlgs.net]
- Lightweightgrid使ってる人いる？
- 824 名前：１３２人目の素数さん [2020/01/01(水) 10:26:07.74 ID:ECfnCY0s.net]
- f が具体的にどんな関数なのかは教えずに、 f が t の関数であることをMathematicaに教えるにはどうすればいいのでしょうか？
- 825 名前：１３２人目の素数さん mailto:sage [2020/01/01(水) 12:56:51 ID:6ySof4cw.net]
- f[t]
- 826 名前：１３２人目の素数さん [2020/02/08(土) 22:12:11 ID:t6F7rXiU.net]
- iとjが添え字で、x[i,j]を総和記号で足し合わせるとき、 j≠iのときは足さないのはどうやって入力するの？
n m
Σ Σ x[i,j]
i=1 j=1
j≠i
みたいなやつ。意味わかります？
- 827 名前：１３２人目の素数さん mailto:sage [2020/02/08(土) 23:18:56.59 ID:ozXSNVVR.net]
- 全部足して対角成分を引けばぁ
- 828 名前：809 [2020/02/09(日) 08:32:21.15 ID:ofw5z8DB.net]
- ＞＞810 レスを有り難うございます。
ただ、それはわかってますが、もうちょっとスマートに書きたいということです。
- 829 名前：１３２人目の素数さん mailto:sage [2020/02/10(月) 01:08:16 ID:esjbRF9d.net]
- ラテン語または�
- 830 名前：Mリシア語の「マテーマタ」(古希:Μαθηματα) は「学ばれるべき事ども(Lesson)」という意味である。 []
- [ここ壊れてます]
- 831 名前：１３２人目の素数さん mailto:sage [2020/02/11(火) 09:40:57 ID:EOv4xygx0.net]
- 対角成分が発散してないなら引くかな
その方が早そうだし
発散するときは
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]
とか？
- 832 名前：809 [2020/02/14(金) 11:36:31 ID:bQ3El55L.net]
- >>817 お返事が遅くなってすみません。有り難うございます。それでうまくいきました。
- 833 名前：１３２人目の素数さん [2020/02/29(土) 09:41:39 ID:hF3qF87H.net]
- Wolfram Languageの実践的な解説をしているサイトを教えてください。
Mathematicaの本を見ても基本的なことしか書いていないですよね。
- 834 名前：１３２人目の素数さん [2020/02/29(土) 09:44:08 ID:hF3qF87H.net]
- 合成関数について質問です。
f[x_, y_] := x^2 + y^3
g[t_] := {t^2, 3*t+1}
とします。
f[g[t]] とはできないですよね？
f[ｖ_] := v[[1]]^2 + v[[2]]^3
g[t_] := {t^2, 3*t+1}
としないといけないんですか？
- 835 名前：１３２人目の素数さん mailto:sage [2020/02/29(土) 16:37:52 ID:/LxUqYKZ0.net]
- 引数がx,yの2個の場合しか定義してないのに
1個しか与えてないのだから未定義扱いになるのは当然では？
- 836 名前：１３２人目の素数さん mailto:sage [2020/03/01(日) 19:38:58 ID:OthADZfA.net]
- >>820
f@@g[t] で如何
- 837 名前：１３２人目の素数さん mailto:sage [2020/03/25(水) 23:44:43 ID:MX3wKSs0.net]
- Mathematicaの年間ライセンスが三月十四日円周率キャンペーンで安くなってるけど買おうかな？
- 838 名前：１３２人目の素数さん [2020/04/15(水) 17:37:36 ID:2L6Mosno.net]
- 斎藤正彦著『斎藤正彦線型代数学』を読んでいます。
正規行列はユニタリ行列によって対角化できるという定理が書いてあります。
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}
]
- 839 名前：１３２人目の素数さん [2020/04/15(水) 17:38:14 ID:2L6Mosno.net]
- 斎藤正彦著『斎藤正彦線型代数学』を読んでいます。
以下の定理が書いてありますが、 A = U*T （T が上三角行列）の場合しか証明されていません。
2.5.17 【定理】
任意の n 次正則行列 A はユニタリ行列 U と、対角成分が正の実数であるような上三角行列（下三角行列でもいい） T の積 U*T （T*U でもいい）として
一意的に表わされる。
- 840 名前：１３２人目の素数さん [2020/04/15(水) 17:38:46 ID:2L6Mosno.net]
- 他のケースについて書いておきます：
(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 が上三角行列）のタイプの分解の一意性から自明ですよね。
- 841 名前：１３２人目の素数さん [2020/04/15(水) 17:39:28 ID:2L6Mosno.net]
- 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]}];
- 842 名前：１３２人目の素数さん mailto:sage [2020/04/15(水) 20:00:57 ID:wsWfw9sb.net]
- >>824,825,826,827
こいつは数学の本スレで連投を繰り返している荒らしです
レスしないでください
- 843 名前：１３２人目の素数さん mailto:sage [2020/05/23(土) 02:52:22 ID:WyYxD60w.net]
- Mathematicaって二重根号に弱い
- 844 名前：１３２人目の素数さん [2020/05/31(日) 01:53:39 ID:T6EnPpfh.net]
- Wolfram Engine太っ腹
- 845 名前：１３２人目の素数さん mailto:sage [2020/05/31(日) 07:05:46.84 ID:Pc2IYhYT.net]
- Wolfram Engine+Jupyter notebookしか使ったことないから、本物のMathematicaがどんな感じなのかわからん
- 846 名前：１３２人目の素数さん [2020/05/31(日) 10:13:42.35 ID:T6EnPpfh.net]
- 正規MathematicaからWolfram Engine呼び出して使えば・・・・・
- 847 名前：１３２人目の素数さん [2020/05/31(日) 13:20:21.86 ID:T6EnPpfh.net]
- フロントエンドカーネルと並列処理のローカルカーネルとしては使える
並列処理のリモートカーネルとして使えないものか
- 848 名前：１３２人目の素数さん [2020/05/31(日) 15:00:12.58 ID:T6EnPpfh.net]
- RpiのMathematicaフロントエンドを高性能のWindowsマシンにインストールしたEngineにリモート接続してフル機能使えるやん
- 849 名前：１３２人目の素数さん [2020/05/31(日) 15:49:49.22 ID:RYSmKOHb.net]
- Jyupiter notebokk+Wolfram Engineより、Mathematica FE on RPi+Wolfram Engineのほうが遙かに幸せになれる
- 850 名前：１３２人目の素数さん mailto:sage [2020/05/31(日) 17:42:02.25 ID:Pc2IYhYT.net]
- ラズパイのMathematicaって仮想マシンでも動くの？
- 851 名前：１３２人目の素数さん [2020/05/31(日) 18:01:02.37 ID:RYSmKOHb.net]
- Raspberry Pi Desktop X86を仮想化作業中
ここに入っていればRaspberry Piも要らないね
Mathematica 6のフロントエンドで動けば助かったけど、6.0 FE+12.1 Engineの組み合わせはちょっと無理っぽかった
Macで動いたって言うブログ見たけれど
- 852 名前：１３２人目の素数さん [2020/05/31(日) 18:05:16.11 ID:RYSmKOHb.net]
- Engineをフロントエンドと並列計算用サブカーネルに使うことは簡単にできる
並列計算用サブカーネルにのみ使って、正規の8カーネル＋Engine 8カーネルで16カーネル並列計算できないか
なんとなくやり方はわかったけど、これは後でやってみる
- 853 名前：１３２人目の素数さん [2020/05/31(日) 20:38:01.89 ID:RYSmKOHb.net]
- やっぱx86版にはMathematicaはないみたいね
- 854 名前：１３２人目の素数さん [2020/06/01(月) 10:42:10 ID:MCEKpUDq.net]
- 5つの8カーネル並列計算プロセスを１台にインストールされたMathematicaフロントエンドで制御して計算する
|
 |
