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 Vue3 大型图片数据拖动排序
  • ¥15 划分vlan后不通了
  • ¥15 GDI处理通道视频时总是带有白色锯齿
  • ¥20 用雷电模拟器安装百达屋apk一直闪退
  • ¥15 算能科技20240506咨询(拒绝大模型回答)
  • ¥15 自适应 AR 模型 参数估计Matlab程序
  • ¥100 角动量包络面如何用MATLAB绘制
  • ¥15 merge函数占用内存过大
  • ¥15 使用EMD去噪处理RML2016数据集时候的原理
  • ¥15 神经网络预测均方误差很小 但是图像上看着差别太大