| 查看: 2596 | 回复: 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 --------------------------------------------------------------------------------------------------------- 谢谢！ |
回复此楼» 猜你喜欢
宿州学院学报
已经有3人回复
-
4,4二甲基联苯干啥用，有懂得吗
已经有3人回复
-
医学类期刊求推荐
已经有6人回复
-
26/27申博自荐
已经有10人回复
-
东北林业大学材料科学与工程学院“一流”A+学科国家级人才团队课题组招收2026级博士生
已经有3人回复
-
生活琐事由它去
已经有4人回复
-
提交了我也来说说感想
已经有12人回复
-
青B发送上会通知了吗
已经有9人回复
-
西安交大新媒学院副院长用撤稿论文结题
已经有6人回复
-
论文撤稿了
已经有8人回复