二次,三次多项式拟合

y=[4, 6.4, 8, 8.4, 9.28, 9.5, 9.7, 9.86, 10, 10.2, 10.32, 10.42, 10.5, 10.55, 10.58, 10.6];

答案： (填写程序语句) 二次多项式拟合 x=1:1:16;

y=[4, 6.4, 8, 8.4, 9.28, 9.5, 9.7, 9.86, 10, 10.2, 10.32, 10.42, 10.5, 10.55, 10.58, 10.6];

a=polyfit(x,y,2) a =

-0.0445 1.0711 4.3252

ezplot('-0.0445*x^2+1.0711*x+4.3252') 三次多项式拟合 x=1:1:16;

y=[4, 6.4, 8, 8.4, 9.28, 9.5, 9.7, 9.86, 10, 10.2, 10.32, 10.42, 10.5, 10.55, 10.58, 10.6];

a=polyfit(x,y,3) a =

0.0060 -0.1963 2.1346 2.5952

ezplot('0.0060*x^3-0.1963*x^2+2.1346*x+2.5952') 1,求下面的优化问题: min -5x1+4x2+2x3 6x1-x2+x3<=8 x1+2x2+4x3<=10 3>=x1>=-1; 2>=x2>=0; x3>=0;

用lingo求解。

7.2.4 解方程 1、代数方程 格式：solve (f,t)

功能：对变量t 解方程f=0，t 缺省时默认为x 或最接近字母x 的符号变量。 例如：求解一元二次方程f=a*x^2+b*x+c的实根，

>> syms a b c x

>> f=a*x^2+b*x+c; >> solve (f,x) ans=

[1/2/a*(-b+(b^2-4*a*c)^ (1/2))]

[1/2/a*(-b-(b^2-4*a*c)^ (1/2))] 2、微分方程

格式：dsolve(‘s’, ’s1’, ’s2’,…, ’x’)

1

其中s为方程；s1,s2,……为初始条件，缺省时给出含任意常数c1,c2,……的通解；x为自变量，缺省时默认为t 。

例如：求微分方程y1y2的通解 >> dsolve(‘Dy=1+y^2’) ans=

tan(t+c1)

1,方程 sin(x)0.3ex0的根. 答案：1.076 (填写程序语句)

可分别用matlab的指令sovle求解，或用lingo求解。

2