查看: 0  |  回复: 0
【悬赏金币】回答本帖问题,作者卡卡艾物理将赠送您 8 个金币

卡卡艾物理

新虫 (初入文坛)

[求助] 已知场分布,如何求出他的坡印廷,也就是二维能流分布

clear all
beta=0.05;
lambda=632.8*10.^-6;
k=2*pi/lambda;
A0=1;
w1=0.1;
w2=0.1;
w0=0.1;
a=0.1;
b=0.1;
c=0.5;
z=0.1;
w=2.*pi./(beta);
aa=cos(z./beta);
bb=sin(z./beta).*beta;
cc=-sin(z./beta)./beta;
dd=cos(z./beta);

[x,y] = meshgrid(-1:.005:1);
M=i.*c./w0.^2+i.*k.*aa./2./bb;
A=-(i.*k.*dd./2./bb+k.^2./4./bb.^2./M).*(x.^2+y.^2)...
   +i*k./8./bb./M.^2.*(x./w1.^3+y./w2.^3)...
   +i*k./2./bb./M.*(a.*x./w1+b.*y./w2)...
   +1./96./M.^3.*(1./w1.^6+1./w2.^6)...
   +1./8./M.^2.*(a./w1.^4+b./w2.^4)...
   +1./4./M.*(a.^2./w1.^2+b.^2./w2.^2);
B=i.*k.*x./2./bb./M./w1+a./2./M./w1.^2+1./16./M.^2./w1.^4;
C=i.*k.*y./2./bb./M./w2+b./2./M./w2.^2+1./16./M.^2./w2.^4;
E=i.*A0.*k./2./bb./M.*exp(A).*airy(B).*airy(C);

II=abs(E).^2;
mesh(x,y,II)
view(2)


%%
figure(2)
Es=conj(E);
c0=3.*10.^8;
Edx2=diff(E,2,1);
[m1 n1]=size(Edx2);
Edy2=diff(E,2,2);
[m2 n2]=size(Edy2);
Esdx2=diff(Es,2,1);
[m3 n3]=size(Esdx2);
Esdy2=diff(Es,2,2);
[m4 n4]=size(Esdy2);
sx=(E(1:m3,1:n3).*Esdx2-Es(1:m1,1:n1).*Edx2).*i.*w./2;
sy=(E(1:m4,1:n4).*Esdy2-Es(1:m2,1:n2).*Edy2).*i.*w./2;
sz=w*k*abs(E).^2;
mm=min([m1 m2 m3 m4  n1 n2 n3 n4]);
s1=sx(1:mm,1:mm)+sy(1:mm,1:mm)+sz(1:mm,1:mm);
[dsx dsy]=gradient(s1);
quiver(x(1:mm,1:mm),y(1:mm,1:mm),dsx,dsy,'linewidth',2)
hold off , axis image
axis([-1 1 -1 1])
图片就是计算坡印廷的公式,不知道这边哪里出问题了 ,就是算不出文献里面的图

已知场分布,如何求出他的坡印廷,也就是二维能流分布
7(6R4%A(9C1TH[LTR0M0VT9.png
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 卡卡艾物理 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 发文 +13 张靖雪 2018-09-20 29/1450 2018-09-21 22:25 by shangfusheng
[硕博家园] 北航和浙大的材料 谁强 +11 红黑色的天空 2018-09-21 14/700 2018-09-21 22:07 by wbwgm
[考研] 考研交流 +4 前面有樱桃 2018-09-20 11/550 2018-09-21 16:51 by 前面有樱桃
[硕博家园] 小弟今年博一,导师补助没有是什么情况?? +28 呜呜拥有你 2018-09-17 28/1400 2018-09-21 15:34 by carangle
[考博] 申请国外博士 +7 13261202922 2018-09-20 7/350 2018-09-20 21:38 by Angela_Gao
[催化] 请问一下 实验室打水的这个泵名字叫什么? +8 卜卜卜卜卜 2018-09-18 16/800 2018-09-20 18:44 by 卜卜卜卜卜
[论文投稿] 论文格式,小白求助 +6 buyu2011 2018-09-20 13/650 2018-09-20 18:05 by 45℃天空下
[论文投稿] IEEE ACCESS重投祈福 +4 tendertime 2018-09-20 4/200 2018-09-20 16:20 by wang31aa
[高分子] 阻燃剂 10+4 8314 2018-09-19 9/450 2018-09-20 13:10 by water_man
[考博] 考博小论文 +3 水木年华欧欧 2018-09-19 10/500 2018-09-19 14:38 by 水木年华欧欧
[硕博家园] 夜深了 +5 蓝色的梦 2018-09-19 6/300 2018-09-19 07:59 by Little Bug
[硕博家园] 可以建议BAT购买国外大部分学术期刊数据库阅读下载权限建一个全球最大的学术网站吗? +9 生信-家 2018-09-18 17/850 2018-09-18 23:59 by lx12633036
[博后之家] 中山大学材料科学与工程学院翟文涛教授团队招聘博士后/特聘研究员/特聘副研究员 +3 wtzhai123 2018-09-15 5/250 2018-09-18 19:38 by wtzhai123
[考博] 有点丧,心不在焉,想某个静静 +18 音轨 2018-09-17 20/1000 2018-09-18 15:33 by 音轨
[博后之家] Postdoctoral Position: Ultrafast Transient Absorption Spectroscopy (15-SEP-2018) +3 ALESMusa 2018-09-15 3/150 2018-09-18 15:07 by 木匠的理想
[硕博家园] 导师这么回复是什么意思,要我? +21 潇潇别离 2018-09-15 36/1800 2018-09-17 08:21 by sunshine???
[论文投稿] SCI in press 算见刊吗 +3 我的名字叫ZQ 2018-09-16 3/150 2018-09-16 20:13 by 纳斯达克ML
[硕博家园] 想法论文。。 +9 微笑思密达5 2018-09-15 11/550 2018-09-16 19:42 by koalm
[考博] 考博 +4 520881 2018-09-15 4/200 2018-09-15 21:12 by 李的妮妮
[考研] 协和基础所好,还是首医基础所好考? +3 北京小爷 2018-09-15 4/200 2018-09-15 18:19 by wzarlll
信息提示
请填处理意见