查看: 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个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[有机交流] Buchwald-Hartwig C-N偶联反应,做了好几次了,还没有产物 20+4 liushuting12 2018-12-13 20/1000 2018-12-17 00:32 by cunzaigan
[有机交流] 噻唑羧酸的tlc点板 +4 18352810677 2018-12-13 6/300 2018-12-17 00:19 by cunzaigan
[文学芳草园] 奈何 +7 mingfengpbpb 2018-12-16 15/750 2018-12-16 21:30 by mingfengpbpb
[博后之家] 跟了神一样的老板,我前世造了什么孽啊 +13 bujunyang 2018-12-14 20/1000 2018-12-16 19:18 by 西瓜雨
[硕博家园] 大家都觉得我应该读博 +23 夏蒟蒻just 2018-12-15 49/2450 2018-12-16 18:29 by 夏蒟蒻just
[找工作] 万华vs中石化上海工程公司 +3 liruchun1993 2018-12-12 5/250 2018-12-16 16:37 by liruchun1993
[分析] 哪里能测锡的穆斯堡尔谱 +3 lithiumcell 2018-12-15 4/200 2018-12-16 16:25 by 斩龙子小飞侠
[硕博家园] 熵焓变均增加能说明啥 +3 17332179410 2018-12-16 3/150 2018-12-16 12:23 by great_zhw
[硕博家园] 现在要自闭了 +10 adolp 2018-12-15 22/1100 2018-12-16 00:04 by 卡罗拉之夜
[论文投稿] 时间都去哪儿了 +3 springhope 2018-12-15 3/150 2018-12-15 22:09 by liouzhan654
[数学] 一个极限的计算与证明 +4 Edstrayer 2018-12-12 10/500 2018-12-15 13:53 by qlf2007
[职场人生] 当我在犹豫要不要辞职的时候,同事先辞职了 +3 xiaoxiao1231 2018-12-11 4/200 2018-12-14 19:11 by gxly1984
[找工作] 年薪三十万还是问题吗 +6 路随心转1 2018-12-11 13/650 2018-12-14 14:34 by 路随心转1
[电化学] 用RHE氢参比电极电位准确,实验数据更可靠,是这样吗? 2+3 scdb123 2018-12-13 4/200 2018-12-14 09:48 by minbest
[公派出国] 请问CSC对英国大学联合培养博士的板凳费是什么政策 +4 天天最动人 2018-12-13 8/400 2018-12-14 09:29 by 望着大河弯弯
[仿真模拟] 充气过程仿真的发散问题 200+3 brucexmc 2018-12-11 9/450 2018-12-14 09:05 by 良人
[有机交流] 重金属镍的去除方法 50+5 hjr123 2018-12-12 9/450 2018-12-13 20:11 by 天山明月
[分析] OpenLAB CDS EZChrom在哪可以下载? +4 yekenyh 2018-12-11 5/250 2018-12-13 15:22 by wslwydgv
[硕博家园] 有机知识学习 +3 123不要动 2018-12-12 5/250 2018-12-13 00:19 by 积跬步
[硕博家园] 全程15万的工程博士有没有必要读 +13 梦翎翊 2018-12-10 22/1100 2018-12-10 20:50 by socitey
信息提示
请填处理意见