- 博客(23)
- 资源 (1)
- 收藏
- 关注
RSS订阅
原创 浅水方程的的精确黎曼求解器——推导浅水方程（一）
一、浅水方程二维浅水方程可根据无粘N-S方程沿水深方向积分可得到，主要运用莱布尼茨公式。运用公式后进行合并，图中给出三维的合并：浅水方程的连续方程到此就推出了将欧拉方程（无粘N-S方程）中的压力梯度项做出假设：压强分布满足静水压强分布进一步得到简化公式：这里的s代表水位s=z将连续方程两边乘以u，动量方程乘以h：合并便可以得到u方向的动量方程：同理v方向至此，二维...
2018-11-16 11:29:11
9122
转载 神经网络学习——第一个代码学习
学习的博客来自：人工神经网络——反向传播算法(BP)以及Python实现一共分为两个py文件1、BPNN"""输入层：3隐藏层：3输出层：2"""import numpy as npimport mathimport matplotlib.pyplot as plt# 激励函数及相应导数def sigmoid(x): return 1.0 / (1 + np.ex...
2020-04-03 11:28:26
480
转载 支持向量机通俗导论——学习（二）
优化后的SMO算法：步骤1：先“扫描”所有乘子，把第一个违反KKT条件的作为更新对象，令为a1；步骤2：在所有不违反KKT条件的乘子中，选择使|E1 −E2|最大的a2进行更新，使得能最大限度增大目标函数的值（类似于梯度下降. 此外，而，求出来的E代表函数ui对输入xi的预测值与真实输出类标记yi之差）。最后，每次更新完两个乘子的优化后，都需要再重新计算b，及对应的Ei值。————————...
2020-03-20 20:27:31
227
原创 Hec-GeoRAS心得
Hec-GeoRAS作为Hec-RAS的前期地理处理软件，和GIS联合使用，方便不少。Geo可以为RAS提供所有的地形数据，直接输出RAS可用的文件，减少制作文件的麻烦以及缺少河道数据时的造假。利用dem可得到相对正确的信息。
2020-03-01 14:47:48
3009
1
原创 CAD等高线转地形高程散点
日常工程设计中经常拿到CAD等高线，可MIKE、SMS水利软件需要地形高程散点书，目前有南方CASS一键提取，也有Arcgis提取法。每每均是叫甲方用南方CASS转化后给我，今第一次用Arcgis操作，收货颇多，这里分享利用Arcgis的提取方法，该方法是由多人努力完成，自己再多次尝试，现记录以便日后所用以及他人所用。1、方法一：打开ArCcatalog工具，连接到cad文件夹打开cad文件，可...
2019-11-30 00:02:48
12376
原创 MIKE与SMS网格的区别——个人感受
在商业软件盛行的水利行业中，以MIKE、SMS、HEC等系列软件最为著名。软件简介就不用描述了，可自行百度和论文查找，这里主要想说明下MIKE和SMS在二维网格生成的特点，仅作为学习过程中的笔记，侵权即删。 本文根据两款软件对连续弯道进行网格划分，弯道是具有比较好的说服力。
2019-06-11 22:12:42
8871
7
原创 一维扩散方程差分格式的数值计算
摘要：采用FTCS显格式、Dufort-Frankel显格式、Lassonen隐格式、Grank-Nicolson隐格式对一维扩散方程进行数值计算，得到不同时间y方向扩散的速度。结果表明了不同差分格式的差别。关键词：扩散方程；差分格式；FORTRAN；module global implicit none real::u(1000,1000)end module prog...
2018-12-24 16:19:12
12963
5
原创 黎曼不变量——特征法的学习
一、特征法参考这个网络讲义1.苏剑林. (2017, Dec 07). 《一阶偏微分方程的特征线法 》[Blog post]. Retrieved from https://kexue.fm/archives/47182.知乎解答：https://www.zhihu.com/question/40518374二、黎曼不变量在非齐次项为零时，变量沿着特征线不变，故称为黎曼不变量:dw1/r1...
2018-12-06 11:43:26
4420
1
转载 一维Burgers方程的学习——来自流沙公众号
对二维Burger方程的认识，该方程import numpy as npimport matplotlib.pyplot as pltfrom mpl_toolkits.mplot3d import Axes3Dnx = 41ny = 41nt = 120c = 1dx = 2 / (nx - 1)dy = 2 / (ny - 1)sigma = 0.0009nu = ...
2018-10-17 22:53:34
2148
转载 对二维扩散方程的学习——来自流沙公众号
代码：import numpy as npimport matplotlib.pyplot as pltfrom mpl_toolkits.mplot3d import Axes3Dnx = 31 # 定义x方向节点数ny = 31 # 定义y方向节点数nt = 17 # 定义时间步长nu = 0.05 # 扩散系数dx = 2 / (nx - 1)dy = 2 / (ny ...
2018-10-16 21:04:40
12009
1
转载 对二维线性对流方程的学习——来自流沙公众号
许久没接触python，又有点忘了继续学习流沙公众号对这个方程，非定常项，以及x和y方向的对流项碰巧今天看安德鲁的计算流体力学偏微分方程所讲，利用克莱姆法则和特征值法来判断方程的属性例子为一阶，在做习题碰到二阶。此为闲话同时，根据文章所写，稍加修改，将for循环和数组放在同一程序中进行对比，以温习python，发现遗忘甚多。from mpl_toolkits.mplot3d imp...
2018-10-16 19:14:52
2955
1
转载 Python网络爬虫——学习心得
网络爬虫，即Web Spider，是一个很形象的名字。把互联网比喻成一个蜘蛛网，那么Spider就是在网上爬来爬去的蜘蛛。网络蜘蛛是通过网页的链接地址来寻找网页的。从网站某一个页面（通常是首页）开始，读取网页的内容，找到在网页中的其它链接地址，然后通过这些链接地址寻找下一个网页，这样一直循环下去，直到把这个网站所有的网页都抓取完为止。如果把整个互联网当成一个网站，那么网络蜘蛛就可以用这个原理...
2018-09-29 21:56:29
8314
转载 对Laplace方程的学习——来自流沙公众号
对Laplace方程的学习代码如下：import numpy as npimport matplotlib.pyplot as pltimport timedx = 0.1dy = 0.1dx2 = dx*dxdy2 = dy*dydef laplace(u): nx,ny = u.shape #查看数组的维数 for i in range(1, nx-1): ...
2018-09-28 22:13:12
1232
转载 一维Burgers方程的学习——来自流沙公众号
伯格斯方程(Burgers equation) 是一个模拟冲击波的传播和反射的非线性偏微分方程。伯格斯方程也称为粘性伯格斯方程，只适用于一个空间维度。
2018-09-27 21:07:16
23646
12
转载 一维扩散方程的学习——来自流沙公众号
来自流沙公众号认识：1、理论：差分离散，较为基础2、第一次接触matplotlib中ion()和ioff()。ion()开启了交互模式,ioff() 感觉像暂停图像同时，查阅其他博客得到：python可视化库matplotlib有两种显示模式：阻塞（block）模式交互（interactive）模式在Python Consol命令行中，默认是交互模式。而在python脚本中，ma...
2018-09-26 21:29:30
11441
5
转载 一维常系数对流方程的学习——来自流沙公众号
一维对流扩散方程的学习——来自流沙公众号认识：1、时间步长越小越靠后移动；网格越小则波的形状越一致,波形失真在减小，引出Courant数：u*dt/dx<sigma；2、对numpy中的ones（）、linspace（）、zeros（）有基础认识；3、再次练习了matplotlib中的pyplot。4、明白了一维方程编程的简单运算，基础import numpy as np# n...
2018-09-25 22:07:35
7975
1
原创 ArcGIS流域提取教程
ArcGIS流域提取教程 ——Hydraulic Science淘宝店 第一步：下载DEM数字高程数据 免费下载地址：http://www.gscloud.cn/ 方法1: 选择数据资源>>DEM数字高程数据>>GDEMV2 30M 分辨率数字高程数据>>根据研究区域的经纬度下载；方法2: 选择高级检索&g...
2018-09-01 13:08:12
23756
8
Exact Riemann Solver for the Euler Equations
2018-09-06
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人