24小时热门版块排行榜     先丰纳米 买纳米材料就点这里

查看: 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个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 投稿IEEE Communication Surveys and Tutorials,给的意见是修改重投,希望大吗? +3 我发SCI 2018-04-22 7/350 2018-04-22 11:17 by zbking
[硕博家园] 发泄发泄 +13 椰子子 2018-04-21 20/1000 2018-04-22 10:33 by 胖子ing
[硕博家园] 人生何其难 +9 初秋的夏天 2018-04-20 14/700 2018-04-22 08:45 by 初秋的夏天
[高分子] 双峰油溶性丙烯酸压敏胶的合成方案 +3 liaoyufeng 2018-04-20 3/150 2018-04-22 08:20 by jiduquegai
[考研] 三跨二战东北师范比较文学拟录取 +9 淡薄K 2018-04-19 14/700 2018-04-22 07:06 by 淡薄K
[硕博家园] 读研期间的学习 +7 叫我杰哥就好 2018-04-21 15/750 2018-04-21 22:29 by tbo3780662
[机械] 隐隐觉得选错了专业,还有救吗? +9 HerrOrange 2018-04-18 18/900 2018-04-21 21:55 by wanghaogea
[职场人生] 找个合伙创业的都找不到 +13 nhhuang1223 2018-04-20 14/700 2018-04-21 20:56 by kingkongking
[考研] 政治考研 +6 一个小小学渣 2018-04-20 8/400 2018-04-21 14:30 by 考研加油12
[留学DIY] UNSW 2018 Semester 2 奖学金有人收到了吗? +4 Foreversuiyi 2018-04-16 9/450 2018-04-20 18:28 by Foreversuiyi
[考研] 2018湖南大学化学专业考研423分经验回馈(703有机+831物理化学理) +4 e_mucher 2018-04-17 9/450 2018-04-20 17:33 by e_mucher
[论文投稿] 《机械设计与制造》期刊,我真的算服了 +6 aladinn 2018-04-20 6/300 2018-04-20 14:10 by wtbccq505
[教师之家] 高校求职交流 +19 bravo1990 2018-04-18 20/1000 2018-04-19 17:01 by SunflowerL
[职场人生] 创业 +3 810843542 2018-04-19 7/350 2018-04-19 16:45 by 810843542
[考博] 武大资环学院英语线怎么还不出 +3 sun小宏宏 2018-04-18 18/900 2018-04-18 23:16 by 梅山赤兔子
[教师之家] 老师们好,求指点 +4 qewhuinqi 2018-04-18 12/600 2018-04-18 19:42 by qewhuinqi
[硕博家园] 人生都压在了科研上,没办法输 +14 杨洋1994 2018-04-18 14/700 2018-04-18 14:08 by 寻水的鱼xzy
[考博] 祈祷有个好结果 +9 ldf二呆 2018-04-16 16/800 2018-04-18 07:03 by ldf二呆
[硕博家园] “如果说你是海上的烟火,我是浪花的泡沫 +14 光源呢呢 2018-04-17 26/1300 2018-04-17 22:48 by 暮饰如纱
[考研] 做人真难 +5 小优优…… 2018-04-16 5/250 2018-04-17 03:52 by 帝企鹅7度
信息提示
请填处理意见