吴一奇-CSDN博客

原创 C++中用宏定义注释

摘要：本文介绍了如何用宏来定义注释，以控制是否需要注释掉cout。1. 起心动念在调试程序的时候，经常需要用cout输出一些内容，以判断程序是否运行正常，以及进行到哪一步了。但是在程序发布的时候，又不需要这些cout的内容，这时候就需要手动把cout内容删除或者注释掉。很繁琐，尤其是你注释完了之后发现还需要调试，就又得解注释。下面提供一种用一条宏定义语句实现注释或解注释所有cout语句的...

2019-04-26 17:32:01 5166 7

原创 3D中的方位与角位移

摘要：本文讲述了3D中描述方位与角位移的方法：矩阵、欧拉角和四元数，以及它们优缺点和适用场景。1. 简介“方位”是指物体的朝向，是一个状态；“旋转”可以把物体从一个方位变到另一个方位，是一个动作；“角位移”是旋转的量。方位和角位移的关系类似于“位置”和“位移”的关系。物体的位置是不能用绝对坐标来描述的，我们必须将其放置在一个坐标系内，通过描述物体相对于参考点的位移来描述物体的位置。例如...

2019-04-21 17:34:27 3661 1

原创几何角度理解叉乘

简介：本文解释了叉乘的计算（行列式）和几何解释之间的关系1.叉乘的计算及几何解释在我学习叉乘（本文只考虑三维向量）的时候，老师教了我们两个东西：叉乘的计算方法：用行列式。设有向量v⃗=(v1,v2,v3),w⃗=(w1,w2,w3),\vec{v}=(v_1, v_2, v_3), \quad \vec{w} = (w_1, w_2, w_3),\quadv=(v1,v2,v3...

2019-04-14 16:46:12 6657

原创精品资源推荐

1. 简介本篇用于记录我在学习过程中发现的比较好的博客、论文、教程等，涉及的范围比较广，也比较杂，我会尽量按类别来分类。本文将持续更新v0.0.1_201810082. 数学类2.1 傅里叶分析傅里叶分析之掐死教程（完整版）更新于2014.06.062.2 MCMC（马尔科夫链蒙特卡洛）MCMC(一)蒙特卡罗方法to be continued…...

2018-10-08 18:20:54 533 2

原创《有限元分析及应用》第5章：变形体力学方程求解的试函数方法的原理

《有限元分析及应用》第5章：变形体力学方程求解的试函数方法的原理5.1 变形体（弹性）力学方程求解方法分类及试函数方法5.1.1 变形体（弹性）力学方程求解方法变形体（弹性）力学方程求解方法主要有两类：直接法：解析法（analytical method）半解析法（semi_inverse method）差分法（finite difference method）间...

2018-04-17 21:36:53 2349

原创第四章 SVD

SVD may be the most important matrix decomposition of all, for both theoretical and computational purposes.4.1 Introduction Theorem 4.1.1 (SVD Theorem) Let A∈Rn×mA∈Rn×mA \in R^{n\times m} ...

2017-07-25 18:10:53 695

原创如何更好的理解矩阵

简介：介绍一些理解矩阵的技巧 1. 把矩阵看成线性变换Ax 就是把向量xx变换成另一个向量。Ax 的结果就是矩阵AA的列向量的线性组合

2017-07-14 11:24:23 635

原创 ATA对称正定的证明

1. 结论 Let A∈Rn×mwithn>m,A∈Rn×mwithn>m,A\in R^{n\times m} \quad with \quad n>m, then matrix ATAATAA^TA is positive semidefinite, If rank(A)=m(i.e.Ahasfullrank),rank(A)=m(i.e.Ahasfullrank),ra...

2017-07-14 10:02:51 13680

原创第三章最小二乘问题

摘要：本篇主要对教材第三章有关最小二乘的重要概念做一总结

2017-07-13 18:15:53 1587

原创 stack around xxx is corrupted

当数组大小定义太小，在使用时越界。就可能出现stack around xxx is corrupted错误。

2017-07-07 21:20:52 942

原创机器学习中训练数据集，交叉验证数据集，测试数据集的作用

1. 简介在Andrew Ng的机器学习教程里，会将给定的数据集分为三部分：训练数据集（training set）、交叉验证数据集（cross validation set）、测试数据集（test set）。三者分别占总数据集的60%、20%、20%。那么这些数据集分别是什么作用呢？

2017-04-14 15:27:18 18718

原创 c语言赋值运算返回值

结论：赋值运算返回左值的引用。见下图：

2017-04-13 21:35:43 12376 11

原创 visual studio技巧总结

本篇文章以visual studio 2017为例，介绍visual studio实用技巧，包括快捷键的使用、一些非常好用的功能、调试技巧等。学会了这些一定会让你的编程事半功倍。包括：Quick Launch、Go To All、Peek Definition等。

2017-04-13 16:27:38 5008

原创 visual studio属性管理器（property manager）上各项的含义

1. 属性管理器以最新版的visual studio 2017为例，新建一个console application。第一步：view -> Property Manager看到的界面如下图所示： 2. 默认项什么含义以Debug | x64的配置为例，1. Microsoft.Cpp.x64.user：系统默认的属性表，全路径为, C:\Users\horswing\Ap

2017-04-11 16:09:51 10448 2

原创矩阵计算一些重要的结论及其证明

对称正定矩阵对角线元素之和等于特征值之和；对称正定矩阵的对角线元素都大于0；对称正定矩阵的特征值都大于0

2017-04-10 19:40:04 4018

原创区分指针常量和常量指针的技巧

一个技巧帮你彻底搞定指针常量和常量指针，再复杂也不怕。

2017-01-21 23:00:25 570

原创如何保证矩阵计算的精确度：backward stable & well conditioned

大纲：介绍范数；介绍条件数;介绍条件数的几何意义；介绍有条件数得出的误差估计；介绍backward stable；介绍backward stable + well conditioned 可以得出求解结果精确的结论；介绍如何判断是否backward stable。介绍实际操作中该怎么做：1）判断是否well conditioned，如果不是该怎么调整；2）通过residual判断是否...

2017-01-21 20:29:42 2298 1

原创 “补码”的数学原理

摘要：本文介绍补码进行运算为什么是有效的。通过求同“同余”将减法变成加法；而求补码就是求同余。

2017-01-20 21:23:23 6667 9

原创 “_MSC_VER” 不匹配项

1. 问题我出现的错误是： error LNK2038: 检测到“_MSC_VER”的不匹配项: 值“1600”不匹配值“1900”2. 含义：visual studio 2010（1600）生成的东西和visual studio 2015生成的东西“不匹配”。3. 解决方案1）清理解决方案2）更改平台工具集3）更换附加库

2017-01-14 10:23:15 7676

原创迭代法01--定常迭代法

本文主要介绍定常迭代法，给出迭代公式，分析其收敛性，包括：- Jacobi 迭代- Gauss-Seidel 迭代- Successive Over-Relaxation

2017-01-07 00:01:02 4059

原创 C语言课程常见问题

1. 上机naive问题输出的格式一定要与题目要求一致：不要多空格、不要少空格不要多换行、不要少换行scanf后面的变量要加&：int a;float b;char c;scanf("%d%f%c", &a, &b, &c);窗口一闪而过，在main函数最后加getchar();#include int main(){ /*your code here

2016-12-23 23:10:05 1970 1

原创矩阵分解总结

1. 简介对于线性系统Ax=bAx=b，其解为x=A−1bx=A^{-1}b。但通常情况下，求矩阵的逆十分耗时，应尽量避免。常用的方法包括直接法（Direct Method）和迭代法。　　直接法的核心就是矩阵分解（Matrix Factorization），包括：　SVD分解QR分解LU分解（等价于高斯消元）Cholesky 分解LDLT分解

2016-12-23 21:27:31 2713

原创详解二进制文件和文本文件的区别暨换行和回车详解

介绍了换行和回车符的历史及区别；C语言处理二进制文件和文本文件的唯一区别在于对换行的处理

2016-12-05 21:15:26 10479 1

原创换行和回车（草稿）

本篇主要介绍换行和回车，以及在C语言中的应用。1. 换行和回车的历史2. 在各个操作系统中的表示3. 按回车会得到什么getchar值获得’\n’；getch()只获得’\r’.

2016-12-05 19:49:13 1151

原创基础牛顿法及其变形（basic Newton's method & its variants）

本篇详细介绍用于求解一元函数零点的牛顿法及其各种形式的变形，并给出收敛性质。1. 简介1.1 问题引出问题：求函数f:I⊆R→Rf:I⊆R→Rf : I \subseteq R\to R的零点αα\alpha（f(α)=0f(α)=0f(\alpha)=0）。这个问题可以使用如下形式的单点迭代法（one-point iteration），得到迭代函数g(x)g(x)...

2016-12-04 11:22:09 5166 4

原创布料仿真论文推荐

本篇主要介绍布料仿真领域里的经典论文，例如，large steps in cloth simulation。如果对这个主题特别感兴趣的可以留言，我会加快我的步伐。 to be continued……

2016-11-30 15:47:55 2569 11

原创 Interactive Simulation of Elastic Deformable Materials论文中的公式纠错

简介了布料仿真的两种动力学方程：基于惩罚力（penalty force）的和基于约束力（constraint force）的。纠正了论文Interactive Simulation of Elastic Deformable Materials里的公式错误，并给出了结果对比。

2016-11-30 15:18:39 1002 5

原创布料仿真中常用积分方法

1. 简介布料仿真中，我们通常将布料剖分为三角形网格（或四边形网格），并用弹簧-质点模型构造动力学系统：质点即三角形的顶点，弹簧即三角形的边。质点在外力（如，重力）和内力（弹簧力）的作用下根据牛顿第二定律可得到一组常微分方程，通常我们很难求其解析解，只能求其数值解。常用方法有：显式/隐式欧拉法，Symplectic Euler，Midpoint method，Leapfrog integrati

2016-11-28 22:36:41 4166 3

原创非线性方程（组）的求解

1. 非线性方程的求解1.1 简介　非线性方程是指含有指数和余弦函数等非线性函数的方程，例如，ex−cos(πx)=0ex−cos(πx)=0e^x-cos(\pi x)=0 与线性方程相比，无论是解的存在性，还是求解的计算公式，非线性方程问题逗逼线性问题要复杂的多，对于一般线性方程f(x)=0f(x)=0f(x)=0，既无直接法可用，也无一定章程可寻。常用的方法有：实根的对分...

2016-11-25 15:44:00 52213 4

原创前向后项差分和显式隐式欧拉法

> 本文主要介绍前向后向差分，显式隐式欧拉法及其稳定性分析。> 前向差分对应显式方法，后向差分对应隐式方法。显式欧拉法是比较流行的显式方法，隐式欧拉法是比较流行的隐式方法。> 显式欧拉法条件稳定，对积分步长有要求，隐式欧拉法无条件稳定，对积分步长无要求（理论上如此，但实际使用中对积分步长仍然有要求，只是比显式欧拉宽松）。

2016-11-15 16:58:24 65359 13

原创最优化方法综述

本篇文章对常用的最优化方法进行简要介绍和比较，主要涉及的方法有：梯度下降法牛顿法、拟牛顿法和高斯牛顿法共轭梯度法Levenberg–Marquardt模拟退火法1. 简介最优化问题简介2. 梯度下降法批量梯度下降法和随机梯度下降法3. 牛顿法4. 共轭梯度法5. Levenberg-Marquardt方法6. 模拟退火法（启发式方法）由于最近比较忙，之后会一一展开。

2016-11-12 11:47:45 1636