2011-03-11 39 views
0

我正在嘗試Mathematica編程,並認爲我會嘗試使用球形極座標手動進行手動計算(磁場相量從M類的E &)。我創建了一個變量,並試圖把它的捲曲:脅迫mathematica象徵性地評估一個球形極性捲曲表達式?

Needs["VectorAnalysis`"] 
SetCoordinates[Spherical] 
SetAttributes[ k, Constant ] 
eE := {0, 0, (Sin[Ttheta]/Rr) (1 - I/(k Rr)) e^{I k Rr}} 
Curl[ eE ] 

這不會實際評估衍生品像我想的那樣,只給:基本上

   I k Rr 
       e  (-I + k Rr) Sin[Ttheta] 
\[Curl]{0, 0, {-------------------------------}} 
         2 
        k Rr 

,它只是隨地吐痰我的輸入反饋給我作爲輸出。 SimplifyFullSimplify不會改變結果。

我猜的一個猜想是,這是因爲我沒有將k指定爲常量,所以我補充說(如上),但這並沒有什麼區別。

回答

3

更改e^{I k Rr}E^(I k Rr)

  1. {}指的載體,以及數學接受其所有的功能載體,但輸出向量,這是不是你想要的。例如,e^{1,2,3}變爲{e^1, e^2, e^3}。因此,您編寫表達式的方式在第一個列表中的第三個位置有一個元素列表,它將Mathematica關閉。

  2. 常數e是在Mathematica中的E

+0

謝謝。我天真地認爲輸入語法像乳膠一樣。 – 2011-03-12 03:19:32

2

兩個問題:

首先,在數學符號E,不e,是指數不斷ê

其次,你養E到列表的功率({...},又名List[...]),在這裏我想你,而不是意在使用的括號:

In[17]:= eE:={0,0,(Sin[Ttheta]/Rr) (1-I/(k Rr)) E^(I k Rr)} 
In[18]:= Curl[eE] 

Out[18]= {(2 E^(I k Rr) (1-I/(k Rr)) Cos[Ttheta])/Rr^2, (Csc[Ttheta] 
    (-I E^(I k Rr) k (1-I/(k Rr)) Sin[Ttheta]^2-(I E^(I k Rr) 
    Sin[Ttheta]^2)/(k Rr^2)))/Rr,0} 

HTH!

+0

繪製該字段的最簡單方法是什麼? – 2011-03-11 12:12:53

+0

我想['VectorPlot3D'](http://reference.wolfram.com/mathematica/ref/VectorPlot3D.html),座標變換是一種選擇:'k = 10; VectorPlot3D [curl /。 Thread [{Rr,Ttheta,Pphi} - > CoordinatesFromCartesian [{x,y,z}]],{x,-1,1},{y,-1,1},{z,-1,1}] ' – 2011-03-12 04:53:27

+0

工程就像一個魅力。從未使用過這樣的座標... []。看起來非常有用。謝謝。 – 2011-03-12 05:28:03