Курсовая работа: Численные методы интегрирования и оптимизации сложных систем

K1=A*(A_X(k+1,:))';

K2=A*(A_X(k+1,:))'+1/3*K1;

K3=A*(A_X(k+1,:))'+1/6*K1+1/6*K2;

K4=A*(A_X(k+1,:))'+1/8*K1+3/8*K2;

K5=A*(A_X(k+1,:))'+1/2*K1-3/2*K3+2*K4;

A_X(k+2,:)=(A_X(k+1,:))+h/6*(K1'+4*K4'+K5');

else

h1=h;

t=t+h1;

H=(eye(length(A_X(1,:)))-(9*h1/24)*A);

G=(eye(length(A_X(1,:)))+19*h1/24*A)*(A_X(k+1,:))'+h1/24*A*(-5*(A_X(k,:))'+(A_X(k-1,:))')

+h1/24*B*(9*1+19*1-5*1);

A_X(k+2,:)=(inv(H)*G)';

end

Otr(k+1)=t;

k=k+1;

h=-0.43496

end

plot(Otr,A_X(1:k,1),'b-');

grid on

9. Приложение 4 (Листинг скриптов для спектрального анализа)

spectr.m

syms t T;

Kx=(638.89/2)*(t-T).^2-56.35*(1./2)*(-2*(t-T))+5.6;

Ky=(595.24/2)*(t-T).^2;

F2=2*t;

L(2)=F2;

F3=4*t.^2-1;

L(3)=F3;