You can not select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
27 lines
903 B
27 lines
903 B
% function ParameterRange()
|
|
aa=0.5:0.1:2;
|
|
bb=0.5:0.1:2;
|
|
cc= -5:0.1:1;
|
|
|
|
[a,b,c]=meshgrid(aa,bb,cc);
|
|
|
|
constraint1 = -a-b-c-1;
|
|
constraint2 = -a-b+c+1;
|
|
constraint3 = -a+b+c-1;
|
|
constraint4 = a-b+c-1;
|
|
|
|
|
|
domain = max(max(max(constraint1,constraint2),constraint3),constraint4)<=0;
|
|
|
|
disp([num2str(min(constraint1.*domain,[],'all')),' <= C1 <= ',num2str(max(constraint1.*domain,[],'all'))])
|
|
|
|
disp([num2str(min(constraint2.*domain,[],'all')),' <= C2 <= ',num2str(max(constraint2.*domain,[],'all'))])
|
|
|
|
disp([num2str(min(constraint3.*domain,[],'all')),' <= C3 <= ',num2str(max(constraint3.*domain,[],'all'))])
|
|
|
|
disp([num2str(min(constraint4.*domain,[],'all')),' <= C4 <= ',num2str(max(constraint4.*domain,[],'all'))])
|
|
|
|
isosurface(a,b,c,-max(max(max(constraint1,constraint2),constraint3),constraint4),0)
|
|
isocaps(a,b,c,-max(max(max(constraint1,constraint2),constraint3),constraint4),0)
|
|
axis equal;axis tight;
|
|
camlight
|