Решение задач оптимального управления с использованием математической системы MATLAB и пакета имитационного моделирования SIMULINK. Сивохин А.В - 183 стр.

UptoLike

183
DiffDiffFi12=(2*(Dy1*(2*x-3*x^2)+y1*(x^2-x^3)+6*x^3-6*x^4))
%-- дифференцирование подинтегрального выражения
% функционала по а1;
IntDiffDiffFi11=int(DiffDiffFi11,x,0,1)
%-- нахождение производной функционала по a0;
IntDiffDiffFi12=int(DiffDiffFi12,x,0,1)
%-- нахождение производной функционала по a1;
[a0,a1]=solve(IntDiffDiffFi11,IntDiffDiffFi12,a0,a1)
%-- нахождение коэффициентов a0 и a1;
y1=subs(subs(y1,a0,'a0'),a1,'a1')
%-- нахождение выражения для второго приближения;
Dy1=subs(subs(Dy1,a0,'a0'),a1,'a1')
%-- нахождение выражения для производной второго
% приближения;
J1=int(Dy1^2+y1^2+12*x*y1,x,0,1)
%-- нахождение выражения для второго приближения
% функционала;
J1=vpa(J1,6)
%-- вычисление второго приближения функционала.
9.4 Построение имитационных моделей
На рис 9.1 представлена имитационная модель для нахождения экстремали
рассматриваемого функционала, а на рис. 9.2 – осциллограммы работы модели
после нахождения экстремали. Имитационная модель для вычисления
функционала и соответствующая осциллограмма приведены на рис. 9.3 и 9.4.
DiffDiffFi12=(2*(Dy1*(2*x-3*x^2)+y1*(x^2-x^3)+6*x^3-6*x^4))
             %-- дифференцирование подинтегрального выражения
             % функционала по а1;
IntDiffDiffFi11=int(DiffDiffFi11,x,0,1)
             %-- нахождение производной функционала по a0;
IntDiffDiffFi12=int(DiffDiffFi12,x,0,1)
             %-- нахождение производной функционала по a1;
[a0,a1]=solve(IntDiffDiffFi11,IntDiffDiffFi12,a0,a1)
             %-- нахождение коэффициентов a0 и a1;
y1=subs(subs(y1,a0,'a0'),a1,'a1')
            %-- нахождение выражения для второго приближения;
Dy1=subs(subs(Dy1,a0,'a0'),a1,'a1')
            %-- нахождение выражения для производной второго
            % приближения;
J1=int(Dy1^2+y1^2+12*x*y1,x,0,1)
            %-- нахождение выражения для второго приближения
            % функционала;
J1=vpa(J1,6)
            %-- вычисление второго приближения функционала.




                  9.4 Построение имитационных моделей

     На рис 9.1 представлена имитационная модель для нахождения экстремали
рассматриваемого функционала, а на рис. 9.2 – осциллограммы работы модели
после нахождения экстремали. Имитационная модель          для вычисления
функционала и соответствующая осциллограмма приведены на рис. 9.3 и 9.4.




                                    183