Question

MATLAB程序问题

clear all
%Constants (all MKS, except energy which is in eV)
hbar=1.055e-34;m=9.110e-31;epsil=8.854e-12;q=1.602e-19;

%Lattice
Np=100;a=1e-10;X=a*[1:1:Np];t0=(hbar^2)/(2*m*(a^2))/q;L=a*(Np+1);
T=(2*t0*diag(ones(1,Np)))-(t0*diag(ones(1,Np-1),1))-(t0*diag(ones(1,Np-1),-1));
x=a*[1:0.001:Np]
[V,D]=eig(T);D=diag(D);[Enum,ind]=sort(D);

E1=D(ind(2));psi1=abs(V(:,ind(2)));P1=psi1.*conj(psi1);
psi2=((2/L)^(1/2))*sin(2*pi*x/L);P2=psi2.*conj(psi2);%这句有问题，我想画出这个函数的图像，出不来

hold on
h=plot(P1,'b');% Part (b)
h1=plot(P2,'b');% Part (b)
set(h,'linewidth',[3.0])
set(h1,'linewidth',[1.0])
set(gca,'Fontsize',[2])

xlabel(' Lattice site # --->');% Part (b)
ylabel(' Probability ---> ');% Part (b)
grid on

Answer 1

不知道你的程序是要什么样的结果，单纯从代码来看，就是全角和半角的问题。我已经修改如下：
建立M文件eq1．m
function dy=eq1(x,y)
dy=zeros(2,1);
dy(1)=y(2);
dy(2)=1/5*sqrt(1+y(1)^2)/(1-x);
取x0=0，xf=0．9999，建立主程序如下：
x0=0;
xf=0.9999;
[x,y]=ode15s('eq1',[x0 xf],[0 0]);
plot(x,y(:,1),'b.')
hold on
y=0:0.01:2;
plot(1,y,'b*')
得到的图像是

追问

额，看不懂你说的什么。我在这方面是个空白。出问题的这句是一个函数的表达式。得到的图像是曲线。

Answer 2