数学建模社区-数学中国
标题:
求各位大神帮忙解非线性方程组
[打印本页]
作者:
失梦人
时间:
2014-4-6 09:32
标题:
求各位大神帮忙解非线性方程组
我要求一个非线性方程组， sqrt((x+a)^2+y^2+z^2)- sqrt(x^2+y^2+z^2)=340*t1
3 B7 V. b/ e* ^3 K5 W
sqrt((x+a)^2+y^2+z^2)-sqrt((x-a)^2+y^2+z^2)=340*t2
" p! W) u5 s( \% t- W. P. W( Y+ R, Z
sqrt(x^2+y^2+z^2)-sqrt((x-a)^2+y^2+z^2)=340*t3
- D, r; K9 v% P
要求输入参数为a，t1，t2，t3，变量。程序返回方程组的解x，y，z。
1 Q! g, x; Q. M& S0 l
function F=myfun(x)
( z/ s* z% U( S4 |
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 Z* {, @3 d/ m
/ l5 b) s. Q& Q* x, y' K% o
x0=[0;0;0];
# N3 [) l( C: i
options=optimset('Display','iter');
5 ]; i' x' h w2 N2 E% w5 d8 H
[x,fv]=fsolve(@myfun,x0,options)
5 ~, O& d Y4 U) _4 X6 X
$ w) v8 I- d6 [4 X: i
x=x(1);y=x(2);z=x(3);
" o h' {" @. Y4 o1 D1 P( \
9 n$ L4 k9 t$ k9 }
假如t1,2,3和a都是已知数了，那么程序可以直接输出x，y，z。但我如何让他们以输入参数的形式出现呢。就比如说，他们是实际测量得到的数据，数据量很大的时候，总不能每次都重新输一遍程序，假如t1,2,3和a是可以matla独处的，那么在哪里赋值，怎样为之赋值才能让他们以方程组自变量的形式出现呢，也就是如何动态的读入并修改程序里的a，t1，t2，t3的值。
3 e" ~/ i. j$ B! p1 g, T4 ?
; A* P/ @6 b/ d! a, r
新手，体力不够了，只有20，刚发了一个悬赏50体力，结果有乱码，帮帮忙把程序编出来，我用七十体力来答谢你，我的全部家当，真的谢谢您
5 n5 j/ y" z5 n
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5