数学建模社区-数学中国
标题:
求各位大神帮忙解非线性方程组
[打印本页]
作者:
失梦人
时间:
2014-4-6 09:32
标题:
求各位大神帮忙解非线性方程组
我要求一个非线性方程组， sqrt((x+a)^2+y^2+z^2)- sqrt(x^2+y^2+z^2)=340*t1
: x- ~! R: c$ P9 o- D2 c
sqrt((x+a)^2+y^2+z^2)-sqrt((x-a)^2+y^2+z^2)=340*t2
$ ?- c- Q. V% i. e( t+ W
sqrt(x^2+y^2+z^2)-sqrt((x-a)^2+y^2+z^2)=340*t3
( ?) z7 U, h% {4 |# K" J# h8 k
要求输入参数为a，t1，t2，t3，变量。程序返回方程组的解x，y，z。
\' {+ _ d- w0 p! y3 _# G* z) K
function F=myfun(x)
! x' x0 p {" S. V
F=[sqrt((x(1)+a)^2+x(2)^2+x(3)^2)-sqrt(x(1)^2+x(2)^2+x(3)^2)-340*t1;sqrt((x(1)+a)^2+x(2)^2+x(3)^2)-sqrt((x(1)-a)^2+x(2)^2+x(3)^2)-340*t2;sqrt(x(1)^2+x(2)^2+x(3)^2)-sqrt(x(1)^2+x(2)^2+x(3)^2)-340*t3;];
8 t- b& B4 |! |+ f# g% r% r- i2 Z# h
- U( U( I2 E' `" r" e4 M( E
x0=[0;0;0];
f+ K: F" s q
options=optimset('Display','iter');
3 q4 a1 e7 Y& Z8 ]9 o0 [. b* |+ W
[x,fv]=fsolve(@myfun,x0,options)
( G/ A7 s' Z) ^" H2 N7 x
% n; {0 K9 h7 F$ p R8 c
x=x(1);y=x(2);z=x(3);
9 c6 s! q6 N( C H6 o) d
/ B# ]5 `9 j3 @/ K' _" g
假如t1,2,3和a都是已知数了，那么程序可以直接输出x，y，z。但我如何让他们以输入参数的形式出现呢。就比如说，他们是实际测量得到的数据，数据量很大的时候，总不能每次都重新输一遍程序，假如t1,2,3和a是可以matla独处的，那么在哪里赋值，怎样为之赋值才能让他们以方程组自变量的形式出现呢，也就是如何动态的读入并修改程序里的a，t1，t2，t3的值。
1 g! _5 E$ P% z* h
{( @/ k, S+ o# X/ g. B- w) `: }) K9 ~
新手，体力不够了，只有20，刚发了一个悬赏50体力，结果有乱码，帮帮忙把程序编出来，我用七十体力来答谢你，我的全部家当，真的谢谢您
/ J9 g/ Y$ R! w
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5