随笔分类 - Mathematical Modeling
摘要：for mm=1:num_m %电影编号是mm的训练集行号 ff= find(train_vec(:,2)==mm); %train_vec(ff,1) 行号对应的用户编号 count(train_vec(ff,1),mm) = train_vec(ff,3);end 桉的是电影数目，每次...
阅读全文
posted @ 2015-10-22 09:39
加拿大小哥哥
摘要：for iter = 1:num_iters %梯度下降 用户向量 for i = 1:m %返回有0有1 是逻辑值 ratedIndex1 = R_training(i,:)~=0 ; %U(i,:) * V' 第i个用户分别对每个电影的评...
阅读全文
posted @ 2015-10-21 22:05
加拿大小哥哥
摘要：对p= 2，这称为弗罗贝尼乌斯范数（Frobenius norm）或希尔伯特-施密特范数（ Hilbert–Schmidt norm），不过后面这个术语通常只用于希尔伯特空间。这个范数可用不同的方式定义： 这里A*表示A的共轭转置，σi是A的奇异值，并使用了迹函数。弗罗贝尼乌斯范数与Kn上欧几...
阅读全文
posted @ 2015-10-20 21:06
加拿大小哥哥
摘要：一、差分与微分 我自己的理解。二、求解 2.1 矩阵 这就是matlab的计算结果.太小的话放大些:c = 4 5 9 7 2 1 5 2 6>> [x,y]=gradient(c)x = 1.0000 2.5...
阅读全文
posted @ 2015-10-15 17:04
加拿大小哥哥
摘要：1.颜色字符串有'c', 'm', 'y', 'r', 'g', 'b', 'w',和'k'。分别表示青，红紫，黄，红，绿，白和黑。 2.线型字符串有：'-' 为实线, '--' 为虚线, ':' 为点线, '-.' 为点虚线, 及'none' 表示不用线型。 3.标记形式有'+', 'o'...
阅读全文
posted @ 2015-10-12 15:04
加拿大小哥哥
摘要：上面是效果图，看着很不错吧，主要的问题在于用XTickLabel设置横坐标时候，横坐标会扩展，就是说如果label是[1 2 3]，咱就做了三组试验，参数分别是 1 2 3，但是显示是1 2 3 1 2 3 1 2.。。。扩展了，不是需要的，如何解决呢？ 同时设置XTickLabel和XTick，只
阅读全文
posted @ 2015-10-12 15:01
加拿大小哥哥
摘要：在给李X写SVD代码的时候注意到的。>> a = magic(3)a = 8 1 6 3 5 7 4 9 2>> diag(a)ans = 8 5 2>> a = [8 5 2]a = 8 ...
阅读全文
posted @ 2015-10-10 11:24
加拿大小哥哥
摘要：一、问题来源 t = 2xy/(x+y);融合相似度和信任度，我需要验证值域是不是[0,1]；二、求解clear all;clc;%linspace(0:1,0.1)这样是错的，第三个参数是段数，不是间距，应该是11，分成11分x = 0:0.1:1;y = x;[x,y] = meshgrid(...
阅读全文
posted @ 2015-10-10 11:20
加拿大小哥哥
摘要：顶点号 顶点号 权值1 2 4001 3 4502 4 3002 8 2302 9 1403 4 6004 ...
阅读全文
posted @ 2015-09-24 15:10
加拿大小哥哥
摘要：一、需求来源 对空间结构聚类，恰好是圆台，找到了上下底面的方程，所以画图。二、需求解决2.1 绘制平面x = linspace(0,5,100);y = linspace(0,4,100);z = 2.*repmat(x,100,1) + 8.*repmat(y,100,1);surf(x,y,z...
阅读全文
posted @ 2015-09-21 00:33
加拿大小哥哥
摘要：1.不精确set(gca,'XMinorTick','on') 这样的话知识x轴显示了细的密度，网格线并没有变。2.精确 set(gca,'xtick',-1:0.1:1); set(gca,'ytick',-1:0.1:1);
阅读全文
posted @ 2015-09-20 18:06
加拿大小哥哥
摘要：clear all;clc;a = []; %不是null，也不能什么都不是for i=1:10 a = [a i]; end
阅读全文
posted @ 2015-09-10 20:13
加拿大小哥哥
摘要：>> a = [1 2 3 4 8 9];>> b = [4 5 6 1]b = 4 5 6 1>> c = intersect(a,b)c = 1 4 判断相同元素个数用intersect函数。 判断不同元素个数用容斥原理。>> d = leng...
阅读全文
posted @ 2015-09-10 19:56
加拿大小哥哥
摘要：我先想到的是for循环........ 然后sum(find(mod(a:b,c)==0))，从10到100得到874，为什么不对呢？ 比如a = [1 2 3 4 2 3 4 2]，find(a==0)得到的是[2 5 8]，并不是[0 1 0 0 1 0 0 1]，直接判断的话得到的是后...
阅读全文
posted @ 2015-07-28 18:06
加拿大小哥哥
摘要：就是转置的意思，和'一个意思，但是并不重复，因为在cellfun中你无法'这样吧，所以有了这个函数，’只是符号。K>> aa = magic(4)aa = 16 2 3 13 5 11 10 8 9 7 6 1...
阅读全文
posted @ 2015-07-20 09:48
加拿大小哥哥
摘要：binIndices是1*50的cell，每个cell是的1*n（n不定），那么cell2mat(binIndices)得到的是1*sum(cellfun(@length,binIndices))的行向量，比如{1 2 3} {4 5 6}2个cell，得到的是[1 2 3 4 5 6]的行向量...
阅读全文
posted @ 2015-07-18 11:33
加拿大小哥哥
摘要：一.问题来源 看到了min_score_pos = find(A0_scores==min(A0_scores), 1); [r,c] = find(X,k)，返回X中第k个非零元素的行列位置。二.试验K>> a = [1 4 5 69 8 3 23]a = 1 4 5 ...
阅读全文
posted @ 2015-07-17 11:10
加拿大小哥哥
摘要：K>> a = 10*ones(1,10);K>> b = [1 56 23 5 6 45 9 7 89 10];K>> c = b> a(c) = b(c)a = 1 10 10 5 6 10 9 7 10 10 从试验可以看...
阅读全文
posted @ 2015-07-14 18:39
加拿大小哥哥
摘要：怎样用mat2cell将一个100*100的矩阵分成10个10*100的矩阵？ 根据帮助中 c = mat2cell(x,m,n)应该这样写 mat2cell(x,[10 10 10 10 10 10 10 10 10 10]，[100]) 能不能有简便的办法？ 假...
阅读全文
posted @ 2015-07-14 17:48
加拿大小哥哥
摘要：Compute the mean of each vector in cell array C.C = {1:10, [2; 4; 6], []};averages = cellfun(@mean, C)This code returnsaverages = 5.5000 4.0000 ...
阅读全文
posted @ 2015-07-14 17:27
加拿大小哥哥
浙公网安备 33010602011771号