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

登录注册

应《网络安全法》要求，自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 ]

回复此楼

» 猜你喜欢

售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急已经有5人回复
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急已经有4人回复
上海工程技术大学激光智能制造课题组 2027级博士研究生招生已经有6人回复
课题组招2027级博士上海工程技术大学激光智能制造方向已经有6人回复
现代”学阀”该如何界定已经有14人回复
求合成方法已经有8人回复
师弟论文见刊大半年才想起来申请专利，还能抢救一下吗？已经有4人回复
有机合成以后会不会被AI改变？做科研的虫友怎么看已经有8人回复
我的奶奶已经有3人回复

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

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的回帖

abingchem

木虫 (著名写手)

应助: 2 (幼儿园)
金币: 4871
帖子: 1672
在线: 117.6小时
虫号: 206219

引用回帖:

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

为什么结果都是零？

赞一下

回复此楼

10楼2010-12-30 14:12:58

已阅回复此楼关注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的回帖