求助：solve函数求解方程

star19831106

已知：
k1=-1.1;k2=8.72;k3=34.272;k4=26.7328; xp>-2
求解方程
k1*ts+1/2*k2*ts^2+1/6*k3*ts^3+1/24*k4*ts^4=xp
我首先将xp离散成若干值
比如
xp=-2:5;
或者
xp=-1:4;
然后编制如下程序进行计算

1. 
   
2. clc;clear;
   
3. xp=-1:1:3;%xp的值大于-2就可以
   
4. np=length(xp);
   
5. %-----------
   
6. k1=-1.1;k2=8.72;k3=34.272;k4=26.7328;
   
7. %------------
   
8. syms ts;
   
9. %------------
   
10. kts=k1*ts+1/2*k2*ts^2+1/6*k3*ts^3+1/24*k4*ts^4;
    
11. Dkts=diff(kts,ts);
    
12. %------------
    
13. for i=1:np
    
14.     fts=Dkts-xp(i);
    
15.     t0=double(eval(solve(fts,ts)))
    
16. end
    

复制代码

但是无论xp的值取为：

xp=-2:5;
还是
xp=-1:4;
或者其他的值
解出来的第一个t0值都是复数值（其实应该是实数值）


请问为什么无论xp如何取值解出来的第一个t0值都会是复数呢？
是程序编错了吗？


谢谢！！

nwcwww

在MATLAB MUPAD里算出来三个实根：



直接在命令行里也一样：
syms t;
f=(8354*t^3)/1875 + (2142*t^2)/125 + (218*t)/25 + 9/10;
s=solve(f,t);
double(s)


ans =

   -0.1407
   -3.2657
   -0.4397

star19831106

十分感谢您的帮助。
我复制了一下您的代码
结果是
ans =
  -0.1407         
  -3.2657 - 0.0000i
  -0.4397 + 0.0000i
虽然3个都是实数解，但是后两个却表现为虚部为零的复数。
用isreal()函数识别不出来。继续请教这是什么原因。谢谢！

qibbxxt

分别判断请用isreal和arrayfun函数

star19831106

继续求助，请高手帮忙分析。
谢谢！！

messenger

应该是计算精度的问题吧，你就用下面这三个解吧
ans =

   -0.1407
   -3.2657
   -0.4397

star19831106

谢谢qibbxxt 斑竹和messenger总斑竹 的帮助。
我是要求一组方程的解
不是单单要求一个方程的解

我遇到的问题是一楼的代码中
第一个方程的解总是复数

取xp=-1:1:3时，第一个方程的解为复数，其余为实数
取xp=-2:1:3时，也是第一个方程的解为复数，其余为实数。

不知道是什么原因。
谢谢！