MATLAB ERROR Feval требует дескриптор функции в качестве первого аргумента

Question

У меня есть этот код (*) и когда я делаю:

»syms x
»newton_raphson({((5400.*(1 + x)*0^360) - (1450000.*x.*(1 + x).^360))}, diff(((5400.*(1 + x)*0^360) - (1450000.*x.*(1 + x).^360)),1), 0.001, eps, 5, 0.1)

эта ошибка появляется:

Ошибка при использовании аргумента feval должна содержать строку или function_handle.

Ошибка в newton_raphson (строка 10) fz = feval(f,z(1));

Как я могу исправить эту ошибку?

(*)

 function [raiz, zn, fz, i] = newton_raphson(f, flinha, x0, eps, iter_max, debug)
   if nargin < 4
      eps = 1e-6;
   end
   if nargin < 5
      iter_max = 1e3;
   end

   z(1) = x0;
   fz = feval(f,z(1));
   fzlinha = feval(flinha,z(1));
   if (nargin > 5 && debug > 0)
      fprintf(2,'i=%d z=%23.18G fz=%G fzlinha=%G\n',0,x0,fz,fzlinha);
   end
   for i = 1:iter_max
      if abs(fzlinha) == 0 % f'(x0) equal zero 
         disp('O valor da derivada em Xi não pode ser zero');
         z(i+1) = x0;
         return
      end
      z(i+1) = x0 - fz / fzlinha;
      fz = feval(f,z(i+1));
      fzlinha = feval(flinha,z(i+1));
      dif = abs(z(i+1) - x0);
      if (nargin > 5 && debug > 0)
         fprintf(2,'i=%d z=%23.18G fz=%G fzlinha=%G dif=%E\n',i,z(i+1),fz,fzlinha,dif);
      end
      if dif < eps
         break;
      elseif i == iter_max
         disp('Foi excedido o número máximo de iterações (iter_max)'); 
         break
      end 
      x0=z(i+1);
   end
   zn = z';
   raiz = z(i+1);
end`

Answer 1

Вы передаете символическое выражение в функцию, предназначенную для оценки анонимной функции, дескриптора функции или функции в пути Matlab, имя которой указано строкой через feval, Если вы хотите, чтобы Matlab сделал дифференцирование для вас, вы можете сначала использовать символические выражения, а затем преобразовать их в анонимные функции с помощью matlabFunction как это

syms x f Df
%
% Symbolic expressions
f  = (5400.*(1 + x)*0^360) - (1450000.*x.*(1 + x).^360); 
Df = diff(f,x); 
%
%  Convert to anonymous functions
f  = matlabFunction(f ,'Vars',x); 
Df = matlabFunction(Df,'Vars',x);