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

热门推荐

查看: 150  |  回复: 3

deletep

新虫 (初入文坛)

[交流] MATALB生成具有指定自相关函数的粗糙表面已有3人参与

参考文献是陈辉、胡元中的《粗擦表面计算机模拟》,,这个是我模拟的图像结果,,文献中的结果见附件
clear
clc

N=512;
L=0.05;
delta=0.05;%表面均方根粗糙度
betax=10;%x方向的相关长度
betay=10;
C=1;%功率谱密度

dx=L/N;dy=dx;
NN=-N/2:N/2-1;
[Nx,Ny]=meshgrid(NN,NN);
taux=dx.*Nx;tauy=dy.*Ny;

%%生成具有指定自相关函数的粗糙表面
eta=randn(N,N);%高斯分布白噪声
A=fft2(eta);%傅里叶变换
R=zeros(N,N);
R=delta^2*exp(-2.3*((taux/betax).^2+(tauy/betay).^2).^0.5);%自相关函数
Gz=1/(2*pi^2).*fft2(R);%功率谱密度函数
H=(Gz/C).^0.5;%传递函数
Z=H.*A;%表面高度的傅里叶变换
z=ifft2(Z);%表面高度分布

figure(1)
meshz(abs(z));
title('rough surface');
axis square

有没有哪位大神能帮忙看看,问题出在哪儿
更多 0

MATALB生成具有指定自相关函数的粗糙表面
1390271227_63743.jpg


MATALB生成具有指定自相关函数的粗糙表面-1
1390271252_165409.jpg
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

刘忆庵

新虫 (初入文坛)


小木虫: 金币+0.5, 给个红包,谢谢回帖
楼主,解决模拟的粗糙表面不光滑这个问题了吗?
2楼2016-06-29 18:40:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

KA_KA_D

新虫 (初入文坛)


小木虫: 金币+0.5, 给个红包,谢谢回帖
楼主,你的问题解决了吗
3楼2017-11-17 11:37:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

木石君

新虫 (初入文坛)


小木虫: 金币+0.5, 给个红包,谢谢回帖
没看出程序有啥问题啊,哪位大神看出来了
4楼2017-12-18 15:45:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 deletep 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 硕士上了一年半,现在想退学。 (金币+2) +150 zhmou 2018-01-18 153/7650 2018-01-24 07:09 by 王三宁
[教师之家] 今天的狗特别特别的乖 +7 么么茶XYZ 2018-01-23 22/1100 2018-01-24 00:34 by paulinho
[论文投稿] 有投稿ICASSP 2018的吗,交流一下 +4 火力全开skg 2018-01-20 6/300 2018-01-23 21:27 by 火力全开skg
[留学DIY] 教授这样回复我…啥意思…不知道该怎么回复 +4 YuzuJuri 2018-01-23 6/300 2018-01-23 19:26 by bigplume
[硕博家园] 身为PHD,老板跟你说你的工作本科生都能做,你会怎么回复呢? +20 blink星 2018-01-22 23/1150 2018-01-23 14:08 by pigfei
[公派出国] 在职公派 +3 caicc123 2018-01-23 5/250 2018-01-23 13:25 by yitianzhi
[基金申请] 关于协同机制的本子的写法 +3 hibluesky 2018-01-21 6/300 2018-01-23 12:20 by zzhshan911
[硕博家园] 研一上学期稀里糊涂读完了,最后啪,挂了一科专业必修课 +32 带着面具生活 2018-01-17 71/3550 2018-01-23 10:55 by 带着面具生活
[有机交流] 关于回流反应的问题 +5 ZJU之梦 2018-01-17 7/350 2018-01-23 10:44 by cgh_518
[考研] 如果以后打算从事高校学术研究,考研考哪比较好 +6 陌上丶花开 2018-01-22 10/500 2018-01-23 03:53 by 水煮鱼好吃
[找工作] 青岛职业技术学院 +3 xkcyj 2018-01-17 5/250 2018-01-21 21:04 by xkcyj
[硕博家园] 研一都快过完了,感觉自己什么怎么长进呀 +6 J天王盖地虎 2018-01-18 24/1200 2018-01-20 18:31 by J天王盖地虎
[公派出国] 大家帮我看看,我有可能申请博士联合培养的出去么? +8 daliaojie205 2018-01-19 18/900 2018-01-20 12:55 by zmh0611
[无机/物化] 各路神人给个建议呗! +9 迷茫啊迷茫 2018-01-19 13/650 2018-01-19 20:50 by starstar_76
[考研] 立志篇 +7 peace_years 2018-01-19 15/750 2018-01-19 19:12 by 学海小舟
[硕博家园] 做个农民也好 +17 虫子咬一口 2018-01-17 25/1250 2018-01-19 15:37 by DA懒虫子
[版块工场] 【早起签'到'贴】2018年1月19日(q q 48626260 进群申请:小木虫——您个人ID) +66 8475 2018-01-19 66/3300 2018-01-19 08:27 by zxf984
[考研] 19物理系考研 +4 E.DongHuang 2018-01-18 4/200 2018-01-19 02:48 by 火柴人233
[高分子] 低聚物可以用挤出机加工吗 15+5 zengwz 2018-01-18 10/500 2018-01-18 11:04 by 半亩花田116
[版块工场] 【早起签'到'贴】2018年1月18日(q q 48626260 进群申请:小木虫——您个人ID) +43 8475 2018-01-18 44/2200 2018-01-18 10:08 by ltyz2001
信息提示
请填处理意见