24小时热门版块排行榜    

CyRhmU.jpeg
南方科技大学公共卫生及应急管理学院2025级博士研究生招生报考通知
查看: 1079  |  回复: 4

zhangpengju

新虫 (小有名气)

[求助] Matlab中的Matrix dimensions must agree.错误求助

命令文件
figure(3);
x=-15:15;
y=x;
[X,Y]=meshgrid(x,y);
z=zeros(size(X));
psi=psi_3dxy(X,Y,Z);————error
psi2=psi.*psi;
subplot(2,2,1);
mesh(X,Y,psi);
axis tight;
xlabel('x'),ylabel('y'),zlabel('3dxy');
subplot(2,2,2);
mesh(X,Y,psi2);
axis tight;
xlabel('x'),ylabel('y'),zlabel('13dxy1^2');
subplot(2,2,3);
contour(psi2,10);
axis tight;
xlabel('x'),ylabel('y'),title('13dxy1^2');
z=x;
[X,Y,Z]=meshgrid(x,y,z);
psi=psi_3dxy(X,Y,Z);
psi2=psi.*psi;
subplot(2,2,4);
slice(X,Y,Z,psi2,[],[],[0]);
axis tight;
xlabel('x'),ylabel('y'),title('13dxy1^2');
view(0,90);
colorbar;
对应的函数文件
function y=psi_3dxy(X,Y,Z);
r=sqrt(X.*X+Y.*Y+Z.*Z);
u=sqrt(X.*X+Y.*Y);
sin_theta=u./r;
sin_phi=Y./u;
cos_phi=X./u;
y=1/81/sqrt(2*pi).*r.*r.*exp(-r/3).*sin_theta.*sin_theta*2.*sin_phi.*cos_phi;
具体原因找不到啊,哪位大虾指点一下,谢啦!
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

做我开始做的!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ruby0531

新虫 (小有名气)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
zhangpengju(金币+1): ★★★很有帮助 2012-03-02 18:05:54
jjdg(金币+1): 感谢参与 2012-03-03 02:30:22
psi=psi_3dxy(X,Y,Z);
psi 和psi_3dxy的维数不一样吧,所以不能这样赋值?
2楼2012-03-02 15:17:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lurencyj

木虫 (著名写手)

【答案】应助回帖


感谢参与,应助指数 +1
zhangpengju: 回帖置顶 2012-03-02 18:05:31
jjdg(金币+1): 感谢参与 2012-03-03 02:30:16
主程序中,第5行的z=ones(size(X));中的z改为大写字母试试。呵呵

[ 发自手机版 http://muchong.com/3g ]
很女子很弓虽大
3楼2012-03-02 15:26:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangpengju

新虫 (小有名气)

引用回帖:
: Originally posted by lurencyj at 2012-03-02 15:26:39:
主程序中,第5行的z=ones(size(X));中的z改为大写字母试试。呵呵

谢谢,主程序和函数部分的大小写弄乱了,现在好了。
做我开始做的!
4楼2012-03-02 16:23:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lurencyj

木虫 (著名写手)

引用回帖:
4楼: Originally posted by zhangpengju at 2012-03-02 16:23:07:
谢谢,主程序和函数部分的大小写弄乱了,现在好了。

不用客气,给金币就好。哇哈哈!!!
很女子很弓虽大
5楼2012-03-02 16:29:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zhangpengju 的主题更新
信息提示
请填处理意见