| 查看: 2575 | 回复: 0 | ||
colored铜虫 (初入文坛)
|
[求助]
关于 fsolve 函数的使用
|
|
楼主现在正在做一个仿真，遇到一些困难，请各位路过的高手帮帮忙！ 问题如下：主函数产生六个参数lamdad, miud, etad，lamdae,miue,etae，function 需要这六个参数中的一些，请问要怎么读取这些值，而且在主函数中还要用fsolve函数来求解一个三元的非线性方程组（方程组为F）。 ------------------------------------------------------------------------------------------------------ 主函数： for ii = 1 : cd +dd [lamdad, miud, etad] = fun_myunified(ii, ad, bd, cd, dd); for jj = 1 : ce + de [lamdae,miue,etae] = fun_myunified(jj, ae, be, ce, de); ------------------------------------------------------------------------------------------------------ function output = fun（*,*,*） F = [x(1)*gamma(x(2)+1/x(3))/(x(2)^(1/x(3))*gamma(x(2))) - 2* miue.*miud .* (betae./etad).^((miud-miue)./2) .* average_gg((miud+miue)/2 , miud-miue, xi); x(1)^2*gamma(x(2)+2/x(3))/(x(2)^(2/x(3))*gamma(x(2))) - 2* miue.*miud .* (betae./etad).^((miud-miue)./2) .* average_gg((miud+miue)/2+1 , miud-miue, xi); x(1)^3*gamma(x(2)+3/x(3))/(x(2)^(3/x(3))*gamma(x(2))) - 2* miue.*miud .* (betae./etad).^((miud-miue)./2) .* average_gg((miud+miue)/2+2 , miud-miue, xi)]; ---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- 在主函数中又需要用fsolve解出来的三个参数。 ------------------------------------------------------------------------------------------------------- 主函数： for ii = 1 : cd +dd [lamdad, miud, etad] = fun_myunified(ii, ad, bd, cd, dd); for jj = 1 : ce + de [lamdae,miue,etae] = fun_myunified(jj, ae, be, ce, de); x0 = [1; 1; 1]; % Make a starting guess at the solution options = optimoptions('fsolve','Display','iter'); % Option to display output [x1,fval1] = fsolve(@fun,x0,options); % Call solver --------------------------------------------------------------------------------------------------------- 谢谢！ |
回复此楼» 猜你喜欢
材料调剂
已经有6人回复
-
化工求调剂
已经有7人回复
-
321求调剂
已经有10人回复
-
材料专硕322分
已经有13人回复
-
材料295
已经有3人回复
-
11408 359分 求调剂
已经有3人回复
-
320分，材料与化工专业，求调剂
已经有24人回复
-
286求调剂
已经有7人回复
-
化学调剂
已经有10人回复
-
357求调剂
已经有3人回复
120