QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 3042|回复: 0
打印 上一主题 下一主题

[代码资源] 求各位大神帮忙解非线性方程组

[复制链接]
字体大小: 正常 放大
失梦人        

4

主题

11

听众

42

积分

升级  38.95%

  • TA的每日心情
    奋斗
    2015-1-16 20:44
  • 签到天数: 10 天

    [LV.3]偶尔看看II

    自我介绍
    华北电力大学 电气工程自动化 12级七班

    邮箱绑定达人 社区QQ达人

    群组学术交流B

    群组Matlab讨论组

    跳转到指定楼层
    1#
    发表于 2014-4-6 09:32 |只看该作者 |倒序浏览
    |招呼Ta 关注Ta |邮箱已经成功绑定
    20体力
    我要求一个非线性方程组, sqrt((x+a)^2+y^2+z^2)- sqrt(x^2+y^2+z^2)=340*t1; Y- J: v* \) S5 {5 K
                                           sqrt((x+a)^2+y^2+z^2)-sqrt((x-a)^2+y^2+z^2)=340*t2
    / K2 e- U* P) [- X1 x                                        sqrt(x^2+y^2+z^2)-sqrt((x-a)^2+y^2+z^2)=340*t3
    5 K3 v6 E4 t0 s$ p0 a( L" i3 c要求输入参数为a,t1,t2,t3,变量。程序返回方程组的解x,y,z。
    0 Q- o% i% I! T" i6 G5 d* H3 ]% mfunction F=myfun(x)/ O" V6 v6 e6 I. `8 }
    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;];
    3 c- S7 j2 C1 @2 |6 q+ k7 s* z0 P; q- x) B" j+ K' Z! `$ c
    x0=[0;0;0];
    ! ~' ?/ n: A( W& m options=optimset('Display','iter');
    ' l+ p: B( [, Y [x,fv]=fsolve(@myfun,x0,options)
    1 U6 b/ F3 P/ N" h9 q. ^& W
    2 _2 v: Y8 }( i; v1 Bx=x(1);y=x(2);z=x(3);
    & A8 u+ u( \& U2 y4 L4 [1 i& k6 r# O. d! l! a3 ]# w% Y
    假如t1,2,3和a都是已知数了,那么程序可以直接输出x,y,z。但我如何让他们以输入参数的形式出现呢。就比如说,他们是实际测量得到的数据,数据量很大的时候,总不能每次都重新输一遍程序,假如t1,2,3和a是可以matla独处的,那么在哪里赋值,怎样为之赋值才能让他们以方程组自变量的形式出现呢,也就是如何动态的读入并修改程序里的a,t1,t2,t3的值。
    # ?1 g9 f. h1 c3 X; P$ \* d: \+ A! @6 ]0 s& e/ k
    新手,体力不够了,只有20,刚发了一个悬赏50体力,结果有乱码,帮帮忙把程序编出来,我用七十体力来答谢你,我的全部家当,真的谢谢您, U" t# m* ]; Z/ b5 w2 O8 P0 {7 o% w

    zan
    转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2026-4-13 10:35 , Processed in 0.440910 second(s), 54 queries .

    回顶部