3Â¥: Originally posted by wtisunday at 2011-11-24 11:09:43:
»¹Çë°æÖ÷¸ø³öÏêϸµÄ½â´ð°¡£¬Ð»Ð»ÁË

function AA
x=[1 100 1000 10000 100000,20000,900,20];
fsolve(@fun,x)

function ff=fun(x)
w=[0.00100 0.02000 0.1300 0.910000 8.7000 76 1001 4301 20032 100000]';
G1=[1000 2034 3000 4700 6002 10032 50345 109820 620000 1200000]';
G2=[6000 7200 9000 11000 13000 23000 38000 70012 209903 632098]';
for i=5:8
    y(mod(i,5)+1)=x(i);
end
   
for i=1:4
    for j=1:4
        A(i,j)=0,B(i,j)=0,b(i)=0,c(i)=0
        for k=1:10
            A(i,j)=A(i,j)+(w(k)*x(j))^2*(w(k)*x(i))^2/(G1(k))^2/(1+(w(k)*...
                x(j))^2)/(1+(w(k)*x(i))^2);
            B(i,j)=B(i,j)+2*y(i)*w(k)^2*x(i)*(w(k)*x(j))^2/(G1(k)^2)/(1+...
                (w(k)*x(i))^2)^2/(1+(w(k)*x(j))^2);
            b(i)=b(i)+(w(k)*x(i))^2/G1(k)/(1+(w(k)*x(i))^2);
            c(i)=c(i)+2*y(i)*x(i)*(w(k))^2/G1(k)/(1+(w(k)*x(i))^2);
        end
    end
end
b=b';
c=c';
X=x(1:4)';
Y=x(5:8)';
ff1=A*X-b;
ff2=B*Y-c;
ff=[ff1;ff2];