24小时热门版块排行榜    

查看: 7325  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +6 yuleib84 2026-08-26 7/350 2026-08-31 19:46 by 鱼翔浅底1
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[基金申请] 中青基了要发朋友圈吗? +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 国社科又开始会评了,不知道这次命运如何 +6 雨打竹帘 2026-08-30 10/500 2026-08-31 09:29 by huixian257
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 基金未中,这种答复是模板吗? +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +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
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
信息提示
请填处理意见