同一维粗糙面模拟一样,公式和基本理论神马的我就不敲了 有兴趣的同学可以参考
http://files.cnblogs.com/xd-jinjian/%E7%AC%AC%E4%BA%8C%E7%AB%A0%EF%BC%88%E5%BB%BA%E6%A8%A1%EF%BC%89.pdf
二维粗糙面的效果图如下所示
具体代码如下
// 2D_GS_RS.cpp : 定义控制台应用程序的入口点。
//*************本程序基于Monte Carlo方法生成二维高斯粗糙面**************
//by changwei
//**********************************************************************
#include "stdafx.h"
#include<iostream>
#include<iomanip> //格式化输出控制符必包含的头文件
#include<fstream>
#include<cmath>
#include<stdlib.h>
#include<complex>
#define M 64 //x方向采样点数
#define N 64 //y方向采样点数
#define Mdx 8 //x方向单位波长内的采样点数
#define Ndy 8 //y方向单位波长内的采样点数
#define pi 3.141592627 //定义常数pi
using namespace std;
complex<double> unit_i(0.0,1.0);
complex<double> unit_r(1.0,0.0);
//产生随机数的子函数
void random2(double start,double end,double seed,double *rand2) //在用二维数组作为形参时,必须指明列数
{
double s=65536.0;
double w=2053.0;
double v=13849.0;
double T=0.0;
int m;
for(int i=0;i<N*M;i++)
{
T=0.0;
for(int k&