0

私はこのコード (*) を持っています。

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

1 に答える 1

1

無名関数、関数ハンドル、または Matlab パス上の関数を評価するように設計された関数にシンボリック式を渡しています。その名前は、文字列によって示されます。Matlab で微分を行いたい場合は、最初にシンボリック式を使用してから、次のようにして無名関数に変換できます。fevalmatlabFunction

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);
于 2016-03-29T00:26:08.503 に答える