| 查看: 412 | 回复: 0 | |||
[交流]
两种方法解的结果不一样，请哪位大侠解释一下，谢谢！
|
|
第一种解法 syms A B C D E F A1 B1 C1 D1 E1 F1 H0 M2 M3 a b c [A B C D E F A1 B1 C1 D1 E1 F1]=solve('-H0+A/a^2=B+C/a^2', 'D+E/b^2=B+C/b^2', 'D+E/c^2=F', 'A1/a^2=B1+C1/a^2', 'D1+E1/b^2=B1+C1/b^2', 'F1=D1+E1/c^2', 'H0+A/a^2=-B+C/a^2+M2-0.000019*(-B1+C1/a^2)', '-D+E/b^2+M3-0.000019*(-D1+E1/b^2)=-B+C/b^2+M2-0.000019*(-B1+C1/b^2)', '-D+E/c^2+M3-0.000019*(-D1+E1/c^2)=-F', 'A1/a^2=100*(-B1+C1/a^2)-2.75*(-B+C/a^2)', '100*(-D1+E1/b^2)-2.75*(-D+E/b^2)=100*(-B1+C1/b^2)-2.75*(-B+C/b^2)', '100*(-D1+E1/c^2)-2.75*(-D+E/c^2)=-F1', 'A, B, C, D, E, F, A1, B1, C1, D1, E1, F1') a=0.03; b=0.02; c=0.01; H0=10000; M2=0 M3=0 eval(A) eval(B) eval(C) eval(D) eval(E) eval(F) eval(A1) eval(B1) eval(C1) eval(D1) eval(E1) eval(F1) ans = -2.064311914376605e-006 ans = -0.24385437050510 ans = -1.000000695151247e+004 ans = -2.716199541293015e+002 ans = 2.637731770953946e-007 ans = 6.036070327523462e-004 ans = -1.000000695151247e+004 ans = -2.716199541293015e+002 ans = 2.637731770953946e-007 ans = 6.036070327523462e-004 ans = -9.999998475757755e+003 ans = -2.655837992560781e+002 第二种解法 syms M=[1/0.03^2 -1 -1/0.03^2 0 0 0 0 0 0 0 0 0;0 -1 -1/0.02^2 1 1/0.02^2 0 0 0 0 0 0 0;0 0 0 1 1/0.01^2 -1 0 0 0 0 0 0;0 0 0 0 0 0 1/0.03^2 -1 -1/0.03^2 0 0 0;0 0 0 0 0 0 0 1 1/0.02^2 -1 -1/0.02^2 0;0 0 0 0 0 0 0 0 0 1 1/0.01^2 -1;1/0.03^2 1 -1/0.03^2 0 0 0 0 -0.000019 0.000019/0.03^2 0 0 0;0 1 -1/0.02^2 -1 1/0.02^2 0 0 -0.000019 0.000019/0.02^2 0.000019 -0.000019/0.02^2 0;0 0 0 -1 1/0.01^2 1 0 0 0 0.000019 -0.000019/0.01^2 0;0 2.75 -2.75/0.03^2 0 0 0 -1/0.03^2 -100 100/0.03^2 0 0 0;0 2.75 -2.75/0.02^2 -2.75 2.75/0.02^2 0 0 -100 100/0.02^2 100 -100/0.02^2 0;0 0 0 2.75 -2.75/0.01^2 0 0 0 0 -100 100/0.01^2 1] M1=inv(M); T=[10000 0 0 0 0 0 -10000 0 0 0 0 0]' ANS=M1*T ansA=[1 0 0 0 0 0 0 0 0 0 0 0]*ANS ansB=[0 1 0 0 0 0 0 0 0 0 0 0]*ANS ansC=[0 0 1 0 0 0 0 0 0 0 0 0]*ANS ansD=[0 0 0 1 0 0 0 0 0 0 0 0]*ANS ansE=[0 0 0 0 1 0 0 0 0 0 0 0]*ANS ansF=[0 0 0 0 0 1 0 0 0 0 0 0]*ANS ansA1=[0 0 0 0 0 0 1 0 0 0 0 0]*ANS ansB1=[0 0 0 0 0 0 0 1 0 0 0 0]*ANS ansC1=[0 0 0 0 0 0 0 0 1 0 0 0]*ANS ansD1=[0 0 0 0 0 0 0 0 0 1 0 0]*ANS ansE1=[0 0 0 0 0 0 0 0 0 0 1 0]*ANS ansF1=[0 0 0 0 0 0 0 0 0 0 0 1]*ANS ansA =-2.064311333633759e-006 ansB =-1.000000258676058e+004 ansC =2.637731821910219e-007 ansD =-1.000000258676058e+004 ansE =2.637731830247537e-007 ansF =-9.999999949028745e+003 ansA1 =-0.24385431419503 ansB1 =-2.716199124050910e+002 ansC1 =6.036069695545399e-004 ansD1 =-2.716199124050910e+002 ansE1 =6.036069695545584e-004 ansF1 =-2.655838427095454e+002 |
回复此楼» 猜你喜欢
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有3人回复
-
E0414, 我的本子有没有希望？
已经有9人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有3人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有3人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有3人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有4人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有3人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有3人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有4人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有5人回复
找到一些相关的精华帖子，希望有用哦~
请问大侠，SCI投稿过程中的due on 是什么意思？
已经有3人回复
- 求助帮忙做个多元线性回归的分析spss，R均可
已经有4人回复
- 核磁氢谱1.25位置处 出现一个很高很高的峰，哪位大侠解释解释啊！！！！！！
已经有32人回复
- gDNA与基因组DNA的区别
已经有8人回复
- 标题相似论文主要结果不一样算一稿多投吗
已经有10人回复
- 请解释一下spin-coating和dip-coating的区别
已经有10人回复
- 请大侠指教matlab的这个错误提示到底是个神马意思？
已经有12人回复
- 求大侠帮我解释下点蚀为什么易发生在下图的两个区域
已经有8人回复
- 各位大侠，请帮我解释一下这个反应，谢谢
已经有6人回复
- 请大侠帮助分析一下，这个裂纹产生的原因
已经有22人回复
- Materials Studio能带图显示问题
已经有3人回复
- 求助数值计算结果和仿真结果有什么不一样的
已经有6人回复
- 学有机的大侠们帮帮忙，谢谢啦^-^ ^-^ ^-^
已经有6人回复
- 哪位大侠给详细解释一下CCA，CA, PCA分析的区别啊，在此谢过啦
已经有6人回复
- 一道很简单的物理化学题，请大家给我解释一下,谢谢~
已经有14人回复
- 在单晶和多晶衬底上生长的薄膜晶粒尺寸不一样，为什么？
已经有4人回复
- 今天用CASTEP计算声子谱，出现下列提示，有能看懂的大侠帮忙解释下
已经有11人回复
- 半定量内参条带的亮度不一样，各位大侠怎么办？？？急呀！！！！！
已经有22人回复
- 请问两种IV曲线的不同之处在哪里，各有什么优点？求大侠指导
已经有6人回复
- 【求助】二甲硫醚氧化
已经有3人回复
- 【求助】解释一下收率
已经有6人回复
- 【求助】求高人解释一下GPC测出的这些数据各代表的意思！
已经有8人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
