小游园-CSDN博客

原创 pyhton_基尼指数计算

1.定义：基尼指数（基尼不纯度）：表示在样本集合中一个随机选中的样本被分错的概率。注意： Gini指数越小表示集合中被选中的样本被分错的概率越小，也就是说集合的纯度越高，反之，集合越不纯。def gini_index_single(a,b): single_gini = 1 - ((a/(a+b))**2) - ((b/(a+b))**2) return round(single_gini,4)# 求取基尼指数def gini_index(a,b,c,d,e,f...

2020-07-01 14:45:29 1356

原创 ANSYS MESHING网格种类及区别

在ANSYS网格划分中，存在多种网格划分方法，如补片独立、补片协调、多区等方法，不同的划分方法有着不同的功能和结果，在一个复杂结构的网格划分中，不同的网格划分方法通过扮演不同的角色，来解决网格质量问题。下面介绍几种网格划分方法：①补片协调算法（Patch Conforming）补片协调算法是四面体网格类型中的一种算法，也是MESHING默认的一种算法，这种方法计算速度快，但网格质量一般。②补片独立算法（Patch Independent）补片独立算法也是四面体算法中的一...

2020-07-01 14:44:21 11965

原创 fluent动网格profile文件的编写

在fluent中指定部件运动主要有两种方式：profile文件和UDF宏文件。当运动方程已知时，利用UDF宏指定运动非常方便，而profile方法则适合当需要指定的运动无任何函数关系的情况。标准瞬态Profile文件可表述为以下格式：（（profile-name transient n periodic？）(field_name_1 a1 a2 a3…an)(field_name_2 b1 b2 b3…bn)...(field_name_r r1 r2 r3 …rn))

2020-07-01 14:40:47 9728 1

原创 Python绘图布局、图例说明、边界

当绘画完成后，会发现X、Y轴的区间是会自动调整的，并不是跟传入的X、Y轴数据中的最值相同。为了调整区间可以使用下面的方式：ax.set_xlim([xmin, xmax]) #设置X轴的区间ax.set_ylim([ymin, ymax]) #Y轴区间ax.axis([xmin, xmax, ymin, ymax]) #X、Y轴区间ax.set_ylim(bottom=-10) #Y轴下限ax.set_xlim(right=25) #X轴上限具体效果如...

2020-07-01 14:39:54 10638

原创 BP神经网络及其Python实现

import numpy as npdef Sigmoid( v ): y = 1 / ( 1 + np.exp( -v ) ) return ydef TrainWByBP( W1, W2, X, D ): alpha = 0.9 [ M, N ] = np.shape( X ) for i in range( M ): #Forword propagation x = np.array( X[i...

2020-07-01 14:37:54 1590

原创使用python对Apriori算法进行学习

Apriori的原理解释：假设我们在经营一家商品种类并不多的杂货店，我们对那些经常在一起被购买的商品非常感兴趣。我们只有4种商品：商品0，商品1，商品2和商品3。那么所有可能被一起购买的商品组合都有哪些？这些商品组合可能只有一种商品，比如商品0，也可能包括两种、三种或者所有四种商品。我们并不关心某人买了两件商品0以及四件商品2的情况，我们只关心他购买了一种或多种商品from numpy import *def loadDataSet(): return [[1, 3, 4], [2, 3.

2020-06-15 15:24:34 597

原创 Python break语句

Python break语句，就像在C语言中，打破了最小封闭for或while循环。break语句用来终止循环语句，即循环条件没有False条件或者序列还没被完全递归完，也会停止执行循环语句。break语句用在while和for循环中。如果使用嵌套循环，break语句将停止执行最深层的循环，并开始执行下一行代码。如下例：for letter in 'Python': # 第一个实例 if letter == 'h': break print '当前字母 :', letter var =.

2020-06-15 15:23:44 756

原创 python list列表添加元素

Python中使用+运算符可以将多个序列连接起来；列表是序列的一种，所以也可以使用+进行连接，这样就相当于在第一个列表的末尾添加了另一个列表。language = ["Python", "C++", "Java"]birthday = [1991, 1998, 1995]info = language + birthdayprint("language =", language)print("birthday =", birthday)print("info =", in..

2020-06-15 15:23:07 3351

原创 solidworks静力学分析概述

在之前的文章中，我们介绍了solidworks的流体分析模块，本期我们将介绍solidworks的静力学分析模块。静力学分析模块是SOLIDWORKS Simulation插件下的一个子分析模块，这个分析插件还可进行热力分析、频率分析、跌落测试和非线性分析等有限元分析。步骤一：新建一个要分析的模型：步骤二：首先启动SOLIDWORKS Simulation插件：启动后如图所示：步骤三：新建静力学分析算例步骤四：为模型添加固定载荷步骤五：为模型添加外部载荷

2020-06-15 15:22:07 12472 1

原创 python卷积代码

import structimport matplotlib.pyplot as pltimport numpy as npdateMat = np.ones((7,7))kernel = np.array([[2,1,1],[3,0,1],[1,1,0]])def convolve(dateMat,kernel): m,n = dateMat.shape km,kn = kernel.shape newMat = np.ones(((m - km + 1),(n - .

2020-06-15 15:20:55 1169

原创模态分析结果的输出实例

/clear,start mass=83500 /prep7 k,1,0,0,0 k,2,10000,0,0 lstr,1,2 et,1,beam44 et,2,mass21 R,1,186160,178980000000,56238000000,1278,680,128980000000, RMODIF,1,21,1278,680, RMODIF,1,25,76160,110000,0,0,7.0678e-009,0, *do,i,1,5 j=i*1000

2020-06-08 17:44:57 3017

原创使用python对决策树算法进行学习

#-*-coding:utf-8-*-from math import logimport operatordef calcShanonEnt(dataSet): ''' 计算给定数据集的香农熵 :param dataSet: :return:shanonEnt ''' numEntries = len(dataSet) labelCounts={} for featVec in dataSet: currentLab.

2020-06-08 17:43:57 314

原创 Python函数定义和调用

Python中函数的应用非常广泛，比如 input() 、print()、range()、len() 函数等等，这些都是 Python 的内置函数，可以直接使用。除了可以直接使用的内置函数外，Python还支持自定义函数，即将一段有规律的、可重复使用的代码定义成函数，从而达到一次编写、多次调用的目的。如len() 函数，通过它我们可以直接获得一个字符串的长度，如果没有 len() 函数，要想获取一个字符串的长度，可用如下代码：n=0for c in "http://c.biancheng.net

2020-06-08 17:42:07 2062

原创 abaqus钻头钻进模拟

在之前的文章中，笔者介绍了二维切削的模拟方法，本次笔者将介绍三维钻头钻进岩石的模拟方法。步骤一：准备模型由于在abaqus中绘图和装配比较繁琐，笔者选择在solidwork中提前建好模型，并进行装配，再以x.t格式存储和导入abaqus。步骤二：材料属性的设置材料属性的设置比较重要，它直接影响到模拟结果的准确性，由于本次是做算例讲解，为了简便仅给出基础数据：岩石属性：密度：2.79e-009；弹性：杨氏模量：70000，泊松比：0.33塑性：钻头属性：传导

2020-06-08 17:40:47 4905 4

原创 TensorFlow2.0代码案例

import tensorflow as tfimport numpy as npinput_xs = np.random.rand(1000)input_ys = 3 * input_xs + 0.217weight = tf.Variable(1., dtype=tf.float32, name="weight")bias = tf.Variable(1., dtype=tf.float32, name="bias")def model(xs): logits = tf.mult.

2020-06-08 17:39:07 1796

原创陶瓷杆撞击铝板的例子

APDL：/filename,barplate /title,Impact Of A Bar Into A Plate /PREP7 /view,,1,2,3 /plopts,info,1 /pbc,all,1 /pnum,mat,1 /num,1 ET,1,SOLID164 ! --- Enter Ceramic Material Properties For Bar mpmod,1,1 mp,ex,1,330e9 mp,nuxy,1,.24 mp,dens,1,8400 .

2020-06-02 00:26:46 444

原创使用python对信息熵进行计算

在机器学习中的决策树类算法中主要依靠信息熵的大小来选择重要的特征属性作为节点换分数据集从而获得训练结果#-*-coding:utf-8-*-from math import logimport operatordef calcShanonEnt(dataSet): ''' 计算给定数据集的香农熵 :param dataSet: :return:shanonEnt ''' numEntries = len(dataSet) labelCou

2020-06-02 00:25:59 3053

原创 Python for循环语句

Python for循环可以遍历任何序列的项目，如一个列表或者一个字符串。语法：for循环的语法格式如下：for iterating_var in sequence: statements(s)实例for letter in 'Python': # 第一个实例print '当前字母 :', letterfruits = ['banana', 'apple', 'mango']for fruit in fruits: # 第二个实例print '当前水果 :'...

2020-06-02 00:25:07 293

原创 Python算术运算符用法

print()函数使用以%开头的转换说明符对各种类型的数据进行格式化输出。+ 加法运算符加法运算符很简单，和数学中的规则一样：m = 10n = 97sum1 = m + nx = 7.2y = 15.3sum2 = x + yprint("sum1=%d, sum2=%.2f" % (sum1, sum2) )运行结果：sum1=107, sum2=22.50拼接字符串当+用于数字时表示加法，但是当+用于字符串时，它还有拼接字符串（将两个字符..

2020-06-02 00:24:16 1126

原创 solidworks流体模拟分析概述

Solid works是一款极为简便的三维建模软件，同时该软件又集合了运动分析、静力学分析和流体力学分析等模块。本次我们将讲解solidworks流体分析模块。首先我们新建一个流道模型：然后加载solidworks流体分析模块：然后将流道两端封堵，创建封盖：通过向导创建流体分析项目：通过向导设置单位、流体性质、传热方式等。设置完成后即可进入流体分析设置，软件会自动判断计算的流体区域，当然也可通过拖动滚动条自行设置计算域。在流体分析设置中，我们主要设计流

2020-06-02 00:23:11 23389 2

原创接触面上的压力总和

单面接触用于当一个物体的外表面与自身接触或和另一个物体的外表面接触时使用。单面接触是 ANSYS/LS-DYNA 中最通用的接触类型，因为程序将搜索模型中的所有外表面，检查其间是否相互发生穿透。由于所有的外表面都在搜索范围内, 不需要定义接触面与目标面。在预先不知接触情况时，单面接触非常有用。相对于 ANSYS 隐式分析, ANSYS/LS-DYNA 的单面接触不会非常耗时。大多数冲击与碰撞问题需要定义单面接触当接触面之间的穿透超过接触单元厚度的 40%时，单面接触自动释放接触，从而对下列问题造...

2020-05-25 11:19:33 1773 1

原创机器学习笔记对sklearn的线性回归器进行学习

#-*-coding:utf-8-*-from sklearn.datasets import load_bostonfrom sklearn.model_selection import train_test_splitimport numpy as npfrom sklearn.preprocessing import StandardScalerfrom sklearn.linear_model import LinearRegressionfrom sklearn.linear_mo.

2020-05-25 11:18:56 467

原创 python字符串拼接、查找、分割、替换

1.字符串定义字符串：以双引号或单引号包围的数据2.字符串的拼接练习：a = 'hello'b = 'world'c = a + bprint(c)运行结果：helloworld3.字符串的查找方法 count 计数功能，返回自定字符在字符串中的个数 find 查找，返回从左第一个指定字符的索引，找不到返回-1 rfind 查找，返回从右第一个指定字符的索引

2020-05-25 11:18:03 417

原创 ABAQUS切削模拟

步骤：1、模型创建（岩石模型上的横线主要是为了划分网格用）2、此处省略材料设置3、装配将零部件装配到如图所示位置，并在刀尖处设置参考点4、网格设置岩石全局种子大小设为2，在岩石被切削部分布置较小单元的种子，刀具与岩石的接触面设置局部种子。指派单元累心为Explicit-温度-位移耦合，单元以四边形为主。5、分析步设置时间长度设置为0.01。场输出请求勾选STATUS和作用力与反作用力。历程输出请求勾选作用力/反作用力下的RF，反作用力和扭矩，...

2020-05-25 11:17:19 6444 1

原创地层三压力预测

# -*- coding: UTF-8 -*-import mathimport numpy as npimport pandas as pdimport matplotlib.pyplot as pltfrom pylab import mplmpl.rcParams['font.sans-serif'] = ['SimHei']mpl.rcParams['axes.unicode_minus'] = False# ==================================...

2020-05-25 11:16:02 3857 6

原创一个圆形水池的静力分析

APDL:/prep7 !------------------------------------------- ET,1,SHELL63 !壳单元 MP,EX,1,2.1e5 !C20 混凝土 MP,NUXY,1,0.2 R,1,200, , , , , , !壁厚 200 R,2,150, , , , , , !肋厚 150 !------------------------------------------- !建模 k,1,0,0,0 k,2,325,0,0 k,3,32.

2020-05-18 14:27:39 865

原创 pyhton_批量文本整合进单一文本

将多个文本整合到单个文本中方便后续进行数据处理#-*-coding:utf-8-*-import os,shutilfor i in range(1,10): file_list=os.listdir("word_count\\"+str(i)) for file in file_list: with open("word_count\\"+str(i)+'\\'+file,'r',encoding='utf-8',errors='ignore') as f:

2020-05-18 14:27:04 145

原创 Python条件语句

Python条件语句是通过一条或多条语句的执行结果（True或者False）来决定执行的代码块。Python程序语言指定任何非0和非空（null）值为true，0 或者 null为false。Python 编程中 if 语句用于控制程序的执行，基本形式为：if 判断条件：执行语句……else：执行语句……其中"判断条件"成立时（非零），则执行后面的语句，而执行内容可以多行，以缩进来区分表示同一范围。else 为可选语句，当需要在条件不成立时执行内容则可以执行相关语句。...

2020-05-18 14:26:05 355

原创 Python符串格式化输出

print()函数使用以%开头的转换说明符对各种类型的数据进行格式化输出，具体请看下表。转换说明符只是一个占位符，它会被后面表达式（变量、常量、数字、字符串、加减乘除等各种形式）的值代替。【实例】输出一个整数：age=8print("C语言中文网已经%d岁了！"%age)运行结果：C语言中文网已经8岁了！在print()函数中，由引号包围的是格式化字符串，它相当于一个字符串模板，可以放置一些转换说明符（占位符）。本例的格式化字符串中包含一个%s说明符，它最终会被后..

2020-05-18 14:25:20 428

原创 matlab判断图中距离几个点总距离最小的点的坐标

本方法可用于寻找距离图中几个点总距离最小的点的坐标，在储运联合站位置选择等方面可以应用。下面附详细代码：clear;clc;x=[39496497 39498814 39495320 39498306 39501067];y=[4218695 4218832 4218558 4217348 4217954];dist=@(var) sum(sqrt((var(1)-x).^2+(var(2)-y).^2));%var(1)=x;var(2)=yvar0=rand(2,1);[var,minD.

2020-05-18 14:24:37 6367

原创 python双增式轨道设计

import pandas as pdimport numpy as npfrom numpy import sinfrom numpy import cosfrom numpy import tanfrom numpy import arctanfrom numpy import arccosimport matplotlib.pyplot as pltfrom pylab import mplfrom mpl_toolkits.mplot3d import Axes3Dmpl....

2020-05-18 14:23:43 536 1

原创循环对称结构模态分析

基本过程： 1、建模 2、define cyclic 3、定义约束 4、求模态解 5、展开并察看结果 APDL:r1=5 !建模 r2=10 d1=2 nsect=24 alpha_deg=360/nsect alpha_rad=2*acos(-1)/nsect /prep7 csys,1 k,1,0,0,0 k,2,0,0,d1 k,3,r1,0,0 k,4,r1,0,d1 l,3,4 arotat,1,,,,,,1,2,alpha_deg/2 k,7,r2,0.

2020-05-10 18:51:14 1261

原创关于聚类分析方法的一些探讨

在之前的文章中，我们介绍了一种利用岩石抗压强度进行PDC钻头选型的软件，软件需要我们输入岩石的单轴或三轴抗压强度以及层段深度。由于层段的深度对结果产生很大的影响，因此采用了聚类分析的方法对得到的数据进行了处理。本文奖简单介绍聚类分析的基本概念和应用方法。以PDC选型软件为例，需要输入层段深度和对应参数，从而得到整个地层的参数。聚类分析的基本思想是在不破坏原有顺序条件下将N个样品分别分割为L段，设Xkj表示第k（k=1，2。。。N）段第j（j=1，2，3。。。N）段样品的指标，表示第k段指标的平均值，.

2020-05-10 18:50:19 302

matlab 俄罗斯方块游戏

matlab 地球 月球运动

空空如也

matlab 地球月球运动