Question

MATLAB解常微分方程组

function di=f(t,i)
global u1 R1 R2 R3 M12 M13 M23 L1 L2 L3
di=zeros(3,1);
di(1)=(u1-R1*i(1)-M12*di(2)-M13*di(3))/L1;
di(2)=(0-R2*i(2)-M12*di(1)-M23*di(3))/L2;
di(3)=(0-R3*i(3)-M13*di(1)-M23*di(2))/L3;
end

Answer 1

你的常微分方程组的自定义函数表述有问题。

问题一：像下列三个中，不能在等式的右边出现di(1)，di(2)，di(3)变量。

di(1)=(u1-R1*i(1)-M12*di(2)-M13*di(3))/L1;
di(2)=(0-R2*i(2)-M12*di(1)-M23*di(3))/L2;
di(3)=(0-R3*i(3)-M13*di(1)-M23*di(2))/L3;
问题二：如用ode（）求方程的数值解，u1、R1、M12、M13、L1、R2、M23、L2、R3、M13、L3必须是具体的数值，还需要有微分方程组的初始条件。

问题三：如用dsolve（）求方程的解析解，也存在变量过多的问题。