function solve_nonl
clear all;clc
format long
global w
a=0:0.1:10;
n=length(a);
x0=[0.1 0.1 -1];
for i=1:n
    w=a(i);
    xresult=fsolve(@nonl,x0);
    xx(i)=xresult(1);
    yy(i)=xresult(2);
    zz(i)=xresult(3);
end
[a' xx' yy' zz']

plot(a,xx,'k',a,yy,'b',a,zz,'r')

function f=nonl(x)
global w
f(1)=3*x(1)-cos(x(2)*x(3))-1/2;
f(2)=x(1)*x(1)-81*(w*x(2)+0.1)^2+sin(x(2))+1.06;
f(3)=exp(-x(1)*x(2))+20*x(3)+(10*pi-3)/3;

2Â¥: Originally posted by ÔÂÖ»À¶ at 2013-08-22 18:48:16
³õÖµ¿ÉÒÔÔÙµ÷ÊÔһϣ¬½á¹û¼û¸½Í¼£¬½ö¹©²Î¿¼¡£

function solve_nonl
clear all;clc
format long
global w
a=0:0.1:10;
n=length(a);
x0=;
for i=1:n
    w=a(i);
    xresult=fsolve(@nonl,x0);
    xx ...

2Â¥: Originally posted by ÔÂÖ»À¶ at 2013-08-22 18:48:16
³õÖµ¿ÉÒÔÔÙµ÷ÊÔһϣ¬½á¹û¼û¸½Í¼£¬½ö¹©²Î¿¼¡£

function solve_nonl
clear all;clc
format long
global w
a=0:0.1:10;
n=length(a);
x0=;
for i=1:n
    w=a(i);
    xresult=fsolve(@nonl,x0);
    xx ...

4Â¥: Originally posted by ÐÒ¸£Ð¡Öí°¡ at 2013-08-22 20:38:16
ÏÈллÄãµÄ°ïÖú¹þ£¬µ«ÊÇ´ÓͼÖÐÀ´¿´£¬²»·ûºÏ·ÇÏßÐÔ·½³ÌµÄ¶àÖµÐÔÌص㣬¶øÇÒÎÒ×î¹ØÐĵĵط½ÊÇÈçºÎÄܹ»½â³ý·½³ÌµÄ²»Îȶ¨½â£¬Ï£ÍûÄãÄÜÔÚ°ïæ˼¿¼Ò»Ï£¬Ð»À²...