Mathematica 问答社区 - 矩阵中最新问题

嵌套函数如何输入矩阵进行计算

Tue, 07 Apr 2020 02:08:34 +0000

lsEigenvalue[H_] := Module[{emax, emaxspace, emin, eminspace, es},

  es = Sort[Eigenvalues[H], Re[#1] < Re[#2] &];

  emin = First[es];

  emax = Last[es];

  eminspace =

   Orthogonalize[NullSpace[H - emin*IdentityMatrix[Length[H]]];

    emaxspace =

     Orthogonalize[NullSpace[H - emax*IdentityMatrix[Length[H]]];

        Return[{{emin, eminspace}, {emax, emaxspace}}]

        ]

       plotNR[A_, n_] :=

      Module[{t = 0., td = 2 \[Pi]/n, Ht, Kt, points = {},

        segment = {} , emax , emaxspace, emin, eminspace, vp, vm, Q,

        R},

       PrintTemporary[ProgressIndicator[Dynamic[t], {0, 2 \[Pi]}]];

       While[t < 2 \[Pi],

        Ht = (Exp[-I t]*A + Exp[I t]*ConjugateTranspose[A])/2;

        {emax, emaxspace} = Last[lsEigenvalue[Ht]];

        which[(* One dimensional eigenspace *)

         Length[emaxspace] == 1,



         vp = First[emaxspace];

         AppendTo[points, Conjugate[vp].A.vp],

         Length[emaxspace] > 1,

         Kt = (Exp[-I t]*A - Exp[I t]*ConjugateTranspose[A])/(2 I);

         Q = Transpose[emaxspace];

         R = ConjugateTranspose[Q].Kt.Q;



         {{emin, eminspace}, {emax, emaxspace}} = lsEigenvalue[R];



         vp = Q.First[emaxspace];

         vm = Q.First[eminspace];



         AppendTo[segments, {Conjugate[vm].A.vm, Conjugate[vp].A.vp}],



         (* Fail *)

         Ture,

         Print["Error"]

         ];



        t = t + td;

        ];

       Return {DeleteDuplicates[points], DeleteDuplicates[segments]}]

    ]

此代码如何输入矩阵进行计算？

怎么在编号的代码中输入矩阵进行计算？

Tue, 07 Apr 2020 02:08:33 +0000

为何matheamtica计算的不对啊

Thu, 10 Oct 2019 03:29:56 +0000

问题1：MaximalBy[{{4.9375, {x -> -(3/4)}}, {2.0792, {x -> 2/Sqrt[

     3]}}}, First]

结果为：{{4.9375, {x -> -(3/4)}}}

问题2:MaximalBy[{{79/

   16, {x -> -(3/4)}}, {-1 + 16/(

    3 Sqrt[3]), {x -> 2/Sqrt[3]}}}, First]

结果为：{{-1 + 16/(3 Sqrt[3]), {x -> 2/Sqrt[3]}}}

请问大家，两个问题一模一样，为啥计算的结果不一样呢？请大家计算看看是哪里出来问题啊。多谢了！

请证明第二个等号，貌似用到泰勒展开且展开到一阶，其中的sigma是PauliMatrix。

Sat, 16 Mar 2019 03:55:21 +0000

泰勒展开到一阶，其中的sigma是PauliMatrix

矩阵 Do循环提示"Objects of unequal length in {Null}\ {4,1,4,1} cannot be combined. "

Tue, 21 Mar 2017 01:43:06 +0000

（*1、提示标题中的错误后求Mf[[3]]也可以得到一个结果，但是似乎数值不对。 2、最后还要以Mi[[i]]为横坐标，Mf[[i]]为纵坐标作图拟合，这一步不知道应该用什么命令。 3、因为对mma了解很少，所以为了减少可能的报错，语句写的复杂，希望懂的前辈指教一二。 谢谢！*)

ClearAll;
Mi=Table[1,44];
g=7.25;
z13=2*10^9;
zi=10^10;
Mf=Table[1,44];
Do[{mi=10^(i+17);Mi[[i]]=mi;A=NDSolveValue[{M'[a]==M[a]^2*(2*10^(-68)+7*10^(-59)*g/a^3,M[1/(1+zi)]==mi},M,{a,(1/(1+zi),1/(1+z13)}];mf=A[1/(1+z13)]-mi;Mf[[i]]=mf},{i,1,44,1};
Plot

大规模系数矩阵的生成

Fri, 10 Feb 2017 11:55:12 +0000

我的问题是这样的，k[x_, t_] := t - x;
f[x_] := x;
t0 = 0.;
tn = 4.;
n = 400;
h = (tn - t0)/n;
b = f[#] & /@ (h*Range[0, n]);
m = SparseArray[{{i_, j_} /; j < i && j != 1 -> -h*
k[h*i, h*j], {i_, 1} /; i != 1 -> -0.5*h*k[h*i, 0], {i_, i_} ->
1 - 0.5*h*k[h*i, h*i]}, {n + 1, n + 1}];
result1 = LinearSolve[m, b];

矩阵的循环相乘有直接的命令吗？

Thu, 09 Feb 2017 11:04:09 +0000

数的循环相乘有命令，如Product[x+i,{i,1,3}]得到(x+1)(x+2)(x+3)。我想要一个结果如a.b.a.b^2.a.b^3.......的矩阵连乘，有直接的命令吗？如果没有怎样写成一个函数？得到f(m)=:a.b.a.b^2.........a.b^m。请大侠们帮忙啊！！！

怎样把Eigenvectors算出的本征矢归一化

Fri, 19 Aug 2016 02:47:03 +0000

怎样把Eigenvectors算出的本征矢归一化

In[1]:= Eigenvectors[{{a, b}, {c, d}}]

Out[1]= {{-((-a + d + Sqrt[a^2 + 4 b c - 2 a d + d^2])/(2 c)),
1}, {-((-a + d - Sqrt[a^2 + 4 b c - 2 a d + d^2])/(2 c)), 1}}

为什么这个函数没有输出？

Thu, 07 Jul 2016 07:50:19 +0000

ComplexContinuedFraction[x_] := Catch[
  For[k = {}; i = 0,
   i < -1,
   i++,
   If[GIRCF[x, i][[2]] === Indeterminate,
    Append[k, Indeterminate]; Throw[k],
    (*如果余值出现未定义，将未定义加入到列表中并输出结果*)
    If[GIRCF[x, i][[2]] === 0, Throw[k], 
     Append[k, GIRCF[x, i][[1]]]]]]]
(*如果余值未出现0，GIRCF[x,i]\[LeftDoubleBracket]1\[RightDoubleBracket]\
加入到列表直到余值为0并输出结果*)

(*可行性测试*)
ComplexContinuedFraction[466611/200 + (1332201 I)/2000]

以上是没有输出的函数。

==========================================================================

前置代码：

GaussianIntegerRegularContinuedFraction[fn___] := GIRCF[fn];
GIRCF[x_, dx_, k_Integer] := GIRCF`Core[Rationalize[x, dx], k];
GIRCF[x_, k_Integer] := GIRCF`Core[Rationalize[x], k];
(*反正最后都是按有理数处理就不需要数值近似.*)
GIRCF`Core[x_, k_] := Block[{GIRCF`m = GIRCF`Core[x, k - 1][[2]]},
    If[FractionalPart[GIRCF`m] === 0, {GIRCF`m, 0},
    (*余值为0意味着规范连分数化的结束*)

     
     If[IntegerPart[GIRCF`m] === 0,
      {GIRCF`m, Indeterminate},
      (*并非每一个复数都可以在高斯整数下被规范连分数化*)
      {IntegerPart[GIRCF`m], 1/FractionalPart[GIRCF`m]}]]] /; k >= 1;
GIRCF`Core[x_, 0] :=
  If[FractionalPart[x] === 0, {x, 0}, {IntegerPart[x], 1/
    FractionalPart[x]}];
GIRCF`Core[x_, k_] := {IntegerPart[x*10^k], IntegerPart[x*10^k] - x} /;
    k <= -1;



(*可行性测试1*)
a = Rationalize[2333.055 + 666.1005 I];
Do[Print[GIRCF`Core[a, n]], {n, 9}]

(*可行性测试2*)
a = Rationalize[2333.12345 + 666.54321 I];
Do[Print[GIRCF`Core[a, n]], {n, 8}]

Mathematica 问答社区 - 矩阵 中最新问题