把上边的代码写成可以在Visual Studio中实现的形式

好的，以下是将Matlab代码转换为C++代码的实现，可以在Visual Studio中执行：

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

const double PI = 3.14159265357898;
const double MU_0 = PI * 0.0000004;

int main()
{
    int numfre = 100;
    double fre[100];
    double z[4] = { 0, 500, 1000, 2000 };
    double rho[4] = { 10, 100, 100, 1000 };
    double rhoa[100];

    //生成一个对数等分的数组
    double start = 4;
    double end = -4;
    double step = (end - start) / (numfre - 1);
    for (int i = 0; i < numfre; i++)
    {
        fre[i] = pow(10, start + step * i);
    }

    //循环计算每个频率对应的电阻率
    for (int nf = 0; nf < numfre; nf++)
    {
        double kn[4];
        double Z0m[4];
        double Zm = 0.0;

        //计算每一层的本征阻抗
        for (int n = 0; n < 4; n++)
        {
            kn[n] = sqrt(-2.0 * PI * fre[nf] * MU_0 / rho[n]);
            Z0m[n] = -1.0 * 1i * 2 * PI * fre[nf] * MU_0 / kn[n];
        }
        Zm = Z0m[3]; //最后一层本征阻抗

        //从最后一层开始向上逐层计算电阻率
        for (int m = 2; m >= 0; m--)
        {
            double hm = z[m + 1] - z[m];
            double Lm = (Z0m[m] - Zm) / (Z0m[m] + Zm);
            Zm = Z0m[m] * (1 - Lm * exp(-2 * kn[m] * hm)) / (1 + Lm * exp(-2 * kn[m] * hm));
        }

        rhoa[nf] = 1.0 / (2 * PI * fre[nf] * MU_0) * pow(abs(Zm), 2);
    }

    //输出每个频率对应的电阻率
    for (int i = 0; i < numfre; i++)
    {
        cout << "Frequency: " << fre[i] << " Hz, Resistivity: " << rhoa[i] << " Ohm-m" << endl;
    }

    return 0;
}

需要注意的是，这里对一些Matlab函数和语法进行了适当的修改和替换，比如使用了C++中的pow函数代替了Matlab中的logspace函数，使用了cmath库中的数学函数等等。另外，由于C++中没有内置的复数类型，这里使用了实部和虚部分别存储复数，并且使用了1i表示虚数单位。