shuhuiz 发表于 2011-6-3 09:32:41

求助解方程组

请高手指点啊,下面的程序解不到想要的结果,请大家帮忙看看!
%不用m文件法;
R0=zeros(1,281);
n=zeros(1,281); ;
Hv=zeros(1,281);;

syms H;         
for k=1:281,   
   H=k+199;
   f=@(x)([((0.08*log(exp(x(1))+0.002))/(x(1)*exp(x(1))-0.998*x(1)))^x(2)-H*(exp(x(1))+0.002)/(630000*x(1));
210000*x(1)*(x(2)/(x(1)*exp(1)))^x(2)-3*H-22.3]);
                  options=optimset('MaxFunEvals',1000);
   x=fsolve(f,,options);

messenger 发表于 2011-6-6 21:39:04

参考http://forum.simwe.com/thread-902236-1-1.html,http://forum.simwe.com/thread-816392-1-1.html

shuhuiz 发表于 2011-6-10 10:03:29

2# messenger
页: [1]
查看完整版本: 求助解方程组