matlab两个问题，急！！在线等

来源：百度知道 编辑：UC知道 时间：2024/06/08 04:38:16
有一组实验数据如附表 1-1 所示。
X 1 2 3 4 5 6 7 8 9 10
Y 16 32 70 142 260 436 682 1010 1432 1960
请分别用拟合(二阶至三阶)和插值(线性和三次样条)的方法来估测 X=9.5 时 Y 的值。以下是实现一阶拟合的语句。 （xy一一对应）

求解当y(0)=2，z(0)=7时，微分方程组的解
dy/dx-z=sinx
dz/dy+y=x+1

第一题：
x=[1 2 3 4 5 6 7 8 9 10]';
y=[16 32 70 142 260 436 682 1010 1432 1960]';

p1=polyfit(x,y,2) %2次拟合
val1=polyval(p1,9.5) %估值

p2=polyfit(x,y,3) %3次拟合
val2=polyval(p1,9.5) %估值

interp1(x,y,9.5) %线性插值
interp1(x,y,9.5,'spline') %三次样条插值
第二题：
求解当时，微分方程组的解

[y,z]=dsolve('Dy-z=sinx,Dz/Dy+y=x+1','y(0)=2,z(0)=7','x')
解为：
y =

(-(-2*RootOf(Int(1/(-2*_f+2*log(_f+sinx)+15-2*log(7+sinx))^(1/2)/(_f+sinx),_f = 7 .. _Z)+x)+2*log(RootOf(Int(1/(-2*_f+2*log(_f+sinx)+15-2*log(7+sinx))^(1/2)/(_f+sinx),_f = 7 .. _Z)+x)+sinx)+15-2*log(7+sinx))^(1/2)*(RootOf(Int(1/(-2*_f+2*log(_f+sinx)+15-2*log(7+sinx))^(1/2)/(_f+sinx),_f = 7 .. _Z)+x)+sinx)-RootOf(Int(1/(-2*_f+2*log(_f+sinx)+15-2*log(7+sinx))^(1/2)/(_f+sinx),_f = 7 .. _Z)+x)-sinx-x*RootOf(Int(1/(-2*_f+2*log(_f+sinx)+15-2*log(7+sinx))^(1/2)/(_f+sinx),_f = 7 .. _Z)+x)-x*sinx)/(-RootOf(Int(1/(-2*_f+2*log(_f+sinx)+1