Re: [问题] 矩阵元素的重新分配疑问 ?

楼主: celestialgod (天)   2016-06-17 16:00:32
※ 引述《MrsJ (MrsJ)》之铭言:
: 各位先进好 自修matlab中 目前遇到一个小瓶颈
: 没法用(我已知的)语法建构我所需(较庞大的)的矩阵
: 问题简化描述如下:
: 令 7x2 矩阵 MatA=
: [ 2.0000 3.4641
: 4.0000 3.7417
: 5.0000 3.8730
: 6.0000 4.0000
: 8.0000 4.2426
: 10.0000 4.4721
: 13.0000 5.0898 ]
: 这里想构造另一个更大的 15x2 矩阵 MatB, 比如说
: 令 MatB= ( 方便起见 MatB(1,1) = 1.000 )
: [ 1.0000 0
: 2.0000 3.4641
: 3.0000 0
: 4.0000 3.7417
: 5.0000 3.8730
: 6.0000 4.0000
: 7.0000 0
: 8.0000 4.2426
: 9.0000 0
: 10.0000 4.4721
: 11.0000 0
: 12.0000 0
: 13.0000 5.0898
: 14.0000 0
: 15.0000 0 ]
: 亦即规则是, 若 MatB(:,1)的某元素 MatB(N,1) 在MatA(:,1)中是不存在的
: 那么就令 MatB(N,2)=0
: 因为实际需要处理的矩阵dimension更大 不可能蠢到手动插值
: 目前也找不到合适的 fuction 可以直接把MatA的元素分解排进 MatB 中
: 想尝试使用再加一个循环来执行这个要求 ... 我不幸卡住。
A = [randi(25, 10, 1), randn(10, 1)];
B = [(1:20)', zeros(20, 1)];
[existence, idx] = ismember(A(:, 1), B(:, 1));
B(idx(existence), 2) = A(existence, 2);
disp(B)
1.0000 0
2.0000 -0.1303
3.0000 -1.4491
4.0000 0
5.0000 0.1837
6.0000 0
7.0000 0
8.0000 -0.1765
9.0000 -1.3320
10.0000 0
11.0000 0
12.0000 0.4517
13.0000 0
14.0000 0
15.0000 0
16.0000 0
17.0000 0
18.0000 0
19.0000 0.3335
20.0000 0.3914
第一个column不是整数应该也没问题
只是可能数值问题(rounding error),7电脑会把他视作不同
existence是确保B的第一行未包含A的第一行的情况
作者: MrsJ (MrsJ)   2016-06-17 16:12:00
谢谢 我试试..
作者: name0625 (lawpy)   2016-06-17 16:20:00
太强啦!

Links booklink

Contact Us: admin [ a t ] ucptt.com