版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

gzqdyouxia

金虫 (著名写手)

应助: 1 (幼儿园)
贵宾: 0.393
金币: 928.9
帖子: 1600
在线: 310.6小时
虫号: 802190

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

y=z+(7999022433351517/2305843009213693952*z+5509600814111415/72057594037927936*z^3)/(1-a)+1/2*b；
z=y+(-4534249478333089/576460752303423488*y-3611543935826555/576460752303423488*y^3+17523/50000*y^5)/a-1/2*b
注：y，z是因变量，a，b是自变量，其中a的变化范围是（0,1），b的变化范围是（0,0.1）。要求分别画出以a，b为自变量轴，y，z为因变量轴的三维图（即y-a-b，z-a-b）。小虫还有一个问题就是让b为固定值，如b=0.01，怎么画出y-a，z-a（其中a的变化范围是（0,1））的二维图？

[ Last edited by gzqdyouxia on 2010-12-27 at 19:37 ]

回复此楼

» 猜你喜欢

Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有6人回复
2025冷门绝学什么时候出结果已经有4人回复
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗已经有8人回复
寻求一种能扛住强氧化性腐蚀性的容器密封件已经有5人回复
论文投稿，期刊推荐已经有6人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有4人回复
孩子确诊有中度注意力缺陷已经有14人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有5人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复

» 本主题相关价值贴推荐，对您同样有帮助:

MATLAB用于三元回归参数已经有3人回复
有没有什么办法能控制uds求解的标量在整个求解域上的范围啊已经有6人回复
高斯09w做过渡金属模拟，请教出现的L103错误如何解决，断电中断后这样处理对不对已经有3人回复
matlab中二维图作图已经有8人回复
关于污水处理的问题已经有32人回复
如何使用matlab求解非线性方程组的所有整数解？已经有9人回复
Matlab求解偏微分方程组已经有22人回复
matlab求解非线性方程组已经有16人回复
matlab数值求解边界条件微分方程组已经有7人回复
求matlab求解一个方程组的问题已经有6人回复
求高人指点用matlab求解非线性方程组，解决了追加100金币；已经有11人回复
matlab的fsove 命令求解非线性方程组已经有6人回复
【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题已经有21人回复
【求助】matlab怎么求解偏微分方程组啊，先谢谢了已经有13人回复
【求助】请教虫友关于matlab差分方程的求解和画图程序已经有4人回复
【求助】求教matlab解非线性方程组已经有9人回复

» 抢金币啦！回帖就可以得到:
查看全部散金贴

1楼2010-12-27 16:51:57

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

xiegangmai

版主 (职业作家)

引用回帖:

Originally posted by blue.wh at 2011-01-10 10:43:03:

是啊，我算了一下二维的比如b=0.01，在a等于0附近，y，z分别对应三组初值，即使用这个三组初值，运行程序，画出图形，只有一组是正解，其余两组都不是。

改了一下，计算了一下误差，源码如下：

CODE:

clear

clc

eps = 0.0001;

N = 100;

[ a, b ] = meshgrid( ( 0 + eps ) : ( 1 - 2 * eps ) / N : ( 1 - eps ), ( 0 + eps ) : ( 0.1 - 2 * eps ) / N : ( 0.1 - eps ) );

[ m, n ] = size( a );

y = zeros( m, n );

z = zeros( m, n );

F0 = [ 0; 0 ];

% load nlsdat1

options=optimset('Display','iter', ...

    'Algorithm', 'trust-region-reflective');

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( 2, 2, 1 )

surf( a, b, y )

axis equal

title( 'y-a-b' )

subplot( 2, 2, 2 )

surf( a, b, z )

axis equal

title( 'z-a-b' )

% 将a、b、y、z代入方程1，计算误差并画图

diff1 = - y  + z + ( 7999022433351517 / 2305843009213693952 * z + 5509600814111415 / 72057594037927936 * z .^ 3 ) ./ ( 1 - a ) + 1 / 2 * b;

max1 = max( max( diff1 ) )

subplot( 2, 2, 3 )

surf( a, b, diff1 )

title( 'diff1' )

% 将a、b、y、z代入方程2，计算误差并画图

diff2 = - z + y + ( - 4534249478333089 / 576460752303423488 *y - 3611543935826555 / 576460752303423488 * y .^ 3 + 17523 / 50000 * y .^ 5 ) ./ a - 1 / 2 * b;

max2 = max( max( diff2 ) )

subplot( 2, 2, 4 )

surf( a, b, diff2 )

title( 'diff2' )

图：

最大误差：

CODE:

max1 =

      0.000172157751275541

max2 =

     8.02443857179647e-005

赞一下

回复此楼

17楼2011-01-10 13:02:23

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 19 个回答

xiegangmai

版主 (职业作家)

★ ★ ★
gzqdyouxia(金币+50):期待您进一步指教 2010-12-28 00:21:14
robert2020(金币+3):辛苦了！ 2010-12-29 09:34:34

a 和 b 之间有什么关系吗？

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

函数：

CODE:

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 ];

求解：

CODE:

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' )

要画二维曲线，b 为固定值，a在区间递增，求解方程，画图。

如果a 、b没有对应关系，求解出的应该是三维的曲面，大概如下：

CODE:

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' )

在options选项中，可以设置求解算法，具体可参考帮助中的fsolve，我试了试用options=optimset('Display','off', 'ScaleProblem', 'Jacobian' );结果略有不同。

[ Last edited by xiegangmai on 2010-12-27 at 20:20 ]

赞一下(1人)

回复此楼

2楼2010-12-27 19:40:23

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

gzqdyouxia

金虫 (著名写手)

应助: 1 (幼儿园)
贵宾: 0.393
金币: 928.9
帖子: 1600
在线: 310.6小时
虫号: 802190

引用回帖:

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

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

赞一下

回复此楼

3楼2010-12-27 22:18:05

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖