解决非对称微分矩阵黎卡提微分方程
%
% dY (t) / dt = AY + YB - YCY + Q (*)
% Y (t0) = Y0
%
%使用向后微分公式的方法。
%
%输入:大小:矩阵(n, n)。
% B:尺寸矩阵(p, p)。
% C:尺寸矩阵(p, n)。
% D:尺寸矩阵(氮、磷)。
% Y0:尺寸矩阵(氮、磷)。
% tf, t0:临时工。
%输出:
% Ytf: inslant tf的解决方案(*)
% Y:解决方案(*)的t \ [t0 tf]
% t
%作者:LAKHLIFA SADEK。
%最后修改:29/09/2019
%的电子邮件:lakhlifasdek@gmail.com;sadek.l@ucd.ac.ma
% ORCID:https://orcid.org/0000 - 0001 - 9780 - 2592
%
%
%的例子
%
% t0 = 0, tf = 1
%
% A =
% 1 2
% 3 5
%
% B =
% 2 2 1
% 2 3 2
% 2 1 1
%
% C =
% 1 4
% 1 3
% 2 3
%
% D =
% 1 5 3
% 2 6 2
%
%
% Y0 =
%
% 1 4 3
% 1 2 3
%
% [t Y Ytf] = BDF_Diff_Riccati (A, B, C, D, Y0, t0, tf)
%
% Y (tf)的解决方案
%
% Ytf =
%
% 0.0790 0.3842 1.3520
% 1.2402 1.8543 0.7197
%
%为示踪剂的解决方案Y_ {1} fonctuion t (t)
%的阴谋(t Y (1,:))
%包含(“临时工t \ [0 1]”)
% ylabel('解决方案Y_ {1 1} (t)”)
%的传说(“快速公车提供方法”)。