Unable to plot inequality region
    4 views (last 30 days)
  
       Show older comments
    
This is the solution I am trying to mimic, but for my case it does not seem to work:
syms ki kd
cond = [ -9*ki, 10*ki - 10*kd + 10, 360*kd - 24*ki + 1320]; % My 3 inequality conditions
r = -5:0.01:5;
[ki, kd] = meshgrid(r);
condition1 = cond(1) > 0;
condition2 = cond(2) > 0;
condition3 = cond(3) > 0;
output = ones(length(r));
output(~(condition1 & condition2 & condition3)) = 0; % this fails
imshow(output, 'xdata', r, 'ydata', r);
axis on;
the line that fails gives me this error:
Error using symengine
Unable to prove '~0 < -9*ki | ~0 < 10*ki - 10*kd + 10 |
~0 < 360*kd - 24*ki + 1320' literally. Use 'isAlways' to
test the statement mathematically.
Error in sym/subsindex (line 841)
                X =
                find(mupadmex('symobj::logical',A.s,9)) -
                1;
How can I solve this? It seems like it should be such a simple thing to do, but I can't make it work.
0 Comments
Accepted Answer
  KSSV
      
      
 on 18 Jun 2022
        
      Edited: KSSV
      
      
 on 18 Jun 2022
  
      r = -5:0.01:5;
[ki, kd] = meshgrid(r);
cond =  -9*ki>0 & 10*ki - 10*kd + 10>0 & 360*kd - 24*ki + 1320>0; % My 3 inequality conditions
% output = ones(length(r));
% output(~(condition1 & condition2 & condition3)) = 0; % this fails
imshow(cond, 'xdata', r, 'ydata', r);
axis on;
3 Comments
More Answers (0)
See Also
Categories
				Find more on Assumptions in Help Center and File Exchange
			
	Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!


