| 24小時(shí)熱門(mén)版塊排行榜 |
| 10 | 1/1 | 返回列表 |
| 查看: 1586 | 回復(fù): 9 | |||
亮亮1283新蟲(chóng) (小有名氣)
|
[交流]
【求助】請(qǐng)教Mathematica程序問(wèn)題 已有5人參與
|
|
為什么我的程序?qū)σ粋(gè)是函數(shù)求導(dǎo)后,得到的函數(shù)帶入數(shù)值后竟然是復(fù)數(shù)! 以下是程序的一部分,M是一個(gè)2*2矩陣,矩陣中每一項(xiàng)均是Omega的函數(shù),而S1,S2是關(guān)于omega的復(fù)函數(shù) S1 = M[[1, 1]]; S2 = M[[2, 1]]; \[Phi]1 = Arg[S1]; \[Phi]2 = Arg[S2]; \[Tau]1 = D[\[Phi]1, \[Omega]]; 得到的Tau代入數(shù)值后卻是復(fù)數(shù)!不知道怎么回事!請(qǐng)好心人幫忙! |
木蟲(chóng) (正式寫(xiě)手)
獨(dú)自為政

至尊木蟲(chóng) (正式寫(xiě)手)
至尊木蟲(chóng) (正式寫(xiě)手)
新蟲(chóng) (小有名氣)
銅蟲(chóng) (正式寫(xiě)手)
至尊木蟲(chóng) (正式寫(xiě)手)
|
我的意思是說(shuō),先對(duì)Arg求導(dǎo)再代入數(shù)值,和,先對(duì)變量賦值再去求Arg的導(dǎo)數(shù),兩個(gè)結(jié)果是不一樣的 估計(jì)這是mathematica不分?jǐn)?shù)據(jù)類(lèi)型導(dǎo)致的(沒(méi)有經(jīng)過(guò)驗(yàn)證) 你可以看看后一種計(jì)算方法是否違背你的本意,若否,那么問(wèn)題也就解決了 |
|
其實(shí)先賦值和后賦值結(jié)果是一樣的,你給的例子4樓里面少了前面的-3 (x^2) ,加上這一項(xiàng)結(jié)果就相同了 D[Arg[-x^3 + 8 I], x] N[% /. x -> 3 + 4 I] x = 3 + 4 I; N[-3 x^2* \!\(\*SuperscriptBox["Arg", "\[Prime]", MultilineFunction->None]\)[-x^3 + 8 I]] Clear[x]; |
至尊木蟲(chóng) (正式寫(xiě)手)
|
好吧,我再解釋一下這個(gè)問(wèn)題 第一個(gè)程序(后賦值的)是按照樓主的程序思想做的。先對(duì)Arg函數(shù)求導(dǎo),結(jié)果前面有了一項(xiàng)-3 (x^2)。所以當(dāng)x是復(fù)數(shù)時(shí),最終結(jié)果當(dāng)然就是復(fù)數(shù)。 第二個(gè)程序(先賦值的)就不用多解釋了,自然是先帶入x的數(shù)值然后求導(dǎo),所以看起來(lái)是Arg‘[-x^3+8I],因此-3 (x^2) 這一項(xiàng)不是隨便加上的。最終結(jié)果當(dāng)然是實(shí)數(shù)了。不過(guò)我沒(méi)有考慮到的是mathematica對(duì)Arg‘[-x^3+8I]的計(jì)算結(jié)果等同于Arg[-x^3+8I],失策失策!見(jiàn)笑了。 那么對(duì)樓主的問(wèn)題再說(shuō)幾句。1、Arg[]給出結(jié)果應(yīng)當(dāng)是實(shí)數(shù),起碼我沒(méi)有想出來(lái)什么情況這個(gè)函數(shù)會(huì)給出復(fù)數(shù)結(jié)果。2、對(duì)Arg求導(dǎo)會(huì)多出來(lái)一些表達(dá)式,建議樓主仔細(xì)看看這些表達(dá)式的結(jié)果是否是復(fù)數(shù)。3、如果樓主珍惜大家的唾沫星子和腦細(xì)胞并且不違反保密的話,還是把S1等的表達(dá)式貼出來(lái)更好一些。 |
| 10 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|