用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');

你的算式有问题,三个方程式
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的项约掉了,且左边≠右边,方程式有问题

温馨提示:内容为网友见解,仅供参考
无其他回答
相似回答