Justpayne的专栏

转载 Matlab去掉右面、上面的小短线（ticks)

先用plot画个图：x = -pi:pi/10:pi;y = tan(sin(x)) - sin(tan(x));plot(x,y,'rs-','LineWidth',1,'MarkerEdgeColor','k',...'MarkerFaceColor','g','MarkerSize',10);看到右边和上边的小短线(ticks)了嘛，有人就是不想要这些短线

转载 奇异值分解

版权声明：    本文由LeftNotEasy发布于http://leftnoteasy.cnblogs.com, 本文可以被全部的转载或者部分使用，但请注明出处，如果有问题，请联系wheeleast@gmail.com前言：    上一次写了关于PCA与LDA的文章，PCA的实现一般有两种，一种是用特征值分解去实现的，一种是用奇异值分解去实现的。在上篇文章中便是基于特征值

转载 一些英文词的标准缩写

有些词可能共用一些缩写。带星号的缩写或词来源于PeopleSoft标准。The following standard word abbreviations should be used in naming records, fields, and SQRs:Word(s)AbbreviationDescriptionAbbrevi

转载 Event Location ode求解器的事件属性

ode求解器的事件属性检测事件matlab微分方程如何设定变量的范围 如dy1=y2dy2=y1+1其中y1的范围为0【解】M文件：function [value,isterminal,direction] = events1(t,y)value = y(1)-4;isterminal= 1;direction = 0;

转载 矩阵求导公式

基本公式：Y = A * X --> DY/DX = A'Y = X * A --> DY/DX = AY = A' * X * B --> DY/DX = A * B'Y = A' * X' * B --> DY/DX = B * A'1. 矩阵Y对标量x求导：相当于每个元素求导数后转置一下，注意M×N矩阵求导后变成N×M了Y = [y(ij)] -->

转载 matlab中的微分方程

第1节  Matlab能够处理什么样的微分方程？ Matlab提供了解决包括解微分方程在内的各种类型问题的函数： 1． 常规微分方程（ODEs）的初始值问题 初值问题是用MATLAB ODE求解器解决的最普遍的问题。初始值问题最典型的是对非刚性度（？nonstiff）问题应用ODE45，对刚性度（？stiff）问题采用ODE15S。（对于stiffness的解释，请参照“什么是Stif

转载 MATLAB 积分算法

在SIMULINK 的仿真过程中选择合适的算法是很重要的，仿真算法是求常微分方程、传递函数、状态方程解的数值计算方法，这些方法主要有欧拉法( Euler) 、阿达姆斯法(Adams) 、龙格·库塔法(Rung-Kutta) ，这些算法都主要建立在泰勒级数的基础上。欧拉法是最早出现的一种数值计算方法，它是数值计算的基础，它用矩形面积来近似积分计算，欧拉法比较简单，但精度不高，现在已经较少使用。阿达姆

转载 Matlab有关不等式画图问题

如何用Matlab画出27.6xy因为没作特殊说明，我们默认楼主是在xoy平面内讨论（排除柱面情形）问题。对于不等式绘图，mathematica有相应的函数——regionplot（2维或者3维），而matlab没有，但是我们有办法。这样的问题（平面不等式问题）主要有两种方法：1. 绘制曲线，给可行域画上hatch；2. 定义空间曲面，对于不满足要求的点设置

转载 分享一些Matlab调试技巧

Matlab的调试总体分为，直接调试和间接调试。1.直接调试：（1）去掉句末的分号；（2）单独调试一个函数：将第一行的函数声明注释掉，并定义输入量，以脚本方式执行 M 文件；（3）适当地方添加输出变量值的语句；（4）添加keyboard命令；2.工具调试：1.）以命令行为主的调试：（1）设置断点：·        dbstop in mfile:在文件名

转载 拉格朗日对偶PPT

转载 拉格朗日对偶（Lagrange duality）

先抛开上面的二次规划问题，先来看看存在等式约束的极值问题求法，比如下面的最优化问题：            目标函数是f(w)，下面是等式约束。通常解法是引入拉格朗日算子，这里使用来表示算子，得到拉格朗日公式为            L是等式约束的个数。    然后分别对w和求偏导，使得偏导数等于0，然后解出w和。至于为什么引入拉格朗日算子可以求出极值，原因是f(w)的dw变化

转载 matlab ode45 函数传自定义参数用法及定步长ode5解算函数

要用的时候总是忘记，这回给把它写在这里！ %%程序1arg1 = 2;arg2 = 1;[T,Y] = ode45('vdp1000',[0 10],[2 0], [], arg1, arg2);plot(T,Y(:,1),'-o'); %%程序2function dy = vdp1000(t, y, flag, arg1, arg2)

转载 n维线性空间上的几何：直线与平面的方程

高维几何是泛函研究的主要对象之一，可对何谓高维几何，笔者没见到明显的定义。笔者认为，所谓高维几何，就是把高维向量空间，同时视为点的集合，研究她的子集的性质和子集之间的关系的学科，就是高维几何(是否如此，请给予点评)。　　那么点与向量有什么关系，他们的表示方法有何区别，在以往的泛函分析中没有明确的概念。一开始用的概念是向量，但当提到距离时，又变成了点，令人摸不到头绪。　　当然笔者要研究的几何

转载 超平面

在数学中，超平面是 n 维欧氏空间中余维度等于一的线性子空间。这是平面中的直线、空间中的平面之推广。设 F 为域（为初等起见，可考虑F=R）。n 维空间 Fn 中的超平面是由方程定义的子集，其中 是不全为零的常数。在线性代数的脉络下，F-矢量空间 V 中的超平面是指形如的子空间，其中 是任一非零的线性映射。在射影几何中

转载 制作PPT进度条

为了美观PPT，然后在前人基础上改了个进度条，用宏命令完成，用OFFICE2010记得用启动宏的pptm模式保存PPT文件，不然宏代码不会保存。按ALT+F8 创建一个宏 填入下面代码 运行下就得到了进度条。注：添加或删减页面，需要手动运行下宏，执行进度条更新效果图：代码：Sub ProgressBar()' by dukenuke@newsmth.ne

转载 Matlab 的 private 文件夹

matlab用多了，发现在代码开发环境方面还真不如VS做的好，而且在项目管理上做的不如其他语言强大，开发小型代码很不错，如果代码工程量大的话，就会发现很难组织好。一个文件里面，第一个函数是主函数，后面的函数对外都不可见，这个就很别扭，所以要把其他需要供其他文件调用的函数都写成一个个文件，文件少了还好，如果多了，就想分分类，怎么分？要知道一个函数或脚本是无法或不方便调用其他文件夹里的函数的，加入ma

转载 Matlab ode45 解微分方程

在用odesolver(ode45, ode15s, …)来解微分方程的时候，最基本的用法是： [t, y] = odesolver(odefun, tspan, y0); 这里的odefun是待求的微分方程。那么odefun中一般会含有多个系统参数，通常要通过改变参数来观察系统动态的变化。那么如何在调用odesolver的时候传递参数呢？ 以前，我都是用全局变量的写法，将参数在主函数和子

转载 拉格朗日乘子法和KKT条件

拉格朗日乘子法(Lagrange Multiplier)和KKT(Karush-Kuhn-Tucker)条件是求解约束优化问题的重要方法，在有等式约束时使用拉格朗日乘子法，在有不等约束时使用KKT条件。前提是：只有当目标函数为凸函数时，使用这两种方法才保证求得的是最优解。对于无约束最优化问题，有很多经典的求解方法，参见无约束最优化方法。拉格朗日乘子法转换为

转载 关于matlab使用函数作为参数的方法

数值计算程序的通用性——函数名作为参数传递 在进行数值算法的程序设计时，往往需要实现程序的通用性，即同一个程序可以对不同的问题进行求解，如：计算函数f(x)的数值积分的程序，我们需要写一个通用的数值积分程序——Integr()，使其可以对不同的f(x)函数均可以进行积分计算，而不需要每计算一个函数的积分就重新写一个Integr()程序，仅需要另写一个f(x)的表达式程序即可。       若使