Question

N=15//no. of basis

b=%pi^2/4;

e2=30^2;//sq. of stiffness costant

a=2;//width of inf. potential well

h=zeros(N,N);//call h matrix

//hnm matrix

for n=1:N;

    h(n,n)=e2*((1-6/(%pi*%pi*n*n))*(%pi^2/48))+n^2;

    for m=n+1:N;

        q=1/((m-n)^2);

        r=1/((m+n)^2);

        t=((-1)^(m+n)+1)/4;

        h(m,n)=e2*t*(q-r);

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

        end

end

disp(h);

//to find eigen values and matrix

s=eye(N,N);//elementary matrix

disp(s);

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

e1=al./be;

c=R;

disp(e1);//eigen values

disp(c);

[E,k]=gsort(real(e1),"g","d");//sorting eigen values

E=E;

disp(E);

d=c;

d(:,1:N)=c(:,k);//sorting eigen vectors

p=1000;

x=linspace(0,a,p);

psi=zeros(N,p);//call psi wave fn.

for m2=1:N;

    for n2=1:N;

        psi(m2,:)=psi(m2,:)+d(n2,m2)*sqrt(2/a)*sin(n2*%pi*x/a);

    end

    if(m2<=15);

        subplot(5,5,m2);

        plot(x,psi(m2,:));

        end

end

我的问题是这是一个谐振子问题，我选择了无限方井的基础，并绘制了其方波函数。我的问题是随着我的基础增加到20（即N = 20），波形函数图只是倒转了，这是不应该发生的。所以我想知道代码中是否存在任何问题，因为物理不可能出错。

波函数图随着基础数的增加而变化

0 个答案: