编写矩阵程序

求这个矩阵的编程，可以用C++或者R语言，急求！！

陶雯岩 发表于 2016-4-1 10:10
好的，我马上试一下

上面代码第5行那里应该是"输入不能为空"啊

结果是个对称矩阵，用matlab我可以较快的实现，c++的话要慢一点，r不是很熟。

天外来客2013 发表于 2016-4-1 00:20
结果是个对称矩阵，用matlab我可以较快的实现，c++的话要慢一点，r不是很熟。

matlab也可以，麻烦您帮忙写一下程序可以吗？有什么其他问题可以直接联系我，QQ464588223，谢谢！

陶雯岩 发表于 2016-4-1 08:27
matlab也可以，麻烦您帮忙写一下程序可以吗？有什么其他问题可以直接联系我，QQ464588223，谢谢！

我自己验证了下，貌似可以，你可以用你的数据验证下。

1. %计算矩阵B(mxm)
   
2. function B = calB(P)
   
3. % 输入：P,维数mxn
   
4.     if(nargin < 1)
   
5.         error('输出不能为空');
   
6.     end
   
7.     m = size(P,1);%m为P的行数
   
8.     B = zeros(m,m);
   
9.     %计算上三角阵(包括对角线上的)
   
10.     for i = 1:m
    
11.         for j = i:m
    
12.             B(i,j) = cal_IJ(P(i,:),P(j,:));
    
13.         end
    
14.     end
    
15.     %计算最终的对称矩阵
    
16.     B = tril(B',-1) + B;
    
17. end
    
18. 
    
19. %% 计算矩阵P第i和第j行的值——结果中位置(i,j)处的值               
    
20. function sumIJ  = cal_IJ(row_I,row_J)
    
21.     if(~isequal(length(row_I),length(row_J)))
    
22.         error('输入向量的长度不一致！请重新输入');
    
23.     end
    
24.     L = length(row_I);%向量row_I的长度
    
25.     sumIJ = 0;
    
26.     for i = 1:L
    
27.         temp1 = row_I(i)-row_I;
    
28.         temp2 = row_J(i)-row_J;
    
29.         sumIJ = sumIJ + temp1*temp2';
    
30.     end
    
31. end        

复制代码

天外来客2013 发表于 2016-4-1 10:04
我自己验证了下，貌似可以，你可以用你的数据验证下。

好的，我马上试一下

天外来客2013 发表于 2016-4-1 10:04
我自己验证了下，貌似可以，你可以用你的数据验证下。

这个是我的验证例子: