24小时热门版块排行榜    

查看: 338  |  回复: 3

xlk75

铁杆木虫 (正式写手)

[求助] 求教什么方法可以求解下面这个方程组??

求解的方程原型是:
x是未知数  6个未知数  6个方程。
f1=2468.3+137200*tan(x1)/3.8756-x5;
f2=471.6650+294758822624*sin(x2)/3.8756-x5;
f3=170*sin(x2)+160-170*sin(x1)-x6;
f4=x6/93.8-x4/((exp(x4)-exp(-x4))/2);
f5=tan(x3)-(exp(x4)-exp(-x4))/2;
f6=15959.03-x5*tan(x3);
代码详见附件,但是出现如下错误
??? Function 'norm' is not defined for values of class 'sym'.

Error in ==> norm at 28
  [varargout{1:nargout}] = builtin('norm', varargin{:});

Error in ==> newton at 11
if norm(x-x0)[ 来自科研家族 流体数值模拟 ]
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : fun.m
  • 2013-03-04 00:07:23, 357 bytes
  • 附件 2 : dfun.m
  • 2013-03-04 00:07:27, 174 bytes
  • 附件 3 : newton.m
  • 2013-03-04 00:07:31, 788 bytes

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

风神时代
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
xlk75: 金币+20, ★★★★★最佳答案, 希望可以给我破解的版本 2013-03-04 16:37:54
用1stOpt吧,其它软件或算法光初值就是个大问题,而1stOpt一切都免了:
CODE:
Function  2468.3+137200*tan(x1)/3.8756-x5;
          471.6650+294758822624*sin(x2)/3.8756-x5;
          170*sin(x2)+160-170*sin(x1)-x6;
          x6/93.8-x4/((exp(x4)-exp(-x4))/2);
          tan(x3)-(exp(x4)-exp(-x4))/2;
          15959.03-x5*tan(x3);

x1: 0.452183147365698
x5: 19664.3382606743
x2: 2.52352495599274E-7
x6: 85.7218910978855
x4: 0.741679132159666
x3: 0.681757429905715
2楼2013-03-04 08:29:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

方程中有三角函数,因此有多组解,根据实际情况加上约束,下面是有负数的几组结果:

x1: -3.61563437861825
x5: -15694.5569502197
x2: -2.12559574069796E-7
x6: 82.3973353024749
x4: -0.893239170259802
x3: -0.793753200574605

x1: -3.61567173531828
x5: -15696.2275580468
x2: -2.12581539850623E-7
x6: 82.3916857275002
x4: -0.893172609260916
x3: 52.6133751791812
3楼2013-03-04 08:32:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xlk75

铁杆木虫 (正式写手)

问一下这个1stOpt到哪里去下载?

直接输入你的代码就可以得到解了吗?
风神时代
4楼2013-03-04 14:31:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xlk75 的主题更新
信息提示
请填处理意见