您正在关注此问题

您将在关注的内容源中看到更新。
您可能会收到电子邮件，具体取决于您的通信预设项。

How to formalize an optimization problem in Matlab?

1 次查看（过去 30 天）

显示更早的评论

reincornator 2022-5-3

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1710820-how-to-formalize-an-optimization-problem-in-matlab

⋮

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1710820-how-to-formalize-an-optimization-problem-in-matlab

评论： Matt J 2022-5-5

Hello.

How to issue a task in Matlab?

Given:

Vectors a1 and a2.

Find a vector of coefficients x such that:

abs((a2.')*x) -> min

abs ((a1.')*x) >= condition

under the conditions:

- Dimensions of vectors a1, a2 and x from 1 to p

- elements a1, a2 and x are complex numbers

- abs(a1_i)>0, abs(a2_i)>0, abs(x_i) = 1, where i =1..p

15 个评论
显示 13更早的评论隐藏 13更早的评论

Torsten 2022-5-3

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1710820-how-to-formalize-an-optimization-problem-in-matlab#comment_2138510

⋮

链接

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1710820-how-to-formalize-an-optimization-problem-in-matlab#comment_2138510

So you want to maximize abs(a1.'*x). And you want to minimize abs(a2.'*x).

How do you want to combine these two objectives in one ?

reincornator 2022-5-4

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1710820-how-to-formalize-an-optimization-problem-in-matlab#comment_2140240

⋮

链接

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1710820-how-to-formalize-an-optimization-problem-in-matlab#comment_2140240

It is necessary to find a vector x such that:

abs ((a1.')*x) -> max - goal

abs((a2.')*x) < tol, wheh tol->0 - limitations

Torsten 2022-5-4

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1710820-how-to-formalize-an-optimization-problem-in-matlab#comment_2140645

⋮

链接

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1710820-how-to-formalize-an-optimization-problem-in-matlab#comment_2140645

Seems you get a quadratic objective with quadratic constraint. Try fmincon.

Matt J 2022-5-4

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1710820-how-to-formalize-an-optimization-problem-in-matlab#comment_2140775

⋮

链接

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1710820-how-to-formalize-an-optimization-problem-in-matlab#comment_2140775

编辑：Matt J 2022-5-4

abs(x_i) = 1, where i =1..p

Are you sure it's not abs(x_i) <= 1 or sum(abs(xi).^2)=1? It would be hard to satisfy the constraint as posted.

For example, if a2=[1;0;0;0] it would be impossible to satisfy abs(a2.'*x)-->0, because abs(a2.'*x)=abs(x1) and abs(x1) is constrained to 1.

reincornator 2022-5-4

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1710820-how-to-formalize-an-optimization-problem-in-matlab#comment_2140950

⋮

链接

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1710820-how-to-formalize-an-optimization-problem-in-matlab#comment_2140950

Torsten, thanks. The "fmincon" function only works with real numbers. I use complex numbers.

reincornator 2022-5-4

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1710820-how-to-formalize-an-optimization-problem-in-matlab#comment_2140960

⋮

链接

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1710820-how-to-formalize-an-optimization-problem-in-matlab#comment_2140960

编辑：reincornator 2022-5-5

Matt J, thanks. You are right, there is another constraint: abs(a1_i)>0, abs(a2_i)>0, i=1..p.

Torsten 2022-5-4

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1710820-how-to-formalize-an-optimization-problem-in-matlab#comment_2140975

⋮

链接

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1710820-how-to-formalize-an-optimization-problem-in-matlab#comment_2140975

@reincornator

With each optimization software you use, you will have to split in real and imaginary part and thus handle two real-valued numbers simultaneously.

Matt J 2022-5-4

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1710820-how-to-formalize-an-optimization-problem-in-matlab#comment_2140990

⋮

链接

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1710820-how-to-formalize-an-optimization-problem-in-matlab#comment_2140990

编辑：Matt J 2022-5-4

在 MATLAB Online 中打开

You are right, there is another constraint: abs(a1_i)>0, abs(a1_i)>0, i=1..p.

I assume you meant abs(a1_i)>0, abs(a2_i)>0.

That's still not enough. Consider a2=[2,1].'. Then if abs(x1)=abs(x2)=1,

abs(a2.'*x)=abs(2*x1+x2)
            >=abs( 2*abs(x1) -abs(x2))
            =1

So, the constraint abs(a2.'*x)=0 can never be satisfied.

reincornator 2022-5-5

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1710820-how-to-formalize-an-optimization-problem-in-matlab#comment_2142070

⋮

链接

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1710820-how-to-formalize-an-optimization-problem-in-matlab#comment_2142070

@Matt J, thank you. Corrected the condition.

abs((a2^T)*x) < tol, tol-> 0

Matt J 2022-5-5

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1710820-how-to-formalize-an-optimization-problem-in-matlab#comment_2142165

⋮

链接

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1710820-how-to-formalize-an-optimization-problem-in-matlab#comment_2142165

My counter-example still applies. There is little guarantee that a solution will exist under those constraints.

Matt J 2022-5-5

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1710820-how-to-formalize-an-optimization-problem-in-matlab#comment_2142180

⋮

链接

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1710820-how-to-formalize-an-optimization-problem-in-matlab#comment_2142180

If you had abs(a2_i)=1 for all i=1...p. You might be able to guarantee existence of a solution for p>1.

reincornator 2022-5-5

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1710820-how-to-formalize-an-optimization-problem-in-matlab#comment_2142645

⋮

链接

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1710820-how-to-formalize-an-optimization-problem-in-matlab#comment_2142645

@Torsten, in real numbers, the condition abs(x_i)=1, can be written as Re(x_i)=sqrt(1-Im(x_i)^2). How to write such a condition for vector x in constraints?

Matt J 2022-5-5

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1710820-how-to-formalize-an-optimization-problem-in-matlab#comment_2142660

⋮

链接

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1710820-how-to-formalize-an-optimization-problem-in-matlab#comment_2142660

编辑：Matt J 2022-5-5

在 MATLAB Online 中打开

How to write such a condition for vector x in constraints?

Your problem wouldn't be written in terms of x_i in what Torsten is recommending. Your objective and constraints would be written in terms of independent variables (u_i, v_i) where u_i stands in for the real component of x_i and v_i stands in for the imaginary component. The constraint, written differentiably, would be,

u_i^2 + v_i^2 = 1

reincornator 2022-5-5

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1710820-how-to-formalize-an-optimization-problem-in-matlab#comment_2142800

⋮

链接

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1710820-how-to-formalize-an-optimization-problem-in-matlab#comment_2142800

@Torsten @Matt J

How to set 2 restrictions in the "fmincon" function?

Matt J 2022-5-5

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1710820-how-to-formalize-an-optimization-problem-in-matlab#comment_2142855

⋮

链接

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1710820-how-to-formalize-an-optimization-problem-in-matlab#comment_2142855

c and ceq can be vectors.

请先登录，再进行评论。

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

回答（1 个）

Matt J 2022-5-4

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1710820-how-to-formalize-an-optimization-problem-in-matlab#answer_957245

⋮

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1710820-how-to-formalize-an-optimization-problem-in-matlab#answer_957245

编辑：Matt J 2022-5-4

在 MATLAB Online 中打开

N=null(a2.');
a3=a1.'*N;
[~,idx]=max(abs(a3));
x=N(:,idx);

6 个评论
显示 4更早的评论隐藏 4更早的评论

reincornator 2022-5-4

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1710820-how-to-formalize-an-optimization-problem-in-matlab#comment_2140945

⋮

链接

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1710820-how-to-formalize-an-optimization-problem-in-matlab#comment_2140945

在 MATLAB Online 中打开

Ran in:

No, it doesn't work. In addition, no restrictions are set: abs(x) == 1.

randcomplex = @(x,y) rand(x,y).*exp(2i*pi*rand(x,y));

a1 = randcomplex(1,10);

a2 = randcomplex(1,10);

N=null(a2.')

N = 1×0 empty double row vector

a3=a1.'*N

a3 = 10×0 empty double matrix

[~,idx]=max(abs(a3))

idx = 1×0 empty double row vector

x=N(:,idx)

x = 1×0 empty double row vector

Matt J 2022-5-4

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1710820-how-to-formalize-an-optimization-problem-in-matlab#comment_2140985

⋮

链接

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1710820-how-to-formalize-an-optimization-problem-in-matlab#comment_2140985

编辑：Matt J 2022-5-4

在 MATLAB Online 中打开

Ran in:

randcomplex = @(x,y) rand(x,y).*exp(2i*pi*rand(x,y));
a1 = randcomplex(10,1);
a2 = randcomplex(10,1);
N=null(a2.');
a3=a1.'*N;
[maxval,idx]=max(abs(a3));
x=N(:,idx)*conj(a3(idx))/abs(a3(idx))
x = 
   0.1546 - 0.3938i
   0.0998 + 0.0391i
   0.0617 - 0.0385i
   0.0246 + 0.1474i
   0.1035 + 0.0910i
   0.0625 + 0.0216i
  -0.5747 - 0.6314i
   0.0177 + 0.0661i
   0.0371 - 0.0283i
   0.0319 + 0.1476i
a2.'*x
ans = -8.3267e-17 + 5.5511e-17i
maxval,
maxval = 0.7671
abs(a1.'*x)
ans = 0.7671