matlab inline函数使用错误问题

大海绵1990

n和t都是随机生成的10000行两列的数据。我就是想把我的那么多变量循环代入我那个复杂的公式当中然后经过牛顿迭代把公式里面的x算出来然后把每次迭代的x放进一个向量里。总是提示错误使用inline函数，自己改了很多次都不对。请大神指点！谢谢谢谢谢谢！
这是主程序
syms x p s1 k1 s2 k2
quan=[0.5:0.001:0.999]; k=length(quan);
for i=1:k
    for j=1:10000
fun=inline('[p-((23/89)*(1-exp(-((x-86.0175)./s1).^k1))+(66/89)*(1-exp(-((x-86.0175)./s2).^k2))), diff(p-((23/89)*(1-exp(-((x-86.0175)./s1).^k1))+(66/89)*(1-exp(-((x-86.0175)./s2).^k2))))]','p','x','s1','k1','s2','k2');
p=quan(i);s1=n(i,1);k1=n(i,2);s2=t(i,1);k2=t(i,2);
[x_star,index,it] = newton(fun,100);
  y(j,i)=x_star;
    end
end
这是调用牛顿迭代的函数
function [x_star,index,it] = newton(fun,x,ep,it_max)
%求解非线性方程的牛顿法
%第一个分量是函数值，第二个分量是导数值
% x为初始点
% ep为精度,当 | x(k)-x(k-1) |<ep时，终止计算,缺省值为1e-5
% it_max为最大迭代次数,缺省值为100
% x_star为当迭代成功时，输出方程的根
%  当迭代失败，输出最后的迭代值
% index为指标变量，当index=1时，表明迭代成功
% 当index=0时，表明迭代失败（迭代次数>=it_max）
% it为迭代次数
if nargin<4 it_max=100;end                              
if nargin<3 ep=1e-5;end
index=0;k=1;
while k<it_max
x1=x;f=feval(fun,x);
x=x-f(1)/f(2);
    if abs(x-x1)<ep
        index=1;break;
    end
    k=k+1;
end
x_star=x;it=k;


错误提示 错误使用 inline/feval (line 22)
内联函数的输入数目不足。

出错 newton (line 16)
x1=x;f=feval(fun,x);

出错 ci_culculation (line 11)
[x_star,index,it] = newton(fun,100);

13top

https://ww2.mathworks.cn/help/matlab/ref/inline.html

jstzysq

inline和符号计算相当低效，你还把它放在二重循环里，就不要期待能给你算出结果。
把符号计算的步骤放在循环外面，只计算一次。用匿名函数的技巧将计算结果传给循环，这样速度快了n个数量级。

1. clc,clear all,close all
   
2. 
   
3. %执行必要的符号计算
   
4. syms x p s1 k1 s2 k2
   
5. fun=p-((23/89)*(1-exp(-((x-86.0175)./s1).^k1))+(66/89)*(1-exp(-((x-86.0175)./s2).^k2)));
   
6. dfun=diff(fun,x); %导数
   
7. 
   
8. %技巧：将符号计算结果转化成匿名函数
   
9. eval(['fun=@(x,p,s1,k1,s2,k2) ',vectorize(fun),';']);
   
10. eval(['dfun=@(x,p,s1,k1,s2,k2) ',vectorize(dfun),';']);
    
11. 
    
12. %{
    
13. 匿名函数求值的示例用法
    
14. x=86;p=3;s1=3;k1=2;s2=5;k2=2;
    
15. fun(x,p,s1,k1,s2,k2)
    
16. dfun(x,p,s1,k1,s2,k2)
    
17. %}
    
18. 
    
19. %你原来的程序。略有改动
    
20. quan=[0.5:0.001:0.999]; k=length(quan);
    
21. for i=1:k
    
22.     for j=1:10000
    
23.     p=quan(i);s1=n(i,1);k1=n(i,2);s2=t(i,1);k2=t(i,2);
    
24.     [x_star,index,it] = newton(fun,dfun,100,p,s1,k1,s2,k2);
    
25.     y(j,i)=x_star;
    
26.     end
    
27. end
    
28. 
    
29. return
    
30. 
    
31. 
    
32. 
    
33. 
    
34. %newton函数，有改动
    
35. function [x_star,index,it] = newton(fun,dfun,x, p,s1,k1,s2,k2)
    
36.     it_max=100 ;                           
    
37.     ep=1e-5;
    
38.     index=0;k=1;
    
39.     while k<it_max
    
40.         x1=x;
    
41.         f1= fun(x1,p,s1,k1,s2,k2); %匿名函数求值，速度相当快
    
42.         f2=dfun(x1,p,s1,k1,s2,k2);
    
43.         x=x1-f1/f2;
    
44.         if abs(x-x1)<ep
    
45.             index=1;break;
    
46.         end
    
47.         k=k+1;
    
48.     end
    
49.     x_star=x;it=k;
    
50. end

复制代码