机器学习的小学生

功能：主函数管理或融合多个子函数的输出，子函数来处理各个事物。% 多个函数相互调用，同时处理相同的数据。% 或者说主函数是一个不做实事的函数，右各个子函数处理每个事物。x = 5;y = 6;% compute x*y + x/y = 30+5/6 = 30.83func = {@add_sum,{@times,@divide}};%　which add_sum;res = feval(f

来自 gpml-matlab-v3.5-2014-12-08 （高斯过程）的startup.m文件：操作流程：读取当前程序的文件名，接着读取当前文件所在的路径。根据相对路径找到其他引用文件添加到Matlab的包含路径中。% startup script to make Octave/Matlab aware of the GPML package%% Copyright (c) by Carl

利用hessian 函数计算hessian矩阵H，然后利用collect获取一次项的系数F： clear all,close all;syms x y ;% 求解该函数的最小值f = (x - 2)^2 + (y - 3)^2 + 5; % s.t. 1 <= x <= 3,2 <= y <= 4%lb = [1 2]';ub = [3 4]'; % lb = [-inf -inf]' ub

转自： http://blog.youkuaiyun.com/bluebelfast/article/details/17612537      现在学术界越来越多人使用Python 语言取代Matlab 作为研究的主要程序编写工具。python语言看似有点类似matlab语言，对于习惯Matlab 的使用者来说，转换到Python语言应该并不困难，但是有些关键要注意的地方，也就是这两种语言的差别。掌握了两种

Python: arr = np.load('ilsvrc_2012_mean.npy')print arr.shape # (3,256,256) channel * height * width########## 数据的转置 ###############arr = np.transpose(arr,[1,2,0]) # height * width *　channelprint arr.sh

% 查找集合中某个元素的位置C = {'xlh','gyl','xyh'};c = strcmp(C,'gyl'); % logical arrayind0 = find(c==1);% Struct结构体某个字段值的位置 S = repmat(struct('name',[],'age',[],'sex',[]),3,1);S(1).name = 'xlh'; S(1).age = '20

Matlab:“`matlab% hdf5_file_2 clearvars;load train_data_2.mat;load mean_file.mat ; % load mean_filetrain_feat = [train_data{1};train_data{2};train_data{3};train_data{4}];test_feat = train_data{5} ;%%%%%

Python Matlab 功能 str.replace(old,new[,count]) dstStr= strrep(origStr,oldSubstr,newSubStr) filepath = 'E:/data/1.jpg'; filepath2 = strrep(filepath,'data','lfw/data');disp(filepath2); 替

Python:import matplotlib.pyplot as pltimport numpy as npplt.ion() ## Note this correctionfig=plt.figure()#plt.axis([0,1000,0,1])i=0x=list()y=list()while i <1000: temp_y=np.random.random();

实验中，我们经常需要对样本进行随机打乱，一般情况下，首先生成一个随机序列，然后通过下标索引的方式得到打乱后的序列。A 表示原样本，B表示打乱的样本，通过 rand_ind，我们可以知道：B ==> A （B到A的映射）。 而现在我们想知道 A ==> B 的映射。clear all;A = [3 4 5 6 1 2 ];disp('A:');disp(A);if ~exist('tst_ran

1. fix 截断fix(num) 删除 number 参数的小数部分并返回以整数表示的结果。 i=1+fix(365*rand(1,1));%产生[1,365]之间的一个随机数2. waitforbuttonpress 等待按钮按下w=waitforbuttonpress;3. tic,toc 计时tic;func(param);tim...

For 2-D plots…The function IMAGE is what you’re looking for. Here’s an example:img = imread('peppers.png'); %# Load a sample imagescatter(rand(1,20)-0.5,rand(1,20)-0.5); %# Plot some rand

% 定义含有三个样本的符号运算% d : nsample * nlabel% lnd : nsample * nlabelnlabel = 5;nsample = 3;feat_len = 4;% lambda ： nlabel nlabelfor i = 1:nlabel for j = 1:nlabel lambda(i,j) = sym(['lambda' n

理解AFLW数据集提供的头部姿态角度。close all;close all;dbpath = '../data/';dbfile = 'aflw.sqlite';model3d = createMeanFace3DModel(fullfile(dbpath,dbfile));% get one face_idmksqlite('open',fullfile(dbpath,d...

clearall;%% titleinfo_=[ '文档的功能: 是实现对图像的标准化' ... ',标准化的方式采用将所有形状对齐的平均形状. \n'... '文档的创建时间为：2016/11/1/18：37 \n' ... '修改文档的时间为：2018/11/01/15:38'];%% contentdatas = load('datas_for_r...

clc;% task :　remove 1,3 ,1,4 ,原则，移除第一行元素中，不满足纵坐标为 7 的元素。% method 1:A = [1,3;1,7;1,4;3,5;1,10;2,9];B = [];for i = 1:size(A,1) if A(i,1) == 1 && A(i,2) ~= 7 continue; end B = [B ;

描述：每一个样本是通过一个结构体的形式进行描述。多个样本通过一个结构体数组进行描述。其中每一个样本都有一个唯一的标识属性index。我们的目的是找到：某个属性满足条件的，所有样本的表示index。例如，在下面的实验中，我们找到年龄大于21的所有人的index,或者所示唯一的id。 clc;clearvars;close all;n=3;id=[12 3 8];name={'x','l','h

clc;clearvars;fid=fopen('data.txt');tline=fgetl(fid);while ischar(tline) %do some work %first method:use strfind(tline,' ');'seperator' %second method:use textscan(tline,'%S'),以cell的形式扫入到字符串 tl

例如在使用Vlfeat库时，plotbox函数与其他库的plotbox冲突，而我们想使用其他库的plotbox，则首先在Matlab命令行窗口使用命令：which plotbox确定当前使用的plotbox的路径，如：E:\matlabworkplace\headpose_with_block\third_part\vlfeat-0.9.20\toolbox\noprefix\plotbox.m在程

%% show Guass Templateclc;clearvars;close all;%G(x,y,\sigma)=frac{1}{2 \pi \sigma^2} e^{-(x^2+y^2)/{2 \sigma^2}}pi=3.1415926;[x,y]=meshgrid(-10:1:10);sigma=1.6;G=1/(2*pi*sigma^2)*exp(-(x.^2+y.^2)/

函数fitgeotrans的用法： tform = fitgeotrans(movingPoints,fixedPoints,transformationType)example tform = fitgeotrans(movingPoints,fixedPoints,’polynomial’,degree) tform = fitgeotrans(movingPoints,fixe

OpenCV is not designed for solving matrix related operation like: multiplication, eigenvalue decomposition or SVD, etc., so your result is not surprise to me. For your problem, there two kinds of choic

如何对特征维度超过20000的数据进行PCA?当传统的PCA算法应用大非常高维的数据集时，会受到限制。例如，在生物信息学中的基因表达数据。但是通常我们仅需要前两个或者三个主蹭饭进行可视化数据。对于仅抽取前k个成分，我们可以使用基于sensible principal components 分析的PPCA（probabilistic PCA ）。例如使用修改的PCA matlab脚本[ppca.m]，

可以通过将imshow函数属性”Border’的值设置为’tight’，而imshow的默认值为’loose’，然后通过figure1-文件-另存为的方式将图片进行保存：imshow(img,'Border','tight');

clearvars;clc;close all;%创建1*10大小的celldata=cell(1,10);for i=1:10 data{i}.name=char(i+96); data{i}.age=i+10; end%% 访问方式1：a=data(2);% a=1×1 cell,%% 访问方式2：% 花括号{}引用cell中的内容，% 其等价于data(

平均形状的计算算法有两种。1. 粗略的获取平均形状经过三步走：原点化–尺度化–平均化(1)原始带有基准点的数据，例如取出一个图像的基准点:landmark=data{i}.landmark;(2)由于对于不同的样本图像，人脸形状坐标是相对于图像的左上角。因此，需要将其进行原点化：% landmark和shape表示是同一个意思landmark=data{i}.landmark;shape_min

预分配内存简介：对于for，while循环，在循环的过程中每次不断的增加数据结构的大小，影响了性能和内存的使用。重复的调整数据的大小需要Matlab花费额外的时间寻找更大的连续内存块，并且将现在的数组移动到连续的内存块中。通常可以通过预分配数组所需的最大内存来来改善代码的执行时间。预分配数组不同类型的数组的预分配函数： 数组类型 函数 例子 数值数组 zeros y=zeros

function y=func(input1,input2)%函数的功能：%函数的描述：%函数的使用：y=func(input1,input2)%输入：% input1:% input2:%输出：% Y:%例子：y=func(1,'type1');%注意事项：利用函数的适用范围。%文档日期：%标签：%创建日期：%最后更新日期：

1. 对结构体数组(struct array)进行排序:即对s_info进行排序n=3;%% 初始化一个cell数组c_info=cell(n,1);id=[12 3 8];name={'x','l','h'};age=[20 22 24];for i=1:n; c_info{i}.id=id(i); c_info{i}.name=name{i}; c_info{i

将cell数组转换为struct数组n=3;%% 初始化一个cell数组c_info=cell(n,1);id=[12 3 8];name={'x','l','h'};age=[20 22 24];for i=1:n; c_info{i}.id=id(i); c_info{i}.name=name{i}; c_info{i}.age=age(i);end%% 转

clearall;%% titleinfo_=[ '文档的功能: 是实现对图像的标准化' ... ',标准化的方式采用将所有形状对齐的平均形状. \n'... '文档的创建时间为：2016/11/1/18：37 \n'];%% contentdata = load('pose_train_data.mat','s_traindata','s_testdata');t

其类似于C/C++中函数名或者函数的取地址操作，在C/C++中函数的作为参数的传递，可以通过函数名func或者&func的形式作为参数。data=load('../train_task_data.mat','train_task_data');train_task_data=data.train_task_data;opt.lossType='square'; % hinge , squareif

% 28 * 28 for conv3% 14 * 14 for conv4% 7 * 7 for conv5close all;datas = load('datas5224_pure3.mat');datas = datas.datas;f = @(x,meanx,nSigma) exp(-1/2*(x-meanx)'*nSigma*(x-meanx));for i = 1:...