24小时热门版块排行榜    

查看: 7316  |  回复: 12
本帖产生 1 个 博学EPI ,点击这里进行查看

lewmihu

金虫 (小有名气)


[交流] 三元二次方程组求解 怎么算啊 急用

急着要结果,各位同学有会解得帮个忙啊

» 猜你喜欢

» 抢金币啦!回帖就可以得到:

查看全部散金贴

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

各成一家之言

银虫 (小有名气)


这个要借助计算机来计算
2楼2010-12-28 17:36:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiegangmai

版主 (职业作家)


引用回帖:
Originally posted by lewmihu at 2010-12-28 16:36:33:
急着要结果,各位同学有会解得帮个忙啊

四个方程,三个变量,可能无解。
3楼2010-12-28 17:53:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiegangmai

版主 (职业作家)


lewmihu(金币+10): 2010-12-30 09:47:12
引用回帖:
Originally posted by 各成一家之言 at 2010-12-28 17:36:14:
这个要借助计算机来计算

用MATLAB大概算了一下,不知对不对。
CODE:
>> syms h k l
f1 = 2.5712 * sqrt( ( 4 * ( h ^ 2 + h * k + k ^2 ) / ( 3 * 3.232 ^ 2 ) + l ^ 2 / 5.147 ^ 2 ) ) - 1;
f2 = sqrt( 1.2957 * ( h ^ 2 + k ^ 2 + h * k + 0.2957 * l ^ 2 ) ) * cos( 7.78 * pi / 180 ) + 0.5 * h - 0.2957 * l + k;
f3 = sqrt( 1.2957 * ( h ^ 2 + k ^ 2 + h * k + 0.2957 * l ^ 2 ) ) * cos( 44.22 * pi / 180 ) - 0.5 * h + 0.5 * k - 0.2957 * l;
f4 = sqrt( h ^ 2 + k ^ 2 + h * k + 0.2957 * l ^ 2 ) * cos( 71.778 * pi / 180 ) + 0.5 * k + h;

F1 = solve( f1, f2, f3, h, k, l );
F2 = solve( f1, f3, f4, h, k, l );
F3 = solve( f1, f2, f4, h, k, l );
F4 = solve( f2, f3, f4, h, k, l );
>> F1.h

ans =

                      [0.1661439598 - 0.007282097086  i]
                      [                                ]
                      [0.1661439598 + 0.007282097086  i]
>> F2.h

ans =

                       [0.1655947219 + 0.01296695154  i]
                       [                               ]
                       [0.1655947219 - 0.01296695154  i]
>> F3.h

ans =

                      [0.1652496131 + 0.005521010225  i]
                      [                                ]
                      [0.1652496131 - 0.005521010225  i]
>> F4.h

ans =

                                      0.

4楼2010-12-28 17:59:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wfwpq

铁杆木虫 (著名写手)


lewmihu(金币+30, 博学EPI+1): 2010-12-30 09:46:37
>> syms l h k;[l h k]=solve('cos(44.22*pi/180)=(0.5*h+0.2957*l-0.5*k)/sqrt(1.2597*(h^2+k^2+h*k+0.2957*l^2))','cos(71.778*pi/180)=(-0.5*k-h)/sqrt(h^2+k^2+h*k+0.2957*l^2)','2.5712*sqrt(4*(h^2+h*k+k^2)/(3*3.232^2)+l^2/5.147^2)=1')

l =

0.1128765955661879120003069241437


h =

-0.90654616138783667156059105710971


k =

1.2373847796375348971671727500807

此为Matlab准解析解
5楼2010-12-28 18:19:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wfwpq

铁杆木虫 (著名写手)


引用回帖:
Originally posted by lewmihu at 2010-12-28 16:36:33:
急着要结果,各位同学有会解得帮个忙啊

补充一下啊,四个方程只需要三个即可,
若只用后面的三个方程是无解析解的
>> syms l h k;[l h k]=solve('cos(7.78*pi/180)=(-0.5*h+0.2958*l-k)/sqrt(1.2597*(h^2+k^2+h*k+0.2957*l^2))','cos(44.22*pi/180)=(0.5*h+0.2957*l-0.5*k)/sqrt(1.2597*(h^2+k^2+h*k+0.2957*l^2))','cos(71.778*pi/180)=(-0.5*k-h)/sqrt(h^2+k^2+h*k+0.2957*l^2)')
Warning: Explicit solution could not be found.
> In solve at 98

l =

[ empty sym ]


h =

     []


k =

     []
必须使用第一个方程才可以
6楼2010-12-28 18:26:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

辰风

铁虫 (初入文坛)


引用回帖:
Originally posted by lewmihu at 2010-12-28 16:36:33:
急着要结果,各位同学有会解得帮个忙啊

无解吧
7楼2010-12-28 22:35:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunshineck

铁杆木虫 (著名写手)


引用回帖:
Originally posted by 辰风 at 2010-12-28 22:35:55:

无解吧

确实无解
8楼2010-12-29 14:53:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lewmihu

金虫 (小有名气)


引用回帖:
Originally posted by wfwpq at 2010-12-28 18:19:25:
>> syms l h k;[l h k]=solve('cos(44.22*pi/180)=(0.5*h+0.2957*l-0.5*k)/sqrt(1.2597*(h^2+k^2+h*k+0.2957*l^2))','cos(71.778*pi/180)=(-0.5*k-h)/sqrt(h^2+k^2+h*k+0.2957*l^2)','2.5712*sqrt(4*(h^2+h ...

但是这个解代入角度方程不符合啊
9楼2010-12-29 19:59:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wfwpq

铁杆木虫 (著名写手)


>> syms l h k;[l h k]=solve('cos(7.78*pi/180)=(-0.5*h+0.2958*l-k)/sqrt(1.2597*(h^2+k^2+h*k+0.2957*l^2))','cos(71.778*pi/180)=(-0.5*k-h)/sqrt(h^2+k^2+h*k+0.2957*l^2)','2.5712*sqrt(4*(h^2+h*k+k^2)/(3*3.232^2)+l^2/5.147^2)=1')

l =

0.10004078040232283602903575751671


h =

-0.88087352251719425909330478731521


k =

1.2835490072381379552799134869021

你再试试这个看看,这一组解是根据第1、2、4个方程算出来的
10楼2010-12-30 13:38:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wfwpq

铁杆木虫 (著名写手)


>> syms l h k;[l h k]=solve('cos(7.78*pi/180)=(-0.5*h+0.2958*l-k)/sqrt(1.2597*(h^2+k^2+h*k+0.2957*l^2))','cos(44.22*pi/180)=(0.5*h+0.2957*l-0.5*k)/sqrt(1.2597*(h^2+k^2+h*k+0.2957*l^2))','2.5712*sqrt(4*(h^2+h*k+k^2)/(3*3.232^2)+l^2/5.147^2)=1')
Warning: possibly spurious solutions [solvelib::solve_sys]

l =

   0.2192212129795302796638513711869
  -0.2192212129795302796638513711869
  0.10849071521000389499775457145023
-0.10849071521000389499775457145023


h =

  -1.1081530884262844769818164731015
   1.1081530884262844769818164731015
-0.88657016017767077026585302583889
  0.88657016017767077026585302583889


k =

  0.71670695842580947821001081264815
-0.71670695842580947821001081264815
   1.2785744505905856362065314907789
  -1.2785744505905856362065314907789
这是根据第1、2、3个方程计算出来的解,可能会有假解,你可以将h、k、l分别取对应的一组解试试看
11楼2010-12-30 13:44:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wfwpq

铁杆木虫 (著名写手)


引用回帖:
Originally posted by wfwpq at 2010-12-30 13:44:39:
>> syms l h k;[l h k]=solve('cos(7.78*pi/180)=(-0.5*h+0.2958*l-k)/sqrt(1.2597*(h^2+k^2+h*k+0.2957*l^2))','cos(44.22*pi/180)=(0.5*h+0.2957*l-0.5*k)/sqrt(1.2597*(h^2+k^2+h*k+0.2957*l^2))','2.57 ...

补充一下,对应的一组解就是若h取第一个解,则k、l也取第一个解;以此类推
12楼2010-12-30 13:46:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lewmihu

金虫 (小有名气)


非常感谢啊,小木虫真是藏龙卧虎之地,看来应该多过来学习!
13楼2010-12-31 10:30:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lewmihu 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的? +12 爱看书的可乐 2026-08-26 14/700 2026-08-26 12:49 by coolysnow
[基金申请] 为什么国自然不能直接公布 +3 bjdxyxy 2026-08-26 3/150 2026-08-26 12:48 by 小丑鱼爱读书
[基金申请] 系统查不到 +9 董八千 2026-08-26 9/450 2026-08-26 12:30 by yanmofang
[基金申请] 科研孤儿太难了 +21 我4大白菜 2026-08-20 22/1100 2026-08-26 11:29 by nlgza
[基金申请] 怎么查啊 +4 huang1991js 2026-08-26 4/200 2026-08-26 11:26 by panjy.cn
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了,中了面上 +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +6 process2012 2026-08-23 8/400 2026-08-26 11:03 by zjlp
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +24 木水思豆 2026-08-25 27/1350 2026-08-26 08:26 by yanchen918
[基金申请] 有没有大神帮我看看基金代码 28+4 1234567wang 2026-08-24 10/500 2026-08-25 19:15 by lfy8008
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 19/950 2026-08-25 17:57 by chick875
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +3 苏东坡二世 2026-08-23 3/150 2026-08-25 10:35 by shisan1313
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
请填处理意见