24小时热门版块排行榜    

Znn3bq.jpeg
查看: 671  |  回复: 2

coinqin

新虫 (初入文坛)

[求助] 如何用matlab画出以坐标轴X,Y为参数的I(X,Y)的图像,即光强和相位图

MATLAB新手因为要画光强和相位图,U(X,Y)是光场复振幅,是一个e的含虚部的次方的公式的叠加,下面程序中有。
光强I=|U(X,Y)|^2,相位phi=arctan(ImU(X,Y)/ReU(X,Y))

贴上自己编的代码,如果有基础性的错误希望大家指出,谢谢
最后 代码画不出来正确图形,不知道是用哪个函数画图,也不知道是不是MATLAB不能够处理这样一个e的含虚部的次方的公式的叠加
诚心谢谢大神帮忙改错!!!

clc;
n=input('Enter the n value:');
m=input('Enter the m value:');
lambda=6.325*10^(-7);
z=0.1;
d=10^(-3);
p=d;
q=0;
alpha=0;
U=exp(-2i*pi*(X*p+Y*q)/lambda/z)*exp(1i*m*alpha);
for t=1:n-1
    alpha=2*t*pi/n;
    d=sqrt(d^2+2*m*lambda*z/n);
    X=d*cos(t*alpha);
    Y=d*sin(t*alpha);
    U=exp(-2i*pi*(X*p+Y*q)/lambda/z)*exp(1i*m*alpha);
    %I=U*conj(U);
    I=abs(U).^2;             %两种方法都可以求出光强
    phi=atan(U);              %相位
    plot(phi,[-7*10^(-5),7*10^(-5),-7*10^(-5),7*10^(-5)]);              该方法可以直接画出以X,Y为坐标的相位图吗?还是用surf方法,但是提示I必须为矩阵,不知道怎么设置成矩阵。
   % surf(X,Y,I);   
   %plot(X,Y,I);
   % mesh(x,y,I);
    axis auto
   
end
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

coinqin

新虫 (初入文坛)

X=d*cos(t*alpha);
    Y=d*sin(t*alpha);写错了
应该为p=d*cos(t*alpha);
    q=d*sin(t*alpha);
2楼2016-03-02 20:11:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

coinqin

新虫 (初入文坛)

有人可以解答吗,拜托了
3楼2016-03-03 20:12:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 coinqin 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 26年申博自荐-计算机视觉 20+3 光亮我自己 2026-05-01 4/200 2026-05-03 09:01 by lfy8008
[博后之家] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 qs8kylt1wr 2026-05-01 3/150 2026-05-03 07:22 by jmr1kibhvw
[公派出国] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 qs8kylt1wr 2026-05-01 4/200 2026-05-03 07:19 by jmr1kibhvw
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 d8tbypao1q 2026-05-01 3/150 2026-05-03 06:04 by jmr1kibhvw
[找工作] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 d8tbypao1q 2026-05-01 3/150 2026-05-03 05:07 by jmr1kibhvw
[博后之家] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 d8tbypao1q 2026-05-01 4/200 2026-05-03 04:19 by w5idbsc3mz
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 d8tbypao1q 2026-05-01 4/200 2026-05-03 03:51 by w5idbsc3mz
[公派出国] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k17zww5nna 2026-05-01 4/200 2026-05-03 03:04 by w5idbsc3mz
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k17zww5nna 2026-05-01 4/200 2026-05-03 02:36 by w5idbsc3mz
[公派出国] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 k17zww5nna 2026-05-01 5/250 2026-05-03 02:04 by w5idbsc3mz
[找工作] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 k17zww5nna 2026-05-01 5/250 2026-05-03 01:49 by w5idbsc3mz
[考博] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k17zww5nna 2026-05-01 4/200 2026-05-03 01:34 by w5idbsc3mz
[考博] 考博 +4 Eternitf 2026-04-29 4/200 2026-05-02 16:10 by yejenny
[公派出国] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 d8tbypao1q 2026-05-01 3/150 2026-05-02 14:21 by zahslv263y
[教师之家] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 d8tbypao1q 2026-05-01 3/150 2026-05-02 14:06 by zahslv263y
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k17zww5nna 2026-05-01 5/250 2026-05-02 09:21 by dz7y1p1i98
[考博] 药化及相关博士的申请 5+3 中药学渣 2026-04-29 3/150 2026-05-01 18:14 by babero
[文学芳草园] 太白金星有点烦 +3 yexuqing 2026-04-27 3/150 2026-04-30 21:56 by myrtle
[基金申请] 收到国自然专家邀请后几年才会有本子送过来评 +3 barley5 2026-04-30 3/150 2026-04-30 14:58 by aspect3000
[论文投稿] 有没有快的中文核心比较快录用的,纳米材料光催化 5+3 2914439618 2026-04-27 4/200 2026-04-29 08:59 by 北京莱茵润色
信息提示
请填处理意见