数学建模社区-数学中国
标题:
求各位大神帮忙解非线性方程组
[打印本页]
作者:
失梦人
时间:
2014-4-6 09:32
标题:
求各位大神帮忙解非线性方程组
我要求一个非线性方程组， sqrt((x+a)^2+y^2+z^2)- sqrt(x^2+y^2+z^2)=340*t1
$ g1 H7 X$ X% r
sqrt((x+a)^2+y^2+z^2)-sqrt((x-a)^2+y^2+z^2)=340*t2
: O M4 `; F2 k/ Q X
sqrt(x^2+y^2+z^2)-sqrt((x-a)^2+y^2+z^2)=340*t3
6 l5 U" V/ V1 t7 q& B# v% I
要求输入参数为a，t1，t2，t3，变量。程序返回方程组的解x，y，z。
1 y- R% V1 _5 q- f4 d# A
function F=myfun(x)
) P! ~: x9 c/ I5 q
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;];
* L) y) a0 @: }9 c8 M, c7 J+ b
6 E' u, \6 h3 L. d# X
x0=[0;0;0];
( J6 m+ M2 R) \ m% A3 N
options=optimset('Display','iter');
8 q |# ^. X0 f+ W" W( ~
[x,fv]=fsolve(@myfun,x0,options)
, B4 F, T' Q' G) S3 d( c: b0 p
0 e( W/ j) ]' H+ ^* l' b7 L- Y$ e5 N8 s
x=x(1);y=x(2);z=x(3);
9 F( ~- O& Q2 H6 k m. z6 g: x
# D W f; h3 W2 B
假如t1,2,3和a都是已知数了，那么程序可以直接输出x，y，z。但我如何让他们以输入参数的形式出现呢。就比如说，他们是实际测量得到的数据，数据量很大的时候，总不能每次都重新输一遍程序，假如t1,2,3和a是可以matla独处的，那么在哪里赋值，怎样为之赋值才能让他们以方程组自变量的形式出现呢，也就是如何动态的读入并修改程序里的a，t1，t2，t3的值。
) O3 c0 Y& X7 {# [8 l& C
) ?" Y& n9 J( Y3 N- C9 G/ ?5 k
新手，体力不够了，只有20，刚发了一个悬赏50体力，结果有乱码，帮帮忙把程序编出来，我用七十体力来答谢你，我的全部家当，真的谢谢您
4 S6 M) g8 P2 E7 Z0 S
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5