Ошибка при использовании *- внутренняя размерная матрица должна совпадать

Question

У меня возникла проблема с задачей, где я должен сравнить методы решения ОДУ: ode45, методы Эйлера и Гаусса-Лежандра. Здесь я должен рассчитать ошибки для разных шагов.

   h=[0.01 0.05 0.1 0.5];
   func = @(t, y) -2*y+t*sin(t);
   opts = odeset('Reltol',1e-13,'AbsTol',1e-14,'Stats','on');
   f1_f2_matrix = [];
   sinus=@(t) t*(sin(t));
   error2nd = zeros(length(h),1);
   errorInf = zeros(length(h),1);
   error2ndEuler = zeros(length(h),1);
   errorInfEuler = zeros(length(h),1);


   for i=1:length(h)
   t=0:h(i):10;
   euler = zeros(length(t),1);
   fun_builtin = zeros(length(t),1);
   t_builtin = zeros(length(t),1);
   results = zeros(length(t),1);


   %matrix for A=-2 and given h=0.01

   matrix=[ 1-1/2*h, (1/2-sqrt(3)/3)*h; 
  (1/2+sqrt(3)/3)*h, 1-1/2*h];


    for n=2:length(t)
    B = [-2*results(n-1) + sinus(t(n-1)+(1/2-sqrt(3)/6)*h);
        -2*results(n-1) + sinus(t(n-1) + (1/2+sqrt(3)/6)*h)];

    %system of equations solution

    f1_f2_matrix = matrix\B; 

    results(n) = (results(n-1) + h*1/2*(f1_f2_matrix(1) + f1_f2_matrix(2))); 

    %euler
    euler(n) = euler(n-1)+h(i)*func(t(n-1),euler(n-1));

   end


   %ODE45 FUNCTION
[t_builtin,fun_builtin] = ode45(func, t, 0, opts);


   %y'=-2y+tsin(t) errors
 error2nd(i)=norm(fun_builtin-results)/norm(fun_builtin); %root mean square error of my function
errorInf(i)=norm((fun_builtin-results), Inf)/ norm((fun_builtin), Inf); %maximum error of my function

error2ndEuler(i)=norm(euler-results)/norm(fun_builtin); %root mean square error of my function (euler)
errorInfEuler(i)=norm((euler-results), Inf)/ norm((fun_builtin), Inf); %maximum error of my function (euler)

   end

   figure semilogy(h, error2nd, h, errorInf,h, error2ndEuler, h, errorInfEuler)
   title([{('Dependence of root mean square and maximum errors')}; {('from h-step and method of solving differential equation')}])
legend('root mean square error: Gauss-Legendre of 4th order','maximum error: Gauss-Legendre of 4th order method','root mean square error: Euler method', 'maximum error: Euler method', 'Location', 'southeast');
xlabel('step');
ylabel('error'); 

Он должен генерировать график сравнения ошибок для разных методов на разных шагах. Но появляются такие ошибки:

  Error using  * 
  Inner matrix dimensions must agree.

  Error in untitled>@(t)t*(sin(t))

  Error in untitled (line 27)
    B = [-2*results(n-1) + sinus(t(n-1)+(1/2-sqrt(3)/6)*h);

Я не могу найти ошибку. Я пытался применить точку ('.'), Но это не сработало. Где я сделал ошибку?

Answer 1

Использование точки .* обязательно решит проблему. После этого возникает вторая проблема: вы не индексировали переменную h, Как только вы это сделаете, код будет работать отлично.

Answer 2

h=[0.01 0.05 0.1 0.5 1];
func = @(t, y) -2*y+t*sin(t);
opts = odeset('Reltol',1e-13,'AbsTol',1e-14,'Stats','on');
f1_f2_matrix = [];
sinus=@(t) t*(sin(t));
error2nd = zeros(length(h),1);
errorInf = zeros(length(h),1);
error2ndEuler = zeros(length(h),1);
errorInfEuler = zeros(length(h),1);


for i=1:length(h)
t=0:h(i):10;
euler = zeros(length(t),1);
fun_builtin = zeros(length(t),1);
t_builtin = zeros(length(t),1);
results = zeros(length(t),1);


%matrix for A=-2 and given h=0.01
matrix=[ 1-1/2*h(i), (1/2-sqrt(3)/3)*h(i); 
        (1/2+sqrt(3)/3)*h(i), 1-1/2*h(i)];


    for n=2:length(t)
    B = [-2*results(n-1) + sinus(t(n-1)+(1/2-sqrt(3)/6)*h(i));
         -2*results(n-1) + sinus(t(n-1) + (1/2+sqrt(3)/6)*h(i))];

%system of equations solution

    f1_f2_matrix = matrix\B; 

    results(n) = (results(n-1) + h(i)*0.5*(f1_f2_matrix(1) + f1_f2_matrix(2))); 

%euler
    euler(n) = euler(n-1)+h(i)*func(t(n-1),euler(n-1));

   end


%ODE45 FUNCTION
   [t_builtin,fun_builtin] = ode45(func, t, 0, opts);


%y'=-2y+tsin(t) errors
error2nd(i)=norm(fun_builtin-results)/norm(fun_builtin); %root mean square error of my function
errorInf(i)=norm((fun_builtin-results), Inf)/ norm((fun_builtin), Inf); %maximum error of my function

error2ndEuler(i)=norm(euler-results)/norm(fun_builtin); %root mean square error of my function (euler)
errorInfEuler(i)=norm((euler-results), Inf)/ norm((fun_builtin), Inf); %maximum error of my function (euler)

end

figure 
semilogy(h, error2nd, h, errorInf,h, error2ndEuler, h, errorInfEuler)
title([{('Dependence of root mean square and maximum errors')}; {('from h-step and method of solving differential equation')}])
legend('root mean square error: Gauss-Legendre of 4th order','maximum error: Gauss-Legendre of 4th order method','root mean square error: Euler method', 'maximum error: Euler method', 'Location', 'southeast');
xlabel('step');
ylabel('error');