| 查看: 7075 | 回复: 12 | ||
| 本帖产生 1 个 博学EPI ,点击这里进行查看 | ||
[交流]
三元二次方程组求解 怎么算啊 急用
|
||
» 抢金币啦!回帖就可以得到:
中国石油大学(华东)吴传德教授团队(国家杰青)2026硕、博招生
+2/222
中国石油大学(华东)吴传德教授团队(国家杰青)2026硕、博招生
+2/204
西南科技大学曹克课题组招收2026级申请考核制有机化学博士研究生
+1/176
南方科技大学物理系夏秀杨课题组招收博士生/博士后 计算与理论软物质/生物物理方向
+1/30
上海交通大学大气环境科学课题组招收2026年入学博士生
+1/30
上海交大药学院侯四化课题组招收2名2026年秋季入学申请-考核制博士生
+1/29
香港中文大学(深圳)陈筱萌 课题组招生公告(博士 / 博后 / 硕士 / RA)
+1/21
深圳大学柔性电子材料方向“申请-考核制”博士生招生
+2/12
中南林业科技大学生物质绿色转化与功能材料课题组2026年博士招生
+1/12
上海工程技术大学张培磊教授团队招收博士生
+1/12
中科院深圳先进院-免疫治疗方向-招收1名博士生(26年9月入学)
+1/9
意大利华人老师University of Padova-全额奖学金博士
+1/9
都放假了嘛?
+1/9
双一流天津工业大学电信学院李鸿强教授招收2026年申请审核制博士3人
+1/8
化学行业,研发出创新的东西是做成项目给公司吃提成,还是自己搞小作坊倒卖?
+1/7
【科研助理招聘-北京理工大学-集成电路与电子学院-国家杰青团队】
+1/6
【经验分享】CRISPR基因敲除细胞系构建全流程踩坑指南——从递送方式选择到克隆筛选
+1/4
海南省医学科学院沈俊辰课题组生信分析技术员招聘启事
+1/4
2026年博士申请考核+福州大学+管理科学与工程
+1/2
求助化学专业科技论文写作的课件及电子版教材
+1/1
2楼2010-12-28 17:36:14
3楼2010-12-28 17:53:52
4楼2010-12-28 17:59:43
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
|
补充一下啊,四个方程只需要三个即可, 若只用后面的三个方程是无解析解的 >> 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
7楼2010-12-28 22:35:55
8楼2010-12-29 14:53:05
9楼2010-12-29 19:59:57
|
>> 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
|
>> 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
12楼2010-12-30 13:46:07
13楼2010-12-31 10:30:17













回复此楼