24小时热门版块排行榜

返回列表

colored

铜虫 (初入文坛)

应助: 0 (幼儿园)
金币: 168.5
帖子: 25
在线: 10.7小时
虫号: 3525071
注册: 2014-11-07
性别: MM
专业: 通信理论与系统

[求助] 关于 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
---------------------------------------------------------------------------------------------------------
谢谢！

回复此楼

» 猜你喜欢

1楼 2016-03-24 20:07:02

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 colored 的主题更新

返回列表