Unable to plot inequality region
1 回表示 (過去 30 日間)
古いコメントを表示
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 件のコメント
採用された回答
KSSV
2022 年 6 月 18 日
編集済み: KSSV
2022 年 6 月 18 日
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 件のコメント
その他の回答 (0 件)
参考
カテゴリ
Help Center および File Exchange で Assumptions についてさらに検索
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!