24小时热门版块排行榜    

查看: 1512  |  回复: 6
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

459354755

新虫 (初入文坛)

[求助] matlab最后一步的作图问题 已有1人参与

程序现在都写好了,就是最后的作图程序了,本人没有学过matlab所以麻烦写出程序。  谢谢!
程序如下:
k=2*pi/(532*10^-9);
z=0.1;
m=1;
l=0;
theta=0;
r0=0.005;
nnx=256;
xx=-0.005:0.01/255:0.005;
yy=-0.005:0.01/255:0.005;
for nn=1:nnx
   for mm=1:nnx
    phi=atan2(yy(mm), xx(nn));
    rho=sqrt(yy(mm).^2+xx(nn).^2);
f(nn,mm)=quad(@(r)r.*exp(-0.1*i*k*r.^1.5)*k*i.^(m+2).*besselj(m,-k*rho*r/z).*sin(m*phi+2*pi*(r/r0).^l+theta),0,r0);
end
end
差最后一个出图像的语言。
最后结果类似于那张图片。
程序调试时间很久,所以只求正确的图像求法。
谢谢各位!

matlab最后一步的作图问题
Z8RST[5IFPCQ_3{UJ[QU[}7.png
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

459354755

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by eidie at 2015-05-07 18:38:17
fmax=max(max(f));
fmin=min(min(f));
f2=(f-fmin)/(fmax-fmin); %数值太小,归一化
imagesc(f2);
colormap jet;
colorbar;

fmax=max(max(f));
fmin=min(min(f));
f2=(f-fmin)/(fmax-fmin); %数值太小,归一化[0, 1]
imagesc(f2);
colormap jet;
colorbar;
Error using image
Image CData can not be complex

Error in imagesc (line 19)
    hh = image(varargin{1},'CDataMapping','scaled');
对不起,运行之后程序错误,你可不可以稍微改下 谢谢
3楼2015-05-08 08:12:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

eidie

木虫之王 (文坛精英)

文献杰出贡献

【答案】应助回帖

感谢参与,应助指数 +1
fmax=max(max(f));
fmin=min(min(f));
f2=(f-fmin)/(fmax-fmin); %数值太小,归一化[0, 1]
imagesc(f2);
colormap jet;
colorbar;
跟一个女博士相亲,我问她:“你学历这么高,嫁的出去吗?”她说:“这个问题,等我回去写篇论文回答你。”
2楼2015-05-07 18:38:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

eidie

木虫之王 (文坛精英)

文献杰出贡献

引用回帖:
3楼: Originally posted by 459354755 at 2015-05-07 12:12:38
fmax=max(max(f));
fmin=min(min(f));
f2=(f-fmin)/(fmax-fmin); %数值太小,归一化
imagesc(f2);
colormap jet;
colorbar;
Error using image
Image CData can not be complex

Error in imagesc (line  ...

根据提示,你自己运算出现复数。不知道你要显示的是模,还是相位。
跟一个女博士相亲,我问她:“你学历这么高,嫁的出去吗?”她说:“这个问题,等我回去写篇论文回答你。”
4楼2015-05-08 08:52:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

459354755

新虫 (初入文坛)

引用回帖:
4楼: Originally posted by eidie at 2015-05-08 08:52:34
根据提示,你自己运算出现复数。不知道你要显示的是模,还是相位。...

相位      谢谢
5楼2015-05-08 09:23:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 建议基金发布提前给出明确的时间点 +10 kulium 2026-08-21 12/600 2026-08-21 19:46 by kulium
[基金申请] 科研孤儿太难了 +16 我4大白菜 2026-08-20 17/850 2026-08-21 19:07 by sakuraai
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +15 医学老男孩 2026-08-20 19/950 2026-08-21 18:54 by 医学老男孩
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] filecode,4个jtjc了 +11 ziyangfang 2026-08-19 14/700 2026-08-21 16:51 by applepetty
[基金申请] 让我中一个面上吧! +11 大萍1987 2026-08-20 13/650 2026-08-21 14:48 by 20120902066
[基金申请] 感觉是下周放榜了 +7 angus9576 2026-08-17 12/600 2026-08-21 13:38 by weiyin
[基金申请] 我面上完蛋了 +7 且听虎啸 2026-08-20 8/400 2026-08-21 12:31 by 酷酷墨镜
[基金申请] 只有每年这种时候来逛逛小木虫 +23 yaoyewhu2008 2026-08-20 25/1250 2026-08-21 12:19 by Poppy1104
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 今天系统多次维护,明天很可能放榜! +10 zju2000 2026-08-16 11/550 2026-08-20 20:02 by cl479861084
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[教师之家] 为什么余额宝的年化利率越来越低?主要原因有哪些? +5 瞬息宇宙 2026-08-15 5/250 2026-08-19 20:23 by super2002521
[基金申请] 今天放榜吗? +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
信息提示
请填处理意见