将符号表达式转换为Matlab中的函数句柄R2013b

Question

我很努力地将符合我的土木工程项目要求的函数从符号表达式转换为函数。我需要使用fzero来查找函数的根。这里H应该是变量，我需要找出H的值。该功能是这样

function x_c = f_x_c(s,H0,VA,Lo,qc,EAo,NF,Sj,Fj) 

if (s < 0) || (s > Lo) 
    disp('The value of s is invalid') 
    disp(['s = ' num2str(s)]); 
    return 
end 


C1 = H/qc; 
if NF == 0 
    n = 0; 
    sn = 0; 
    sum_Fj = 0; 
end 


if NF >= 1 
    Sj_Q = [0; Sj; Lo]; 
    %Determine n and sn if 0 <= s < Lo: 
    if s < Lo 
     STOP = 0; 
     k = 0; 
     while STOP == 0 
      k = k + 1; 
      if (s >= Sj_Q(k,1)) && (s < Sj_Q((k + 1),1)) 
       STOP = 1; 
      end 
     end 
     n = k - 1; 
     sn = Sj_Q(k,1); 
    end 
    %Determine n and sn if s = Lo: 
    if s == Lo 
     n = NF; 
     sn = Sj(NF,1); 
    end 
    sum_Fj = sum(Fj(1:n,1)); 
end 


x_c = (H/EAo)*s; 
x_c = x_c + C1*asinh((qc*s - VA + sum_Fj)/H) + ... 
    - C1*asinh((qc*sn - VA + sum_Fj)/H); 


for j = 1:n 
    sk = Sj_Q((j + 1),1); 
    sk_1 = Sj_Q(j,1); 
    sum_Fj = sum(Fj(1:(j - 1))); 


    x_c = x_c + ... 
     + C1*asinh((qc*sk - VA + sum_Fj)/H) + ... 
     - C1*asinh((qc*sk_1 - VA + sum_Fj)/H); 

end 

我想用这个f_x_c.m文件中的主要文件，其中我会找到这个方程的根。 有人可以指导我如何做到这一点？

我已经尝试使用下面的代码，但我没有成功。

if (s < 0) || (s > Lo) 
    disp('The value of s is invalid') 
    disp(['s = ' num2str(s)]); 
    return 
end 

C1 = @(H) (H/qc); 
if NF == 0 
    n = 0; 
    sn = 0; 
    sum_Fj = 0; 
end 


if NF >= 1 
    Sj_Q = [0; Sj; Lo]; 
    %Determine n and sn if 0 <= s < Lo: 
    if s < Lo 
     STOP = 0; 
     k = 0; 
     while STOP == 0 
      k = k + 1; 
      if (s >= Sj_Q(k,1)) && (s < Sj_Q((k + 1),1)) 
       STOP = 1; 
      end 
     end 
     n = k - 1; 
     sn = Sj_Q(k,1); 
    end 
    %Determine n and sn if s = Lo: 
    if s == Lo 
     n = NF; 
     sn = Sj(NF,1); 
    end 
    sum_Fj = sum(Fj(1:n,1)); 
end 


x_c =@(H) (H/EAo)*s; 
x_c =@(H) (x_c(H) + (C1(H))*asinh((qc*s - VA + sum_Fj)/H) + ... 
    - (C1(H))*asinh((qc*sn - VA + sum_Fj)/H)); 


for j = 1:n 
    sk = Sj_Q((j + 1),1); 
    sk_1 = Sj_Q(j,1); 
    sum_Fj = sum(Fj(1:(j - 1))); 


    x_c =@(H) (x_c(H) + ... 
     + C1(H)*asinh((qc*sk - VA + sum_Fj)/H) + ... 
     - C1(H)*asinh((qc*sk_1 - VA + sum_Fj)/H)); 

end 

编辑：

我要解决的主要文件中的以下公式：

equation = f_x_c(inext_length, H0, vertical_reaction, inext_length, qc, EAo, NF, hanger_arc_length, point_hanger_force) + 1400; 
% Whatever equation f_x_c returns, I have to add another number to it(like here it is 1400), then solve this equation using fzero. 

因此，在主文件中，我写这样的：

equation = @(H) f_x_c(inext_length, H0, vertical_reaction, inext_length, qc, EAo, NF, hanger_arc_length, point_hanger_force); 

equation = @(H) (equation(H) + 1400); 
answer = fsolve(equation, H0); 

Answer 1

您的问题的模拟答案可能看起来像

function x_c = f_x_c(H,A,B,C,D) 
    x_c = H*A; 
    x_c = x_c + B*asinh(C/H) - B*asinh(D/H); 
end 

，并求解一个呼叫

H = fzero(@(H)(f_x_c(H,1,1,1,1)+1400),1);