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文章，我:8.O.55.1.O54,科目全,可伽急已经有3人回复
浙江师范大学是怎么坑我的已经有7人回复
求教各位大神：2026教育部人文社科青年基金项目何时公示呀？已经有6人回复
面上有专家说实验设备不是我们单位的已经有8人回复
科研求助已经有4人回复
博士去二本高校当辅导员，如何调整心态？已经有11人回复
江苏省自然基金什么时候出结果已经有6人回复
闲聊已经有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的回帖

blue.wh

铁杆木虫 (知名作家)

应助: 972 (博后)
金币: 8637.6
帖子: 7822
在线: 3306.3小时
虫号: 551495

引用回帖:

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不从0开始算起，比如：a = 1.0000e-010 : 1 / N : 1;会出现这个结果：
??? Error using ==> plot
Vectors must be the same lengths.

Error in ==> Untitled at 19
plot( a( 2: end ), y, 'r' )
谢谢~

[ Last edited by blue.wh on 2010-12-29 at 11:24 ]

赞一下

回复此楼

9楼2010-12-29 11:20:35

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