- 博客(30)
- 资源 (2)
- 问答 (3)
- 收藏
- 关注
RSS订阅原创 PID控制
PV：process variableSP：setpointdifference：也称error，是PV与SP之间的差别。目的：给定SP，使得PV很快调整到SP且保持稳定。绝对误差：SP - PV，即当前误差，对应比例单元所有时间的误差和：即当前以及之前的所有误差之和，对应积分单元死区时间：指系统的延迟，给烤箱指定温度，烤箱需要一段时间才能到达指定的温度比例单元(P)、积分单元(I)和微分单元(D)分别对应目前误差、过去累计误差及未来误差。https://web.archive.org.
2021-01-28 11:00:11
2374
1
原创 一稿多投
一稿多投一稿多发https://www.zhihu.com/question/404961473http://muchong.com/html/201204/4371264.html中英文都投
2021-01-14 19:56:29
571
1
原创 如何写论文
https://wenku.baidu.com/view/c62761a40029bd64783e2cfe.html写论文、投稿、改稿新：新问题、新方法、新结果提纲outline标题title——摘要abstract——引言introduction——正文main body——结论conclusion——致谢acknowledgement——文献refernces——附录appendix思想上的洗礼读书学问，本欲开心明月，利于引耳奥古斯特Henry Augustus .
2021-01-14 11:05:40
172
原创 数值分析库
https://en.wikipedia.org/wiki/List_of_numerical_libraries#Fortran每种语言都有对应的数值分析库FortranMINPACK：函数最小化、最小二乘问题，2001年最后更新EISPACK：特征值与特征向量、奇异值分解，2001年最后更新LINPACK：线性方程组、最小二乘问题、奇异值分解，2001年最后更新LAPACK：linear algebra package，线性方程组、最小二乘问题、奇异值分解、特征值与特征向量，20
2021-01-13 15:02:05
202
2
原创 梯度投影算法
也称GLP投影法投影向量的线性变换。原向量乘以投影矩阵变换成新向量。物体v，物体的影子v'，投影矩阵PPv=v'PPv=Pv'=v'投影矩阵满足 PPv=Pv，投影两次与投影一次的效果一样U，V为P的像空间与零空间1. 任意x∈U，P(x) = x，也就是投影两次与投影一次效果一样。2. 整个向量空间W为U与V两个子向量空间的直和，直和就是矩阵对应元素相加。x∈W，x可写作两个向量之和 x = u + v，满足u = Px，v = (E - P)x正交投影正交投影：..
2021-01-12 21:48:25
12614
4
原创 Quadratic programming
Quadratic programming原问题：minf(x)如果Q是半正定矩阵，那么f(x)是一个凸函数。相应的二次规划为凸二次规划问题；此时若约束条件定义的可行域不为空，且目标函数在此可行域有下界，则该问题有全局最小值。如果Q是正定矩阵，则该问题有唯一的全局最小值。若Q为非正定矩阵，则目标函数是有多个平稳点和局部极小点的NP难问题。如果Q=0，二次规划问题就变成线性规划问题。当二次规划只有等式约束时可以用线性方程求解。否则，常用的解法有：有效集法（active set me
2021-01-12 16:02:47
392
2
原创 Givens旋转与Householder反射
参考wiki如何理解Householder变换和Givens旋转? - 荆哲的回答 - 知乎 https://www.zhihu.com/question/51432776/answer/125819331Givens旋转矩阵Gii = cosθ，Gjj = cosθ，Gij = sinθ，Gji = -sinθ，Gkk = 1（k != i,j）乘积 G(i, j, θ)x 表示向量 x 在 (i,j)平面中的逆时针旋转 θ 弧度。旋转：旋转矩阵乘以一个向量，不改变向量的大小，改
2021-01-12 10:32:49
1423
3
原创 奇异值
https://www.cnblogs.com/LeftNotEasy/archive/2011/01/19/svd-and-applications.htmlhttps://blog.csdn.net/dss875914213/article/details/84678957定义特征值分解奇异值分解Σ对角元素即为奇异值酉矩阵与正交矩阵U为方阵，前者为正交矩阵在实数域，后者为酉矩阵在复数域...
2021-01-11 21:25:55
1015
原创 governing equation
质量守恒方程动量守恒方程角动量守恒方程能量守恒方程熵守恒方程麦克斯韦方程组麦克斯韦——法拉第感应电场方程，说明感应电场产生磁场安培——麦克斯韦感应磁场方程，说明感应磁场产生电场高斯电通量方程，说明电荷产生电场高斯磁通量方程，说明不存在磁荷连续介质力学中的控制方程板理论（Kirchhoff—Love板理论，Mindlin—Reissner板理论）铁木辛柯梁理论涡流脱落非定常流的有限体积法声学理论大涡模拟terminology explanation状态方程state
2021-01-10 21:16:22
874
原创 显式与隐式
显式与隐式explicit与implicit区别：当前的量由之前的量计算得到；当前的量由之前和当前的量计算得到。表达式——显式——隐式从中可看出，隐式还多了一步，解方程。还有两者结合的显式-隐式法举例——前向欧拉法、反向欧拉法、CN法 前向欧拉法即式1反向欧拉法即可以看出有两个根，但因为初始条件为正，取正根式2Crank Nicolson法即举例——前向-反向欧拉法前向-反向欧拉法即CFL条件偏微分方程收敛的条件，现在多用于差分方程收敛的条件.
2021-01-10 20:05:00
4906
2
原创 ansys求解瞬态动力学
想结识做仿真（固体力学、流体力学、电磁学、热等等）的朋友！瞬态动力学方程求解方法：向前差分时间积分法、NewMark时间积分法(改进后的HHT时间积分法）。前者用于显式瞬态动力学分析，LS-DYNA模块。后者用于隐式瞬态动力学分析。只说NewMark法。显式与隐式NewMark法使用NewMark法求解当前位移Un+1将NewMark法的两个方程改写为由三个方程得出式1用前一次求解的位移、速度、加速度可以推出当前的位移。从方程求解得到位移，然后用NewMark法的两
2021-01-10 17:35:48
3088
1
原创 ansys求解特征值问题
来自理论手册15.14 缩减 子空间 分块兰索斯 非对称 阻尼 QR阻尼 命令 modopt,reduc modopt,subsp modopt,lanb modopt,unsym modopt,damp modopt,qrdamp 不推荐求解屈曲 对称矩阵 对称矩阵，不能用于屈曲 非对称矩阵 对称或者非对称的阻尼矩阵 对称或者非对称的阻尼矩阵 M C K ..
2021-01-09 19:02:31
1680
原创 CPU与GPU
CPU GPU 中文名称 中央处理器 图形处理器 专注于 低延迟 高吞吐量 擅长 处理串行指令 处理并行指令 核心 强大的少量核心 弱的很多核心 特点 控制器强大可以预测执行与乱序执行 架构允许内存的延迟 速度 高 比CPU都高 内存消耗 高 低 CPU由算术逻辑单元ALU与控制逻辑单元CU组成。ALU进行计算与临时存储，CU进行控制可以计算与控制CPU只有几个Die...
2021-01-09 14:28:10
187
原创 拟牛顿法
牛顿法及其扩展在当前迭代点用目标函数二阶展开近似目标函数拟牛顿法海森矩阵近似矩阵的B的求取方1方2参考博客《Numerical Optimization Ch7. Large-Scale Unconstrained Optimization》https://www.cnblogs.com/kisetsu/p/9157371.html...
2021-01-06 21:35:01
101
1
原创 线性共轭方向法
能用最简单的话描绘出你所说的问题才是好样的，现实中往往怕嘲笑或为了装高级，而将描绘复杂化。工科生不需要学理论？大错特错。线性共轭方向法求解线性方程组问题，在数值计算中有对应的解法，在最优化中将其转化为一个最优问题。求解转化为求解初始点x0，初始方向d01.计算步长αk，αk = argmin{f(xk+α*dk)}，α>=02.xk+1 = xk + αk*dk3.计算方向dk+1，使其与dk,...d0关于A共轭4.k=k+1算法当然有终止计算的规则，我没写。共轭.
2021-01-06 20:14:27
702
原创 非线性规划
整数规划决策变量只能取整数的规划问题动态规划按时间或空间将多阶段决策问题划分为多个阶段，对每个阶段进行决策，前一次决策结果影响当前决策，当前决策影响下一次决策。非线性规划目标函数、约束条件至少一个是决策变量的非线性函数可行解：满足约束条件的解可行域：约束条件构成的区域非线性规划问题包括无约束问题与有约束问题。非线性规划问题引入了不等式约束，标志现代规划理论的开始。无约束最优化变量轮换法最速下降法牛顿法共轭梯度法变尺度法有约束最优化Zoutend
2021-01-05 21:46:08
1528
1
原创 线性规划
定义：目标函数、约束条件为决策变量的线性函数最简单的求解方式就是图解法最优问题用数学形式表达出来为规划问题标准形式的数学表示求和形式矩阵形式向量形式将非标准形式转化为标准形式求最小，加负号变为求最大小于等于不等式约束，左端加松弛变量变为等式约束大于等于不等式约束，左端加剩余变量变为等式约束b小于0，Ax=b两端同时乘-1x可正可负，令x=x'-x''，两个新变量x',x''>=0解的概念可行解：满足约束条件Ax=b与非负约束条件x>=0的解
2021-01-05 19:29:54
925
原创 拉格朗日乘子法、惩罚函数法
卡罗需－库恩－塔克条件Karush-Kuhn-Tucker Conditions，KKT条件，KT条件一个非线性规划（Nonlinear Programming）问题能有最优化解法的条件。最优解要满足以下条件——必要引入了乘子如果满足以下条件则是最优解——充分目标函数、约束条件都是凸函数仿射变换向量间的线性加平移变换拉格朗日乘数法也称拉格朗日乘子法约束问题只有等式约束条件，则可通过乘子将约束问题转化为无约束问题乘子的求解增广拉格朗日惩罚函数法
2021-01-05 15:49:53
3186
原创 big O与small o，以及复杂度
P问题、NP问题、NP完全问题、NP困难问题这些问题就是从时间复杂度为多项式时间出发定义的时间复杂度时间复杂度分类：多项式时间、常数时间、对数时间等等。时间复杂度用big O表示空间复杂度big O又称为渐进符号，是用于描述函数渐进行为的数学符号无穷大渐进T(n)趋于无穷大时，n^2占主导地位分析算法效率时，n为输入的数据大小，T为计算时间，当n越来越大时，T≈n^2，称该算法具有n^2阶的时间复杂度无穷小渐进samll o无穷小的比阶：高阶无穷小、
2021-01-05 12:13:57
1219
原创 特征向量的研究
一个特征值对应几个特征向量？特征值非重根时，该特征值对应一个特征向量特征值为n重根时，该特征值最多对应n个特征向量不能超过n，所以也有可能只对应一个特征向量特征向量间的研究特征向量结果不唯一，那么在你和我计算都对的情况下，我们的特征向量之间成什么关系？研究3x3矩阵，秩为1、2、3的情况时，求解的特征向量+++++++++++++++++++++++++++++++++++++++++++++++++++++++此例中rank(A)=1++++++++++++++++++++++
2021-01-05 10:56:13
3484
原创 代理模型
代理模型仿真计算时不使用流体力学、固体力学的知识进行有限元、有限体积仿真，而是通过事先的训练数据（输入与输出）建立起代理模型（黑匣子），给定代理模型输入则会得到输出。最重用的代理模型有多项式响应曲面法，克里金法，梯度增强克里金法（GEK），支持向量机，空间映射，和人工神经网络。响应曲面法独立变数X与响应Y1.找一个最优模型——响应曲面设计一阶回归模型二阶回归模型2.通过最优化方法求解——响应曲面最佳化...
2021-01-04 20:50:57
4794
原创 线程、进程的概念（二）
在看代理模型就看到了代理模式内存管理引用计数：是指将资源的被引用次数保存起来，当被引用次数变为零时就将其释放的过程。垃圾回收：当某个程序占用的一部分内存空间不再被这个程序访问时，这个程序会借助垃圾回收算法向操作系统归还这部分内存空间。垃圾回收可以通过引用计数来实现，垃圾回收就是一种内存管理方式。内存管理，是指软件运行时对计算机内存资源的分配和使用的技术。其最主要的目的是如何高效、快速的分配，并且在适当的时候释放和回收内存资源。 可以想像一个很大的程序，当他运行时变成进程，而大部分的存储器空间
2021-01-04 20:36:18
79
原创 structure dynamic modify
结构动力学修改在design process时，考虑结构强度、刚度等等，不考虑振动与噪声。在prototype stage时，要考虑振动与噪声。（我没搞懂这个设计与原型是什么意思。）SDM问题：更正模型，使得振动与噪声很小。动力学方程 MX''+CX'+KX=F，修改M C K。动力学方程X=H*F，修改频响函数H。通过修改M C K 后，会得到新模型的模态。可计算响应，从中选出最优的M、C、K，即为最终的模型。只要知道原模型的模态与ΔM和ΔK，就可以求出新模型的模态。利用模态叠..
2021-01-03 12:51:14
313
原创 threading模块
推荐使用KiteThread多线程会有一个主线程与多个子线程。子线程就是通过Thread创建的线程。默认没有守护线程，当主线程执行结束后退出，不影响子线程的执行。设置子线程为守护线程时，主线程一旦执行结束，不论子线程是否执行完，直接终止子线程执行。在这种情况下，如果设置子线程同步，主线程会被阻塞，等到所有子线程执行结束，主线程才会执行结束。start启动线程joint同步线程，会阻塞主线程直到所有子线程执行完，如果设置了时间，超时后停止阻塞。setDaemon守护线程如果有需求，还可以
2021-01-02 21:47:23
123
原创 突破GIL
IMOOC，讲师链接https://www.imooc.com/t/5683545如何实现单进程多线程真正的并行1. 使用其他解释器2. 用C重写python线程创建的过程，就可以实现单进程多线程真正达到并行反之目前来看，Cpython的GIL仍将存在。双核CPU进行三种对比老师在虚拟机中的Ubuntu系统下运行三种情形：单进程多线程多进程调用c语言写的可执行文件具体说一下第三个：编译好c语言后可得到可执行文件，通过ctypes运行该可执行文件。可执行文件没有GIL。结
2021-01-02 11:14:58
132
原创 GIL的意义
只要是单进程不论是单还是多线程，无论是单核还是多核CPU，都是每次执行一个线程，前者是串行，后者是并发。只要是多进程就可以在多核CPU上并行，当然如果在单核CPU还是并发。为何不取消GIL？CPU bound：计算密集，CPU一直运作。使用多进程。IO bound：输入输出密集，CPU大部分时间是闲置，比如内存IO，网络IO，磁盘IO，sleep函数。使用单进程多线程与多进程单线程相比更省资源。补充1计算机有五部分：控制器、运算器，对应CPU存储器，对应内存条与硬盘输入设备，对
2021-01-02 10:17:23
108
原创 GIL
anaconda：包含conda与python以及一堆包。conda是一个IDE。miniconda：包含conda与python。pycharm：是一个IDE，所以安装pycharm后还需要安装python。
2021-01-01 21:34:48
114
4
原创 编译器与解释器
anaconda与pycharmanaconda：包含conda与python以及一堆包。conda是一个IDE。miniconda：包含conda与python。pycharm：是一个IDE，所以安装pycharm后还需要安装python。虚拟机python虚拟机：PVMjava虚拟机：JVM和VMware这种软件类似。编译器与解释器（我认为）原代码→编译器→ 机器码（可执行文件）→运行→结果原代码→解释器→ 字节码（pyc文件）→ 虚拟机→结果编译器优点：编译.
2021-01-01 18:31:41
228
原创 线程、进程的概念（一）
线程thread操作系统能够运算调度的最小单位。由操作系统内核调度的内核线程与由用户进程自行调度的用户线程或者混合调度。进程process进程是程序的实例，若干进程对应一个程序，进程被加载到存储器中。每个进程皆可以同步（循序）或异步（平行）的方式独立运行。在一个处理器上通过分时复用表现出平行运行的感觉，在多个处理器上真正平行运行。时分多路复用也称时间共享，分时复用。两个以上的信号轮流占用通道。多线程、多进程运行一个程序会有若干进程，每个进程会有若干线程。软件多线程：只有
2021-01-01 14:51:13
155
原创 碰撞
碰撞在物理学中表现为两粒子或物体间极短的相互作用。 碰撞前后参与物发生速度、动量、能量改变。分类弹性碰撞是碰撞前后整个系统的动能不变的碰撞。弹性碰撞的必要条件是动能没有转成其他形式的能量。非弹性碰撞是碰撞后整个系统的部分动能转换成至少其中一碰撞物的内能，使整个系统的动能无法守恒。弹性碰撞物体A与B相同质量A速度V，B速度0，碰撞后，A速度0，B速度V。A速度V，B速度-V/2，碰撞后，A速度-V/2，B速度V。A速度V，B速度V/2，碰撞后，A速度V/2，B速度V。可以看出相同
2021-01-01 10:09:32
461
masm link debug dosbox
2018-06-03
深度学习英伟达显卡配置
2023-09-07
问一道关于信号处理的问题，帮我解答一下？bai'tuo
2019-10-27
pyinstaller打包后exe不能运行
2019-10-06
不知道如何通过浏览器访问ubuntu sever
2017-08-16
TA创建的收藏夹 TA关注的收藏夹
TA关注的人