〓Mathematica〓
..
369:通りがかり
03/02/06 00:03
>362
>367
>368 の方法で良いけど、ここも参考になるかも。
URLリンク(www53.tok2.com)
370:132人目の素数さん
03/02/06 01:21
>362(=>367だとは思うが)
det(A)≠0なら
x=(0,0,0,0,0,0,0,0,0,0,0)'
det(A)=0なら
NullSpace[A]で基本解の導出
371:132人目の素数さん
03/02/07 18:13
mathematica でPlotしたグラフをWordで表示する方法を教えてもらえませんか
372:132人目の素数さん
03/02/07 19:05
プロットしたグラフを右クリックして
copy as より bitmap で保存して word に貼り付ければ?
373:132人目の素数さん
03/02/07 20:20
372さん、ありがとうございます
374:132人目の素数さん
03/02/08 01:15
ボイコット・ザ・Mathematica。
もっとフリーソフトをつかう、作るべし。
MathematicaにせよMapleにせよ囲い込まれるだけだし、
ある日突然サポートがとまったりするリスクを考えよ。
バージョンアップで毎年金を毟り取られるのに快感を覚えるな。
他と著しく違う変なシンタックスに染まっても、システムが
替われば、そのスキルは無用の長物。
もっと全数式処理ユーザーが結束して、どうでもよいところは
共通仕様化することを要求していくべきだし、なるべくフリーな
ソースのある、内部で採用で採用されているアルゴリズムが
公開されているものをつかうようにしないと、本当に中世の
暗黒時代が再来するだけだよ。
375:132人目の素数さん
03/02/08 09:54
どなたか
Mathematicaの入門書を教えてください
どうしても自力じゃ使いこなせません
よろしくお願いします
あるていど使いこなしてみたいんです
376:132人目の素数さん
03/02/08 13:57
はやわかり Mathematica 第2版 榊原 進
は、どうですか
377:132人目の素数さん
03/02/08 14:09
メイプル マスキャド
と比較してどうなんでしょ?
378:132人目の素数さん
03/02/08 14:16
ますます出来るMathmatica =標準版=って
どうですか
わかりやすいですか?
でも14000は高いですよね?
379:132人目の素数さん
03/02/08 22:59
>374 じゃ、 MAXIMA ね。
まあ、正直僕も Mathematica の税金的な面は痛い所。
昔は3万円!でバージョンアップしてくれたのにね。。。
特に日本だと代理店経由しないとだめで、個人の気楽な
用途からどんどん遠ざかります。
>375
何したいのかによりますよ。
あ)とにかく操作入門レベル。コンピューター使用暦ゼロに近い
い)直面した問題に対処したい。コンピューター使用暦数年。
う)情報処理系的側面重視。コンピュータープログラミング要望。
でいろいろありますね。
まあ自分が知っている分野で、 Mathematica という言葉が
表題にあれば、大抵簡単な入門はのっているよ
380:132人目の素数さん
03/02/08 23:32
理工系の教養としての情報科学―Mathematicaによる入門
ローマン メーダー (著), Roman Maeder (原著), 時田 節 (翻訳)
URLリンク(www.amazon.co.jp)
381:132人目の素数さん
03/02/08 23:33
>>380
Mathematicaでのプログラミングにはお勧め
382:132人目の素数さん
03/02/09 11:42
Mathematicaって
普通のアプリにあるようなバグフィックスのアップデータって無いの?
バグなんか無いって事なんだろうか。
383:132人目の素数さん
03/02/09 12:20
>>382
バージョンアップ時に小出しにバグフィックス致します。
バージョンアップ料金は10万円以上かかります。
384:132人目の素数さん
03/02/09 12:40
俺の持ってるのは3.0.1.1xってなってる
3.0.0.0持ってる人もいるの?
385:132人目の素数さん
03/02/09 15:02
MuPADってだめなの?
386:132人目の素数さん
03/02/09 19:36
まてぃまてぃかで
世界征服したいんですけど?
387:132人目の素数さん
03/02/11 18:22
Plot した図を Showで表示します。
2つの違う関数を表示するのだけど、その線がどっちかって説明とA4に印刷するた
めの方法を教えてもらえませんか。
簡単に言えば、レポートに提出するためのものです
388:132人目の素数さん
03/02/11 22:38
>>385
フリーなのは学生と教員だけだろ。
389:132人目の素数さん
03/02/11 23:24
>387
help に線の太さやパターンを変える方法がのってるよ。
wolfram のサイトにいけば、日本語で読める可能性有り。
390:132人目の素数さん
03/02/15 23:16
>>388
個人利用はフリーと聞いたが
391:132人目の素数さん
03/02/21 15:50
URLリンク(risa.cs.ehime-u.ac.jp)
392:132人目の素数さん
03/03/01 00:06
mathematicaのstudent版って卒業してしまうと使えなくなってしまうのですか?
買ってから4年間、みたいな制約があるんでしょうか?
393:Lisper
03/03/08 15:57
>>392 GNU Maximaなら永遠に無料でつ.シミュレーション板Maxima擦れを参照.
どうせMathematicaをそれ程クリティカルな業務に使ってないのでしょうから,Maximaで十分でつよ.
3Dグラフがマウスでぐりぐり回せますぜ.どうでつか.
394:山崎渉
03/03/13 13:23
(^^)
395:132人目の素数さん
03/03/16 08:28
396:とある馬鹿 ◆BAKAB.w.so
03/03/16 12:33
Mathematica を早く使いこなせるように私はなりたい ううう ( つДT)
大学でしかできないからなぁ・・・
やっぱり自腹で買ったほうがいいのかな
それとも研究室の金で買ったほうがいいのかな・・・
397:132人目の素数さん
03/03/16 13:00
大学でできれば十分ジャン。
たまには家に帰りたいもんだ。
398:とある馬鹿 ◆BAKAB.w.so
03/03/16 14:16
>たまには家に帰りたいもんだ
なんか(・∀・)カコイイ!!!!!
399:質問!
03/03/16 14:27
(関数を定義する)
仮に・・・a_{i}:実数
s[a_{1},a_{2},....,a_{n}]:=a_{n-1} : 与えられた数のn-1番目を取り出す
という関数sを定義したいですが,変数をリストとしてではなく,
与えたいと思ってます.どうすれば良いのでしょうか?
ちなみに,一般のnについてできるようにしたいです.
ちょうど,Max[a_{1},a_{2},....,a_{n}]の入力の仕方で
答えを得るようにしたいということです.
s[{a_{1},a_{2},....,a_{n}}]から,a_{n-1}を得るのは,簡単ですけど,
s[a_{1},a_{2},....,a_{n}]:これのときは?#n:一般の自然数#
400:399
03/03/16 18:21
誰かわかる方いませんか?
401:132人目の素数さん
03/03/16 20:27
>>399-400
これでいいだろ。
In[10]:=
s[x__]:=List[x][[-2]]
In[11]:=
s[1,2,3,4,5,6]
Out[11]=
5
402:132人目の素数さん
03/03/16 20:40
これでもいい。
In[16]:=
s[x__]:={x}[[-2]]
In[15]:=
s[1,2,3,4,5,6]
Out[15]=
{1,2,3,4,5,6}
403:402訂正
03/03/16 20:53
In[16]:=
s[x__]:={x}[[-2]]
In[17]:=
s[1,2,3,4,5,6]
Out[17]=
5
404:399-400
03/03/16 21:43
なるほど!!!
そういうことで回避できるんだー!
401-403 : ありがとうございました!
うーんこれは使える・・・
405:132人目の素数さん
03/03/18 22:30
>>183
先生, 質問です.
StudentからProfessionalにUpgradeすると,
パスワードの発行形態は後者に準じたものに変わるのでしょうか?
406:132人目の素数さん
03/03/18 22:39
>>393
使えるのなら、使っていいんじゃないの?
version3.0 for Students を持っているけど、いまだに使っているよ。
インストールしなおす時には、もちろん、きーじぇ・・・スミマセン。
407:kawabata君の友人
03/03/19 00:48
>>405
はい、発行形態がかわります。
ただし、Proffessionalへアップグレードしても
(1)ホームユース(名前は確かじゃないですが「一本」のライセンスで、
いつも使用している研究室や職場だけでなく、自宅でも使用することが
できるシステム)(この場合は自宅の端末専用に1年間有効のMathIDベースの
パスワードが発行されます。当然CDもMathID仕様のものが追加で交付
されます。1年ごとに追加料金がかかります
(2)そのほかにも北米ライセンスを、一時的に他の地域で使用する場合、
または、その逆で、インターナショナルライセンスを北米で使う場合
などにも、時限パスワードが追加交付されるはずです
408:132人目の素数さん
03/03/19 01:24
マス 書いてもいいか
409:132人目の素数さん
03/03/19 09:47
>>407
ということは,
For StudentからProffessionalへアップグレードしても"MathIDは常に必要"ということですか?
最初からProffessionalの場合みたく, 実質的には(自分の所有する)複数のPCに入れ放題ということにはならないのか……
410:kawabata君の友人
03/03/20 00:43
>>409
プロフェッショナルはMathID依存でない18桁〜25桁のパスワードがでる
はずです。しかしホームユースや地域外での利用に対しては、別途料金がかかり
その料金徴収、有効期限管理、追加で使用を認めるマシンを限定し特定する必要
があるため、MathID形態のCD並びにパスワードを追加交付しているのでしょう。
411:baggio
03/03/20 02:06
質問です。mathematicaを使って、非線形連立方程式を解きたいのですが、どうしたら
いいのでしょうか?初心者の質問で申し訳ないですが、教えてください。
412:132人目の素数さん
03/03/20 22:41
>410
うーん。通常の Professional ではなく、プレミアサポート
(定価30%、PC版なら約10万)でホームユースが
出来る、と考えた方がいいでしょう。
413:132人目の素数さん
03/03/20 23:44
>411
代数方程式なら
Solve[方程式のリスト,未知数のリスト]
とします。
それ以外の場合はFindRootを使います。
その場合は解の初期近似を与える必要があります。
FindRoot[方程式のリスト,{未知数,その初期値}のリスト]
414:132人目の素数さん
03/03/23 03:14
Mathematicaで整数計画問題解けますか?
415:baggio
03/03/28 01:10
>413 132人目の素数さんありがとうございます。
一回やってみます。わざわざ解答ありがとうでした。
416:132人目の素数さん
03/04/15 01:53
Linux版のMathematicaで、スタンドアローンのプロフェッショナル版を
買おうかと思っているのですが、Linuxのバージョンをアップデート
したり、あるいはハードを変更したり、マザボードやCPUをとりかえたり
ネットワークカードを変えたり、IPアドレスを変えたりしたら、
使えなくなってしまうのではないかと恐れています。
何が、保存されていれば(たとえばMACアドレスとかIPアドレス)、
環境を変化させても、ライセンスを取り直さなくても使えるのでしょうか?
417:132人目の素数さん
03/04/15 16:13
MathematicaではLogListPlotやLogLogListPlotのような片対数グラフや両対数グラフを
表示できて便利なんですけど、エラーバー付きのLogListPlotやLogLogListPlotのような
ものはないんでしょうか?
対数表示でない通常のエラーバー付きグラフならErrorListPlotがあるんですが
ErrorLogListPlotやErrorLogLogListPlotのようなものがないのか知りたいです。
418:132人目の素数さん
03/04/15 16:25
>>417
臭うぞ! Q.manの臭いがプンプンするッ!
419:mathmania ◆uvIGneQQBs
03/04/15 16:45
>> 417
Exp,Log,ErrorListPlotを組み合わせて自分で関数を作れ。
>> 418
417がQ.manならば、Q.manは変数変換の公式なぞろくに使えないだろう。
420:132人目の素数さん
03/04/15 18:01
Q.manってだれ
421:Q.man
03/04/15 18:23
ゎたすがQ.manです。
422:132人目の素数さん
03/04/16 10:14
あっそ
423:P.man
03/04/16 23:09
ゎたすがP.manです。
424:132人目の素数さん
03/04/19 18:31
(・∀・)ゲハハハハ
425:山崎渉
03/04/20 03:59
∧_∧
( ^^ )< ぬるぽ(^^)
426:apico
03/04/23 20:52
Mathematicaで3次元のグラフを表示させる方法を教えて下さい。
以前販売元へ質問状をFAXとメール送ったのですが、その後連絡がありません。
日本人には冷たいのかな??
閑話休題、
[1]曲面z=Sqrt[x*y]と3平面z=0,x=2,y=3によって囲まれた立体を表示させたい
のと
[2]曲面z=Sqrt[a^2-y^2](y≧a/2)と4平面x=0,x=2a,y=a/2,z=0によって囲まれる
立体を表示させたいのです。
表示させようと思ってすでに半年。自分でもどうやって聞いたらよいのか分からな
いまま経ってしまいました。もしこれが無理なら日本語でこういった質問の出来る
メールアドレスを教えて下さい。御願いします。
427:mathmania ◆uvIGneQQBs
03/04/25 18:11
>> 426
表示しようと思って半年もかかるくらいなら、手書きの方が正確で早い。
だが、それでもコンピュータで表示したいなら、
? Plot3D
を使うか、
あるいは
<< Graphics`をやって
? ParametricPlot3D
でヘルプを出して関数を理解してくれ。
428:132人目の素数さん
03/04/25 18:50
>>426
誰でも、すぐできることだから。
429:132人目の素数さん
03/04/25 19:45
はじめまして
今日mathematica 4.2 for students買ってきました
今後お世話になるかもしれませんがよろぴく
430:132人目の素数さん
03/04/27 20:40
>>427, 428
「囲まれた立体」を表示させるのは
Plot3Dでは不可能だし
ParametricPlot3Dでも[1], [2]を表示するのはすくなくとも困難だ。
まぁ428には不可能だな。
431:132人目の素数さん
03/04/28 22:16
>>426
? ClipFill
432:132人目の素数さん
03/04/29 07:49
ClipFillねー
これでいいのかな?
Plot3D[Sqrt[x*y],{x,0,2},{y,0,3},
PlotRange->{0,3},
PlotPoints->30,
ClipFill->Automatic,
BoxRatios->{2,3,3}]
433:132人目の素数さん
03/04/29 22:25
>>432
Clipしてないのだが?
434:132人目の素数さん
03/05/01 18:17
mathematicaって将来性あるの?
435:132人目の素数さん
03/05/02 08:41
MaximaやMuPADなど無料のソフトがもっと
グレードアップしてくると
今の商売は難しくなるね。
我々はその方が有り難いけど。
436:132人目の素数さん
03/05/02 08:42
優良スレに認定されてます
◇ 「知」の欺瞞
スレリンク(philo板)
◇ アーベル賞にセール
スレリンク(newsplus板)
◇ 数学の本
スレリンク(math板)
◇ 代数学
スレリンク(math板)
◇ 数学セミナー
スレリンク(math板)
◇ Poincare 予想 解決
スレリンク(math板)
◇ Lie群・Lie環
スレリンク(math板)
◇ グロタンディック
スレリンク(math板)
◇ 基礎論
スレリンク(math板)
◇ アーベル賞にセール
スレリンク(math板)
◇ 素数判定は「決定的」多項式時間
スレリンク(math板)
◇ スペクトル系列
スレリンク(math板)
◇ unix と数学
スレリンク(unix板)
437:132人目の素数さん
03/05/03 20:39
いまいち盛り上がらないので。
Mathematica programmingコンテスト。
お題は、
トランプ(Range[52])をシャッフルする関数を書け。
とか言ってみる。
438:132人目の素数さん
03/05/03 20:43
>437
ドナルド・トランプですか?
439:437
03/05/03 20:53
>>438
普通のトランプです。
ドナルドさんが52人おられるのでしたらそれでも結構です。
440:132人目の素数さん
03/05/15 01:19
>>437
こんなのとか。
Shuffle[x_List] :=
Prepend[Shuffle[Rest[#]], First[#]] &@
RotateLeft[x, Random[Integer, Length[x] - 1]]
Shuffle[{}] = {}
Shuffle[Range[52]]
441:山崎渉
03/05/21 22:59
━―━―━―━―━―━―━―━―━[JR山崎駅(^^)]━―━―━―━―━―━―━―━―━―
442:132人目の素数さん
03/05/25 05:59
29
443:山崎渉
03/05/28 14:44
∧_∧
ピュ.ー ( ^^ ) <これからも僕を応援して下さいね(^^)。
=〔~∪ ̄ ̄〕
= ◎―◎ 山崎渉
444:132人目の素数さん
03/06/02 16:04
Mathematicaが原因でMATHEMATICA.EXEにエラーが発生しました。
Mahtmeticaは終了します。
問題が解決しない場合は、コンピュータを再起動してください。
と出たので再起動したのですが、何回再起動しても同じエラーが出て
起動してくれません。
フォント関係の設定をいじった後このエラーが出たので
元に戻したいと思っているのですが
設定ファイルはどこにあるかわかりませんか?
レジストリは探してみたのですが、わかりませんでした。
Windows MeでMathemtica4.2のシングルライセンスのを使っています。
445:132人目の素数さん
03/06/02 18:40
再インストールしる!
446:132人目の素数さん
03/06/02 20:41
Mathematicaって使いにくくない?ライブラリとか凄いとは思うけど
447:mathmania ◆uvIGneQQBs
03/06/02 20:46
コマンドが長いのは確かに使いづらい。
もっとも、私は最近Mathematicaをやっていない。
C言語(ForTranも少し)を使ってる。
理由はというと、、、Mathematicaが使えなくなったから。
まぁ、使いやすさの面でいえば、大きなプログラムには、C言語の方がいいだろう。
448:444
03/06/02 23:50
>>445
再インストールしても同じエラーが出るんです。
設定ファイルが残ってるんじゃないかと疑ってるのですが
その場所が分かりません。
449:132人目の素数さん
03/06/03 22:42
>>446
原始的なプログラミング言語と不統一な仕様の
mapleよりは使いやすいと思う。
言語としてつかえてなんぼだから
素人にはお薦めできない、って感じはあるかもね。
>>447
補完機能をつかうべし。
450:132人目の素数さん
03/06/07 17:45
初心者にお勧めのMathemticaの
本は何でしょうか?
451:132人目の素数さん
03/06/07 21:09
>448
OSを再インストールしる!
それでだいじょうぶだ! 安心しろ
>450
白石さんのがいいぞ!
安心して買え
452:132人目の素数さん
03/06/07 21:25
あなたが探してるのってこれだよね?二日間無料だしね♪
URLリンク(alink3.uic.to)
453:132人目の素数さん
03/06/07 22:15
mathmaniaスレかと思った
と言ってみる
454:132人目の素数さん
03/06/08 00:02
>>450
Mathematicaプログラミング J.ゲイロード・N.カーミン・R.ウエリン
を読んで、その後で
プログラミングMATHEMATICA. R・メーダー 著、時田 節
を読むとよくわかると思われ。
特にメーダーの本は(全くの初心者以外には)一読の価値あり。
非常に下手な翻訳だけど。
455:132人目の素数さん
03/06/08 00:09
>>248
Windows Me はやめた方がいいですよ。
2KかXPにしましょう。
(マジレス)
456:132人目の素数さん
03/06/08 02:09
Mathematicaって有効数字が6桁の表示で、例えば
1.23456789とか1234.56789って直接入力すると
1.23456や1234.56と表示されます。
# 表示だけがそうなるだけで内部的には入力した数のままです
これを以降の入力からすべて指定した有効数字で表示させたいのですが
どうやったら良いのでしょうか?NumberFormではそれで指定した場合だけ
表示されるのでこれを以降の入力すべてで有効にしたいのでよろしく
お願いいたします。
457:132人目の素数さん
03/06/08 10:38
>>451 >>454
どうもありがとう!
458:132人目の素数さん
03/06/08 10:49
>>456
ちょっと邪魔くさいけど、関数を定義すればどう?たとえば
F[x_]:=N[x, 30];
F[E^3]
とするとe^3が30桁表示されるはず。
459:mathmania ◆uvIGneQQBs
03/06/08 13:02
456の方法でも有効数字は確かに表示される。
だが、SetPrecision[E^3,Infinity]か、SetAccuracy[E^3,Infinity]はどうだろう?
まぁ、分数表示になるが、その数に2^52(?)を掛ければ小数点表示がわかるだろう。
あぁ、それと、私はMathematica復活したわけではないので。
460:444
03/06/08 21:28
>>451
結局OS再インストールしましたよ。。。
したくなかったけど。
>>455
大学のPCだから自由にできません。
経費で2000買ってくれなかった。
461:132人目の素数さん
03/06/09 17:46
>>458-459
Mathematicaでの以降の入力すべてで指定した有効数字で表示させたいということ
でしょうから、その方法ではその関数内の引数にだけ有効であって意味がありません。
何かシステムの設定でできそうな気もしますがそれらしいものが見当たりませんね。
どうやるんだろうか?
462:132人目の素数さん
03/06/09 18:00
>>461
SetPrecisionは計算精度を変更させるだけであってその表示は依然として
有効数字6桁です。
463:132人目の素数さん
03/06/10 19:55
>>456
編集→環境設定→書式設定→式の書式→表示設定→PrintPrecision
464:132人目の素数さん
03/06/10 23:57
>>463
書式設定でしたか。見落としてました。大変助かりました。ありがとうございます。
465:mathmania
03/06/17 23:17
駄スレ漏れスレage
466:132人目の素数さん
03/06/18 00:04
あなたが探してるのってこれだよね?この中にあったよ♪
URLリンク(endou.kir.jp)
467:132人目の素数さん
03/06/21 17:16
4次元データを視覚化したいのですが、mathematicaに限らず
何かよい方法(ツール)はありませんか?
等高線図を複数枚、平行になるよう並べて
ナナメ上から見下ろすようにすればよいと思うのですが、
よいアプリケーションを探しています。
他の方法でも、可能ならば何でもいいです。
468:132人目の素数さん
03/06/22 17:01
>>467
いってることが良くわかりません。等高線図を斜めから見下ろす?
mathematicaでの実現ならば、三次元にスライスしたときの
データをどうやって表示するかという話になるのだろうけど、
それはデータの形式に依存するからなんともいえないし。
469:467
03/06/23 13:26
>>468
間違いがありました。すみません。
まず画面に出したいのは4次元データ中の3次元分なので、
「3次元データを表示したい」と訂正します。
データは格子状データです。
古い人間なのでfortran形式で言うと
do k=1,nz
do j=1,ny
do i=1,nx
write(6,*) x(i),y(j),z(k),data(i,j,k)
end do
end do
end do
の形式の等間隔のデータと思って下さい。
やりたいのは、これまた解かりにくい表現かもしれませんが
gnuplotで
set contour base
unset surface
set view 60,10
くらいの時にプロットした等高線図を複数枚、一つの画像にプロットしたいのです。
470:132人目の素数さん
03/06/23 20:03
>>469
gnuplot云々はわかりませんが、
ListContourPlot3Dかな?
471:132人目の素数さん
03/06/24 22:42
Ver.5 ですね。
それに $1880 5.0 Eng vs \448800 日本語版 4.2only
US ってこれでも Premier Service 付きだって。。。。
潮時かね
472:467
03/06/26 18:39
忙しかったので、見るの遅くなって済みません。
>>470
>ListContourPlot3Dかな?
手持ちのマニュアルには
ListContourPlotとListPlot3Dしか見付かりません。。
(ver5からの新機能ですか?)
それら二つから推測すると、
ListContourPlot3D[{{{},{}},{{},{}}}]
のような3重ネストの配列を入れてプロットするのでしょうか?
値段高い…
473:132人目の素数さん
03/06/26 23:19
>>472
4.2にはあります。
Onlinehelpで出てくるはずです。
または、
<<Graphics`
?ListContourPlot3D
474:467
03/06/27 14:28
みなさま、御親切に
本当にありがとうございます!
ちょっといろいろ試してみます。
475:132人目の素数さん
03/06/27 15:22
足してnになる1以上の整数の組み合わせを並べる時、並べ方が最大になるような
整数の組み合わせをMathematica(あるいはC言語でも)で求める方法を教えてください。
わかりにくいので例を挙げますと n=5 のとき
足して5になる1以上の整数の組み合わせは
5, 41, 32, 311, 221, 2111, 11111 であり
並べ方が何通りあるかはそれぞれ 1, 2, 2, 3, 3, 4, 1
で最大値は 2111 の場合で 4 に、
n=6の場合は 321 と 2211 の場合で最大値は 6 になります。
一般に、数字の個数(311なら3個、11111なら5個)がm個、1がk_1個、2がk_2個、…nがk_n個の
場合、並べ方はm!/(k_1!・k_2!・・・k_n!)になります。
476:mathmania ◆uvIGneQQBs
03/06/27 15:36
Re:>475
最大値をどうやって計算するかは分かるだろう。
あとは、k_1,…,k_nの与え方だが、
これはk_n,k_{n-1},…,k_2,k_1の順に再帰的に与えればいいのではないか?
n=5のときを例に挙げると、
まず、k_5=1ととり、1*k_1+2*k_2+3*k_3+4*k_4=0になるk_1,k_2,k_3,k_4の組合わせは0,0,0,0しかない。
次にk_5=0ととるとき、1*k_1+2*k_2+3*k_3+4*k_4=5となるk_1,k_2,k_3,k_4を選ぶことを考える。(a)
(a)の中でk_4=1ととり、1*k_1+2*k_2+3*k_3=1となるようにk_1,k_2,k_3を選ぶようにする。(k_1=0,k_2=k_3=0しかない。)
(a)の中で、今度はk_4=0ととり、1*k_1+2*k_2+3*k_3=5となるようにk_1,k_2,k_3を選ぶ。(b)
(b)の中で、k_3=1ととると、(k_2,k_1)=(1,0)と(k_2,k_1)=(0,2)となる。
(b)の中で、k_3=0ととると、k_2=2,1,0の3通りにわかれる。
k_2=2ととるとき、k_1=1となり、k_2=1ととるとき、k_1=3となり、k_2=0ととると、k_1=5となる。
この文では分かりにくいかもしれないが、とにかくやってみてくれ。
477:474
03/06/27 16:10
>>476
その考え方はわかりました。
初心者のものでプログラムの書き方を勉強してきます。
478:132人目の素数さん
03/06/27 16:11
477=475です
479:132人目の素数さん
03/06/27 17:06
>>475
どうぞ。
<< DiscreteMath`Combinatorica`
<< Statistics`DataManipulation`
Max @@ ((Length[#]!)/(Times @@ (#[[1]]! &) /@ Frequencies[#]) &) /@
Partitions[5]
480:132人目の素数さん
03/06/27 20:12
Mathematica Ver.5の情報はどこにあるんですか?
何が変わるのでしようか
おしえてね!
481:132人目の素数さん
03/06/27 20:24
過去にレスがついていなかったのですが、
ますます出来るMATHEMATICA-標準偏-
URLリンク(www.catkk.co.jp)
どうでしょうか?
482:132人目の素数さん
03/06/27 20:25
>>480
URLリンク(www.wolfram.com)
483:475
03/06/27 23:37
>>479
よくわからないのですが、Partitions[5]の 5 を n に変えれば
足してnになる1以上の整数の組み合わせの並べ方が最大に
なる整数の組み合わせが求められるんでしょうか?
484:483
03/06/27 23:39
すいません。” 5 を n に” は ” 5 を任意の整数に”という意味です
485:132人目の素数さん
03/06/27 23:43
>>483
f[n_Integer] := Max @@ ((Length[#]!)/(Times @@ (#[[1]]! &) /@ Frequencies[#]) &) /@
Partitions[n]
ってこと?
もしかして、数学っぽいnについての式を欲しがってるんですか?
それならそれは数学の問題。
プログラムで求めることはできません。
Cでなんて論外。
486:479 == 485
03/06/27 23:45
>>484
See f in >>485.
487:475
03/06/27 23:54
>>486
ありがとうございます。手元にMathematicaがないので
少し勉強して、学校ので試してみたいと思います。
数学的に(nを使って)解くことが出来たらそれが一番いいのですが、
それはここで聞くことではないので他のスレで聞いてみます。
488:132人目の素数さん
03/07/04 13:36
再セットアップしたらアクティベーションできなくなった
489:132人目の素数さん
03/07/04 19:39
URLリンク(www.wolfram.com)
うわーん!
Ver.4.2の情報しかないじゃないか。
しかも日本語だー
490:132人目の素数さん
03/07/05 01:36
>489
どこかに English ってあるからそこから英語ページへ。
トップにあるよ。
まあ大規模行列と数値計算の速度改良が主ですね。
試行錯誤でアルゴリズム作るときは、全作業時間において
コンパイル系とためはるつもりなのかな
なお日本では日本語版のみしか売れない。
誰だよな日本語ヘルプキット強制バンドル計画したやつは。
491:132人目の素数さん
03/07/05 21:11
URLリンク(www.wolfram.com)
に五版の概要がありますね。
英語版でいいから新しいやつを動かしてみたいねえ。
ほんと。
492:132人目の素数さん
03/07/05 23:02
私が作りました♪
URLリンク(nuts.free-city.net)
493:488
03/07/07 15:52
誰か助けてくれ
494:132人目の素数さん
03/07/07 19:47
>>493
サポートに連絡しろ。
495:132人目の素数さん
03/07/08 01:13
>>494
わかったやってみる
dくす
496:475
03/07/08 16:22
>>479
Max @@ ((Length[#]!)/(Times @@ (#[[1]]! &) /@ Frequencies[#]) &) /@
Partitions[5] を実行したところ
Part::"partd": "パート指定5[[1]]の長さはオブジェクトの深さを超えています."
1/5[[1]]!
となったのですが、直すべきところがわからないので教えてください。
497:132人目の素数さん
03/07/08 21:15
>>496
ちゃんと上の二行を入力してください。
# そのエラーはPartitions[5]が評価されていないため。
498:132人目の素数さん
03/07/08 22:23
最大尤度法で推定をしていて、結構複雑なモデルなんですが、推定したい変数が5個
あるんです。それで、尤度関数を作って、偏微分してFindRootつかって求めようした
んですが、出来ないんです。
そこで、正規分布使って、普通に正規分布を推定してみようって思って、同じような
方法を使ったんですが、これでもうまくいきません。エラーとしては
{{初期値の尤度}},{{初期値の尤度}}は、{a,b}={1,2}では長さ2のリストになって
いません
というような感じです。最大尤度法をやってるプログラムのサンプルどっかないかなあ
って探したんですが、見つからないし・・・・。泣きそうなんですけど、どなたかお助け
ください。
499:132人目の素数さん
03/07/08 22:43
>>498
その説明で他人がなにかわかると思うか?
> んですが、出来ないんです。
ふーん出来ないんだ。(どこがどう出来ないのか全くわからん)
> {{初期値の尤度}},{{初期値の尤度}}は、{a,b}={1,2}では長さ2のリストになって
文法間違いだな。 (syntax errorなんだからソース見ないとどうしようもない)
> って探したんですが、見つからないし・・・・。泣きそうなんですけど、どなたかお助け
まぁ冷静になれ。
500:132人目の素数さん
03/07/08 23:08
>>497
わかりました。
501:498
03/07/08 23:11
そうですか・・・・すいません。どっかサンプルプログラム探すしかない
ですかね・・・
正規分布を推定する場合、乱数を出して、その後、尤度関数を出せたとして
K1[a_,b_]=D[L[a,b],a]
K2[a_,b_]=D[L[a,b],b]
として、
FindRoot[{K1=0,K2=0},{a,-1},{b,2}]
としてるんですけど。。この部分では、間違いがないとは思ってるんですが
尤度関数が間違ってるのかな・・・・。はあー泣きそう。。。
502:132人目の素数さん
03/07/08 23:12
☆頑張ってまーす!!☆女の子が作ったサイトです☆
☆見て見て!!
URLリンク(yahooo.s2.x-beat.com)
503:132人目の素数さん
03/07/08 23:20
>>501
Mathematicaの基礎からやりなおしってかんじ。
とりあえずダメなところは指摘しとくが、
マニュアルか教科書(上でも挙がってる)読み直せ。
> K1[a_,b_]=D[L[a,b],a]
= と := の違いが良くわかってるのなら何もいわないが。
> FindRoot[{K1=0,K2=0},{a,-1},{b,2}]
K1 = 0 は K1 という変数を 0にする。
書くなら {K1[a, b] == 0, K2[a,b] == 0} 。
504:498
03/07/08 23:20
あー勿論
FindRoot[{K1[a,b]=0,K2[a,b]=0},{a,-1},{b,2}]
としてます。。。。。
サンプルプログラム載ってるサイトってないですかねえ・・・・。。。
505:503
03/07/08 23:22
>>504
その2行目をもういちど読み直せ。
506:498
03/07/08 23:23
:=と=の違いってえーと、すみません、良く分ってないです。。。
確か:=は、なんか固定するみたいな意味だった気が・・・・
507:498
03/07/08 23:25
>>505
==ですね・・・すみません、それは、やってます・・・・。
508:503
03/07/08 23:28
>>506
で、聞いてる方としてはあいかわらずどこがうまく行かないのか
わからないんだが。
で、どううまくいかいの?
エラーメッセージは?
509:498
03/07/08 23:37
エラーメッセージは、FindRootのところなんですよー。
>>498に書いたようなエラーメッセージが出ます。
FindRootの文法は、合ってると思うんですが。。。
510:503
03/07/08 23:45
>>509
498を読んだところでは第一パラグラフと第二パラグラフは
別の手段を試してるようにしか読めない(FindRootで失敗したから別の手段)。
どの関数でエラーが出力されたかもわからない。
エラーがあなたの言語に翻訳されているがそれが適切かどうかわからない。
すでに書いたがsyntax errorはソースをみないと何もわからない。
技術的なことを質問する方法をしらない、
シンタックスエラーが自分で解決できない、
レベルだったらソース張るべきだよ。
変にプライドがあるのかそういう「初心者」は多いけどね。
情報を出さずに、自分のやってることは間違ってない!、みたいな。
たんに誰かが作ったライブラリが欲しいのなら、
まず、googleなりwolframのサイトなりを検索し、
もっと適切な質問のしかたをすべき。
おれは今日はもう寝ます。
511:503
03/07/09 00:02
>>509
寝ようと思ってたが気づいたので。
多分 L が リストを返してるからだ。
512:132人目の素数さん
03/07/09 14:24
ndist2=NormalDistribution[0,4];
p1=Table[ Random[ndist2],{100}];
L[a_,b_] = Product[(1/(Sqrt[2*π]*a))*Exp[(-({p1[[i]] - b}^2/(2 *a^2)))],{i,1,100}]
K1[a_,b_]=D[L[a,b],a]
K2[a_,b_]=D[L[a,b],b]
FindRoot[{R1[a,b]==0,R2[a,b]==0},{a,{1,2}},{b,{-2,4}}]
ソースは、これです。
>>511
あーそれは、どうすれば、解消できますかね。リストで返さずに関数として認識
してもらわないとつらいっすね・・・。
513:498
03/07/09 15:12
あっ最後の式は、
FindRoot[{K1[a,b]==0,K2[a,b]==0},{a,{1,2}},{b,{-2,4}}]
です。
514:132人目の素数さん
03/07/09 16:05
今課題をやっています。
誰か偏微分係数の出し方を教えてください。
ちなみにMathematika 4.1を使用しています。
ヘルプの使い方がわからなくて・・・
515:514
03/07/09 16:19
解決すますた。
失礼。
516:503
03/07/09 17:37
>>498
ソースはコピペしてそのまま張れ。
タイプミスもあるしな。
πってPiだよな?
<<Statistics`NormalDistribution`
してるんだよな?
原因はp1のまえにある'{'だ。'{'はリスト構築子。
L[a_,b_] = Product[(1/(Sqrt[2*π]*a))*Exp[(-({p1[[i]] - b}^2/(2 *a^2)))],{i,1,100}]
↓におきかえろ。
L[a_,b_] = Product[(1/(Sqrt[2*π]*a))*Exp[(-((p1[[i]] - b)^2/(2 *a^2)))],{i,1,100}]
517:132人目の素数さん
03/07/09 21:13
このスレッドの人から見てS言語はどうですか?
社会学系だとmatlabやmathematicaより使われていますが…
518:498
03/07/10 16:44
>>516
ありがとうございました。それで正規分布の推定は、できました。ただ、本当に
求めたいのは、
L[a1_, a2_, b1_, b2_, b3_] = Sum[y1[[i]]*Log[F[a1 - (b1*
x1[[i]] + b2*x2[[i]] + b3* x3[[i]])]] +
y2[[i]]*
Log[F[a2 - (b1*x1[[i]] + b2*x2[[i]] + b3*x3[[i]])] - F[
a1 - (b1*x1[[i]] + b2*x2[[i]] +
b3*x3[[i]])]] +
y3[[i]]*Log[
1 - F[a2 - (b1*x1[[i]] + b2*x2[[i]] + b3*x3[[i]])]], {i, 1,
100}];
という結構複雑なもので、これも{を(にしてみたんですが、やはりリストで返してる
みたいです。Sumの扱い方が悪いんでしょうか・・・・
てか、見にくくてすみません。
519:132人目の素数さん
03/07/10 19:28
>>518
a1,a2,b1,b2,b3,F,y1,... のどれかにリストが入ってんだろ…
同じことを何度も書くが、部分だけ張るなよ…
プログラミング自体初めてなのかも知れないが、レベル低すぎる。
一冊Mathematicaの本読め。
520:498
03/07/10 21:10
すみません、一遍に貼れないんで2回に分けますけど、これが全部です。。。
<<Statistics`ContinuousDistributions`
ndist=NormalDistribution[0,1];
x1=Table[ Random[ndist],{100}];
x2=Table[ Random[ndist],{100}];
x3=Table[ Random[ndist],{100}];
y=Table[2*x1[[i]]+3*x2[[i]]-x3[[i]],{i,100}];
y1=Table[If[y[[i]]<=-5,1,0],{i,100}];
y2=Table[If[-5<y[[i]]<=5,1,0],{i,100}];
y3=Table[If[5<y[[i]],1,0],{i,100}];
521:498
03/07/10 21:11
F[x_] = Integrate[Exp[z]/{1 + Exp[z]}^2, {z, -Infinity, x}];
L[a1_,a2_,b1_,b2_,b3_]=Sum[y1[[i]]*Log[F[a1-(
b1*x1[[i]]+b2*x2[[i]]+b3*
x3[[i]])]]+
y2[[i]]*Log[F[a2-(b1*x1[[i]]+b2*x2[[i]]+b3*x3[[i]])]-
F[a1-(b1*x1[[i]]+b2*x2[[i]]+
b3*x3[[i]])]]+y3[[i]]*Log[1-
F[a2-(b1*x1[[i]]+b2*x2[[i]]+b3*x3[[i]])]],{i,1,100}];
K1[a1_,a2_,b1_,b2_,b3_]=D[L[a1,a2,b1,b2,b3],a1];
K2[a1_,a2_,b1_,b2_,b3_]=D[L[a1,a2,b1,b2,b3],a2];
K3[a1_,a2_,b1_,b2_,b3_]=D[L[a1,a2,b1,b2,b3],b1];
K4[a1_,a2_,b1_,b2_,b3_]=D[L[a1,a2,b1,b2,b3],b2];
K5[a1_,a2_,b1_,b2_,b3_]=D[L[a1,a2,b1,b2,b3],b3];
FindRoot[
{K1[a1,a2,b1,b2,b3]\[Equal]0,
K2[a1,a2,b1,b2,b3]\[Equal]0,
K3[a1,a2,b1,b2,b3]\[Equal]0,
K4[a1,a2,b1,b2,b3]\[Equal]0,
K5[a1,a2,b1,b2,b3]\[Equal]0},
{a1,0},{a2,1},{b1,1},{b2,3},{b3,1}
]
522:503
03/07/10 21:20
>>521
F[x_]の定義を良く読め。
それからMathematicaでプログラミングするときは
出力を確認しながらやれ。
F[x_] の行をセミコロン無しで打つか、
F[x]とか試しに打ってみながらやってれば間違えないだろ?
あとは自分でやれ。
523:132人目の素数さん
03/07/11 14:51
Solve[]や、NSolve[]で、求めた複数の解のうち、
ひとつだけを使いたいのですが、
たとえばt = %[[-1]], t = Last[%]などとしても、
{x ->0.234hogehoge}と t に入ってしまいます。
この x->を取り除くにはどうすればいいのでしょうか?
524:523
03/07/11 15:00
x /. %
でできました。トンクス
525:132人目の素数さん
03/07/12 00:43
>521
えー、 Sum[<<>>,{i,1,100}] ってあるけど
数値やるなら NSum じゃ駄目なん?
どうせFindRoot なんだし
526:132人目の素数さん
03/07/14 23:27
あのー四捨五入するような関数ってないですか?切り捨てでもいいんですけど。
マニュアルで調べてみたんですが見つからなくて・・・・
527:132人目の素数さん
03/07/14 23:28
Roundじゃないの?
528:132人目の素数さん
03/07/14 23:36
0.5足して切り捨てろ。
529:山崎 渉
03/07/15 12:38
__∧_∧_
|( ^^ )| <寝るぽ(^^)
|\⌒⌒⌒\
\ |⌒⌒⌒~| 山崎渉
~ ̄ ̄ ̄ ̄
530:132人目の素数さん
03/07/16 23:06
URLリンク(shop.ponycanyon.co.jp)
531:132人目の素数さん
03/07/17 07:26
age
532:132人目の素数さん
03/07/17 14:11
>>530
これなら3,800円か。安いな。
533:132人目の素数さん
03/07/17 16:26
質問なんですが、例えば
x=k/10
y=k+1
For[k=1,k<11,k++,
r=x+y]
みたいに繰り返し計算をさせる場合に、このrは、毎回更新されてしまうんで
その結果を格納してリストの形にしたいとするとどーすればいいですか?
勿論、これぐらい単純だとTable使えばいいんですが、やりたいのはもっと
ややこしいんで、結果をリストにする格納する方法を教えてください。
よろしくお願いします。
534:132人目の素数さん
03/07/17 19:13
>>533
単純であろうがなかろうがTableでいいはずだ。
どうしてもというなら
l ={}
For[k=1,k<11,k++,r=x+y;l = Prepend[r,l]]
等。
しかしmathematicaではFor文を使う必要はない(その上
遅いので使うべきでもない)ので再考することを勧める。
535:132人目の素数さん
03/07/18 14:39
とぇ
536:132人目の素数さん
03/07/22 12:25
FindMinimumを使って最小値とそのときの変数の値を求めると、
{234,{a -> 3,b -> 2}}
みたいな感じだと思うんですが、これをいっぱい出して、a,bの平均値を出すには
どうすればいいでしょうか?
537:132人目の素数さん
03/07/23 01:12
>>546
漠然としすぎ。
こういうことか?
<< Statistics`DescriptiveStatistics`
data = {x, y} /. (#[[2]] &) /@
Flatten[Table[
FindMinimum[Sin[x + y], {x, i}, {y, j}], {i, 0, 3}, {j, 0, 3}],
1];
Mean /@ Transpose @ data
538:132人目の素数さん
03/07/23 01:12
>>537
s/546/537
539:537=538
03/07/23 01:13
536駄・・・まちがえすぎだ。
540:536
03/07/24 13:07
すみません、537のプログラム見て解読してたんですけど、ちょっと良くわかんなくて。
たとえば、
FindMinimum[f[a,b],{a,1},{b,2}]
とした場合、{最小値,{a -> 0.89898},{b->0.1}]
みたいにした場合、a->の部分を消すのはどうすればいいですか?
541:536
03/07/24 14:36
なんか読みにくいですね・・えーと置換なんかでうまくa->の部分を消せない
かなあと・・・無理なんですかねえ。。。
542:132人目の素数さん
03/07/24 14:37
開店セール見に来てーーーーーーーーーー!DVDを見ようー!GO!
▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼
開店セール1枚 500円より→→→URLリンク(www.dvd-exp.com)<)
アダルトDVDの価格はなーんと1枚900円→→→URLリンク(www.dvd-exp.com)<)
セルビデオ店の人気商品が24時間お買い物が出来ます。
見やすい画面・注文簡単・時間をかけてゆっくり選んでください!
URLリンク(www.dvd-exp.com)←へGO!
DVDエクスプレス
543:132人目の素数さん
03/07/24 15:18
>>541
簡単に消せるように"->"つきで出力されている。
s = {a -> 1, b -> 2}
{a,b} /. s
544:132人目の素数さん
03/07/24 22:08
s はリストみたいになっているから
s[[1]] や s[[2]] とかで取り出しても良いんじゃない?
545:132人目の素数さん
03/07/25 04:10
>>544
543で解決してます。”a->”みたいなのを取り除きたいってのが541の質問だろ。
ってか、房は黙ってろ(w
546:536
03/07/25 12:02
>>543
ありがとうございました。出来ました。
>>544
それだと、a -> 0.89898などと出ます。
547:132人目の素数さん
03/07/25 13:33
s[[2,2]]
とやるんじゃ
548:132人目の素数さん
03/07/25 14:11
【クレジットカードのショッピング枠を現金化!】
クレジットカードで買い物をすると、商品代金の80〜90%のキャッシュバック!!
10万円の買い物をすれば8〜9万円の現金がお手元に!!
■■朝日カード■■
URLリンク(www2.pekori.to)
549:132人目の素数さん
03/07/25 21:29
>>547
自分で一回やってみれば?違うのよくわかるから。
ほとんど使ったことないでしょ?分りもしないのにアドバイスしようとしない
方がいいじゃないのか?
550:132人目の素数さん
03/07/25 22:37
>549
残念ですが{a->1,b->2}[[2,2]] で 2 が出てくるよ。
[[]] は別に行列以外でも使える、数式構造の場所指定でした。
まあこのやり方は、-> の値を取り出すのにふさわしくないとは
僕も思いますがね。
551:132人目の素数さん
03/07/25 22:40
In[4]:=
{a -> 1, b -> 2} [[{1,2},2]]
Out[4]=
{1,2}
List[Rule[a, 1], Rule[b, 2]]
In[6]:=
FullForm[{{1,2},{a,b}}]
Out[6]//FullForm=
List[List[1,2],List[a,b]]
552:132人目の素数さん
03/07/25 22:52
>>550
質問の意図と違うから。その意味での違うなんだけど?なにがでしただよw
553:550
03/07/25 23:02
>552
確かに s[[2,2]] だと、本当に欲しい {1,2} は
出てきませんね。失礼スマソ。
ただ -> に対して、通常考える /. 以外の [[]] 指定方法に
目を付けた点は悪くはないのでは、とちょっと擁護
してみたかったのでした。
554:132人目の素数さん
03/07/26 10:30
>>553
> 目を付けた点は悪くはないのでは、とちょっと擁護
何がやりたいのかわかりにくい、バグに弱い、
一般的な状況(a^2 + aを求めるとか)では使えない…
センス悪すぎるよ。
555:132人目の素数さん
03/07/26 20:47
リストから第一行目、第十一行目、第二十一行目・・・というふうに取り出したい
んですけど、どーすればいいですか?
556:132人目の素数さん
03/07/26 22:46
>554
Mathematicaのルールどおりだからバグに弱くない
おまえがセンス悪、
それとも正確ワルか
557:132人目の素数さん
03/07/26 23:07
_,.. ---- .._
,. '" `丶、
/ ` 、
,..-‐/ ...: ,ィ ,.i .∧ , ヽ.
. ,:' .l .::;',. :::;/..://:: /,':/ ', l、 .i ヽ
. ,' ..::| .::;',' :;:','フ'7フ''7/ ',.ト',_|, , ',.',
,' .::::::!'''l/!:;'/ /'゙ / '! ゙;:|:、.|、| 'l
. ,'. .:::::::{ l'.l/ 、_ _,. 'l/',|.';|
l :::::::::::';、ヾ  ̄ `‐-‐'/! ';. '
. ! :::::::::::/ `‐、 ゝ |'゙ |
| ::::::::/ \ 、_, _.,.,_ ノ::: !
|::::/. _rl`': 、_ ///;ト,゙;:::::./ マターリと良スレ。
.. `´ /\\ `i;┬:////゙l゙l ヾ/ 数ヲタがうなぎを食べられますように・・・
,.:く::::::::`:、\ 〉l゙:l / !.|
. /:.:.:.:\:.:.:.:.`:、ソ/:.:| | |
/.:.:.:.:.:.:.:.:.:\:.:.:.:У:.:;l /./
. /:.:.:.:.:.:.:.r'´`‐,`、:/.,.:‐{ | !`:、
,'.:.:.:.:.:.:.:.:.';_,゚.,ノ.:./,:':.:.:.:', | |`、:|
!:.:.:.:.:.:.:.:.:.:.゙、:.::/:.:.:.:.:.:.ヽ, / ,!:.:`、
558:132人目の素数さん
03/07/27 00:52
鰻食いてーー
559:132人目の素数さん
03/07/27 00:53
でも貧乏だからとても買えねーよ(;´Д`)
560:132人目の素数さん
03/07/27 01:07
>>556
> Mathematicaのルールどおりだからバグに弱くない
あきれますな。悪いこと言わんからもっと謙虚になれ。
プログラムくんだ経験自体無いだろ?
たとえばFindRootで解が見付からなくてそのまま式自体が
かえったらどうなる?
[[]]ってのはかなりError proneだぞ。
次ページ最新レス表示スレッドの検索類似スレ一覧話題のニュースおまかせリスト▼オプションを表示暇つぶし2ch
5174日前に更新/264 KB
担当:undef