查看: 443|回复: 0|关注: 0
打印 上一主题 下一主题

[未答复] Matlab模拟氢原子电子云,显示:错误使用 symengine Array sizes must match.?

[复制链接]

新手

5 麦片

财富积分


050


2

主题

2

帖子

0

最佳答案
跳转到指定楼层
1#
发表于 2019-4-25 14:59:52 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式
a0=1;     %归一化的波尔半径
syms x r;n=3,l=2,m=1;
N=-1*sqrt((2/n/a0)^3*factorial(n-l-1)/2/n/(factorial(n+l))^3);
t=laguerreL(n,l);
q=subs(t,x,2*r/n/a0);
[y,z]=meshgrid([-25:0.33:25]);
k=sqrt(y.^2+z.^2);
f=(N*exp(-1*r/n/a0))*(2*r/n/a0)^l*q;
f1=subs(f,r,k);    %径向波函数
theta=z.*(k.^(-1));
M=(-1)^m*sqrt((2*l+1)*factorial(l-abs(m))/4/pi/factorial(l+abs(m)));
k1=legendre(l,m);
L1=subs(k1,x,cos(theta));
theta2=pi/2-theta;
w=M.*subs(k1,x,cos(pi/2-theta2));   %球谐波函数
W=f1.*w;    %电子在整个空间的波函数
surf(z,y,abs(W).^2);   %画图输出电子在整个空间的彩色立体曲面电子云图
shading interp;
ylabel('y/a0');
xlabel('z/a0');
title('氢原子电子云');


您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

站长推荐上一条 /3 下一条

快速回复 返回顶部 返回列表