| 查看: 353 | 回复: 1 | |||
[交流]
两种方法解的结果不一样,请哪位大侠解释一下,谢谢!
|
|
第一种解法 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 第二种解法 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 |
» 猜你喜欢
最近几年招的学生写论文不引自己组发的文章
已经有5人回复
职称评审没过,求安慰
已经有54人回复
26申博自荐
已经有3人回复
A期刊撤稿
已经有4人回复
» 本主题相关价值贴推荐,对您同样有帮助:
求大神们帮我翻译一下下面的内容,谢谢了!
已经有2人回复
求助帮忙做个多元线性回归的分析spss,R均可
已经有4人回复
核磁氢谱1.25位置处 出现一个很高很高的峰,哪位大侠解释解释啊!!!!!!
已经有32人回复
gDNA与基因组DNA的区别
已经有8人回复
标题相似论文主要结果不一样算一稿多投吗
已经有10人回复
请解释一下spin-coating和dip-coating的区别
已经有10人回复
请大侠指教matlab的这个错误提示到底是个神马意思?
已经有12人回复
求大侠帮我解释下点蚀为什么易发生在下图的两个区域
已经有8人回复
各位大侠,请帮我解释一下这个反应,谢谢
已经有6人回复
请大侠帮助分析一下,这个裂纹产生的原因
已经有22人回复
大侠们,求救!看不懂这信上说的是什么意思……
已经有6人回复
求助数值计算结果和仿真结果有什么不一样的
已经有6人回复
哪位大侠给详细解释一下CCA,CA, PCA分析的区别啊,在此谢过啦
已经有6人回复
一道很简单的物理化学题,请大家给我解释一下,谢谢~
已经有14人回复
在单晶和多晶衬底上生长的薄膜晶粒尺寸不一样,为什么?
已经有4人回复
求助:请高手帮忙解释一下自扩散系数,互扩散系数之间的关系,送金币!
已经有3人回复
今天用CASTEP计算声子谱,出现下列提示,有能看懂的大侠帮忙解释下
已经有11人回复
半定量内参条带的亮度不一样,各位大侠怎么办???急呀!!!!!
已经有22人回复
请问两种IV曲线的不同之处在哪里,各有什么优点?求大侠指导
已经有6人回复
【求助】解释一下收率
已经有6人回复
【求助】求高人解释一下GPC测出的这些数据各代表的意思!
已经有8人回复
|
交流]两种方法解的结果不一样,请哪位大侠解释一下,谢谢! 第一种解法 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 第二种解法 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 |
2楼2015-02-04 16:03:48













回复此楼