版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

卡卡爱物理

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 189.5
帖子: 57
在线: 6.7小时
虫号: 5350362
注册: 2016-12-13
专业: 凝聚态物性 II ：电子结构

[求助] 求帮找找小错误已有1人参与

x=(-10:0.01:10)*w0;z=0:0.01:10;y=5;%zp*0.5*pi;[x,z]=meshgrid(x,z);a=A0*kw/zp./(kw./zp.*cos(z/zp)+2*1i.*sin(z/zp));b=exp(omg^2*w0*w0.*tan(z/zp)./(tan(z/zp)-i*kw./zp));c=cosh(kw./zp*omg*x./(kw./zp.*(cos(z/zp))+2*i.*sin(z./zp)));d=cosh(kw./zp*omg*y./(kw./zp.*(cos(z/zp))+2*i*sin(z./zp)));e=exp((kw^2*gamma1^2*p0.*tan(z/zp)-2*i*kw./zp).*(x.^2+y.^2)... /w0*w0./(2*i*kw./zp-4*tan(z/zp)));phi=a.*b.*c.*d.*e;mesh(x,z,abs(phi.^2))view(2)%shading interp
想得到这样的图

发自小木虫Android客户端

回复此楼

» 猜你喜欢

1楼 2017-04-19 10:23:01

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

专家经验: +1059
应助: 1712 (讲师)
贵宾: 8.888
金币: 68161.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3292.2小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
jjdg: 金币+2, 辛苦了 2017-04-19 23:08:56
卡卡爱物理(jjdg代发): 金币+20, 代发金币 2017-04-28 00:42:21
jjdg: 金币-30, 求助者不兑现承诺的，按1.5-20倍扣除 2017-04-28 00:43:04
jjdg: 金币+30, 不好意思，点错了 2017-04-28 00:43:43

没有定义的变量，都按等于1来赋值，否则原代码无法运行，你根据实际数值自行修改吧。

CODE:

clear;clc

w0=1;

zp=1;

A0=1;

kw=1;

omg=1;

gamma1=1;

p0=1;

x=(-10:0.01:10)*w0;

z=0:0.01:10;

y=5;%zp*0.5*pi;

[x,z]=meshgrid(x,z);

a=A0*kw/zp./(kw./zp.*cos(z/zp)+2*1i.*sin(z/zp));

b=exp(omg^2*w0*w0.*tan(z/zp)./(tan(z/zp)-1i*kw./zp));

c=cosh(kw./zp*omg*x./(kw./zp.*(cos(z/zp))+2*1i.*sin(z./zp)));

d=cosh(kw./zp*omg*y./(kw./zp.*(cos(z/zp))+2*1i*sin(z./zp)));

e=exp((kw^2*gamma1^2*p0.*tan(z/zp)-2*1i*kw./zp).*(x.^2+y.^2)...

    /w0*w0./(2*1i*kw./zp-4*tan(z/zp)));

phi=a.*b.*c.*d.*e;

mesh(x,z,abs(phi.^2))

view(2)%shading interp

untitled.png

赞一下

回复此楼

MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。

2楼2017-04-19 14:53:18