Eigen(C++)复数操作技术问询：矩阵向量与复数常量相乘及元素赋值

嘿，刚好我对Eigen的复数操作挺熟的，来帮你解决这个赋值的问题～

你提到kx.real()(0, 0) = 1.0;能正常运行，但kz_r.imag()(0, ii) =1.0;编译报错，大概率是两种原因：要么kz_r是const限定的矩阵（此时imag()返回的是不可修改的const表达式），要么kz_r是个临时对象（比如某个函数直接返回的矩阵，你没存到变量里就链式调用imag()）。

下面给你几种靠谱的复数矩阵元素赋值方式，按需选用：

1. 直接赋值完整的复数对象 #

这是最直观不容易出错的方式，直接给矩阵元素赋一个std::complex<T>类型的值：

#include <complex>
#include <Eigen/Dense>

// 定义一个2x2的double类型复数矩阵
Eigen::MatrixXcd kz_r(2, 2);
int ii = 0;

// 方式1：构造复数对象赋值
kz_r(0, ii) = std::complex<double>(0.0, 1.0);

// 方式2：C++14及以上可用复数字面量（需要引入命名空间）
using namespace std::literals::complex_literals;
kz_r(0, ii) = 0.0 + 1.0i;

2. 单独修改已有元素的实部/虚部 #

如果你只想修改某个元素的虚部（或实部），直接定位到元素再修改更清晰：

// 方式a：用imag()成员函数直接赋值
kz_r(0, ii).imag(1.0);

// 方式b：访问虚部的引用赋值
kz_r(0, ii).imag() = 1.0;

对比你之前的写法，这种是先定位到单个元素再修改属性，而非先取整个矩阵的虚部矩阵再修改元素，能避开很多const或临时对象的坑。

3. 批量修改整个矩阵的虚部 #

如果需要给整个矩阵的虚部统一赋值一个实矩阵，可以这么写：

Eigen::MatrixXd imag_matrix(2, 2);
imag_matrix.setConstant(1.0); // 把所有元素设为1.0
kz_r.imag() = imag_matrix; // 给kz_r的整个虚部赋值

注意这种写法要求kz_r是可修改的（非const），且不是临时对象。

总的来说，优先用前两种方式，代码更清晰也不容易踩坑～

内容的提问来源于stack exchange，提问作者user3533030