ОШИБКА MATLAB Feval требует дескриптор функции в качестве первого аргумента

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

»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`

compiler-errors arguments matlab newtons-method function-handle
person UsherUser077    schedule 28.03.2016    source источник

Ответы (1)


arrow_upward
1
arrow_downward

Вы передаете символическое выражение в функция, предназначенная для оценки анонимной функции, дескриптор функции или функция на пути Matlab с ее < href="http://www.mathworks.com/help/matlab/ref/feval.html#buyb9wa-2" rel="nofollow">имя, указанное строкой через 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);
person TroyHaskin    schedule 29.03.2016