试题(有效约束回归分析区间估计迭代法收敛域).docVIP

数学实验试题 2003.6.22 上午 班级 姓名 学号 得分 说明： （1）第一、二、三题的答案直接填在试题纸上； （2）第四题将数学模型、简要解题过程和结果写在试题纸上；卷面空间不够时，可写在背面； （3）考试时间为90分钟。 一．（10分，每空2分）（计算结果小数点后保留4位有效数字） 某两个地区上半年6个月的降雨量数据如下（单位：mm）： 月份 1 2 3 4 5 6 地区A 25 99 46 33 70 54 地区B 10 50 30 20 45 30 在90%的置信水平下，给出A地区的月降雨量的置信区间： [32.35 76.65] x=[25 99 46 33 70 54]; [mu, sigma, muci, sigmaci]=normfit(x,0.1) 输出结果： muci = 32.348211120103912 76.651788879896088 在90%的置信水平下，A地区的月降雨量是否不小于70（mm）？ 是 H0：≥70 x=[25 99 46 33 70 54]; [h,sig,ci] =ttest(x,70,0.1,-1) 输出结果： h =0 在90%的置信水平下，A、B地区的月降雨量是否相同？ 否 x=[25 99 46 33 70 54]; y=[10 50 30 20 45 30]; [h,sig,ci]=ttest2(x,y,0.1) 输出结果： h =1 A地区某条河流上半年6个月对应的径流量数据如下（单位：m3）：110，184，145，122，165，143。该河流的径流量y与当地的降雨量x的线性回归方程为 y=91.12+0.9857*x ；若当地降雨量为55mm，该河流的径流量的预测区间为 [130.9 159.7] （置信水平取90%）。 y=[110 184 145 122 165 143]; x=[25 99 46 33 70 54]; n=6 T=[ones(n,1),x]; [b,bint,r,rint,s]=regress(y,T); b,bint,s, rcoplot(r,rint) 输出结果： b =91.115018618121596 0.985657150737830 少量数据回归分析区间预测：y=b0+b1*x %预测区间函数M文件源程序： function y30=guessless(x0,x,s2,b0,b1,alpha) s=sqrt(s2); %将s2开方即得 xx=mean(x); %x的均值 sxx=(length(x)-1)*var(x); y0=b0+b1*x0; y30=[y0-tinv(1-alpha/2,length(x)-2)*s*sqrt((x0-xx)^2/sxx+1/length(x)+1), y0+tinv(1-alpha/2,length(x)-2)*s*sqrt((x0-xx)^2/sxx+1/length(x)+1)]; %利用（30）式预测区间 %预测55mm降雨量对应河流径流量所在区间源程序(需要先运行回归程序！) x=[25 99 46 33 70 54]; y55=guessless(55,x, s(4),b(1),b(2),0.1) 输出结果： y55 =1.0e+002 * 1.309185906962643 1.597337331211401 还可以应用POLYTOOL命令： y=[110 184 145 122 165 143]; x=[25 99 46 33 70 54]; Polytool(x,