| 查看: 7109 | 回复: 12 | ||
| 本帖产生 1 个 博学EPI ,点击这里进行查看 | ||
[交流]
三元二次方程组求解 怎么算啊 急用
|
||
» 抢金币啦!回帖就可以得到:
Call for paper: Frontiers in Plant Science 中科院2区
+2/100
【技术服务】研究所横向技术开发及服务
+1/99
加拿大阿尔伯塔大学招收电磁、无线通信、机器学习方向全奖硕士/博士/博士后/访问学者
+1/56
【武汉高校】【新能源材料与器件课题组】 招收优秀硕士调剂生!!!
+5/45
上海交通大学电催化方向博后招聘-有锂电、高分子、燃料电池背景者优先
+1/35
结节穿刺等待中,求良性
+1/34
南京林业大学-国家级青年人才团队 招2026级申请考核制博士(合成化学、植物化学方向)
+1/19
国家级领军人才团队高分子材料方向2026年博士研究生招生
+1/16
上海交通大学药学院高杨课题组诚聘博士后
+1/14
江苏大学环境学院课题组长期接收学硕专硕
+1/11
加拿大阿尔伯塔大学招收电磁、无线通信、机器学习方向全奖硕士/博士/博士后/访问学者
+2/10
福州大学国家级人才汤育欣教授课题组2026年招收”申请考核制“储能电池方向博士研究生
+1/9
新加坡国立大学刘小钢课题组招聘博士后
+1/9
英国埃克塞特大学 & 法国巴黎萨克雷大学联合培养博士
+1/8
课题组招收环境及相关专业调剂硕士研究生(欢迎优秀学生加入)
+1/7
重庆大学药学院闫海龙课题组拟招收2026年申请考核制博士研究生
+1/5
请问为什么我的贴子自动到“已删除的主题”里而没有被发布?
+1/5
连发两篇 Nature!大样本单细胞测序解析认知、衰老与灵长类进化
+1/3
2026年浙江清华长三角研究院智慧能源中心招聘博士后-电气方向
+1/1
招收环工、环科、化工、应化、化学、生工、生物、材化等相关专业研究生,要求过国家线
+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













回复此楼