MATLAB求解隐式微分方程（组）—符号计算函数solve

一般遇到的微分方程均可以转化为标准形式，然后用ode函数求解。实际上，有部分微分方程无法转化为标准形式。这里介绍一种 隐式微分方程（组）的求解方法——利用solve函数求出微分方程的显示表达。
需要注意：有时候有些隐式方程并不能通过符号计算函数solve求出显示表达式，遇到这种情况则需要其他方法求解。
隐式方程(组)如下：
y2'*(y2*cos(4*t)-y1^3)-t/5*y1'=0;
t*sin(y2)/8-2*y2*y2'+sqrt(t)*y1'=0;
求解范围 tspan = [1 20];
初始条件 y1(0)=1,y2(0)=1.
示例代码如下：
**********************************
完整代码请加文末企鹅号
function implicit_ODE_solve
%用符号积分求解dy的解析表达式
**************
tspan = [1 20];
y0 = [1;1];
[t,y] = ode45(DyDtAnony,tspan,y0);%调用ode45求解
figure;
plot(t,y(:,1),'k-');
hold on
plot(t,y(:,2),'k:');
%图例,位置自动选择最佳位置
L = legend('{\ity}_1(t)','{\ity}_2(t)','Location','best');
set(L,'fontname','Times New Roman');
xlabel('\itt','fontsize',16);
% --------------------------------------------
请不要做伸手党，示例完整代码收费10 RMB
——————