北师MATLAB课件实验课三.pptVIP

数学模型实验课（三） 矩阵、解方程组、拟合 （一）解方程组——矩阵形式 康尔乃奶粉 32.4元 400g; 67.1元 900g. 400 k1 + 4002/3 k2 = 32.4 900 k1 + 9002/3 k2 = 67.1 a1=[400,900];a2=a1.^(2/3);A=[a1’,a2’]; b=[32.4,67.1]’;k=A\b 解得: k1 = 0.0538, k2 = 0.2005 模型: C(W)=0.0538 W + 0.2005 W2/3. c(w)=k(1).*w+k(2).*(w.^(2/3)); 预测: W=1800, W=2500, C(W) w=[1800, 2500];c(w) C(1800)=126.48, C(2500)=171.41 实际: C(1800)=115.90, C(2500)=146.85 (二)数据中总体信息的初步描述 在齿轮加工中，齿轮的径向综合误差是个随机 变量，对40件同样的齿轮进行的测量，误差数据 (mm)如下， 16 20 18 13 14 19 19 18 8 17 24 21 28 13 16 28 18 33 24 18 17 24 13 11 21 28 28 18 22 28 20 11 19 15 10 12 22 21 27 33 其中最大值为33，最小值为8。 误差服从什么样的分布？ Matlab程序 x=[16 20 18 13 14 19 19 18 8 17 24 21 28 13 16 28 18 33 24 18 17 24 13 11 21 28 28 18 22 28 20 11 19 15 10 12 22 21 27 33]; g=[7.5 9.5 11.5 13.5 15.5 17.5 19.5 21.5 23.5 25.5 27.5 29.5 33.5]; f=histc(x,g); %计算频数 h=f/length(x); %计算频率 figure(1); bar(g,f,histc) %画频数直方图 h=lillietest(x) %用Lilliefore检验,当h=0时，确定x服从正态分布。 figure(2);h=normplot(x) %正态概率纸检验 . （三）解超定方程：A x = y 如果 设：Ax+ε= y，求 x 使得误差 ||ε|| 最小。 A’A x = A’y x = (A’A)-1A’y 正规方程组 人的身高与体重 W=bL3=L3b 方法I——最小二乘 w=[12,17,22,35,48,54,66,75]; l=[0.86,1.08,1.16,1.35,1.55,1.67,1.78,1.85]; L=l.^3; b=L’\w’; w1=b.*L; plot(L,w,’*’,L,w1,’-r’) 模型二：W＝rL3， 其中 a=w./l.^3; ll=length(a); r=sum(a)/ll; plot(l,w,’o’,l,r.*L,’-b’) 练习： 模型III ：W＝rL3 其中 (三) 拟合 1.直线拟和 数据（xi, yi）,拟合直线 y = a + bx, 误差最小 模型：yi = a + bxi +ε 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,n); 4.线性回归 [b,bint,R,Rint,stats]=regress(y’,X); b：回归系数，bint：置信区间， R：剩余（误差），Rint 置信区间， Stats：统计量：r2，F，p，σ2。 例：根据1977年的旧格局中级别抓举的世界纪录 组建模型描述举重成绩与运动员体重间的关系 重量级 53 56 60 67.5 75 82.5 90 110 110 举重量 109 120.5 130 141.5 157.5 170