| 查看: 2632 | 回复: 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 --------------------------------------------------------------------------------------------------------- 谢谢！ |
回复此楼» 猜你喜欢
澳大利亚 Murdoch University 全奖博士招生（3个名额）地质化工冶金领域
已经有36人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有28人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有10人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有10人回复
-
阿姨
已经有4人回复
-
一作与独作在应聘高校教师时区别大吗
已经有4人回复
-
静等基金结果
已经有21人回复
-
FileCode能看出啥？
已经有28人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有6人回复
-
filecode
已经有16人回复
