24小时热门版块排行榜    

查看: 1160  |  回复: 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的回帖
回帖置顶 ( 共有1个 )

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的回帖
普通回帖

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的回帖

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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +9 yunziaaaaa 2026-03-01 10/500 2026-03-02 15:18 by 半截的诗0927
[考博] 诚招农业博士 +3 心欣向荣 2026-02-28 3/150 2026-03-02 13:33 by 时间不狗
[考研] 292求调剂 +7 yhk_819 2026-02-28 7/350 2026-03-02 12:43 by 无际的草原
[考研] 求调剂 +3 熬夜的猫头鹰 2026-03-02 3/150 2026-03-02 11:45 by 刘兵
[考研] 276求调剂 +4 路lyh123 2026-02-28 5/250 2026-03-02 11:20 by yuchj
[考研] 284求调剂 +10 天下熯 2026-02-28 11/550 2026-03-02 11:03 by 无际的草原
[考研] 264求调剂 +4 巴拉巴拉根556 2026-02-28 4/200 2026-03-02 10:48 by yuchj
[考研] 279求调剂 +3 dua1 2026-03-01 4/200 2026-03-02 00:23 by 大脸蛋子
[基金申请] 成果系统访问量大,请一小时后再尝试。---NSFC啥时候好哦,已经两天这样了 +4 NSFC2026我来了 2026-02-28 4/200 2026-03-01 22:37 by 铁门栓
[考研] 一志愿中南大学理学化学 +4 15779376950 2026-03-01 5/250 2026-03-01 19:00 by Fff-1
[考研] 321求调剂一志愿东北林业大学材料与化工英二数二 +4 虫虫虫虫虫7 2026-03-01 7/350 2026-03-01 16:52 by caszguilin
[考研] 0856材料求调剂 +4 麻辣鱿鱼 2026-02-28 4/200 2026-03-01 16:51 by caszguilin
[考研] 285求调剂 +8 满头大汗的学生 2026-02-28 8/400 2026-03-01 16:47 by caszguilin
[基金申请] 刚录用,没有期刊号,但是在线可看的论文可以放为代表作吗 10+3 arang1 2026-03-01 3/150 2026-03-01 16:43 by babero
[考研] 313求调剂 +3 水流年lc 2026-02-28 3/150 2026-03-01 16:01 by 新能源达人
[考研] 298求调剂 +9 人间唯你是清欢 2026-02-28 12/600 2026-03-01 14:23 by Ducount.Y
[考研] 调剂 +3 简木ChuFront 2026-02-28 3/150 2026-03-01 11:46 by 王伟要上岸啊
[考博] 博士自荐 +4 kkluvs 2026-02-28 4/200 2026-03-01 10:19 by 馥安馥安
[考研] 304求调剂 +3 52hz~~ 2026-02-28 5/250 2026-03-01 00:00 by 52hz~~
[高分子] 求环氧树脂研发1名 +3 孙xc 2026-02-25 11/550 2026-02-28 16:57 by ichall
信息提示
请填处理意见