查看: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 基金激励了一批人,同时也扼杀了一批人! +15 平凡冰雪花 2018-07-15 15/750 2018-07-16 06:49 by huangyq79
[基金申请] 工材青基涵审意见ABC能上会吗? +11 tianjia901 2018-07-15 12/600 2018-07-16 05:05 by pigthree
[硕博家园] 跟了个水货导师,真心醉了 +38 ddc7451 2018-07-10 63/3150 2018-07-16 03:47 by Wei_ren
[论文投稿] 专硕考博 +5 被拒20次 2018-07-15 10/500 2018-07-15 21:57 by 能源大师
[基金申请] 祈福基金中 (金币+2) +14 liuxin590 2018-07-14 14/700 2018-07-15 21:42 by ocbc
[硕博家园] 导师不让读博士 +36 haiyan1217 2018-07-11 53/2650 2018-07-15 18:16 by haiyan1217
[留学DIY] 现在套磁成功了应该做什么 +7 Soy123 2018-07-15 13/650 2018-07-15 17:25 by Carrie_g
[硕博家园] 考博英语推荐 +4 qmf561242 2018-07-15 6/300 2018-07-15 17:10 by qmf561242
[硕博家园] 高晓松竟然成了哈佛大学研究员! +21 a155644 2018-07-14 58/2900 2018-07-15 15:50 by cscshu
[找工作] 面试如何穿搭 +10 CassielYuer 2018-07-11 10/500 2018-07-15 11:25 by 小耿子
[有机交流] 关于C-N偶联,很简单的反应做不出来 +4 佟耳朵 2018-07-14 4/200 2018-07-14 20:28 by gtywugang
[基金申请] 今晚保密工作做的真好! +29 堂吉可德 2018-07-10 37/1850 2018-07-14 10:05 by a2004280229
[找工作] 什么样的工作安稳轻松 +8 归船系浅渡 2018-07-13 9/450 2018-07-13 21:42 by cln116
[硕博家园] 投SCI时,有没有Supplementary data 会影响文章质量和评价吗? +4 18354212358 2018-07-12 4/200 2018-07-12 22:37 by 029539zbb
[硕博家园] 求XRD分析软件 +8 芝芝芝 2018-07-10 16/800 2018-07-12 09:01 by li1102289962
[基金申请] 一般单位自评,反映好的,反而中不了。大家有这种感觉吗? +19 natureding 2018-07-10 24/1200 2018-07-12 07:32 by 天天向上3035
[考研] 【招生】华中科技大学计算机学院陈敏教授招收硕博士研究生 +3 hustepic_ 2018-07-09 3/150 2018-07-11 17:10 by qiulaiqi
[分析] 邻菲罗林测铁时空白为什么也有颜色 +3 akdjd 2018-07-11 3/150 2018-07-11 15:28 by luoyuan121
[硕博家园] 做院士的学生 +16 888888就是觉 2018-07-09 18/900 2018-07-09 14:02 by cococ_可可西
[硕博家园] 国外的大学博士,排名很靠后的值得去吗? +8 luokereng 2018-07-09 12/600 2018-07-09 13:18 by luokereng
信息提示
请填处理意见