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

查看: 300  |  回复: 6

deletep

新虫 (初入文坛)

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

参考文献是陈辉、胡元中的《粗擦表面计算机模拟》,,这个是我模拟的图像结果,,文献中的结果见附件
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的回帖

w_tlvz_r

新虫 (初入文坛)


小木虫: 金币+0.5, 给个红包,谢谢回帖
楼主请问问题解决了吗?求指导一下,可有偿回报
5楼2018-02-02 15:51:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yanglongren

新虫 (小有名气)

6楼2018-02-03 10:17:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yanglongren

新虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
可以三个方面解决,一是N值变大,二是采取平滑编码shading  interp或者其他的平滑命令,三是扩大Z的值,使得视角深一点

发自小木虫IOS客户端
7楼2018-02-03 10:22:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 deletep 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[材料综合] 扩散偶的制备? +4 lenko 2018-04-21 8/400 2018-04-22 00:24 by 西厢REPTILE
[考博] 有梦想谁都了不起,有勇气就会有奇迹 +9 ltz 2018-04-21 9/450 2018-04-21 21:53 by wanghaogea
[硕博家园] 南京大龄在职男博士,还没对象,工作过5年多,博士怎样完成下去 +42 ay8r92 2018-04-15 44/2200 2018-04-20 21:20 by cm14704019
[有机交流] 有没有查文献的免费网址 +8 小白鼠110 2018-04-19 9/450 2018-04-20 13:36 by sj513
[硕博家园] 放养模式研一 该如何努力 +7 Stella- 2018-04-19 7/350 2018-04-19 15:55 by 萌萌哒小兔子
[无机非金属] 600℃时候将样品在马弗炉中拿出来对炉子有多大伤害? +6 徐成世 2018-04-18 8/400 2018-04-19 08:55 by 爱吃汤圆
[文学芳草园] 时间赠人阅历,世事尽可原谅 +7 孤鸿飘逸 2018-04-18 7/350 2018-04-19 00:14 by 我是本笨蛋
[考研] 考研二战失败 +23 谢绝采访 2018-04-17 28/1400 2018-04-18 21:03 by 淡淡的幸福
[硕博家园] 求一份本科毕业答辩ppt模板 谢谢学长学姐 +4 啊披pa批啵 2018-04-18 5/250 2018-04-18 16:34 by 啊披pa批啵
[公派出国] 8和11只是两个状态,与录取关系不大 +12 CQU小猪儿 2018-04-16 14/700 2018-04-18 12:07 by evanbjfu
[考博] 鉴定完毕,浙大就是肉 +3 andwp3574 2018-04-18 4/200 2018-04-18 11:37 by kawala
[精细化工] 透明质酸 +3 一缕清风888 2018-04-16 6/300 2018-04-18 11:31 by ngdzx
[考博] 2018年上海大学博士英语考试 +17 哼哼哈哈68 2018-04-15 40/2000 2018-04-18 10:55 by qiongxiaox
[考研] 学校选择 +13 刘晨涛 2018-04-15 22/1100 2018-04-18 07:43 by ctgucquli
[论文投稿] 摘要写的不好,被拒绝,还要不要修改一下投全文 +5 静电电 2018-04-17 13/650 2018-04-17 23:48 by 静电电
[论文投稿] International Journal of rock mechanics and mining sciences 的论文录用了 +9 沧浪之水_ 2018-04-16 10/500 2018-04-17 23:00 by 逝水无痕1990
[考研] 博士求调剂 +6 slp326 2018-04-17 7/350 2018-04-17 17:17 by 拜拜学术
[考研] 考研科目 +7 子潭 2018-04-17 14/700 2018-04-17 13:04 by static~state
[找工作] 二战失败找工作 +6 团子不圆123 2018-04-15 13/650 2018-04-16 21:29 by 眾下葉の甯
[硕博家园] 论文开题~如此的猝不及防,4天准备,,,全无头绪 +10 功夫小白 2018-04-15 11/550 2018-04-16 11:14 by popping周
信息提示
请填处理意见