TigerTai98-优快云博客

原创 Numpy 向量化操作加速计算，解决向量外积求和慢的问题

在计算函数二阶导的时候，我们通常需要计算Hessian 矩阵，此时会遇到向量外积求和的问题∑i=1nuivi,ui∈Rn×1,vi∈R1×n\sum_{i=1}^n u_i v_i, u_i \in \mathbb{R}^{n\times 1}, v_i \in \mathbb{R}^{1\times n}i=1∑nuivi,ui∈Rn×1,vi∈R1×n如果我们一个一个地计算 uiviTu_i v_i^TuiviT 并将它们用for循环加起来，那么我们需要计算 nnn 个矩阵相乘，是

2021-05-08 02:13:35 731 8

原创 Numpy 通过矩阵操作避免for循环之 [None, :, :]运用

我们知道python中 for 循环的效率是不高的，相比之下，直接运用矩阵向量操作（如点乘）可以让代码运行更快。因为在复现paper时候遇到了这样的问题，在这里总结一下我的解决方法。[None, :, :] 介绍通过None indexing，可以将一个二维矩阵扩展成一个三维矩阵，算是一种 broadcasting吧。一个大小为[a, b]的矩阵M，在A=M[None, :, :]操作下，A的...

2019-01-21 14:55:15 5198

原创一个行列式求导

一个行列式求导公式d∣A∣dt=∣A∣tr(A−1∗dAdt),&nbsp;A∈Rn×n\frac{d|A|}{dt} = |A|tr(A^{-1}*\frac{dA}{dt}),\ A\in R^{n \times n}dtd∣A∣=∣A∣tr(A−1∗dtdA),&nbsp;A∈Rn×n证明如下首先我们有∣A(aij+ϵ)∣−∣A(aij)∣=ϵAijd∣A∣daij=Aij...

2018-11-01 17:02:57 11652 2

原创 Tensorflow: Model parallelism 模型并行计算

在tensorflow官方tutorial上给出了多GPU的用法，但那是基于data-parallelism的计算，主要思想是将数据划分成不同部分，用同一个模型进行计算但是我在写代码中发现，会出现单个模型过大无法再单个GPU上运行，这时候就需要model-parallelism上网查找了很多资料后，发现这个博主写的不错，附带了github代码，How to Use Distributed ...

2018-07-20 04:58:25 5533 7

原创编程作业三

[算法导论] 编程作业三:矩阵链乘法与最大子数组PB15000105 肖映泰一、作业任务1.采⽤用动态规划算法编程实现矩阵连乘的最优括号化方案，并输出最优括号化方案2.采用分治法和贪心算法实现最大子数组二、算法原理矩阵链乘法给定n个矩阵的链<A1,A2,...,An><A1,A2,...,An>，矩阵的规模为pi−1×pi(1≤i≤n...

2018-04-29 11:41:57 487

原创 Latex中文编辑

Latex中文编辑只需要将汉字放在CJK环境中就可以使用中文编辑，记得导入CJK宏包\documentclass{11pt}{article}\usepackage{CJK}\begin{document}\begin{CJK}{GBK}{song}恭喜编译成功\end{CJK}\end{document}...

2018-04-20 11:50:10 1605

原创极小曲面

极小曲面肖映泰 PB15000105一、算法定义：平均曲率处处为0 H(vi)=0, ∀iH(vi)=0, ∀iH(v_i)=0,\ \forall i 平均曲率 limlen(γ)→01len(γ)∫v∈γ(vi−v)ds=H(vi)nilimlen(γ)→01len(γ)∫v∈γ(vi−v)ds=H(vi)ni\lim_{len(\gamma)\to0}\...

2018-04-14 13:09:06 7692 6

原创数学建模作业一

人口模型PB15000105 肖映泰深圳在改革开放后迎来了巨大的发展，本人从深圳统计局获取了深圳市1979-2016年来的人口数，希望通过建立适当的模型预测深圳市接下来的人口增长情况，用以给出合理的政策建议。这里我们采用第三列数据，即分析总人口变化规律下面主要从两个模型入手分析和预测人口变化情况一、指数增长模型——Malthus模型 1.模型假设基本假设：人口的自...

2018-03-20 17:19:18 1532 1

翻译 NormFace

NormFaceNormFace: L2 HyperSphere Embedding for Face Verificationhttp://www.cs.jhu.edu/~alanlab/Pubs17/wang2017normface.pdf Baseline Model Original Accuracy Finetune Using Normalization...

2018-02-20 10:22:09 919

原创 1.Tensorflow安装--基于Windows

Tensorflow安装–基于Windows现在tensorflow支持windows下3.5以上的python安装，这里我们采用python3.5.2安装1.安装Python3.5.2 python3.5.2安装地址在上面的网站上找到python3.5.2版本，选一个链接进行安装，我选的是Download Windows x86-64 executable installer2...

2018-02-09 12:04:15 271

原创 1.Tensorflow 安装 --基于 Ubuntu

Tensorflow安装–基于Ubuntu最开始阶段，我们需要安装Tensorflow，而在Ubntu系统下，安装Tensorflow是很简单的 1.安装pip $ sudo apt-get install python-pip python-dev2.安装virtualenv虚拟环境$ sudo pip install virtualenv --upgrade3.创建新...

2018-02-08 22:35:14 287

原创 eigenface

eigenface实现人脸识别采用PCA方法，利用YALE的人脸数据库做人脸识别，15个人，165张图片，每张图片大小为100*100这里yale文件中是源图片，通过changename.py将名字更改，存放到picture文件下eigenface具体步骤 1.读取训练所需要的图片，这里使用subjectxx.normal.bmp文件，共15张图片，每张图片拉直为一个（1，1

2018-01-27 16:22:15 1662 1

原创 windows10 dlib19 andaconda python3.6

非常简洁的安装dlib19.07环境：WINDOWS10 anaconda py3.6 准备材料：需要将anaconda 的pip 路径加入环境变量（注意需要排除windows本身含有的python 的pip的影响，这个可以通过改名等方式实现，总之，注意将pip安装到anaconda目录下)资源下载：这里是wheel文件的下载，直接pip install就可以安装，非常方便。点击下载whl文件

2017-10-30 14:37:37 731

翻译 python 列表操作

Think Python 列表相关操作# -*- coding: utf-8 -*-"""Created on Fri Oct 27 22:05:42 2017@author: 123"""cheeses=['Cheddar','Edam','Gouda']numbers=[42,123]empty=[]print(cheeses,numbers,empty)#list支持赋值，但是字符

2017-10-27 22:54:13 298

原创 python pandas基础

读取文档#导入csv函数food_info=pandas.read_csv("food_info.csv")print(type(food_info))print(food_info.dtypes)基本操作#帮助文档print(help(pandas.read_csv))food_info.head(3)#前三行数据food_info.tail(5)#后5行数据food_info.c

2017-08-07 00:23:54 401

原创 python numpy基础

基本类型import numpy#向量vector=numpy.array([5,10,15,20])vector==10#矩阵matrix=numpy.array([[5,10,15],[20,25,30],[35,40,45]])matrix==25类型转化astype函数 dtype函数vector1=numpy.array(["1","2","3"])print(vector

2017-08-07 00:14:00 465

翻译 python 实现 kNN 算法

kNN算法1.准备带有标签的数据（需要转换成矩阵格式） 2.提供一个待分类的数据 3.求待分类数据与每一个带标签的数据的”距离“ 4.选取“距离”最小的k个带标签数据 5.求k个数据中出现频率最高的数据得出分类结果python 代码实现from numpy import *import operator#inX为待分类数据，dataSet为数据，labels为标签def classify

2017-08-03 19:15:24 666

原创 Python 第5章条件和递归

向下取整和求余数hour=105//60 ##hour=1minute=105%60 ##minute=45键盘输入text=input();name=input("What is your name?\n");两种递归##先后顺序不同，输出顺序不同def print_up(n): if (n < 0): return print_up(n-1) pr

2017-06-24 18:20:45 735

转载机器学习：python 实现一个linear regression

1.原理介绍linear regression步骤： 1.导入数据 2.将数据分为训练集合测试集（linear regression 分为x_train, x_text, y_train, y_test） 3.导入线性回归算法利用训练集计算出模型参数 4.模型检验利用测试集测试真实值和预测值的差异（用x_test计算出y_predict，与y_test做比较，计算误差） 5

2017-06-09 12:14:16 10860

翻译 MFC基本绘图函数（1）

简单数据类型Windows POINT,RECT,SIZE结构体封装：typedef struct tagPOINT{ LONG x; LONG y;}POINT,*PPOINT;typedef struct tagSIZE{ LONG cx; LONG cy;}SIZE,*PSIZE;typedef struct _RECT{ LONG le

2017-06-08 15:46:56 1336

翻译 2.1面向对象程序设计基础（3）

1.派生类的定义格式：class 派生类名：[继承方式]基类名{ 派生类新增的成员;}继承方式：private public protected 注意：派生类继承基类的成员函数，但不继承构造函数2.派生类构造函数的定义派生类构造函数名（总参数列表）:基类构造函数名（参数列表）{派生类中新增成员变量初始化语句}注意：若没有基类构造函数，则按默认构造函数初始化基类的变量例子：#inclu

2017-06-08 08:35:18 486

翻译 2.1面向对象程序设计基础（2）

2.1.3对象的动态建立与释放全局数据区：全局变量、静态数据、常量代码区：类成员函数、非成员函数代码栈区：局部变量、函数参数、返回数据、返回地址堆区：动态存储分配，建立后一定要释放（否侧会出现内存泄漏）建立对象格式： CRectangle *pRect=new CRectangle; 或者： CRectangle *pRect=new CRectangle(40,10); ——-

2017-06-07 21:43:09 403

翻译 2.1面向对象程序设计基础

2.1.1类和对象1.类的定义格式： class 类名{ public: ... private: ... protected: ...};2.类成员函数定义格式：void Picture::Chart(int i){ ...}//::是作用域运算符3.对象格式：Picture circle(0.5,(0,0));2.1.2构

2017-06-07 18:08:59 590