查看: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 找女友 +21 晴朗的天 2018-10-16 49/2450 2018-10-19 08:35 by 晴朗的天
[环境] 抖腿男 +9 zxl960121 2018-10-17 10/500 2018-10-19 08:04 by quan3145
[硕博家园] 脱单比脱贫都难 +14 加油August 2018-10-18 18/900 2018-10-19 07:27 by lkk.1985
[找工作] 找工作的心情。有点低落。 +10 hemiao811025 2018-10-18 10/500 2018-10-19 07:25 by lkk.1985
[硕博家园] 大神们 +4 一只蝉 2018-10-18 5/250 2018-10-18 20:05 by 一只蝉
[考博] 有人想要来日本读硕或博的吗? +23 ddv8888 2018-10-15 50/2500 2018-10-18 16:18 by 落月gaea
[考博] 师兄师姐 +5 944081651 2018-10-17 7/350 2018-10-18 11:52 by 召唤师10359
[硕博家园] 刚入学,很迷茫 +16 rgt1996 2018-10-17 37/1850 2018-10-18 08:58 by 情天恨海
[教师之家] 民办本科当辅导员 +4 wushaghr 2018-10-17 11/550 2018-10-18 07:17 by 琪琪你好
[博后之家] 浙江大学医学院博士后招聘启事 +3 ricky49 2018-10-17 3/150 2018-10-17 22:39 by liuhongying1984
[找工作] 工作决定 +27 百草醇酸 2018-10-16 41/2050 2018-10-17 22:16 by zlmqsj
[硕博家园] 想去西藏高校工作一年,一般硕士可以过去吗 +6 第三个 2018-10-16 8/400 2018-10-17 11:12 by 第三个
[硕博家园] 面试 +3 nianhua~ 2018-10-16 3/150 2018-10-17 00:13 by 快到锅里来
[论文投稿] 关于电力系统自动化学报投稿 5+3 Dr-Steven 2018-10-16 4/200 2018-10-16 23:13 by nono2009
[硕博家园] 怎么摆脱失恋 +45 我就是我1 2018-10-15 71/3550 2018-10-16 22:10 by 红楼望月12
[育儿交流] 谁带娃有分歧,怎么办?? (金币+2) +20 朵儿儿 2018-10-13 23/1150 2018-10-16 16:20 by 猪猪唤醒
[考博] 一起交流! +3 bcwah717 2018-10-16 3/150 2018-10-16 12:28 by 梦时无意
[硕博家园] 求女生的观点~也欢迎男生说说自己的真实看法 +57 梦想成真2014 2018-10-12 131/6550 2018-10-16 10:48 by 梦想成真2014
[硕博家园] 博士半年了,基础差,求安慰 +8 就你还想毕业 2018-10-15 11/550 2018-10-15 19:25 by 20092290
[找工作] 焦躁 +10 1114130692 2018-10-12 13/650 2018-10-14 17:48 by purpleP
信息提示
请填处理意见