tDatakietailambdaBi
x0065009801470187024302950344039804480499
y1237232223215205196187178169159
y20154022703270403049705760644070707570795
y326775268452714727549280522838928645287562864528243
figure1pf仁polyfitxy12
px1polyvalpf1x
plotxpx1k
grid
xlabel转速比i
ylabel变矩比K
title二阶多项式拟合k曲线

pause
figure2
pf2polyfitxy23
px2polyvalpf2x
plotxpx2b
grid
xlabel转速比i
ylabel效率eta
title三阶多项式拟合eta曲线

pause
ffigure3pf3polyfitxy34px3polyvalpf3xplotxpx3rgridxlabel转速比iylabel泵轮转矩系数lambdaBtitle四阶多项式拟合lambdaB曲线figure4pf仁polyfitxy12px1polyvalpf1xplotxy1orxpx1kgridxlabel转速比iylabel变矩比Ktitle二阶多项式拟合k曲线Lege
d原始数据拟合曲线各的原始特性参数数据点和二阶拟合曲线绘制在同一张图形中pausefigure5pf2polyfitxy23px2polyvalpf2x
fplotxy2mxpx2bgridxlabel转速比iylabel效率etatitle三阶多项式拟合eta曲线Lege
d原始数据拟合曲线0各的原始特性参数数据点和三阶拟合曲线绘制在同一张图形中pausefigure6pf3polyfitxy34px3polyvalpf3xplotxy3pkxpx3rgridxlabel转速比iylabel泵轮转矩系数lambdaBtitle四阶多项式拟合lambdaB曲线Lege
d原始数据拟合曲线0各的原始特性参数数据点和四阶拟合曲线绘制在同一张图形中
y1poly2strpf1x
Ki曲线的二阶拟合公式
y2poly2strpf2x
y3poly2strpf3x运行后的结果如下
i曲线的三阶拟合公式Bi曲线的四阶拟合公式
运行后的二阶，三阶，四阶拟合曲线函数为y1001325xA218035x2491y2012713人316598xA224499x00025474y31067407xA41999852人3958404人287272x269754
f四、微分方程求解。（25分）自己选择确定一个三阶微分方程，自己设置初始条件，用
ode45方法求微分方程的解。
3
2
要求：例如：屮2嗨4如8yt1y00刨0ki，
dt
dt
dt
dt
d2y0dt20
1）仿真时间t30秒
2）结果绘制在一张图中包括yt曲线，一阶y1曲线，二阶y