数学模型实验课二.pptVIP

y=[17,16,17,23,26,27,41,49]; x=[1953,2015,2015,2821,3049,3049,5133,5592]; X=[ones(8,1),x’];b=X\y’; y1=b(1)+b(2).*x; plot(x,y1,’r’,x,y,’*r’) 2. 直线拟合的效果（相关系数） X1=[x’,y’];m=size(X1);n=m(1);p=m(2); X=[ones(n,1),x’];B=X\y’;a=B(1);b=B(2:p); y1=a+b*x;plot(x,y1,’r’,x,y,’*r’); Cov=X1’*(eye(n)-ones(n)./n)*X1; Lyy=Cov(p,p);Lxx=Cov(1:(p-1),1:(p-1)); R=(b’*Lxx*b)/Lyy;r=sqrt(R); 3. 多项式拟合 b=polyfit(x,y,1); 4.线性回归 [b,r,j,k,l]=regress(y’,X);