| 查看: 261 | 回复: 2 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】请帮忙看看问题出在那里了 已有1人参与
|
|||
|
clear; clc; L=0.632; w0=400; k=2*pi/L; z=linspace(0,50,100); r=linspace(0,0.1*w0,200); eta=linspace(0,2*pi,200); [rho,theta]=meshgrid(r,eta); [x,y]=pol2cart(theta,rho); E=1./(L.*z).*(k*w0./z).^4.*besselj(4,k*w0.*rho./z)./rho.*exp(1i*(3*theta+k.*(rho.^2/(2.*z)+z))); I=(1/(L*z)*(k*w0/z)^4*besselj(4,k*w0.*rho/z)./rho).^2; p=patch(isosurface(x,y,z,E,0)); set(p, 'FaceColor', 'red', 'EdgeColor', 'none'); daspect([1 1 1]) view(3) camlight; lighting phong 结果提醒： Error using ==> rdivide Matrix dimensions must agree. 请帮忙，谢谢！ |
回复此楼» 猜你喜欢
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
基金中了
已经有14人回复
-
关于filecode
已经有7人回复
-
娱乐
已经有6人回复
-
化学口download_prp&fileCode的固定段好像这几天一直没变，有变的大神么？
已经有4人回复
-
关于代码变化问题，想知道的进来
已经有12人回复
-
固定端突然变了，今天
已经有10人回复
-
听说今天filecode变了
已经有47人回复
lijinfeng042
木虫 (小有名气)
Matlab
- 仿真EPI: 2
- 应助: 1 (幼儿园)
- 金币: 2156.1
- 散金: 115
- 帖子: 291
- 在线: 31.5小时
- 虫号: 1019062
- 注册: 2010-05-15
- 性别: GG
- 专业: 通信理论与系统
★ ★
adu886886(金币+2):鼓励交流 2010-06-03 08:27:40
shenailin(金币+5):谢谢，高手！ 2010-06-03 15:23:44
adu886886(金币+2):鼓励交流 2010-06-03 08:27:40
shenailin(金币+5):谢谢，高手！ 2010-06-03 15:23:44
|
令z=k*w0.*rho 运算后 whos可知 whos Name Size Bytes Class Attributes E 1x100 800 double L 1x1 8 double ans 200x200 320000 double eta 1x200 1600 double k 1x1 8 double r 1x200 1600 double rho 200x200 320000 double s 200x200 320000 double theta 200x200 320000 double w0 1x1 8 double x 200x200 320000 double y 200x200 320000 double z 1x100 800 double z是一个1*100矩阵 而 s是200*200矩阵 形不同 ./(点除)当然出错了啊 |
3楼2010-06-02 20:29:58
lijinfeng042
木虫 (小有名气)
Matlab
- 仿真EPI: 2
- 应助: 1 (幼儿园)
- 金币: 2156.1
- 散金: 115
- 帖子: 291
- 在线: 31.5小时
- 虫号: 1019062
- 注册: 2010-05-15
- 性别: GG
- 专业: 通信理论与系统
2楼2010-06-02 18:03:43
100