該樓層疑似違規已被系統折疊?隱藏此樓查看此樓

有了這個，算一些常用張量就可以方便很多了。

這里給出的程序中的默認情況是計算帶點球對稱靜態解的，可以直接給出二階協變度規(就是常用的那個度規的倒數)，Christoff聯絡、Ricci張量與Ricci標量。

別的形式的度規嘛，改改函數就好了。

syms?M;

syms?u;

syms?sita;

syms?t;

syms?phi;

syms?Gamma;

syms?Ricci;

syms?R;

syms?Q;

r=u;

delta(1,1)=1;

delta(2,2)=1;

delta(3,3)=1;

delta(4,4)=1;

v(1)=t;

v(2)=u;

v(3)=sita;

v(4)=phi;

g(1,1)=simplify(-(1-2*M/r+Q/r/r));

g(2,2)=simplify(1/(1-2*M/r+Q/r/r)*diff(r,u)^2);

g(3,3)=simplify(r^2);

g(4,4)=simplify(r^2*cos(sita)^2);

gb=delta/g

for?i=1:4

for?j=1:4

for?k=1:4

tmp=0;

for?l=1:4

tmp=tmp+gb(k,l)/2*(diff(g(l,i),v(j))+diff(g(l,j),v(i))-diff(g(i,j),v(l)));

end

Gamma(i,j,k)=simplify(tmp);

end

end

end

Gamma

for?i=1:4

for?j=1:4

tmp=0;

for?k=1:4

tp=0;

for?l=1:4

tp=tp+Gamma(i,j,l)*Gamma(l,k,k)-Gamma(i,k,l)*Gamma(j,l,k);

end

tmp=tmp+diff(Gamma(i,j,k),v(k))-diff(Gamma(i,k,k),v(j))+tp;

end

Ricci(i,j)=simplify(tmp);

end

end

Ricci

R=0;

for?i=1:4

for?j=1:4

R=R+Ricci(i,j)*gb(i,j);

end

end

R=simplify(R)

總結

以上是生活随笔為你收集整理的矩张量计算 matlab,用来计算广义相对论常用张量计算的Matlab程序的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。