QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 2869|回复: 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, I0 E7 D+ g* B" Q3 N
                                           sqrt((x+a)^2+y^2+z^2)-sqrt((x-a)^2+y^2+z^2)=340*t2
    5 }( s/ p3 U/ E                                        sqrt(x^2+y^2+z^2)-sqrt((x-a)^2+y^2+z^2)=340*t3
    / [8 f7 e  G! N' O# V要求输入参数为a,t1,t2,t3,变量。程序返回方程组的解x,y,z。9 z* `- E1 Q; h4 e! A8 ^: L( G, S
    function F=myfun(x)/ Z  B! o" j) y7 a
    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;];
    " g/ R) q9 F' _3 G) d5 j  `' }+ D# \+ O& a. V. ?6 E
    x0=[0;0;0]; ( C: o, P* ]% s0 |) X0 E
    options=optimset('Display','iter'); - c5 e, U4 C2 c$ r! S9 D& o: Y% a
    [x,fv]=fsolve(@myfun,x0,options) 8 H& g: f- J0 z( `! s6 P$ d* W
      x7 l5 {. l. r  F, {: E
    x=x(1);y=x(2);z=x(3);$ N$ J" C4 f# T9 o$ Q2 o, P1 `* R; O
    1 W" q0 R  i8 K- d. Z2 D
    假如t1,2,3和a都是已知数了,那么程序可以直接输出x,y,z。但我如何让他们以输入参数的形式出现呢。就比如说,他们是实际测量得到的数据,数据量很大的时候,总不能每次都重新输一遍程序,假如t1,2,3和a是可以matla独处的,那么在哪里赋值,怎样为之赋值才能让他们以方程组自变量的形式出现呢,也就是如何动态的读入并修改程序里的a,t1,t2,t3的值。
    - }! l' P6 o# H: X$ e
    9 @! r3 k7 s) W8 D; j' ]& D5 I- F! ?新手,体力不够了,只有20,刚发了一个悬赏50体力,结果有乱码,帮帮忙把程序编出来,我用七十体力来答谢你,我的全部家当,真的谢谢您, y2 K* |$ a+ V9 _# G6 G

    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, 2025-10-30 15:34 , Processed in 2.172622 second(s), 54 queries .

    回顶部