| 查看: 7079 | 回复: 12 | ||
| 本帖产生 1 个 博学EPI ,点击这里进行查看 | ||
[交流]
三元二次方程组求解 怎么算啊 急用
|
||
» 抢金币啦!回帖就可以得到:
湖南大学材料院陶益杰老师招收2026年联合培养硕士学生一名
+1/182
人间烟火,实则就是追求最简单的快乐
+1/97
香港科技大学(广州)黄加强课题组智能电池方向博士招聘
+1/81
哈工大医康学院材料模拟计算方向人才招聘
+1/79
欢迎报考南京农业大学植物环境适应课题组课题组2026级博士生。
+1/76
好玩的不敢搞,能搞的不挣钱,能挣钱的我不会做
+1/70
山东征女友,坐标济南
+1/61
同济大学脑机智能团队脑机接口方向招生招聘
+1/33
南方科技大学物理系夏秀杨课题组招收博士生/博士后 计算与理论软物质/生物物理方向
+1/30
西交利物浦大学招收26年【全奖】博士生1名(空间智能沉浸式手术导航)
+1/26
深圳大学信息功能电子材料方向“申请-考核制”博士生招生
+2/20
深圳大学柔性电子材料方向“申请-考核制”博士生招生
+2/12
上海工程技术大学张培磊教授团队招收博士生
+1/11
岭南大学(香港)诚招材料方向研究助理/RA/博士后
+1/7
澳科大药诚招2026年秋季药剂学/生物材料硕士研究生
+1/5
双一流天津工业大学电信学院李鸿强教授招收2026年申请审核制博士3人
+1/3
求 Habashi F主编的《handbook of extractive metallurgy》一书
+3/3
哈工大(深圳)物理招收2026年9月入学博士生1个名额
+1/1
求助化学专业科技论文写作的课件及电子版教材
+1/1
广东工业大学-化学工程专业博士生招生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













回复此楼