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 ...

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 ...

4楼: Originally posted by 幸福小猪啊 at 2013-08-22 20:38:16
先谢谢你的帮助哈，但是从图中来看，不符合非线性方程的多值性特点，而且我最关心的地方是如何能够解除方程的不稳定解，希望你能在帮忙思考一下，谢啦...