QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 3098|回复: 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*t16 u9 a# T0 V: h# r8 b% r
                                           sqrt((x+a)^2+y^2+z^2)-sqrt((x-a)^2+y^2+z^2)=340*t24 c' B( g- @6 P8 P! a) a" ~: [
                                            sqrt(x^2+y^2+z^2)-sqrt((x-a)^2+y^2+z^2)=340*t3
    ' v0 z* F; h  H6 X( f/ D+ h! E* Q要求输入参数为a,t1,t2,t3,变量。程序返回方程组的解x,y,z。& L" k9 T% B7 N( t6 j2 A) E
    function F=myfun(x)( y, q  I: g! R7 r+ t$ S" |# s
    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;];
    / u1 O' Y1 Y* m4 q4 e+ b$ ~8 m* l5 y8 E1 v$ K
    x0=[0;0;0]; . [3 f: L' G" B2 ~
    options=optimset('Display','iter'); - I8 S* ?, O; i* w1 }/ v7 f
    [x,fv]=fsolve(@myfun,x0,options)
    ( V' D3 M( z+ n% O. j  p  I  }7 j
    x=x(1);y=x(2);z=x(3);
    3 Z# R3 s0 P+ o1 n$ ]1 p+ j5 o" P6 D+ a# N; Q7 A5 e3 `& Y
    假如t1,2,3和a都是已知数了,那么程序可以直接输出x,y,z。但我如何让他们以输入参数的形式出现呢。就比如说,他们是实际测量得到的数据,数据量很大的时候,总不能每次都重新输一遍程序,假如t1,2,3和a是可以matla独处的,那么在哪里赋值,怎样为之赋值才能让他们以方程组自变量的形式出现呢,也就是如何动态的读入并修改程序里的a,t1,t2,t3的值。' M7 U1 D/ X, z5 o
    9 n6 k7 c/ H' M" `2 t+ K
    新手,体力不够了,只有20,刚发了一个悬赏50体力,结果有乱码,帮帮忙把程序编出来,我用七十体力来答谢你,我的全部家当,真的谢谢您
    6 H2 a! j- ?" L, C- P* S1 O

    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-7-10 15:51 , Processed in 0.287046 second(s), 54 queries .

    回顶部