这段代码运行显示：no math for operator 第26行无法运行

可能是因为在第26行使用了`exp`函数，但是程序中没有包含`math.h`头文件，导致编译器无法识别`exp`函数。你可以在程序开头添加`#include <cmath>`，来包含数学函数库，或者使用`std::exp`来调用`exp`函数。具体代码如下：

#include <iostream>
#include <complex>
#include <cmath>
using namespace std;

int main() {
    const int numfre = 100;
    double fre[numfre];
    complex<double> Z0m[4], Zm;
    const int nc = 4;
    double z[nc] = {0, 500, 1000, 2000};
    double rho[nc] = {10, 100, 100, 1000};
    const double wPI = 3.14159265357898;
    const double MU_0 = wPI * 0.0000004;
    complex<double> kn[nc];

    // 计算频率数组
    for (int i = 0; i < numfre; i++) {
        fre[i] = exp(log(10) * (4 - (-4)) / (numfre - 1) * i - 4);
    }

    // 计算每一层的本征阻抗
    for (int nf = 0; nf < numfre; nf++) {
        for (int n = 0; n < nc; n++) {
            kn[n] = sqrt(-complex<double>(0, 1) * 2 * wPI * fre[nf] * MU_0 / rho[n]);
            Z0m[n] = -complex<double>(0, 1) * 2 * wPI * fre[nf] * MU_0 / kn[n];
        }

        // 计算每一层的传输矩阵
        Zm = Z0m[nc - 1];
        for (int m = nc - 2; m >= 0; m--) {
            double hm = z[m + 1] - z[m];
            complex<double> Lm = (Z0m[m] - Zm) / (Z0m[m] + Zm);
            Zm = Z0m[m] * (1 - Lm * std::exp(-2 * kn[m] * hm)) / (1 + Lm * std::exp(-2 * kn[m] * hm));
        }

        // 计算每个频率下的电阻率
        double rhoa = (1 / (2 * wPI * fre[nf] * MU_0)) * norm(Zm);
        cout << fre[nf] << " " << rhoa << endl;
    }
    return 0;
}

这样修改后，就可以正确地运行了。