Question

用matlab解方程组，无解析解，求助！

方程没有解析解，怎么解出来啊。
vm1=0.0149; vm2=0.0136; vm3=0.0130;

vt1=0.0044;vt2=0.0103;vt3=0.0162;
b=0.125;

求解y0,a0,t0
程序如下：
[y0,a0,t0]=solve('(vm1+y0)=a0*exp(-b*(vt1+t0))','(vm2+y0)=a0*exp(-b*(vt2+t0))','(vm3+y0)=a0*exp(-b*(vt3+t0))','y0','a0','t0')；

Answer 1

你的算式有问题，三个方程式
vm1+y0=a0*exp(-b*(vt1+t0)),vm2+y0=a0*exp(-b*(vt2+t0)),vm3+y0=a0*exp(-b*(vt3+t0))
把已知代入后，0.0149+y0=a0*exp(-0.125*(0.0044+t0)),0.0136+y0=a0*exp(-0.125*(0.0103+t0)),0.0130+y0=a0*exp(-0.125*(0.0162+t0))
1式减2式，2式减3式得：0.0013=a0*[exp(-0.125*(0.0044+t0))-exp(-0.125*(0.0103+t0))],0.0006=a0*[exp(-0.125*(0.0103+t0))-exp(-0.125*(0.0162+t0))]
两式相除，左边=2.1667，
右边=[exp(-0.125*(0.0044+t0))-exp(-0.125*(0.0103+t0))]/[exp(-0.125*(0.0103+t0))-exp(-0.125*(0.0162+t0))]含t0项约掉
=[exp(-0.125*0.0044)-exp(-0.125*0.0103)]/[exp(-0.125*0.0103)-exp(-0.125*0.0162)]=1.007
左边≠右边

追问

方程没有问题，因为曲线都有了，应该是小数位数引起的误差问题。你看那个方程为什么用matlab解不出来呢

追答

我上面分析了，解的过程中把含t0的项约掉了，且左边≠右边，方程式有问题