Производная Matlab возвращает константу

Question

Я определил функцию, которая возвращает производную от символьного выражения или дескриптора функции (плюс еще несколько).

function df = der(f) 
    if isa(f, 'cfit') || isa(f, 'sfit') 
     df = @(x) differentiate(f, x); 
    elseif isa(f, 'sym') || isa(f, 'function_handle') 
     syms r 
     F = sym(f); 
     df = matlabFunction(diff(F), 'Vars', r); 
    else 
     error('Not a derivative of a known type') 
    end 
end 

Это хорошо работает, но когда я питаюсь в строке, я получаю то, что не создает массивы.

>> df = der(@(r) r) 

df = 

    @(r)1.0 

>> df([1,2,3]) 

ans = 

    1 

Есть в любом случае, чтобы поймать ли выходная функция является постоянным и изменить выход, так что поведение,

>> df = der(@(r) r) 

df = 

    @(r)1.0*ones(size(r)) 

>> df([1,2,3]) 

ans = 

    1 1 1 

Answer 1

Вот одно решение.

function df = der(f) 
    if isa(f, 'cfit') || isa(f, 'sfit') 
     df = @(x) differentiate(f, x); 
    elseif isa(f, 'sym') || isa(f, 'function_handle') 
     syms r 
     F = sym(f); 
     df = matlabFunction(diff(F), 'Vars', r); 
    % These next four lines have been added: 
     c = df(ones(1,2));    
     if length(c) == 1 
      df = @(x) c*ones(size(x)); 
     end        
    else 
     error('Not a derivative of a known type') 
    end 
end 

Оценка df = der(@(r) 20*r); df([1,2,3]) теперь дает выходной сигнал

ans = 

    20 20 20 

Answer 2

Чтобы объяснить ошибку в коде, syms r создает скалярной символическую переменную. Насколько мне известно, нет надежного решения для получения векторизованной версии.

Простой "исправить", оценить значение каждого в отдельности:

>> arrayfun(df,[1,2,3]) 

ans = 

    1  1  1