二维粗糙面的模拟

同一维粗糙面模拟一样，公式和基本理论神马的我就不敲了 有兴趣的同学可以参考

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&