大概是这样，我不是专业的程序猿，有些专业的东西可能不太了解

这回是用python写的一个SQGQ法生成多孔介质的程序，如果人为地更改参数降低循环次数的话，是能跑出结果的

但按照正常的来跑，根本没有头，不清楚跟数据溢出有没有关

IDE用的是VScode，大概这样

一下是我的代码，献丑了

import numpy as np
import random
from matplotlib import pyplot as plt

max_i = max_j = 100 
Dp = 0.1                                        #生长核分布概率
Vr = 0.3                                        #孔隙度
PP = [0.2,0.2,0.2,0.2,0.2,0.2,0.2,0.2]          #生长概率

#初始化孔隙数量
NumVoid = 0

#初始孔隙坐标列表
ListVoid = []

#初始化空间矩阵
SpaceArray = np.random.rand(max_i,max_j)
for i in range(max_i):
    for j in range(max_j):
        if SpaceArray[i,j] < Dp:
            #设孔隙格点为1
            SpaceArray[i,j] = 1
            #记录孔隙格点的坐标
            ListVoid.append([i,j])
        else:
            #设实体格点为0
            SpaceArray[i,j] = 0
    
#计算已确定的孔隙格点数
NumVoid = np.sum(SpaceArray)

#------------------------------------------------------------------

#计算满足孔隙度要求的孔隙格点个数
NumVoid_Need = Vr*(max_i*max_j)

#开始生长
while NumVoid<NumVoid_Need:
    for PointVoid in ListVoid:
        Point_i = PointVoid[0]
        Point_j = PointVoid[1]

        #孔隙向上发展
        #排除格点在上边界的情况
        if Point_i > 0:
            i = Point_i - 1
            j = Point_j
            if SpaceArray[i,j]==0 and random.random()<PP[0]:
                SpaceArray[i,j] = 1

        #孔隙向下发展
        #排除格点在下边界的情况
        if Point_i < max_i - 1:
            i = Point_i + 1
            j = Point_j
            if SpaceArray[i,j]==0 and random.random()<PP[1]:
                SpaceArray[i,j] = 1

        #孔隙向左发展
        #排除格点在左边界的情况
        if Point_j > 0:
            i = Point_i
            j = Point_j - 1
            if SpaceArray[i,j]==0 and random.random()<PP[2]:
                SpaceArray[i,j] = 1

        #孔隙向右发展
        #排除格点在右边界的情况
        if Point_j < max_j - 1:
            i = Point_i
            j = Point_j + 1
            if SpaceArray[i,j]==0 and random.random()<PP[3]:
                SpaceArray[i,j] = 1

        #孔隙向左上发展
        #排除格点在左边界或上边界的情况
        if Point_j>0 and Point_i>0:
            i = Point_i - 1
            j = Point_j - 1
            if SpaceArray[i,j]==0 and random.random()<PP[4]:
                SpaceArray[i,j] = 1

        #孔隙向右上发展
        #排除格点在右边界或上边界的情况
        if Point_j<max_j-1 and Point_i>0:
            i = Point_i - 1
            j = Point_j + 1
            if SpaceArray[i,j]==0 and random.random()<PP[5]:
                SpaceArray[i,j] = 1

        #孔隙向左下发展
        #排除格点在左边界或下边界的情况
        if Point_j>0 and Point_i<max_i-1:
            i = Point_i + 1
            j = Point_j - 1
            if SpaceArray[i,j]==0 and random.random()<PP[6]:
                SpaceArray[i,j] = 1

        #孔隙向右下发展
        #排除格点在右边界或下边界的情况
        if Point_j<max_j-1 and Point_i<max_i-1:
            i = Point_i + 1
            j = Point_j + 1
            if SpaceArray[i,j]==0 and random.random()<PP[7]:
                SpaceArray[i,j] = 1
        

    #重新计算已确定的孔隙格点数
    NumVoid = np.sum(SpaceArray)
print(NumVoid/(max_i*max_j))

plt.imshow(SpaceArray)
plt.show()

c语言中关于末尾空格的问题

qq_36722997

4个回复