hongyu.. 2022-02-25 15:53 采纳率: 100%
浏览 91
已结题

matlab关于cvx的问题

问题遇到的现象和发生背景

cvx的SDP问题,我在约束条件中想写凸函数小于等于一个实数

问题相关代码,请勿粘贴截图

rk=1e-5;
M1=h_BI'*theta0'*(h_I1'*h_I1)*theta0*h_BI;
M2=h_BI'*theta0'*(h_I2'*h_I2)*theta0*h_BI;
W1t=eye(Nt,Nt);
W2t=eye(Nt,Nt);
Zxt=eye(Nt,Nt);
cvx_begin sdp
    variables W1x(Nt,Nt) W2x(Nt,Nt) 
    variable Zx(Nt,Nt) complex semidefinite
    variable p13  nonnegative
    variable p23  nonnegative
    variable p14  nonnegative
    variable p24  nonnegative
obj=-log(real(trace(Zx*M1))+sigmak^2+real(trace(M1*W1x))+real(trace(M1*W2x)))-log(real(trace(Zx*M2))+sigmak^2+real(trace(M2*W1x)))+real(trace(M2*W2x))-real(trace((-(1/log(2))*(M1/real((trace(Zxt*M1))+sigmak^2+real(trace(M1*W2t)))+M2/(real(trace(Zxt*M2))+sigmak^2+real(trace(M2*W1t)))))'*Zx))-real(trace((-(1/log(2))*M2/(real(trace(Zxt*M2))+sigmak^2+real(trace(M2*W1t))))'*W1x))-real(trace((-(1/log(2))*M1/real((trace(Zxt*M1))+sigmak^2+real(trace(M1*W2t))))'*W2x))-re;
minimize obj
    subject to 
    trace(W1x)+trace(Zx)<=P;
    trace(W2x)+trace(Zx)<=P;
-log(real(trace(Zx*M1))+sigmak^2+real(trace(M1*W1x))+real(trace(M1*W2x)))-log(real(trace(Zx*M2))+sigmak^2+real(trace(M2*W1x))+real(trace(M2*W2x)))-real(trace((-(1/log(2))*(M1/(real(trace(Zxt*M1))+sigmak^2+real(trace(M1*W2t)))+M2/real((trace(Zxt*M2))+sigmak^2+real(trace(M2*W1t))))))'*(Zx-Zxt))-real(trace((-(1/log(2))*M2/real((trace(Zxt*M2))+sigmak^2+real(trace(M2*W1t)))))'*(W1x-W1t))-real(trace((-(1/log(2))*M1/real((trace(Zxt*M1))+sigmak^2+real(trace(M1*W2t)))))'*(W2x-W2t))<=-rk;
[sigmak^2*(2^(re)-1)-p13,zeros(1,M_irs);zeros(M_irs,1),p13*epsilong_e^(-2)*eye(M_irs)]+(S3*theta0*h_BI)*((2^(re)-1)*Zx-W1x)*((S3*theta0*h_BI)')>=0;
 [sigmak^2*(2^(re)-1)-p23,zeros(1,M_irs);zeros(M_irs,1),p23*epsilong_e^(-2)*eye(M_irs)]+(S3*theta0*h_BI)*((2^(re)-1)*Zx-W2x)*((S3*theta0*h_BI)')>=0;
 [sigmak^2*(2^(re)-1)-p14,zeros(1,M_irs);zeros(M_irs,1),p14*epsilong_e^(-2)*eye(M_irs)]+(S4*theta0*h_BI)*((2^(re)-1)*Zx-W1x)*((S4*theta0*h_BI)')>=0;
 [sigmak^2*(2^(re)-1)-p24,zeros(1,M_irs);zeros(M_irs,1),p24*epsilong_e^(-2)*eye(M_irs)]+(S4*theta0*h_BI)*((2^(re)-1)*Zx-W2x)*((S4*theta0*h_BI)')>=0;
Zx>=0;
W1x>=0;
W2x>=0;
cvx_end
运行结果及报错内容

出错 main (第 105 行)


-log(real(trace(Zx*M1))+sigmak^2+real(trace(M1*W1x))+real(trace(M1*W2x)))-log(real(trace(Zx*M2))+sigmak^2+real(trace(M2*W1x))+real(trace(M2*W2x)))-real(trace((-(1/log(2))*(M1/(real(trace(Zxt*M1))+sigmak^2+real(trace(M1*W2t)))+M2/real((trace(Zxt*M2))+sigmak^2+real(trace(M2*W1t))))))'*(Zx-Zxt))-real(trace((-(1/log(2))*M2/real((trace(Zxt*M2))+sigmak^2+real(trace(M2*W1t)))))'*(W1x-W1t))-real(trace((-(1/log(2))*M1/real((trace(Zxt*M1))+sigmak^2+real(trace(M1*W2t)))))'*(W2x-W2t))<=-rk;

出错内容:
SDP constraint {matrix} <= {scalar} valid only if the scalar is zero.

我的解答思路和尝试过的方法

我把实数移到不等式左边,让整体小于等于0,但是报错为:Both sides of an SDP constraint must be affine.

我想要达到的结果

不报错

  • 写回答

2条回答 默认 最新

  • 有问必答小助手 2022-02-28 09:32
    关注

    你好,我是有问必答小助手,非常抱歉,本次您提出的有问必答问题,技术专家团超时未为您做出解答


    本次提问扣除的有问必答次数,已经为您补发到账户,我们后续会持续优化,扩大我们的服务范围,为您带来更好地服务。

    本回答被题主选为最佳回答 , 对您是否有帮助呢?
    评论
查看更多回答(1条)

报告相同问题?

问题事件

  • 系统已结题 3月8日
  • 已采纳回答 2月28日
  • 创建了问题 2月25日

悬赏问题

  • ¥15 对于相关问题的求解与代码
  • ¥15 ubuntu子系统密码忘记
  • ¥15 信号傅里叶变换在matlab上遇到的小问题请求帮助
  • ¥15 保护模式-系统加载-段寄存器
  • ¥15 电脑桌面设定一个区域禁止鼠标操作
  • ¥15 求NPF226060磁芯的详细资料
  • ¥15 使用R语言marginaleffects包进行边际效应图绘制
  • ¥20 usb设备兼容性问题
  • ¥15 错误(10048): “调用exui内部功能”库命令的参数“参数4”不能接受空数据。怎么解决啊
  • ¥15 安装svn网络有问题怎么办