MCR-WPT系统，保持其它参数不变，在改变发射侧电感情况下，利用 MATLAB 仿真阻抗角与频率、发射线圈电感的变化的关系图像 matlab三维

以下是一个简单的示例代码，可以用于仿真MCR-WPT系统的阻抗角与频率、发射线圈电感的变化关系：

% 参数设置
C1 = 10e-9; % 电容值
R1 = 10; % 电阻值
L1 = 10e-6; % 发射线圈电感值
L2 = [5e-6, 7.5e-6, 10e-6, 12.5e-6, 15e-6]; % 接收线圈电感值
f = logspace(0, 6, 1000); % 频率范围

% 计算阻抗角
[Z, theta] = MCR_WPT_impedance(C1, R1, L1, L2, f);

% 绘制三维图像
figure;
surf(L2, f, theta);
xlabel('发射线圈电感值 (H)');
ylabel('频率 (Hz)');
zlabel('阻抗角 (rad)');

其中，`MCR_WPT_impedance` 函数是一个自定义函数，用于计算MCR-WPT系统的阻抗角。以下是一个简单的实现：

function [Z, theta] = MCR_WPT_impedance(C1, R1, L1, L2, f)
% 计算MCR-WPT系统的阻抗角

Z = zeros(length(L2), length(f));
theta = zeros(length(L2), length(f));

for i = 1:length(L2)
    for j = 1:length(f)
        omega = 2 * pi * f(j);
        Z1 = R1 + 1i * omega * L1;
        Z2 = 1i * omega * L2(i);
        ZC = -1i / (omega * C1);
        Zeq = Z1 + Z2 + ZC;
        Z(i, j) = abs(Zeq);
        theta(i, j) = angle(Zeq);
    end
end

end

这段代码会生成一个三维图像，横轴为发射线圈电感值，纵轴为频率，纵轴为阻抗角。你可以通过修改参数来探索不同情况下的阻抗角与频率、发射线圈电感的变化关系。