matlab编译和运行结果不一样,谁能告诉我为什么Matlab两次运行结果不一样？

最新推荐文章于 2024-04-29 13:03:03 发布

转载最新推荐文章于 2024-04-29 13:03:03 发布 · 1.4k 阅读

标签

#matlab编译和运行结果不一样

博主在Matlab中遇到一个问题，即同一段代码在不同时间运行得到的结果不一致。具体涉及复杂运算和符号求解过程。在解决一个高阶方程组后，再次运行之前的简单求解程序时，出现了无法找到显式解的警告。寻求经验丰富的用户帮助理解这一现象。

syms x y

d=x^2 + x*y + y - 3;

c=x^2 - 4*y + 3 ;

[x,y] = solve(d,c)

%运行以上程序，刚开始的时候能够出现结果。

当我从事一次较复杂的运算

eq1 =

23/50-(2-(1-3022314549036573/151115727451828646838272/rp)^2)^2*(1+81602492823987471/75557863725914323419136000/rp-2256193152629581474387615421329263/5708990770823839524233143877797980545530986496000/rp^2+12174680555738555285211392885523367861254587051997/3450873173395281893717377931138512726225554486085193277581262111899648000/rp^3)*(1-3022314549036573/151115727451828646838272/rp)^2/(1-exp(-23*(2-(1-3022314549036573/151115727451828646838272/rp)^2)^2*(1+81602492823987471/75557863725914323419136000/rp-2256193152629581474387615421329263/5708990770823839524233143877797980545530986496000/rp^2+12174680555738555285211392885523367861254587051997/3450873173395281893717377931138512726225554486085193277581262111899648000/rp^3)*D/(1-69513234627841179/1511157274518286468382720/rp+5270540279624568869318437644157833/11417981541647679048466287755595961091061972992000/rp^2+193248897710135798177958617230529648591342651619/107839786668602559178668060348078522694548577690162289924414440996864000/rp^3)/(2050338190066411/38685626227668133590597632-142525639672509649983195570538569/58460065493236116728147393308651320786237301719040/rp+10806390017597554265558859216130742149928435847363/441711766194596082395824375185729628956870974218904739530401550323154944000/rp^2+396225595163328829636228797640143382163006246482900686446669409/4171849679533027504677776769862406473833407270227837441302815640277772901915313574263597826048000/rp^3))*(1-(2-(1-3022314549036573/151115727451828646838272/rp)^2)^2*(1+81602492823987471/75557863725914323419136000/rp-2256193152629581474387615421329263/5708990770823839524233143877797980545530986496000/rp^2+12174680555738555285211392885523367861254587051997/3450873173395281893717377931138512726225554486085193277581262111899648000/rp^3)*(1-3022314549036573/151115727451828646838272/rp)^2))

eq2 =

7/24*rp^2/D-12

[D,rp = solve(eq1,eq2)

之后再运行先前的程序

其运行结果就变成一下这种情况了

Warning: Explicit solution could not be found.

> In solve at 140

In sym.solve at 49

In Untitled at 5

x =

[ empty sym ]

y =

[]

由于本人功力尚浅，还望见多识广的朋友们给与指点，将不甚感谢！

[本帖最后由 qi_matlab 于 2010-1-11 10:51 编辑]