我编了下面一个程序,是二叉树结构划分的,为什么每次调用就出现下面的错误 One or more output arguments not assigned during call to 'parsubsample'?请求高手指点!谢谢
function [Y1,f,Y2,g,Y3,h,Y4,l,Y5,o,Y6,p,Y7,q,Y8,r,P,c]=parsubsample(N,z1,z2,z3,m)
%the order of the splitting variable both in the left and right branch is z1,z2,z3
%to find the corresponding partion of the tree and the corresponding
%subsamples when svar(1,8)=z1,svar(2,4)=svar(2,,12)=z2,svar(3,2)=svar(3,6)=svar(3,10)=svar(3,14)=z3
m=length(z1);
%m is the number of elements in z1, z2 and z3
c(1)=0;
%to find the corresponding partion of the tree and the corresponding
%subsamples when svar(1,8)=z1,svar(2,4)=svar(2,,12)=z2,svar(3,2)=svar(3,6)=svar(3,10)=svar(3,14)=z3
if N(1,8)=='terminal'
%if the root node is the terminal node
P=sym('[ ]');
f=0;g=0;h=0;l=0;o=0;p=0;q=0;r=0;
for t=2:m
f=f+1;
Y1(f,:)=[z1(t),z2(t),z3(t)];
c(t)=1;
end
Y2=zeros(0,3);Y3=zeros(0,3);Y4=zeros(0,3);Y5=zeros(0,3);
Y6=zeros(0,3);Y7=zeros(0,3);Y8=zeros(0,3);
end
%there is no terminal node in the first row of the tree
if N(2,4)=='terminal'&N(2,12)=='terminal'
%if the left and right nodes in the second row are both terminal nodes
P1=sym('[z1<=0]');
P2=sym('[z1>0]');
P=[P1,P2];
f=0;g=0;h=0;l=0;o=0;p=0;q=0;r=0;
for t=2:m
if z1(t-1)<=0
f=f+1;
Y1(f,:)=[z1(t),z2(t),z3(t)];
c(t)=1;
else
g=g+1;
Y2(g,:)=[z1(t),z2(t),z3(t)];
c(t)=2;
end
end
Y3=zeros(0,3);Y4=zeros(0,3);Y5=zeros(0,3);
Y6=zeros(0,3);Y7=zeros(0,3);Y8=zeros(0,3);
end
if N(2,4)=='terminal'&N(3,10)=='terminal'&N(3,14)=='terminal'
%if the left in the second row and the third and fourth node in the third row are all terminal nodes
P1=sym('[z1<=0]');
P2=sym('[z1>0z2<=0]');
P3=sym('[z1>0z2>0]');
P=[P1,P2,P3];
f=0;g=0;h=0;l=0;o=0;p=0;q=0;r=0;
for t=2:m
if z1(t-1)<=0
f=f+1;
Y1(f,:)=[z1(t),z2(t),z3(t)];
c(t)=1;
elseif z1(t-1)>0&z2(t-1)<=0
g=g+1;
Y2(g,:)=[z1(t),z2(t),z3(t)];
c(t)=2;
elseif z1(t-1)>0&z2(t-1)>0
h=h+1;
Y3(h,:)=[z1(t),z2(t),z3(t)];
c(t)=3;
end
end
Y4=zeros(0,3);Y5=zeros(0,3);Y6=zeros(0,3);
Y7=zeros(0,3);Y8=zeros(0,3);
end
if N(2,4)=='terminal'&N(4,9)=='terminal'&N(3,14)=='terminal'
%if the left in the second row ,the fourth node in the third row and the fifth and sixth node in the fourth row are all terminal nodes
P1=sym('[z1<=0]');
P2=sym('[z1>0z2<=0z3<=0]');
P3=sym('[z1>0z2<=0z3>0]');
P4=sym('[z1>0z2>0]');
P=[P1,P2,P3,P4];
f=0;g=0;h=0;l=0;o=0;p=0;q=0;r=0;
for t=2:m
if z1(t-1)<=0
f=f+1;
Y1(f,:)=[z1(t),z2(t),z3(t)];
c(t)=1;
elseif z1(t-1)>0&z2(t-1)<=0&z3(t-1)<=0
g=g+1;
Y2(g,:)=[z1(t),z2(t),z3(t)];
c(t)=2;
elseif z1(t-1)>0&z2(t-1)<=0&z3(t-1)>0
h=h+1;
Y3(h,:)=[z1(t),z2(t),z3(t)];
c(t)=3;
elseif z1(t-1)>0&z2(t-1)>0
l=l+1;
Y4(l,:)=[z1(t),z2(t),z3(t)];
c(t)=4;
end
end
Y5=zeros(0,3);Y6=zeros(0,3);Y7=zeros(0,3);Y8=zeros(0,3);
end
if N(2,4)=='terminal'&N(3,10)=='terminal'&N(4,13)=='terminal'
%if the left in the second row ,the third node in the third row and the seventh and eighth node in the fourth row are all terminal nodes
P1=sym('[z1<=0]');
P2=sym('[z1>0z2<=0]');
P3=sym('[z1>0z2>0z3<=0]');
P4=sym('[z1>0z2>0z3>0]');
P=[P1,P2,P3,P4];
f=0;g=0;h=0;l=0;o=0;p=0;q=0;r=0;
for t=2:m
if z1(t-1)<=0
f=f+1;
Y1(f,:)=[z1(t),z2(t),z3(t)];
c(t)=1;
elseif z1(t-1)>0&z2(t-1)<=0
g=g+1;
Y2(g,:)=[z1(t),z2(t),z3(t)];
c(t)=2;
elseif z1(t-1)>0&z2(t-1)>0&z3(t-1)<=0
h=h+1;
Y3(h,:)=[z1(t),z2(t),z3(t)];
c(t)=3;
elseif z1(t-1)>0&z2(t-1)>0&z3(t-1)>0
l=l+1;
Y4(l,:)=[z1(t),z2(t),z3(t)];
c(t)=4;
end
end
Y5=zeros(0,3);Y6=zeros(0,3);Y7=zeros(0,3);Y8=zeros(0,3);
end
if N(2,4)=='terminal'&N(4,9)=='terminal'&N(4,13)=='terminal'
%if the left in the second row ,the fifth,sixth,seventh and eighth node in the fourth row are all terminal nodes
P1=sym('[z1<=0]');
P2=sym('[z1>0z2<=0z3<=0]');
P3=sym('[z1>0z2<=0z3>0]');
P4=sym('[z1>0z2>0z3<=0]');
P5=sym('[z1>0z2>0z3>0]');
P=[P1,P2,P3,P4,P5];
f=0;g=0;h=0;l=0;o=0;p=0;q=0;r=0;
for t=2:m
if z1(t-1)<=0
f=f+1;
Y1(f,:)=[z1(t),z2(t),z3(t)];
c(t)=1;
elseif z1(t-1)>0&z2(t-1)<=0&z3(t-1)<=0
g=g+1;
Y2(g,:)=[z1(t),z2(t),z3(t)];
c(t)=2;
elseif z1(t-1)>0&z2(t-1)<=0&z3(t-1)>0
h=h+1;
Y3(h,:)=[z1(t),z2(t),z3(t)];
c(t)=3;
elseif z1(t-1)>0&z2(t-1)>0&z3(t-1)<=0
l=l+1;
Y4(l,:)=[z1(t),z2(t),z3(t)];
c(t)=4;
elseif z1(t-1)>0&z2(t-1)>0&z3(t-1)>0
o=o+1;
Y5(o,:)=[z1(t),z2(t),z3(t)];
c(t)=5;
end
end
Y6=zeros(0,3);Y7=zeros(0,3);Y8=zeros(0,3);
end
if N(2,12)=='terminal'&N(3,2)=='terminal'&N(3,6)=='terminal'
%if the right in the second row and the first and second node in the third row are all terminal nodes
P1=sym('[z1<=0z2<=0]');
P2=sym('[z1<=0z2>0]');
P3=sym('[z1>0]');
P=[P1,P2,P3];
f=0;g=0;h=0;l=0;o=0;p=0;q=0;r=0;
for t=2:m
if z1(t-1)<=0&z2(t-1)<=0
f=f+1;
Y1(f,:)=[z1(t),z2(t),z3(t)];
c(t)=1;
elseif z1(t-1)<=0&z2(t-1)>0
g=g+1;
Y2(g,:)=[z1(t),z2(t),z3(t)];
c(t)=2;
elseif z1(t-1)>0
h=h+1;
Y3(h,:)=[z1(t),z2(t),z3(t)];
c(t)=3;
end
end
Y4=zeros(0,3);Y5=zeros(0,3);Y6=zeros(0,3);
Y7=zeros(0,3);Y8=zeros(0,3);
end
if N(2,12)=='terminal'&N(4,1)=='terminal'&N(3,6)=='terminal'
%if the right in the second row ,the second node in the third row and the first and second node in the fourth row are all terminal nodes
P1=sym('[z1<=0z2<=0z3<=0]');
P2=sym('[z1<=0z2<=0z3>0]');
P3=sym('[z1<=0z2>0]');
P4=sym('[z1>0]');
P=[P1,P2,P3,P4];
f=0;g=0;h=0;l=0;o=0;p=0;q=0;r=0;
for t=2:m
if z1(t-1)<=0&z2(t-1)<=0&z3(t-1)<=0
f=f+1;
Y1(f,:)=[z1(t),z2(t),z3(t)];
c(t)=1;
elseif z1(t-1)<=0&z2(t-1)<=0&z3(t-1)>0
g=g+1;
Y2(g,:)=[z1(t),z2(t),z3(t)];
c(t)=2;
elseif z1(t-1)<=0&z2(t-1)>0
h=h+1;
Y3(h,:)=[z1(t),z2(t),z3(t)];
c(t)=3;
elseif z1(t-1)>0
l=l+1;
Y4(l,:)=[z1(t),z2(t),z3(t)];
c(t)=4;
end
end
Y5=zeros(0,3);Y6=zeros(0,3);Y7=zeros(0,3);Y8=zeros(0,3);
end
if N(2,12)=='terminal'&N(4,5)=='terminal'&N(3,2)=='terminal'
%if the right in the second row ,the first node in the third row and the third and fourth node in the fourth row are all terminal nodes
P1=sym('[z1<=0z2<=0]');
P2=sym('[z1<=0z2>0z3<=0]');
P3=sym('[z1<=0z2>0z3>0]');
P4=sym('[z1>0]');
P=[P1,P2,P3,P4];
f=0;g=0;h=0;l=0;o=0;p=0;q=0;r=0;
for t=2:m
if z1(t-1)<=0&z2(t-1)<=0
f=f+1;
Y1(f,:)=[z1(t),z2(t),z3(t)];
c(t)=1;
elseif z1(t-1)<=0&z2(t-1)>0&z3(t-1)<=0
g=g+1;
Y2(g,:)=[z1(t),z2(t),z3(t)];
c(t)=2;
elseif z1(t-1)<=0&z2(t-1)>0&z3(t-1)>0
h=h+1;
Y3(h,:)=[z1(t),z2(t),z3(t)];
c(t)=3;
elseif z1(t-1)>0
l=l+1;
Y4(l,:)=[z1(t),z2(t),z3(t)];
c(t)=4;
end
end
Y5=zeros(0,3);Y6=zeros(0,3);Y7=zeros(0,3);Y8=zeros(0,3);
end
if N(2,12)=='terminal'&N(4,1)=='terminal'&N(4,5)=='terminal'
%if the right in the second row ,the first,second,third and fourth node in the fourth row are all terminal nodes
P1=sym('[z1<=0z2<=0z3<=0]');
P2=sym('[z1<=0z2<=0z3>0]');
P3=sym('[z1<=0z2>0z3<=0]');
P4=sym('[z1<=0z2>0z3>0]');
P5=sym('[z1>0]');
P=[P1,P2,P3,P4,P5];
f=0;g=0;h=0;l=0;o=0;p=0;q=0;r=0;
for t=2:m
if z1(t-1)<=0&z2(t-1)<=0&z3(t-1)<=0
f=f+1;
Y1(f,:)=[z1(t),z2(t),z3(t)];
c(t)=1;
elseif z1(t-1)<=0&z2(t-1)<=0&z3(t-1)>0
g=g+1;
Y2(g,:)=[z1(t),z2(t),z3(t)];
c(t)=2;
elseif z1(t-1)<=0&z2(t-1)>0&z3(t-1)<=0
h=h+1;
Y3(h,:)=[z1(t),z2(t),z3(t)];
c(t)=3;
elseif z1(t-1)<=0&z2(t-1)>0&z3(t-1)>0
l=l+1;
Y4(l,:)=[z1(t),z2(t),z3(t)];
c(t)=4;
elseif z1(t-1)>0
o=o+1;
Y5(o,:)=[z1(t),z2(t),z3(t)];
c(t)=5;
end
end
Y6=zeros(0,3);Y7=zeros(0,3);Y8=zeros(0,3);
end
% there is no terminal node in the first and second row of the tree
if N(3,2)=='terminal'&N(3,6)=='terminal'&N(3,10)=='terminal'&N(3,14)=='terminal'
%if the first,second,third and fourth node in the third row are
%terminal node
P1=sym('[z1<=0z2<=0]');
P2=sym('[z1<=0z2>0]');
P3=sym('[z1>0z2<=0]');
P4=sym('[z1>0z2>0]');
P=[P1,P2,P3,P4];
f=0;g=0;h=0;l=0;o=0;p=0;q=0;r=0;
for t=2:m
if z1(t-1)<=0&z2(t-1)<=0
f=f+1;
Y1(f,:)=[z1(t),z2(t),z3(t)];
c(t)=1;
elseif z1(t-1)<=0&z2(t-1)>0
g=g+1;
Y2(g,:)=[z1(t),z2(t),z3(t)];
c(t)=2;
elseif z1(t-1)>0&z2(t-1)<=0
h=h+1;
Y3(h,:)=[z1(t),z2(t),z3(t)];
c(t)=3;
elseif z1(t-1)>0&z2(t-1)>0
l=l+1;
Y4(l,:)=[z1(t),z2(t),z3(t)];
c(t)=4;
end
end
Y5=zeros(0,3);Y6=zeros(0,3);Y7=zeros(0,3);Y8=zeros(0,3);
end
if N(4,1)=='terminal'&N(3,6)=='terminal'&N(3,10)=='terminal'&N(3,14)=='terminal'
%if the second,third and fourth node in the third row and the first and second node in the fourth row are
%terminal node
P1=sym('[z1<=0z2<=0z3<=0]');
P2=sym('[z1<=0z2<=0z3>0]');
P3=sym('[z1<=0z2>0]');
P4=sym('[z1>0z2<=0]');
P5=sym('[z1>0z2>0]');
P=[P1,P2,P3,P4,P5];
f=0;g=0;h=0;l=0;o=0;p=0;q=0;r=0;
for t=2:m
if z1(t-1)<=0&z2(t-1)<=0&z3(t-1)<=0
f=f+1;
Y1(f,:)=[z1(t),z2(t),z3(t)];
c(t)=1;
elseif z1(t-1)<=0&z2(t-1)<=0&z3(t-1)>0
g=g+1;
Y2(g,:)=[z1(t),z2(t),z3(t)];
c(t)=2;
elseif z1(t-1)<=0&z2(t-1)>0
h=h+1;
Y3(h,:)=[z1(t),z2(t),z3(t)];
c(t)=3;
elseif z1(t-1)>0&z2(t-1)<=0
l=l+1;
Y4(l,:)=[z1(t),z2(t),z3(t)];
c(t)=4;
elseif z1(t-1)>0&z2(t-1)>0
o=o+1;
Y5(o,:)=[z1(t),z2(t),z3(t)];
c(t)=5;
end
end
Y6=zeros(0,3);Y7=zeros(0,3);Y8=zeros(0,3);
end
if N(4,5)=='terminal'&N(3,2)=='terminal'&N(3,10)=='terminal'&N(3,14)=='terminal'
%if the first,third and fourth node in the third row and the third and fourth node in the fourth row are
%terminal node
P1=sym('[z1<=0z2<=0]');
P2=sym('[z1<=0z2>0z3<=0]');
P3=sym('[z1<=0z2>0z3>0]');
P4=sym('[z1>0z2<=0]');
P5=sym('[z1>0z2>0]');
P=[P1,P2,P3,P4,P5];
f=0;g=0;h=0;l=0;o=0;p=0;q=0;r=0;
for t=2:m
if z1(t-1)<=0&z2(t-1)<=0
f=f+1;
Y1(f,:)=[z1(t),z2(t),z3(t)];
c(t)=1;
elseif z1(t-1)<=0&z2(t-1)>0&z3(t-1)<=0
g=g+1;
Y2(g,:)=[z1(t),z2(t),z3(t)];
c(t)=2;
elseif z1(t-1)<=0&z2(t-1)>0&z3(t-1)>0
h=h+1;
Y3(h,:)=[z1(t),z2(t),z3(t)];
c(t)=3;
elseif z1(t-1)>0&z2(t-1)<=0
l=l+1;
Y4(l,:)=[z1(t),z2(t),z3(t)];
c(t)=4;
elseif z1(t-1)>0&z2(t-1)>0
o=o+1;
Y5(o,:)=[z1(t),z2(t),z3(t)];
c(t)=5;
end
end
Y6=zeros(0,3);Y7=zeros(0,3);Y8=zeros(0,3);
end
if N(4,9)=='terminal'&N(3,2)=='terminal'&N(3,6)=='terminal'&N(3,14)=='terminal'
%if the first,second and fourth node in the third row and the fifth and sixth node in the fourth row are
%terminal node
P1=sym('[z1<=0z2<=0]');
P2=sym('[z1<=0z2>0]');
P3=sym('[z1>0z2<=0z3<=0]');
P4=sym('[z1>0z2<=0z3>0]');
P5=sym('[z1>0z2>0]');
P=[P1,P2,P3,P4,P5];
f=0;g=0;h=0;l=0;o=0;p=0;q=0;r=0;
for t=2:m
if z1(t-1)<=0&z2(t-1)<=0
f=f+1;
Y1(f,:)=[z1(t),z2(t),z3(t)];
c(t)=1;
elseif z1(t-1)<=0&z2(t-1)>0
g=g+1;
Y2(g,:)=[z1(t),z2(t),z3(t)];
c(t)=2;
elseif z1(t-1)>0&z2(t-1)<=0&z3(t-1)<=0
h=h+1;
Y3(h,:)=[z1(t),z2(t),z3(t)];
c(t)=3;
elseif z1(t-1)>0&z2(t-1)<=0&z3(t-1)>0
l=l+1;
Y4(l,:)=[z1(t),z2(t),z3(t)];
c(t)=4;
elseif z1(t-1)>0&z2(t-1)>0
o=o+1;
Y5(o,:)=[z1(t),z2(t),z3(t)];
c(t)=5;
end
end
Y6=zeros(0,3);Y7=zeros(0,3);Y8=zeros(0,3);
end
if N(4,13)=='terminal'&N(3,2)=='terminal'&N(3,6)=='terminal'&N(3,10)=='terminal'
%if the first,second and third node in the third row and the seventh and eighth node in the fourth row are
%terminal node
P1=sym('[z1<=0z2<=0]');
P2=sym('[z1<=0z2>0]');
P3=sym('[z1>0z2<=0]');
P4=sym('[z1>0z2>0z3<=0]');
P5=sym('[z1>0z2>0z3>0]');
P=[P1,P2,P3,P4,P5];
f=0;g=0;h=0;l=0;o=0;p=0;q=0;r=0;
for t=2:m
if z1(t-1)<=0&z2(t-1)<=0
f=f+1;
Y1(f,:)=[z1(t),z2(t),z3(t)];
c(t)=1;
elseif z1(t-1)<=0&z2(t-1)>0
g=g+1;
Y2(g,:)=[z1(t),z2(t),z3(t)];
c(t)=2;
elseif z1(t-1)>0&z2(t-1)<=0
h=h+1;
Y3(h,:)=[z1(t),z2(t),z3(t)];
c(t)=3;
elseif z1(t-1)>0&z2(t-1)>0&z3(t-1)<=0
l=l+1;
Y4(l,:)=[z1(t),z2(t),z3(t)];
c(t)=4;
elseif z1(t-1)>0&z2(t-1)>0&z3(t-1)>0
o=o+1;
Y5(o,:)=[z1(t),z2(t),z3(t)];
c(t)=5;
end
end
Y6=zeros(0,3);Y7=zeros(0,3);Y8=zeros(0,3);
end
if N(4,1)=='terminal'&N(4,5)=='terminal'&N(3,10)=='terminal'&N(3,14)=='terminal'
%if the third and fourth node in the third row and the first,second,third and fourth node in the fourth row are
%terminal node
P1=sym('[z1<=0z2<=0z3<=0]');
P2=sym('[z1<=0z2<=0z3>0]');
P3=sym('[z1<=0z2>0z3<=0]');
P4=sym('[z1<=0z2>0z3>0]');
P5=sym('[z1>0z2<=0]');
P6=sym('[z1>0z2>0]');
P=[P1,P2,P3,P4,P5,P6];
f=0;g=0;h=0;l=0;o=0;p=0;q=0;r=0;
for t=2:m
if z1(t-1)<=0&z2(t-1)<=0&z3(t-1)<=0
f=f+1;
Y1(f,:)=[z1(t),z2(t),z3(t)];
c(t)=1;
elseif z1(t-1)<=0&z2(t-1)<=0&z3(t-1)>0
g=g+1;
Y2(g,:)=[z1(t),z2(t),z3(t)];
c(t)=2;
elseif z1(t-1)<=0&z2(t-1)>0&z3(t-1)<=0
h=h+1;
Y3(h,:)=[z1(t),z2(t),z3(t)];
c(t)=3;
elseif z1(t-1)<=0&z2(t-1)>0&z3(t-1)>0
l=l+1;
Y4(l,:)=[z1(t),z2(t),z3(t)];
c(t)=4;
elseif z1(t-1)>0&z2(t-1)<=0
o=o+1;
Y5(o,:)=[z1(t),z2(t),z3(t)];
c(t)=5;
elseif z1(t-1)>0&z2(t-1)>0
p=p+1;
Y6(p,:)=[z1(t),z2(t),z3(t)];
c(t)=6;
end
end
Y7=zeros(0,3);Y8=zeros(0,3);
end
if N(4,1)=='terminal'&N(4,9)=='terminal'&N(3,6)=='terminal'&N(3,14)=='terminal'
%if the second and fourth node in the third row and the first,second,fifth and sixth node in the fourth row are
%terminal node
P1=sym('[z1<=0z2<=0z3<=0]');
P2=sym('[z1<=0z2<=0z3>0]');
P3=sym('[z1<=0z2>0]');
P4=sym('[z1>0z2<=0z3<=0]');
P5=sym('[z1>0z2<=0z3>0]');
P6=sym('[z1>0z2>0]');
P=[P1,P2,P3,P4,P5,P6];
f=0;g=0;h=0;l=0;o=0;p=0;q=0;r=0;
for t=2:m
if z1(t-1)<=0&z2(t-1)<=0&z3(t-1)<=0
f=f+1;
Y1(f,:)=[z1(t),z2(t),z3(t)];
c(t)=1;
elseif z1(t-1)<=0&z2(t-1)<=0&z3(t-1)>0
g=g+1;
Y2(g,:)=[z1(t),z2(t),z3(t)];
c(t)=2;
elseif z1(t-1)<=0&z2(t-1)>0
h=h+1;
Y3(h,:)=[z1(t),z2(t),z3(t)];
c(t)=3;
elseif z1(t-1)>0&z2(t-1)<=0&z3(t-1)<=0
l=l+1;
Y4(l,:)=[z1(t),z2(t),z3(t)];
c(t)=4;
elseif z1(t-1)>0&z2(t-1)<=0&z3(t-1)>0
o=o+1;
Y5(o,:)=[z1(t),z2(t),z3(t)];
c(t)=5;
elseif z1(t-1)>0&z2(t-1)>0
p=p+1;
Y6(p,:)=[z1(t),z2(t),z3(t)];
c(t)=6;
end
end
Y7=zeros(0,3);Y8=zeros(0,3);
end
if N(4,1)=='terminal'&N(4,13)=='terminal'&N(3,6)=='terminal'&N(3,10)=='terminal'
%if the second and third node in the third row and the first,second,seventh and eighth node in the fourth row are
%terminal node
P1=sym('[z1<=0z2<=0z3<=0]');
P2=sym('[z1<=0z2<=0z3>0]');
P3=sym('[z1<=0z2>0]');
P4=sym('[z1>0z2<=0]');
P5=sym('[z1>0z2>0z3<=0]');
P6=sym('[z1>0z2>0z3>0]');
P=[P1,P2,P3,P4,P5,P6];
f=0;g=0;h=0;l=0;o=0;p=0;q=0;r=0;
for t=2:m
if z1(t-1)<=0&z2(t-1)<=0&z3(t-1)<=0
f=f+1;
Y1(f,:)=[z1(t),z2(t),z3(t)];
c(t)=1;
elseif z1(t-1)<=0&z2(t-1)<=0&z3(t-1)>0
g=g+1;
Y2(g,:)=[z1(t),z2(t),z3(t)];
c(t)=2;
elseif z1(t-1)<=0&z2(t-1)>0
h=h+1;
Y3(h,:)=[z1(t),z2(t),z3(t)];
c(t)=3;
elseif z1(t-1)>0&z2(t-1)<=0
l=l+1;
Y4(l,:)=[z1(t),z2(t),z3(t)];
c(t)=4;
elseif z1(t-1)>0&z2(t-1)>0&z3(t-1)<=0
o=o+1;
Y5(o,:)=[z1(t),z2(t),z3(t)];
c(t)=5;
elseif z1(t-1)>0&z2(t-1)>0&z3(t-1)>0
p=p+1;
Y6(p,:)=[z1(t),z2(t),z3(t)];
c(t)=6;
end
end
Y7=zeros(0,3);Y8=zeros(0,3);
end
if N(3,2)=='terminal'&N(4,5)=='terminal'&N(4,9)=='terminal'&N(3,14)=='terminal'
%if the first and fourth node in the third row and the third,fourth,fifth and sixth node in the fourth row are
%terminal node
P1=sym('[z1<=0z2<=0]');
P2=sym('[z1<=0z2>0z3<=0]');
P3=sym('[z1<=0z2>0z3>0]');
P4=sym('[z1>0z2<=0z3<=0]');
P5=sym('[z1>0z2<=0z3>0]');
P6=sym('[z1>0z2>0]');
P=[P1,P2,P3,P4,P5,P6];
f=0;g=0;h=0;l=0;o=0;p=0;q=0;r=0;
for t=2:m
if z1(t-1)<=0&z2(t-1)<=0
f=f+1;
Y1(f,:)=[z1(t),z2(t),z3(t)];
c(t)=1;
elseif z1(t-1)<=0&z2(t-1)>0&z3(t-1)<=0
g=g+1;
Y2(g,:)=[z1(t),z2(t),z3(t)];
c(t)=2;
elseif z1(t-1)<=0&z2(t-1)>0&z3(t-1)>0
h=h+1;
Y3(h,:)=[z1(t),z2(t),z3(t)];
c(t)=3;
elseif z1(t-1)>0&z2(t-1)<=0&z3(t-1)<=0
l=l+1;
Y4(l,:)=[z1(t),z2(t),z3(t)];
c(t)=4;
elseif z1(t-1)>0&z2(t-1)<=0&z3(t-1)>0
o=o+1;
Y5(o,:)=[z1(t),z2(t),z3(t)];
c(t)=5;
elseif z1(t-1)>0&z2(t-1)>0
p=p+1;
Y6(p,:)=[z1(t),z2(t),z3(t)];
c(t)=6;
end
end
Y7=zeros(0,3);Y8=zeros(0,3);
end
if N(3,2)=='terminal'&N(4,5)=='terminal'&N(3,10)=='terminal'&N(4,13)=='terminal'
%if the first and third node in the third row and the third,fourth,seventh and eighth node in the fourth row are
%terminal node
P1=sym('[z1<=0z2<=0]');
P2=sym('[z1<=0z2>0z3<=0]');
P3=sym('[z1<=0z2>0z3>0]');
P4=sym('[z1>0z2<=0]');
P5=sym('[z1>0z2>0z3<=0]');
P6=sym('[z1>0z2>0z3>0]');
P=[P1,P2,P3,P4,P5,P6];
f=0;g=0;h=0;l=0;o=0;p=0;q=0;r=0;
for t=2:m
if z1(t-1)<=0&z2(t-1)<=0
f=f+1;
Y1(f,:)=[z1(t),z2(t),z3(t)];
c(t)=1;
elseif z1(t-1)<=0&z2(t-1)>0&z3(t-1)<=0
g=g+1;
Y2(g,:)=[z1(t),z2(t),z3(t)];
c(t)=2;
elseif z1(t-1)<=0&z2(t-1)>0&z3(t-1)>0
h=h+1;
Y3(h,:)=[z1(t),z2(t),z3(t)];
c(t)=3;
elseif z1(t-1)>0&z2(t-1)<=0
l=l+1;
Y4(l,:)=[z1(t),z2(t),z3(t)];
c(t)=4;
elseif z1(t-1)>0&z2(t-1)>0&z3(t-1)<=0
o=o+1;
Y5(o,:)=[z1(t),z2(t),z3(t)];
c(t)=5;
elseif z1(t-1)>0&z2(t-1)>0&z3(t-1)>0
p=p+1;
Y6(p,:)=[z1(t),z2(t),z3(t)];
c(t)=6;
end
end
Y7=zeros(0,3);Y8=zeros(0,3);
end
if N(3,2)=='terminal'&N(3,6)=='terminal'&N(4,9)=='terminal'&N(4,13)=='terminal'
%if the first and second node in the third row and the fifth,sixth,seventh and eighth node in the fourth row are
%terminal node
P1=sym('[z1<=0z2<=0]');
P2=sym('[z1<=0z2>0]');
P3=sym('[z1>0z2<=0z3<=0]');
P4=sym('[z1>0z2<=0z3>0]');
P5=sym('[z1>0z2>0z3<=0]');
P6=sym('[z1>0z2>0z3>0]');
P=[P1,P2,P3,P4,P5,P6];
f=0;g=0;h=0;l=0;o=0;p=0;q=0;r=0;
for t=2:m
if z1(t-1)<=0&z2(t-1)<=0
f=f+1;
Y1(f,:)=[z1(t),z2(t),z3(t)];
c(t)=1;
elseif z1(t-1)<=0&z2(t-1)>0
g=g+1;
Y2(g,:)=[z1(t),z2(t),z3(t)];
c(t)=2;
elseif z1(t-1)>0&z2(t-1)<=0&z3(t-1)<=0
h=h+1;
Y3(h,:)=[z1(t),z2(t),z3(t)];
c(t)=3;
elseif z1(t-1)>0&z2(t-1)<=0&z3(t-1)>0
l=l+1;
Y4(l,:)=[z1(t),z2(t),z3(t)];
c(t)=4;
elseif z1(t-1)>0&z2(t-1)>0&z3(t-1)<=0
o=o+1;
Y5(o,:)=[z1(t),z2(t),z3(t)];
c(t)=5;
elseif z1(t-1)>0&z2(t-1)>0&z3(t-1)>0
p=p+1;
Y6(p,:)=[z1(t),z2(t),z3(t)];
c(t)=6;
end
end
Y7=zeros(0,3);Y8=zeros(0,3);
end
if N(3,2)=='terminal'&N(4,5)=='terminal'&N(4,9)=='terminal'&N(4,13)=='terminal'
%if the first node in the third row and the third,fourth,fifth,sixth,seventh and eighth node in the fourth row are
%terminal node
P1=sym('[z1<=0z2<=0]');
P2=sym('[z1<=0z2>0z3<=0]');
P3=sym('[z1<=0z2>0z3>0]');
P4=sym('[z1>0z2<=0z3<=0]');
P5=sym('[z1>0z2<=0z3>0]');
P6=sym('[z1>0z2>0z3<=0]');
P7=sym('[z1>0z2>0z3>0]');
P=[P1,P2,P3,P4,P5,P6,P7];
f=0;g=0;h=0;l=0;o=0;p=0;q=0;r=0;
for t=2:m
if z1(t-1)<=0&z2(t-1)<=0
f=f+1;
Y1(f,:)=[z1(t),z2(t),z3(t)];
c(t)=1;
elseif z1(t-1)<=0&z2(t-1)>0&z3(t-1)<=0
g=g+1;
Y2(g,:)=[z1(t),z2(t),z3(t)];
c(t)=2;
elseif z1(t-1)<=0&z2(t-1)>0&z3(t-1)>0
h=h+1;
Y3(h,:)=[z1(t),z2(t),z3(t)];
c(t)=3;
elseif z1(t-1)>0&z2(t-1)<=0&z3(t-1)<=0
l=l+1;
Y4(l,:)=[z1(t),z2(t),z3(t)];
c(t)=4;
elseif z1(t-1)>0&z2(t-1)<=0&z3(t-1)>0
o=o+1;
Y5(o,:)=[z1(t),z2(t),z3(t)];
c(t)=5;
elseif z1(t-1)>0&z2(t-1)>0&z3(t-1)<=0
p=p+1;
Y6(p,:)=[z1(t),z2(t),z3(t)];
c(t)=6;
elseif z1(t-1)>0&z2(t-1)>0&z3(t-1)>0
q=q+1;
Y7(q,:)=[z1(t),z2(t),z3(t)];
c(t)=7;
end
end
Y8=zeros(0,3);
end
if N(4,1)=='terminal'&N(3,6)=='terminal'&N(4,9)=='terminal'&N(4,13)=='terminal'
%if the second node in the third row and the first,second,fifth,sixth,seventh and eighth node in the fourth row are
%terminal node
P1=sym('[z1<=0z2<=0z3<=0]');
P2=sym('[z1<=0z2<=0z3>0]');
P3=sym('[z1<=0z2>0]');
P4=sym('[z1>0z2<=0z3<=0]');
P5=sym('[z1>0z2<=0z3>0]');
P6=sym('[z1>0z2>0z3<=0]');
P7=sym('[z1>0z2>0z3>0]');
P=[P1,P2,P3,P4,P5,P6,P7];
f=0;g=0;h=0;l=0;o=0;p=0;q=0;r=0;
for t=2:m
if z1(t-1)<=0&z2(t-1)<=0&z3(t-1)<=0
f=f+1;
Y1(f,:)=[z1(t),z2(t),z3(t)];
c(t)=1;
elseif z1(t-1)<=0&z2(t-1)<=0&z3(t-1)>0
g=g+1;
Y2(g,:)=[z1(t),z2(t),z3(t)];
c(t)=2;
elseif z1(t-1)<=0&z2(t-1)>0
h=h+1;
Y3(h,:)=[z1(t),z2(t),z3(t)];
c(t)=3;
elseif z1(t-1)>0&z2(t-1)<=0&z3(t-1)<=0
l=l+1;
Y4(l,:)=[z1(t),z2(t),z3(t)];
c(t)=4;
elseif z1(t-1)>0&z2(t-1)<=0&z3(t-1)>0
o=o+1;
Y5(o,:)=[z1(t),z2(t),z3(t)];
c(t)=5;
elseif z1(t-1)>0&z2(t-1)>0&z3(t-1)<=0
p=p+1;
Y6(p,:)=[z1(t),z2(t),z3(t)];
c(t)=6;
elseif z1(t-1)>0&z2(t-1)>0&z3(t-1)>0
q=q+1;
Y7(q,:)=[z1(t),z2(t),z3(t)];
c(t)=7;
end
end
Y8=zeros(0,3);
end
if N(4,1)=='terminal'&N(4,5)=='terminal'&N(3,10)=='terminal'&N(4,13)=='terminal'
%if the third node in the third row and the first,second,third,fourth,seventh and eighth node in the fourth row are
%terminal node
P1=sym('[z1<=0z2<=0z3<=0]');
P2=sym('[z1<=0z2<=0z3>0]');
P3=sym('[z1<=0z2>0z3<=0]');
P4=sym('[z1<=0z2>0z3>0]');
P5=sym('[z1>0z2<=0]');
P6=sym('[z1>0z2>0z3<=0]');
P7=sym('[z1>0z2>0z3>0]');
P=[P1,P2,P3,P4,P5,P6,P7];
f=0;g=0;h=0;l=0;o=0;p=0;q=0;r=0;
for t=2:m
if z1(t-1)<=0&z2(t-1)<=0&z3(t-1)<=0
f=f+1;
Y1(f,:)=[z1(t),z2(t),z3(t)];
c(t)=1;
elseif z1(t-1)<=0&z2(t-1)<=0&z3(t-1)>0
g=g+1;
Y2(g,:)=[z1(t),z2(t),z3(t)];
c(t)=2;
elseif z1(t-1)<=0&z2(t-1)>0&z3(t-1)<=0
h=h+1;
Y3(h,:)=[z1(t),z2(t),z3(t)];
c(t)=3;
elseif z1(t-1)<=0&z2(t-1)>0&z3(t-1)>0
l=l+1;
Y4(l,:)=[z1(t),z2(t),z3(t)];
c(t)=4;
elseif z1(t-1)>0&z2(t-1)<=0
o=o+1;
Y5(o,:)=[z1(t),z2(t),z3(t)];
c(t)=5;
elseif z1(t-1)>0&z2(t-1)>0&z3(t-1)<=0
p=p+1;
Y6(p,:)=[z1(t),z2(t),z3(t)];
c(t)=6;
elseif z1(t-1)>0&z2(t-1)>0&z3(t-1)>0
q=q+1;
Y7(q,:)=[z1(t),z2(t),z3(t)];
c(t)=7;
end
end
Y8=zeros(0,3);
end
if N(4,1)=='terminal'&N(4,5)=='terminal'&N(4,9)=='terminal'&N(3,14)=='terminal'
%if the fourth node in the third row and the first,second,third,fourth,fifth and sixth node in the fourth row are
%terminal node
P1=sym('[z1<=0z2<=0z3<=0]');
P2=sym('[z1<=0z2<=0z3>0]');
P3=sym('[z1<=0z2>0z3<=0]');
P4=sym('[z1<=0z2>0z3>0]');
P5=sym('[z1>0z2<=0z3<=0]');
P6=sym('[z1>0z2<=0z3>0]');
P7=sym('[z1>0z2>0]');
P=[P1,P2,P3,P4,P5,P6,P7];
f=0;g=0;h=0;l=0;o=0;p=0;q=0;r=0;
for t=2:m
if z1(t-1)<=0&z2(t-1)<=0&z3(t-1)<=0
f=f+1;
Y1(f,:)=[z1(t),z2(t),z3(t)];
c(t)=1;
elseif z1(t-1)<=0&z2(t-1)<=0&z3(t-1)>0
g=g+1;
Y2(g,:)=[z1(t),z2(t),z3(t)];
c(t)=2;
elseif z1(t-1)<=0&z2(t-1)>0&z3(t-1)<=0
h=h+1;
Y3(h,:)=[z1(t),z2(t),z3(t)];
c(t)=3;
elseif z1(t-1)<=0&z2(t-1)>0&z3(t-1)>0
l=l+1;
Y4(l,:)=[z1(t),z2(t),z3(t)];
c(t)=4;
elseif z1(t-1)>0&z2(t-1)<=0&z3(t-1)<=0
o=o+1;
Y5(o,:)=[z1(t),z2(t),z3(t)];
c(t)=5;
elseif z1(t-1)>0&z2(t-1)<=0&z3(t-1)>0
p=p+1;
Y6(p,:)=[z1(t),z2(t),z3(t)];
c(t)=6;
elseif z1(t-1)>0&z2(t-1)>0
q=q+1;
Y7(q,:)=[z1(t),z2(t),z3(t)];
c(t)=7;
end
end
Y8=zeros(0,3);
end
%there is no terminal node in the first,second and third row in the
%tree
if N(4,1)=='terminal'&N(4,5)=='terminal'&N(4,9)=='terminal'&N(4,13)=='terminal'
%if the first,second,third,fourth,fifth,sixth,seventh and eighth node in the fourth row are
%terminal node
P1=sym('[z1<=0z2<=0z3<=0]');
P2=sym('[z1<=0z2<=0z3>0]');
P3=sym('[z1<=0z2>0z3<=0]');
P4=sym('[z1<=0z2>0z3>0]');
P5=sym('[z1>0z2<=0z3<=0]');
P6=sym('[z1>0z2<=0z3>0]');
P7=sym('[z1>0z2>0z3<=0]');
P8=sym('[z1>0z2>0z3>0]');
P=[P1,P2,P3,P4,P5,P6,P7,P8];
f=0;g=0;h=0;l=0;o=0;p=0;q=0;r=0;
for t=2:m
if z1(t-1)<=0&z2(t-1)<=0&z3(t-1)<=0
f=f+1;
Y1(f,:)=[z1(t),z2(t),z3(t)];
c(t)=1;
elseif z1(t-1)<=0&z2(t-1)<=0&z3(t-1)>0
g=g+1;
Y2(g,:)=[z1(t),z2(t),z3(t)];
c(t)=2;
elseif z1(t-1)<=0&z2(t-1)>0&z3(t-1)<=0
h=h+1;
Y3(h,:)=[z1(t),z2(t),z3(t)];
c(t)=3;
elseif z1(t-1)<=0&z2(t-1)>0&z3(t-1)>0
l=l+1;
Y4(l,:)=[z1(t),z2(t),z3(t)];
c(t)=4;
elseif z1(t-1)>0&z2(t-1)<=0&z3(t-1)<=0
o=o+1;
Y5(o,:)=[z1(t),z2(t),z3(t)];
c(t)=5;
elseif z1(t-1)>0&z2(t-1)<=0&z3(t-1)>0
p=p+1;
Y6(p,:)=[z1(t),z2(t),z3(t)];
c(t)=6;
elseif z1(t-1)>0&z2(t-1)>0&z3(t-1)<=0
q=q+1;
Y7(q,:)=[z1(t),z2(t),z3(t)];
c(t)=7;
elseif z1(t-1)>0&z2(t-1)>0&z3(t-1)>0
r=r+1;
Y8(r,:)=[z1(t),z2(t),z3(t)];
c(t)=8;
end
end
end