求助：使用MATLAB求解耦合数值微分方程问题

对应的微分方程图示：

原代码存在的关键问题 #

• 符号函数与数值求解接口错误：通过matlabFunction转换符号ODE后，直接取M(1)传入ode45是错误的，转换后的函数应返回两个方程的导数向量，而非单个方程。
• 初始条件矛盾：代码中定义A1(0)=1、A2(0)=0，但调用ode45时传入的初始值是[2 0]，两者完全不符。
• 冗余的数组运算符号：符号运算中使用.*和.^是多余的，符号环境下直接用*和^即可实现对应运算。
• 复数处理逻辑不清晰：未明确处理复数方程的传递逻辑，可能导致ode45无法正确识别复数导数。

修正后的MATLAB代码 #

clc;
clear;

% 定义物理常数
lambda1 = 1560e-9;
c = 3e8;
d_eff = 27e-12;
omega1 = 2*pi*c/lambda1;
omega2 = 2*omega1;
n = 2.2;
k1 = n*omega1/c;
k2 = n*omega2/c;
delta_k = 2*k1 - k2; % 计算相位失配项

% 预计算方程中的常数系数，简化后续表达式
coeff1 = (2*1i*omega1^2*d_eff)/(k1*c^2);
coeff2 = (1i*omega2^2*d_eff)/(k2*c^2);

% 定义ODE数值函数：状态向量y = [A1; A2]，返回导数向量
dAdz = @(z, y) [coeff1 * y(2) * conj(y(1)) * exp(-1i*delta_k*z);
                coeff2 * y(1)^2 * exp(1i*delta_k*z)];

% 设置初始条件与求解区间
y0 = [1; 0]; % A1(0)=1, A2(0)=0
zspan = [0, 20];

% 使用ode45求解复数耦合ODE
sol = ode45(dAdz, zspan, y0);

% 生成采样点并计算数值解
z_eval = linspace(zspan(1), zspan(2), 1000);
y_eval = deval(sol, z_eval);
A1_eval = y_eval(1,:);
A2_eval = y_eval(2,:);

% 绘制解的幅值曲线（复数解的幅值通常为核心关注量）
figure;
plot(z_eval, abs(A1_eval), 'LineWidth', 1.5);
hold on;
plot(z_eval, abs(A2_eval), 'LineWidth', 1.5);
xlabel('z');
ylabel('|A(z)|');
legend('A1(z)','A2(z)');
grid on;
title('耦合微分方程解的幅值分布');

修正说明 #

1. 直接定义数值ODE函数：跳过符号转换步骤，避免接口错误，同时提升计算效率。
2. 统一初始条件：确保初始值与问题定义一致，传入[1; 0]。
3. 预计算常数项：将所有固定常数提前计算，简化ODE表达式，增强可读性。
4. 可视化结果：通过deval获取求解区间内的数值解，绘制幅值曲线直观展示能量转换过程。

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