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的会议 要求补上实验结果 +4 红肠大列巴 2018-07-19 10/500 2018-07-20 14:43 by zhangbarton
[硕博家园] 又回到最初的起点 +13 牛奶抹茶 2018-07-13 16/800 2018-07-20 11:18 by 期待成长
[有机交流] 找寻实验室短工 +14 mygrain 2018-07-19 15/750 2018-07-20 10:52 by 浩浩汤汤
[留学DIY] 2020出国读博士求队友 +6 lgj2392 2018-07-18 6/300 2018-07-19 23:46 by 刘佳狗
[签证指南] J1面签 +7 abby0821 2018-07-18 9/450 2018-07-19 19:40 by any_mxm
[硕博家园] 读研成果的不可控因素 +3 灿灿ac 2018-07-17 4/200 2018-07-19 16:46 by jyj4947
[找工作] 有在渭南师范学院工作或者上过学的虫友吗? +17 混世小魔王l 2018-07-14 36/1800 2018-07-19 16:08 by 饭没了吃吃
[硕博家园] 寻找交往对象(本人男) +11 xiarileixin 2018-07-16 12/600 2018-07-18 23:39 by 我变胖求
[教师之家] 心情不爽~~牢骚一下 +13 lambert2014 2018-07-17 13/650 2018-07-18 21:48 by 门外汉2018
[有机交流] 重氮化反应关环关不上,求助. +3 药化狗 2018-07-13 4/200 2018-07-18 20:50 by hyacinth1944
[留学生活] 刮胡子被种族歧视怎么办? +16 霸王防脱 2018-07-14 21/1050 2018-07-18 19:39 by FLIPED
[考研] 有机问题求助 +5 氢氧化铝AL 2018-07-16 7/350 2018-07-18 10:43 by 蝶变ing
[基金申请] 有哪位知道中科院大科学装置科学研究联合基金什么时候开始会评? +7 zilan616 2018-07-16 9/450 2018-07-18 09:53 by wenyanyin
[职场人生] 天下第一奇葩 +8 我爱我家66 2018-07-16 19/950 2018-07-17 16:35 by tla02
[考研] 考研 +5 yxynl 2018-07-17 5/250 2018-07-17 15:34 by supreme健
[基金申请] 信息学部是不是该公布评审专家了? +3 ty_duck 2018-07-16 5/250 2018-07-17 11:03 by brigie2016
[基金申请] 五个评审专家,上会要达到什么条件?什么是a类上会? +7 gaofz 2018-07-14 10/500 2018-07-16 08:34 by xjtu-dogface
[硕博家园] 机械硕士(偏系统开发)考博疑问 +4 LIBO呆鱼 2018-07-15 10/500 2018-07-16 06:04 by LIBO呆鱼
[硕博家园] 求建议 +22 计算qi 2018-07-14 24/1200 2018-07-14 23:36 by 1584715932
[考研] 工业设计跨考理论物理难度真的很大么? +4 我-胖虎 2018-07-13 9/450 2018-07-14 12:13 by zyl雨田
信息提示
请填处理意见