24小时热门版块排行榜    

查看: 450  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 QUjhNVAcOSff 2026-09-13 5/250 2026-09-16 14:39 by i3KBRD19ERkG
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 QUjhNVAcOSff 2026-09-13 6/300 2026-09-16 14:27 by i3KBRD19ERkG
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 vZfe6xYu34yj 2026-09-14 3/150 2026-09-16 08:51 by SjQobnC04j06
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 Aj1rhIDL5ixY 2026-09-14 4/200 2026-09-16 07:03 by LH5NkK5Ud4bh
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 40ms4Wlo4umh 2026-09-14 4/200 2026-09-16 06:27 by LH5NkK5Ud4bh
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-16 05:35 by LH5NkK5Ud4bh
[博后之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 4/200 2026-09-16 04:11 by DgNGHc3h5tPl
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-13 6/300 2026-09-16 03:50 by DgNGHc3h5tPl
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-13 4/200 2026-09-16 03:38 by DgNGHc3h5tPl
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 3n8v2C8RimXI 2026-09-13 5/250 2026-09-16 02:50 by DgNGHc3h5tPl
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 QUjhNVAcOSff 2026-09-13 5/250 2026-09-16 02:14 by DgNGHc3h5tPl
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 Aj1rhIDL5ixY 2026-09-14 3/150 2026-09-15 20:14 by zNcEhTcH7ihq
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:39 by CgyNCDVNhGVg
[找工作] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 15:02 by CgyNCDVNhGVg
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +6 QUjhNVAcOSff 2026-09-13 7/350 2026-09-15 14:22 by VMDqgnlKimtZ
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:50 by 5BDX0d0WFp7t
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:38 by 5BDX0d0WFp7t
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-14 16:13 by Aj1rhIDL5ixY
[基金申请] 求教各位大神:2026教育部人文社科青年基金项目何时公示呀? +5 云宴山人 2026-09-10 7/350 2026-09-14 12:11 by ecnu2013
[考博] 27秋季材料博士申请求助 10+3 Cahal1 2026-09-11 4/200 2026-09-14 08:50 by 北京莱茵编辑
信息提示
请填处理意见