Question

我有一个莫尔斯电势，可以通过scilab求解并找到其能量本征值和向量。但是我遇到了一个问题，因为scilab表示子矩阵定义不正确。

我试图了解为什么我的矩阵无法正确。我是这个平台的新手。

De=4.618;

b=18.118;

x=0.1275;

a=0.35;

hc=197.327;

mu=912847000;

N=45;

c=(De/(2*a*b));

d=exp(2*b*(x-a));

e=exp(2*b*x);

i=exp(b*x);

f=(2*De/(a*b));

g=((a*b)^(2));

h=(De*a*b/(2));

  h=zeros(N,N);

  s=eye(N,N);

  for n=1:N

      h(n,n)=(((n*%pi*hc)^2)/(2*mu*a^(2)))+De-(c*d)+(c*e)+(f*d)-(f*e)-((h*(e-d))/(g+(%pi^(2)*n^(2))))+4*((h*(e-d))/(g+(4*(%pi^(2)*n^(2)))));

      for m=n+1:N

          h(m,n)=4*h*(1-(exp(-2*b*a)*((-1)^(m+n))))*e*((1/(4*g+(((m-n)^(2))*%pi^2)))-(1/(4*g+(((m+n)^(2))*%pi^2))))+4*h*i*(1-(exp(-1*b*a)*((-1)^(m+n))))*e*((1/(g+(((m+n)^(2))*%pi^2)))-(1/(g+(((m-n)^(2))*%pi^2))))

          h(n,m)=h(m,n);

          end

  end

  [al,bl,R]=spec(h,s);

  el=al./bl;

  e=R;

  [el,k]=gsort(el)

  disp(h);

  disp(el)

Answer 1

您将名称h用于两个不同的变量（标量和矩阵），这就是计算失败的原因。顺便说一句，您正在使用许多无用的括号（这不利于阅读您的代码）。每次引用矩阵并在其中定义行之后，都可以用大写的h来替换H，

  H=zeros(N,N);
  s=eye(N,N);

  for n=1:N

      H(n,n)=(((n*%pi*hc)^2)/(2*mu*a^(2)))+De-(c*d)+(c*e)+(f*d)-(f*e)-((h*(e-d))/(g+(%pi^(2)*n^(2))))+4*((h*(e-d))/(g+(4*(%pi^(2)*n^(2)))));

      for m=n+1:N

          H(m,n)=4*h*(1-(exp(-2*b*a)*((-1)^(m+n))))*e*((1/(4*g+(((m-n)^(2))*%pi^2)))-(1/(4*g+(((m+n)^(2))*%pi^2))))+4*h*i*(1-(exp(-1*b*a)*((-1)^(m+n))))*e*((1/(g+(((m+n)^(2))*%pi^2)))-(1/(g+(((m-n)^(2))*%pi^2))))

          H(n,m)=H(m,n);

      end

  end

  [al,bl,R]=spec(H,s);
  el=al./bl;
  e=R;
  [el,k]=gsort(el)

  disp(H);

  disp(el)

在Scilab中错误定义了子矩阵

1 个答案: