怎么用matlab程序画出文中图4、图5、和图11中需要matlab程序画的图，请给出完整的matlab程序

要使用MATLAB绘制文中的图4、图5 和 图11，我们需要根据提供的数据和模型来编写相应的代码。以下是每个图的完整MATLAB程序示例：

### 图4：Junction Capacitance and SPICE Simulated Qgd

#### 数据准备
假设我们有以下数据文件：
- `junction_capacitance_data.mat` 包含变量 `Vds` 和 `C_xss`（Ciss, Coss, Crss）
- `spice_qgd_data.mat` 包含变量 `Vds` 和 `Qgd_dy`

% 加载数据
load('junction_capacitance_data.mat');
load('spice_qgd_data.mat');

% 计算 Cgd_dy
Cgd_dy = gradient(Qgd_dy, mean(diff(Vds)));

% 绘制图4(a)
figure;
plot(Vds, C_xss, 'LineWidth', 2);
xlabel('V_{DS} (V)');
ylabel('Capacitance (F)');
legend('C_{iss}', 'C_{oss}', 'C_{rss}');
title('Junction Capacitances of SiC MOSFET C3M0075120K');
grid on;

% 绘制图4(b)
figure;
plot(Vds, Qgd_dy, 'LineWidth', 2);
xlabel('V_{DS} (V)');
ylabel('Charge (C)');
title('SPICE Simulated Q_{gd}-dy as a Function of V_{DS}');
grid on;

### 图5：Fit Threshold Voltage and Transfer Characteristics

#### 数据准备
假设我们有以下数据文件：
- `threshold_voltage_data.mat` 包含变量 `Tj` 和 `Vth_fit`
- `transfer_characteristics_data.mat` 包含变量 `Vgs` 和 `Ich_fit`, `Ich_spice`

% 加载数据
load('threshold_voltage_data.mat');
load('transfer_characteristics_data.mat');

% 绘制图5(a)
figure;
plot(Tj, Vth_fit, 'LineWidth', 2);
hold on;
plot(Tj, Vth_datasheet, 'r--', 'LineWidth', 2);
xlabel('Temperature (°C)');
ylabel('Threshold Voltage (V)');
legend('Fit Curve', 'Datasheet');
title('Fit Threshold Voltage as a Function of T_j');
grid on;

% 绘制图5(b)
figure;
plot(Vgs, Ich_fit, 'LineWidth', 2);
hold on;
plot(Vgs, Ich_spice, 'r--', 'LineWidth', 2);
xlabel('Gate-Source Voltage (V)');
ylabel('Channel Current (A)');
legend('Fit Curve', 'SPICE Model');
title('Transfer Characteristics of the DUTs');
grid on;

### 图11：Energy Losses as a Function of Gate Resistor and Drain Current

#### 数据准备
假设我们有以下数据文件：
- `energy_loss_data.mat` 包含变量 `Rg`, `Id`, `E_on_exp`, `E_on_model`, `E_on_spice`, `E_on_model_6`

% 加载数据
load('energy_loss_data.mat');

% 绘制图11(a) E_on vs. Rg
figure;
subplot(2, 1, 1);
plot(Rg, E_on_exp, 'ko-', 'LineWidth', 2, 'MarkerFaceColor', 'k');
hold on;
plot(Rg, E_on_model, 'b^-', 'LineWidth', 2, 'MarkerFaceColor', 'b');
plot(Rg, E_on_spice, 'rs-', 'LineWidth', 2, 'MarkerFaceColor', 'r');
plot(Rg, E_on_model_6, 'g*-', 'LineWidth', 2, 'MarkerFaceColor', 'g');
xlabel('Gate Resistor (\Omega)');
ylabel('Turn-On Energy Loss (mJ)');
legend('Experimental', 'Proposed Model', 'SPICE Model', 'Model [6]');
title('E_{on} vs. R_g for IMZ120R030M1H');
grid on;

% 绘制图11(b) E_on vs. Id
subplot(2, 1, 2);
plot(Id, E_on_exp, 'ko-', 'LineWidth', 2, 'MarkerFaceColor', 'k');
hold on;
plot(Id, E_on_model, 'b^-', 'LineWidth', 2, 'MarkerFaceColor', 'b');
plot(Id, E_on_spice, 'rs-', 'LineWidth', 2, 'MarkerFaceColor', 'r');
plot(Id, E_on_model_6, 'g*-', 'LineWidth', 2, 'MarkerFaceColor', 'g');
xlabel('Drain Current (A)');
ylabel('Turn-On Energy Loss (mJ)');
legend('Experimental', 'Proposed Model', 'SPICE Model', 'Model [6]');
title('E_{on} vs. I_d for IMZ120R030M1H');
grid on;

### 注意事项
1. **数据文件**：确保所有数据文件都已准备好并保存在当前工作目录中。
2. **变量命名**：确保数据文件中的变量名称与代码中使用的变量名称一致。
3. **绘图样式**：可以根据需要调整线条颜色、线型、标记等以符合具体需求。

这些代码将生成类似于文章中所示的图形。如果数据格式或变量名称不同，可能需要相应地调整代码。