2014-11-29 3 views
0

'fsolve'로 연습하려고하는데 다음 코드로 어떤 일이 일어나고 있는지 완전히 파악하지 못했습니다. 누구든지 그것에 빛을 비춰 줄 수 있습니까? 여기 MATLAB : fsolve 함수

function N=productivity1(N,Ac,Aw) 
global Thetac Thetaw tau a b 
N=[N(1),N(2)]; 
N=[N(1)-(Thetac/a)^(1/b)*(1+tau)*Ac; 
    N(2)-(Thetaw/a)^(1/b)*(1+tau)*Aw; 
    N(1)+N(2)-1]; %this meant to be a constraint... 
end 

N0=[0.7,0.3]; %initial guess for x 
option=optimset('Display','iter'); 
result=fsolve(@(N)productivity1(N,Ac0,Aw0),N0,option); 

가 상세하다 :

Ac=2 
Aw=20 
AC 및 아은 N (1) + N (2) = 1, AC 및 아의 시작 값

고정되지되도록 값을 취해야

고정 파라미터 값

Thetac=0.6 
Thetaw=0.6 
tau=0.07 
a=0.01 
b=0.8 

출력 t의 값에 관계없이 초기 값의 벡터와 동일 그는 위의 매개 변수. 이것이 제가 직면 한 문제입니다.

답변