How to replace the sdpvar vector x with meshgrid variables to plot it? (YALMIP)

Question

DB 2022-3-30

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1684259-how-to-replace-the-sdpvar-vector-x-with-meshgrid-variables-to-plot-it-yalmip

回答： Saarthak Gupta 2023-12-24

在 MATLAB Online 中打开

I have polynomial f of min degree 2, max degree 6 with variables x1 and x2, which are defined in Matlab by:

x = sdpvar(2,1)

According to this comment, I set the following lines to prepare the polynomial for meshing in order to plot it:

f = f{1};
f = replace(f,'*','.*');                    %replace to dot product
f = replace(f,'^','.^');                    %replace to dot product
f = eval(['@(x)' f])                        %function handle @(x,y)
[X1,X2] = meshgrid(-1:0.01:1,-1:0.01:1);    %create meshgrid
f2 = f([X1;X2]);                            %define f2 as f with x(1)=X1 and x(2)=X2
mesh(X1,X2,f2);                             %plot f2 in X1 and X2

However I get an error because f2 becomes a scalar, instead of a matrix with all the possible solutions along X1 and X2. Therefore, I cannot plot it with mesh. I tried to follow the steps as given in this answer (link), but although both cases are with two variables, I cannot grasp it.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Saarthak Gupta 2023-12-24

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1684259-how-to-replace-the-sdpvar-vector-x-with-meshgrid-variables-to-plot-it-yalmip#answer_1377022

在 MATLAB Online 中打开

Hi,

Looks like you are trying to plot a “sdpvar” expression (YALMIP) over a given domain.

Your code does not correctly parametrize and convert the “sdpvar” expression to a function handle. Based on the given example, "p" represents a bi-variate polynomial in variables “x” and “y”, with a degree of 4. After completing the sum-of-squares decomposition and performing the required adjustments, the function handle is expected to accept two arguments – “x” and “y”. However, the error arises because you have provided only x as the input.

Refer to the following code:

sdpvar x y
[p,c] = polynomial([x y],4);    % generate a parameterized polynomial
% solve sum-of-squares decomposition problem using constraints:
% 1. sum-of-squares constraint defined by sos(p)
% 2. p(pi,2) = 3
solvesos([sos(p), replace(p,[x y],[pi 2]) == 3],[],[],c);
% display evaluated polynomial in symbolic MATLAB form.
f = sdisplay(replace(p,c,value(c)));
f = f{1};
f = replace(f,'^','.^');    % element-wise exponentiation
f = replace(f,'*','.*');    % element-wise multiplication
f = str2func(['@(x,y)' f]); % convert function string to function handle
[X,Y] = meshgrid(-1:0.01:1,-1:0.01:1);
F = f(X,Y);
mesh(X,Y,F);