24小时热门版块排行榜

返回列表

zhufu777

新虫 (初入文坛)

应助: 0 (幼儿园)
帖子: 9
在线: 3.3小时
虫号: 1549388
注册: 2011-12-23
专业: 金属结构材料

[求助] matlab求解非线性方程组

先建立M文件L501如下
function ms=L501(X)
ms(1)=400-106670*(X(2)^2+0.33*X(1)^2*X(3)^2)^0.5;
ms(2)=17542*X(1)^2-13740*X(2)*((X(2)^2+5.3*X(1)^2)^0.57-(X(2)^2+0.33*X(1)^2*X(3)^2)^0.57);
ms(3)=35125*X(1)^3-4971.5*((X(2)^2+5.3*X(1)^2)^1.57-(X(2)^2+0.33*X(1)^2*X(3)^2)^1.57)+13740*X(2)^2*((X(2)^2+5.3*X(1)^2)^0.57-(X(2)^2+0.33*X(1)^2*X(3)^2)^0.57);

然后是运行程序
fsolve('L501',[1 1 1])
ans =
0.0140 0.0037 0.0347

fsolve('L501',[0.1 0.1 0.1])

ans =
0.0122 0.0037 0.0009
不同初值，结果相差很大，怎么办啊，到底应该取什么初值啊？？
谁能帮我求一下解。谢谢

只有6个金币。给5个吧

回复此楼

» 猜你喜欢

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

matlab求解非线性方程组已经有3人回复
MATLAB求解五元非线性方程组，在线等已经有7人回复
非线性二阶微分方程组求解，matlab 已经有12人回复
求助解复杂非线性方程组的好的方法已经有24人回复
matlab非线性方程组该如何编程。求助已经有10人回复
用matlab解非线性方程组已经有4人回复
用matlab求解一个非线性方程组的解已经有4人回复
matlab 求解非线性方程组未果已经有10人回复
用matlab求解非线性方程组说无解，一定是方程组本身无解，还是有可能程序有问题呢？已经有11人回复
matlab求解高阶方程组已经有3人回复
求Matlab解非线性方程代码已经有12人回复
关于matlab求解非线性指数方程组出现问题已经有3人回复
matlab解非线性方程组解不开已经有3人回复
求助matlab 非线性方程组全部解！已经有6人回复
matlab求解非线性方程组，错误提示怎么解决已经有5人回复
求助 MATLAB解方程组-fslove 已经有7人回复
matlab求解非线性方程组，求助！已经有6人回复
如何使用matlab求解非线性方程组的所有整数解？已经有9人回复
matlab求解非线性方程组已经有16人回复
求高人指点用matlab求解非线性方程组，解决了追加100金币；已经有11人回复
matlab的fsove 命令求解非线性方程组已经有6人回复
【求助】matlab求解非线性方程组，并画图处理。要求y,z是实数解！已经有18人回复
【求助】求教matlab解非线性方程组已经有9人回复

qq839767816

1楼 2012-10-29 22:06:21

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

zsy312

银虫 (小有名气)

数学EPI: 1
应助: 8 (幼儿园)
金币: 495.5
红花: 1
帖子: 145
在线: 20.7小时
虫号: 756330
注册: 2009-04-24
性别: GG
专业: 数学

【答案】应助回帖

你算的结果不对，三个函数分别是5.32043, 1.82239, -0.00181376，明显偏大了。
我用mathematica算的结果差异比较大的有四组解

CODE:

{x -> -0.000708432, y -> -0.00289383,  z -> 5.86005}, 

{x -> -6.3323*10^-6, y -> 0.00374988,  z -> 1.10622}, 

{x -> -6.3323*10^-6, y -> 0.00374988, z -> -1.10622},

{x -> -0.000708432, y -> -0.00289383, z -> -5.86005}

代码(初始值随机生成，大量搜寻)：

CODE:

n = 300; err = 10^-10;

s = Table[{}, {i, n}];

For[i = 1, i <= n, i++, rand = RandomReal[{-10, 10}, 3];

  (s[[i]] = 

     FindRoot[{400 - 106670 (0.33 x^2 z^2 + y^2)^0.5 == 0, 

       17542 x^2 - 

         13740 y ((5.3 x^2 + y^2)^0.57 - (0.33 x^2 z^2 + y^2)^0.57) ==

         0, 35125 x^3 + 

         13740 y^2 ((5.3 x^2 + y^2)^0.57 - (0.33 x^2 z^2 + 

              y^2)^0.57) - 

         4971.5 ((5.3 x^2 + y^2)^1.57 - (0.33 x^2 z^2 + y^2)^1.57) == 

        0}, Transpose[{{x, y, z}, rand}]]) // Quiet;];

s1 = {};

For[i = 1, i <= n, i++, 

  If[Norm[{400 - 106670 (0.33 x^2 z^2 + y^2)^0.5, 

       17542 x^2 - 

        13740 y ((5.3 x^2 + y^2)^0.57 - (0.33 x^2 z^2 + y^2)^0.57), 

       35125 x^3 + 

        13740 y^2 ((5.3 x^2 + y^2)^0.57 - (0.33 x^2 z^2 + 

             y^2)^0.57) - 

        4971.5 ((5.3 x^2 + y^2)^1.57 - (0.33 x^2 z^2 + y^2)^1.57)} /. 

      s[[i]]] <= err, s1 = PrependTo[s1, s[[i]]]]];

{Length[s1], s1}

Table[{400 - 106670 (0.33 x^2 z^2 + y^2)^0.5, 

   17542 x^2 - 

    13740 y ((5.3 x^2 + y^2)^0.57 - (0.33 x^2 z^2 + y^2)^0.57), 

   35125 x^3 + 

    13740 y^2 ((5.3 x^2 + y^2)^0.57 - (0.33 x^2 z^2 + y^2)^0.57) - 

    4971.5 ((5.3 x^2 + y^2)^1.57 - (0.33 x^2 z^2 + y^2)^1.57)} /. 

  s1[[i]], {i, Length[s1]}]