| 查看: 260 | 回复: 1 | ||
[求助]
MATLAB求解多元多次方程 求得的解竟然不能是原方程为0
|
|
程序如下： NE = 9; Delta00 = 0.5; Deltat = 0.2; steps = 101; Classicalcorrel00=zeros(1,steps); Timelist00=zeros(1,steps); for i=1:steps Currt=(i-1)*Deltat; Timelist00(1,i)=Currt; p=(2/3+(1-NE*(Delta00^2)*(Currt^2)/4)*exp(-NE*(Delta00^2)*(Currt^2)/8)/3); xx3=-2/3+4*p/3; T003=1/3-2*p/3; T006=-2/3+4*p/3; T007=-2/3; T306=1/3; T315=1-2*p; syms a3 b3 c3 T11 T33 real Eqn100 = a3 - (c3*T007 - 2*T006*T11 + b3*T306 + T315*T33 + xx3)/(1 + 2*b3^2 + 2*(T11^2) + T33^2); Eqn101 = b3 - (T006 + a3*T306)/(1 + a3^2); Eqn102 = c3 - (T003 + a3*T007)/(1 + a3^2); Eqn103 = T11 - ( T007 - a3*T006)/(1 + a3^2); Eqn104 = T33 - ( T007 + a3*T315)/(1 + a3^2); % [a3 b3 c3 T11 T33]=solve(Eqn100,Eqn101,Eqn102,Eqn103,Eqn104,a3,b3,c3,T11, T33); [a3 b3 c3 T11 T33]=solve(Eqn100,Eqn101,Eqn102,Eqn103,Eqn104); Classicalcorrel00(1,i)=((1 + a3^2)*(1 + b3^2 + c3^2 + 2*T11^2 + T33^2) + (26/9)*(1 - 2*p)^2 + ... 8/9 - 2*(a3*xx3 + (2*T11 + T33 + a3*c3)*T007 + b3*a3*T306 + ... c3*T003 + (b3 - 2*a3*T11)*T006 + a3*T33*T315))/8; end MATLAB虽然可以很快得出结果，但是把解代入原方程，Eqn100，Eqn101，Eqn102，Eqn103，Eqn104不为0 比如，以i=1为例 MATLAB得到的解为： K>> [a3 b3 c3 T11 T33] ans = [ -0.69723179291582163445262963600762, -0.86312308374866176586915777567782, 0.90787603708517215995857949120792, 0.53134050208298150303610149633743, -0.51450718729075094864284295767401] 代入Eqn100 K>> Eqn100 Eqn100 = a3 + ((4*T11)/3 + T33 - b3/3 + (2*c3)/3 - 2/3)/(2*T11^2 + T33^2 + 2*b3^2 + 1) K>> a3 + ((4*T11)/3 + T33 - b3/3 + (2*c3)/3 - 2/3)/(2*T11^2 + T33^2 + 2*b3^2 + 1) ans = -0.57062821968242266185604900221681 |
回复此楼» 猜你喜欢
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
听说今天filecode变了
已经有34人回复
-
基金中了
已经有6人回复
-
filecode
已经有14人回复
-
固定端突然变了，今天
已经有7人回复
-
filecode
已经有11人回复
-
大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑
已经有6人回复
-
一个有机合成实验室都需要哪些设备？
已经有12人回复
-
求各位大神看下
已经有33人回复
-
咨询面上基金
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- MATLAB求解非线性方程组
已经有5人回复
- 求Matlab解非线性方程代码
已经有12人回复
- 求助一个关于超越方程用matlab的解法！！
已经有3人回复
- 求高手解答Matlab求解微分方程的初始值和边界值问题
已经有5人回复
- 求Matlab高手解决线性方程组的迭代求解问题
已经有17人回复
- 如何用matlab求解矩阵系数的二阶微分方程
已经有13人回复
- matlab求解非线性方程组，错误提示怎么解决
已经有5人回复
- matlab求解偏微分方程组
已经有3人回复
- 各位朋友，如何用matlab编写这个方程，并求解？
已经有15人回复
- 求助 MATLAB解方程组-fslove
已经有7人回复
- 求助matlab---fsolve解非线性方程组
已经有6人回复
- 1stopt或matlab如何求解以下的非线性方程并拟合出相应曲线？
已经有13人回复
- 高手进来用matlab求解常微分方程，十分感谢！
已经有17人回复
- matlab微分方程组参数拟合，以周为单位求解，汇总后以年为单位进行数值比较
已经有4人回复
- Matlab求解偏微分方程组
已经有22人回复
- matlab求解非线性方程组
已经有16人回复
- matlab数值求解边界条件微分方程组
已经有7人回复
- 求matlab求解一个方程组的问题
已经有6人回复
- 求高人指点用matlab求解非线性方程组，解决了追加100金币；
已经有11人回复
- MATLAB解方程出现??? Subscripted assignment dimension mismatch.
已经有3人回复
- 【求助】matlab 求解微分方程中的未知参数
已经有20人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
- 【求助】matlab怎么求解偏微分方程组啊，先谢谢了
已经有13人回复
- 【求助】利用matlab求解微分方程得出解函数中含0的意思
已经有3人回复
- 【求助】求教matlab解非线性方程组
已经有9人回复
2楼2013-04-19 20:40:51
50