Vector decomposition with Matlab

Question

Gogu 2019-7-11

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/471268-vector-decomposition-with-matlab

编辑： Bruno Luong 2019-7-12

采纳的回答： Bruno Luong

在 MATLAB Online 中打开

Hello !

I want to split a vector b, in two parts, b1 and b2, like that:

b= [2; 1; 2; 1; 0; 1; 2; 1; 2]

b1=[2; 1; 1; 1; 0; 0; 1; 0; 0]

b2=[0; 0; 1; 0; 0; 1; 1; 1; 2]

But my code, work only for n=4; for n>4, don't work fine.

n=4;
b=zeros((n-1)^2,1);
b(1)=2;
b(2:(n-2))=ones(n-3,1);
b(n-1)=2;
for i=1:n-3
    b((n-1)*i+1)=1;
    b((n-1)*(i+1))=1;
end
b((n-1)*(n-2)+1)=2;
b((n-1)*(n-2)+2:(n-1)*(n-2)+n-2)=ones(n-3,1);
b((n-1)^2)=2;
b1(1)=2;
b1(2:(n-2))=1;
b1(n-1)=1;
for i=1:n-3
    b1((n-1)*i+1)=1;
    b1((n-1)*(i+1))=0;
end
b1((n-1)*(n-2)+1)=1;
b1(n*2:(n-1)*(n-2))=0;
  
b1((n-1)*(n-2)+2:(n-1)*(n-2)+n-2)=zeros(n-3,1);
b1((n-1)^2)=0;
b2(1)=0;
b2(2:(n-2))=0;
b2(n-1)=1;
b2((n-1)*(n-2))=1;
b2((n-1)*(n-2)+1)=1;
b2((n-1)*(n-2)+2:(n-1)*(n-2)+n-2)=ones(n-3,1);
b2((n-1)^2)=2;

2 个评论
显示无隐藏无

Torsten 2019-7-11

One can split a vector in two vectors in many ways. What are the rules ?

Gogu 2019-7-11

b=b1+b2

For the first half of elements, the rules are:

the first 2 keep in b1; 1 keep in b1; 2 split in two values of 1; 1 keep in b1

For the second half of elements, the rules are:

1 move to b2; 2 split in two; the last 2, move to b2

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Bruno Luong 2019-7-11

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/471268-vector-decomposition-with-matlab#answer_382844

编辑：Bruno Luong 2019-7-12

在 MATLAB Online 中打开

b= [2; 1; 2; 1; 0; 1; 2; 1; 2]
b1=min(b,1);
b1([find(b==2,1,'first'),find(b==2,1,'last')])=2;
b2=b-b1;
n=length(b);
tail=floor(n/2)+1:n; %or ceil?
[b1(tail),b2(tail)]=deal(b2(tail),b1(tail));
b1
b2

2 个评论
显示无隐藏无

Gogu 2019-7-11

在 MATLAB Online 中打开

Unfortunately, one rule for the second half of elements, are not met (1 move to b2) :

b1 =
   2   1   1   1   0   1   1   1   0
b2 =
   0   0   1   0   0   0   1   0   2

Gogu 2019-7-12

在 MATLAB Online 中打开

Solved !!!

lx = (length(b));
half = ceil(lx/2);
b1(1)=2;
b2(1)=0;
b1((n-1)^2)=0;
b2((n-1)^2)=2;
j=(n-1)^2;
for i=2:j-1
  if b(i)==2
    b1(i)=1;
    b2(i)=1;
  end
end
for i=2:half
  
  if find(b(i)==1)
    b1(i)=1;
    b2(i)=0;
  end
  
end
for i= half + 1 : j-1
  if find(b(i)==1)
    b1(i)=0;
    b2(i)=1;
  end
end
b1'
b2'