| 查看: 7347 | 回复: 12 | ||
| 本帖产生 1 个 博学EPI ,点击这里进行查看 | ||
[交流]
三元二次方程组求解 怎么算啊 急用
|
||
» 抢金币啦!回帖就可以得到:
江苏师范大学计算化学方向专任教师招聘(事业编)
+3/223
大连理工大学大连市复杂工业场景具身智能重点实验室科研助理招聘启事
+2/178
急招科研助理/博士后
+2/166
招收 哈工大深圳﹣大湾区大学联合培养博士生1名,2027年入学
+1/94
美国肯塔基大学Dr. Sheng Tong博士后招聘(空间可控的基因组编辑与肿瘤免疫治疗)
+1/94
高校教师 科研教学岗
+1/80
【推免】】中科院国家级人才团队招收27级推免硕士生 (材料+机械)
+1/71
哈尔滨工业大学(深圳)吴晓军教授课题组招收2027级“申请-考核”博士研究生
+2/18
来,到深技大读研究生!2027 期待与你 共创不凡!
+1/16
中国科大-合肥国家实验室冷原子量子网络团队招聘启事
+2/14
中国石油大学(华东)生物质绿色化学转化和生物基功能高分子团队诚招27级推免生
+1/10
怀柔实验室与中国矿业大学(北京)氢能与燃料电池课题组招收2027年联培博士生
+1/7
宁波诺丁汉大学招收2027年春季入学CO2电催化方向博士生
+1/4
【生物材料与生物打印专题征稿】Scopus期刊,免APC
+1/4
无人机方向 博士后与青年教师招聘公告
+1/3
北理工集成电路杰青团队 | 诚招科助理
+1/3
浙江大学Said Mikki教授课题组招生2027秋博士生-电子科学与技术
+1/3
浙江大学Kemal Celebi教授课题组招收2027年秋季博士生
+1/3
口腔特刊客座招募客编Biofunctional Materials牙科生物功能材料专题,Scopus期刊
+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










回复此楼