数学建模与数学实验三.doc

集美大学实验报告用纸 课程 数学建模与数学实验 日期 年 月 日 班级 学号 姓名 实验项目 Matlab工具箱 实验目的 通过本实验，掌握用matlab优化工具箱和统计工具箱的基本用法，能够求解简单的优化问题，能够对数据进行统计处理，能够对回归模型进行求解。 实验平台 Matlab7.0及以上版本 实验要求 分别建立M-文本文件205页第4题； 完成238页第1题，第2题和第4题； 完成2691题； 对90天报纸的需求量做正态分布的假设检验，对90天报纸的需求量做参数估计，确定报纸的最佳订购量； 完成317页第3题。 完成情况 实验成绩 实验代码及结果截图 1、.实验代码： function f = F_205_4(x,t,y) f = x(1)+x(2)*exp(-x(4)*t)+x(3)*exp(-x(5)*t)-y; clear,clc i = 1:33; t = 10*(i-1); y=[0.844 0.908 0.932 0.936 0.925 0.908 0.881 0.85 0.818 0.784 0.751 0.718 0.685 0.658 0.628 0.603 0.58 0.558 0.538 0.522 0.506 0.49 0.478 0.467 0.457 0.448 0.438 0.431 0.424 0.42 0.414 0.411 0.406]; x0 = [0.5 1.5 -1 0.01 0.02]; [x,norm,r,e,out] = lsqnonlin(@F_205_4,x0,[],[],[],t,y); x,norm 2．1、 clear,clc c = [1 -2 -1]; A1 = [-1 1 2; 1 1 -1;1 2 0;0 1 1]; b1 = [2 4 4 4]; A2 = [0 1 -1]; b2 = 1; v1 = [0 0 0]; [x,f,e,out,lag] = linprog(c,A1,b1,A2,b2,v1); f = -f lag.ineqlin 2.2、clear,clc H = [2 0 0 0 0;0 2 -2 0 0;0 -2 2 0 0;0 0 0 2 -2;0 0 0 -2 2]; c = [-2 0 0 0 0]; A2 = [1 1 1 1 1;0 0 1 -2 -2]; b2 = [5 -3]; [x,f,e,out,lag] = quadprog(H,c,[],[],A2,b2); x f 2.4、（1）function f = p238_4_1(x) f = exp(x(1)*x(2)*x(3)*x(4)*x(5)); if nargout1 g(1) = f*x(2)*x(3)*x(4)*x(5); g(2) = f*x(1)*x(3)*x(4)*x(5); g(3) = f*x(1)*x(2)*x(4)*x(5); g(4) = f*x(1)*x(2)*x(3)*x(5); g(5) = f*x(1)*x(2)*x(3)*x(4); end （2）function [c,ceq,g,geq] = p238_4_2(x) c = []; ceq = [x(1)^2+x(2)^2+x(3)^2+x(4)^2+x(5)^2-10 x(2)*x(3)-5*x(4)*x(5) x(1)^3+x(2)^3+1]; if nargout2 g = []; geq = [2*x(1) 2*x(2) 2*x(3) 2*x(4) 2*x(5); 0 x(3) x(2) 5*x(5) 5*x(4);3*x(1)^2 3*x(2)^2 0 0 0]; end （3）clear,clc v1 = [-2.3,-2.3,-2.3,-2.3,-2.3]; v2 = [2.3,2.3,3.2,3.2,3.2]; x0 = [1 1 1 1 1]; [x,fv,ef,out,lag] = fmincon(@p238_4_1,x0,[],[],[],[],v1,v2,@p238_4_2); x fv 3. 代码如下： clear,clc y1=normpdf([1.5 2 3.5],2,0.5) y2=chi2cdf([15 20 35],20) y3=tinv([0.9 0.8],10) y4=poisspdf(2:5,5) 结果截图: 4.实验代码： clear,clc A = [245 164 215 189 195 168 199 196 200 254 228 160 231 191 170;