风翼冰舟的博客

官方点的叫法，在网上被称为函数句柄用途在于，比如你定义了两个函数：strupper()用于寻找字符串的大写字母数目；strlower()用于寻找字符串的小写字母数目现在使用str在这里面调换使用如下:func='str';%统一前缀myfunc=str2func(func 'lower');%需要的函数后缀部分，该谁就调用谁num=myfunc(s)%s为你的字符串，直接就调用了s

先来看自带函数zscore的使用>> A=[1 2 3;4 5 6];>> [B,A_mean,A_std]=zscore(A)B = -0.7071 -0.7071 -0.7071 0.7071 0.7071 0.7071A_mean = 2.5000 3.5000 4.5000A_std = 2.1213

matlab GPU编程基础

前言推荐两个比较好的教程:BLAS (Basic Linear Algebra Subprograms) LAPACK for Windows命名规范BLAS基本线性代数子程序的函数命令都有一定规范，便于记忆<character> <name> <mod> ()character定义的是数据类型 s 实数域，单精度 c 复数域，单精度 d 实数域，双精度 z 复数域，双

其实也就是用到了Ruslan Salakhutdinov and Geoff Hinton提供的工具包% Version 1.000%% Code provided by Ruslan Salakhutdinov and Geoff Hinton%% Permission is granted for anyone to copy, use, modify, or distribute

分段线性插值，临近插值，球面插值，三次多项式插值！！

前言工程中经常将参数文件存储为bin格式, 但是实际中为了分析其参数, 也不好用C++去读取调试它, 所以可以用matlab或者python去读取它, 但是还是蛮坑的Matlab中的读取和写入写入文件比较坑的是, 一定要注意自己的文件存储的类型, 比如数值是float还是double之类的, 不然很容易在读取的时候出问题.关于数字的写入方法如下:fileID=fopen('MatrixTest.bi

前言最近motion graph相关实验，发现实现运动过渡需要构建运动图，而为了避免运动过渡陷入死胡同，需要对图结构进行裁剪，方法就是计算图模型的极大强联通分量，但是自己懒得去实现，所以就去搜了一下matlab中有无相关实现，果然找到了。国际惯例，参考网址：Find strongly or weakly connected components in graph创建稀疏矩阵...

前言最近复现音乐驱动舞蹈的文章《Dancing-to-Music Character Animation》，用到了与傅里叶变换很相似的称为常Q变换的方法去分割音乐，所以对傅里叶变换做了一个小了解，本文不深入各种乱糟糟的理论，比如什么蝶形算法啥的，单纯讨论离散傅里叶变换(DFT)，我们常见的是快速傅里叶变换(FFT)，其实FFT是对DFT的一个计算优化，主要是剔除DFT里面有周期性之类的被冗余计...

前言上一篇博客讲了离散傅里叶变换，里面的实例是对整个信号进行计算，虽然理论上有N点傅里叶变换(本博客就不区分FFT和DFT了，因为它俩就是一个东东，只不过复杂度不同)，但是我个人理解是这个N点是信号前面连续的N个数值，即N点FFT意思就是截取前面N个信号进行FFT，这样就要求我们的前N个采样点必须包含当前信号的一个周期，不然提取的余弦波参数与正确的叠加波的参数相差很大。如果在N点FFT的时...

简单来说，reshape(A,m,n)就是用来把原矩阵的元素，按照列读取拍成一行，然后按照指定的m*n矩阵按行放好。比如原矩阵为

简单的卷积就不说了，向量卷积用此函数与用conv效果相同，矩阵卷积用此函数与conv2的二维卷积效果相同。此函数的方便之处在于支持三维卷积：其实相对于conv2来说就是省了一个for循环。对于三维卷积，比如A矩阵大小为[2,3,3]，B矩阵大小为[2,3]，计算A与B的卷积结果就是过程就是对于A矩阵的第三维的每一维分别卷积，等价于利用conv2利用卷积核B对A每一维卷积三次。看代码：>>

图像基本的矩阵转换，存储和读取

matlab之norm计算范数

网上的说法：要计算一个系统相频特性，就要用到反正切函数，计算机中反正切函数规定，在一、二象限中的角度为0～pi，三四象限的角度为0～-pi。若一个角度从0变到2pi，但实际得到的结果是0～pi，再由-pi～0，在w=pi处发生跳变，跳变幅度为2pi，这就叫相位的卷绕。unwrap（w）就是解卷绕，使相位在pi处不发生跳变，从而反应出真实的相位变化 "查看 unwrap 的帮助文

【注】给定时间精确到小时的情况下的转换函数：function [ date ] = ConvertDate( x )%将unix时间戳转换为标准时间% date = datestr(1426406400/86400 + datenum(1970,1,1));date = datestr((x-3600*24+8*3600)/86400 + 70*365+19+datanum(190...

A=[1 2 3;4 5 6]A = 1 2 3 4 5 6>> [B,PS]=mapminmax(A,0,1)B = 0 0.5000 1.0000 0 0.5000 1.0000PS = name: 'mapminmax'

简单的调用方法：bsxfun(@plus，A，B)，其中@plus代表的是加法，还可以换成减法@minus，乘法@times等，可以参考matlab里面的帮助文档，对着bsxfun按F1即可。有何用处？   我目前用来就是实现对矩阵A和B的加减乘除操作，它的优势在于，当A和B的维度不相等，并且A和B各自有一个维度必须为1时，简单的说就是使用bsxfun的前提条件是A和B的维度分别为1*m和

原文地址：http://www.ilovematlab.cn/thread-293710-1-1.html-------------------------------------conv2函数----------------------------------------1、用法C=conv2(A,B,shape);        %卷积滤波复制代码A

diff函数是用来求导数的目前用到的调用格式为diff(A，m，n)意思是对矩阵A求m阶导数，n=1时按列求，n=2时按行求举例：>> A=[1 2 3;4 5 6;7 8 9; 10 11 12]A = 1 2 3 4 5 6 7 8 9 10 11 12>> B=dif

主要借鉴此博客代码：http://blog.youkuaiyun.com/sherry_gp/article/details/50560003但是这个博主的代码达不到我想要的效果，所以修改了一下我想要实现的效果是：给定一列的预测标签，以及这一列标签的哪一部分理应属于哪一部分标签。此代码实现的功能是：给定条件：给定预测标签为A=[1 2 1 1 2 2 3 2 3 3 3 3 4 4 1

前言刚体运动旋转一般用：欧拉角、四元数、轴角对等表示，在对某个坐标旋转的时候，只需将欧拉角或四元数转换为旋转矩阵，并与原始坐标相乘，便可得到旋转以后的坐标。这里主要看看欧拉角、四元数和旋转矩阵。国际惯例，参考博客：四元数与旋转【Unity技巧】四元数（Quaternion）和旋转三维转动的欧拉角和转轴转角参数相互转换的谢国芳公式 欧拉角百度百科欧拉角维基百科欧拉角百...