[交流] 【求助】matlab求解非线性方程组，并画图处理。要求y,z是实数解！

Originally posted by blue.wh at 2011-01-06 20:51:41:

你用的方法有问题啊，得到结果不是正解！！！！

clear
clc

N = 100;
[ a, b ] = meshgrid( 0 : 1 / N : 1, 0 : 0.1 / N : 0.1 );
[ m, n ] = size( a );

y = zeros( m, n );
z = zeros( m, n );
F0 = [ 0.1; 0.1 ];
options=optimset('Display','off');

for i = 1 : m
    for j = 1 : n
        F0 = fsolve( @( x )gzqdyouxia( x, a( i, j ), b( i, j ) ), F0, options );
        y( i, j ) = F0( 1 );
        z( i, j ) = F0( 2 );
    end
end
figure( 'name', '三维图')
subplot( 1, 2, 1 )
surf( a, b, y )
title( 'y-a-b' )

subplot( 1, 2, 2 )
surf( a, b, z )
title( 'z-a-b' )

function F = gzqdyouxia( x, a, b )

F = [ - x( 1 )  + x( 2 ) + ( 7999022433351517 / 2305843009213693952 * x( 2 ) + 5509600814111415 / 72057594037927936 * x( 2 ) ^ 3 ) / ( 1 - a ) + 1 / 2 * b;
    - x( 2 ) + x( 1 ) + ( - 4534249478333089 / 576460752303423488 * x( 1 ) - 3611543935826555 / 576460752303423488 * x( 1 ) ^ 3 + 17523 / 50000 * x( 1 ) ^ 5 ) / a - 1 / 2 * b ];

clear
clc

N = 1000;
a = 0 : 1 / N : 1;
b = 0 : 0.1 / N : 0.1;
y = zeros( 1, N );
z = zeros( 1, N );
F0 = [ 0; 0 ];
options=optimset('Display','off');

for i = 2 : N    % a 不取0， 也可以从0开始
    F0 = fsolve( @( x )gzqdyouxia( x, a( i ), b( i ) ), F0, options );
    y( i ) = F0( 1 );
    z( i ) = F0( 2 );
end

plot3( a( 2: end ), b( 2: end ), y, 'r' )
hold on
plot3( a( 2: end ), b( 2: end ), z, 'b' )
xlabel( 'a' )
ylabel( 'b' )
legend( 'y', 'z' )

clear
clc

N = 100;
[ a, b ] = meshgrid( 0 : 1 / N : 1, 0 : 0.1 / N : 0.1);
[ m, n ] = size( a );

y = zeros( m, n );
z = zeros( m, n );
F0 = [ 0; 0 ];
options=optimset('Display','off');

for i = 1 : m
    for j = 1 : n
        F0 = fsolve( @( x )gzqdyouxia( x, a( i, j ), b( i, j ) ), F0, options );
        y( i, j ) = F0( 1 );
        z( i, j ) = F0( 2 );
    end
end

surf( a, b, y )
hold on
surf( a, b, z )
xlabel( 'a' )
ylabel( 'b' )

Originally posted by xiegangmai at 2010-12-27 19:40:23:
a 和 b 之间有什么关系吗？

如果 a 和 b 一一对应，可以用fsolve求解，大概如下：

函数：
[code]function F = gzqdyouxia( x, a, b )

F = [ - x( 1 )  + x( 2 ) + ( 7999022433351517 / 23058430092136 ...

Originally posted by gzqdyouxia at 2010-12-27 22:18:05:

谢谢你的回答，不过问题我没有解决，a，b都是自变量，y，z是关于a和b的二元函数，我要求的是实数对应关系，您能不能棒我再看一下？