| 查看: 7327 | 回复: 12 | ||
| 本帖产生 1 个 博学EPI ,点击这里进行查看 | ||
[交流]
三元二次方程组求解 怎么算啊 急用
|
||
» 抢金币啦!回帖就可以得到:
开心就散金
+5/1635
天平合规性问题
+1/99
原子层沉积(ALD)磁控溅射PECVD等微纳代工服务:18817872921
+1/86
东南大学能源与环境学院-林日琛教授课题组-27级博士招生
+1/82
东南大学能源与环境学院-林日琛教授课题组-27级博士招生
+1/82
东南大学能源与环境学院-林日琛教授课题组-27级博士招生
+1/80
同一份本子,四份夸技术,一份骂意义——国自然评审的坐标系差异
+1/34
中南大学化学化工学院喻桂朋教授团队诚聘博士后(复合材料方向)
+1/33
深圳大学刘翼振课题组诚聘博士后(CRISPR分子诊断)
+1/31
中南大学化学化工学院喻桂朋教授团队诚聘博士后(复合材料方向)
+1/19
中南大学化学化工学院喻桂朋教授团队诚聘博士后(复合材料方向)
+1/13
四川大学高分子学院、空天学院招收研究助理
+1/12
【科研诚信】武汉轻工大学何江玲3篇论文:涉嫌抄袭、成组数据重复、制造虚假显著
+1/11
湖南大学张旺课题组招收推免硕士
+1/9
湖南大学张旺课题组招收研究助理
+1/5
中山大学智能传感仪器方向招2027届考核制博士生
+1/4
西安交通大学储能电催化团队诚招硕士/博士研究生
+1/3
中山大学具身智能方向招2027届考核制博士生
+1/3
多重蛋白检测定量技术赋能心血管疾病基础与转化研究
+1/3
香港大学智能制造实验室招收27年春秋季入学全奖博士和硕士研究生
+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










回复此楼