数学建模作业4.docxVIP

姓名 黄耀东 学号 21480010222 评分 理学院实验报告 课程名称： 数学建模 实验名称： 回归分析 专业班级： 信息与计算科学21-2班 指导教师： 蔡军伟 完成日期： 2023年3月26日 3.1 已知矩阵A=145426563， 解：对于二次型f(x)的取值集合来看，单位球面属于凸集，所以这是一个凸规划问题，局部最优解即为全局最优解，建立非线性规划模型。 min s.t. 基于问题求解的Matlab程序如下： clc,clear, prob = optimproblem; a=[1 4 5;4 2 6;5 6 3]; x=optimvar(x,3,LowerBound,-1,UpperBound,1); prob.Objective=x*a*x; prob.Constraints.con=sum(x.^2)==1; x0.x=rand(3,1); [sol,fval,flag,out]=solve(prob,x0),sol.x 求得的最优解为=0.3130，=0.5774，=-0.7541，对应的最小值z=-3.6687. 3.2. 一个塑料大筐中装满了鸡蛋，两个两个数，余1个鸡蛋；三个三个地数，正好数完；四个四个地数，余1个鸡蛋；五个五个数，余4个鸡蛋；六个六个数，余3个鸡蛋；七个七个数，余四个鸡蛋；八个八个数，余1个鸡蛋；九个九个数，正好数完。建立数学规划模型，求大筐中鸡蛋个数的最小值。 解：引入变量，令y表示鸡蛋的总个数，建立整数线性规划模型： min s.t. 基于问题求解的MATLAB程序如下： clc,clear,prob = optimproblem; a=[2,3,4,5,6,7,8,9]; b=[1,0,1,4,3,4,1,0]; y=optimvar(y,Type,integer,LowerBound,0); x=optimvar(x,8,Type,integer,LowerBound,0); prob.Constraints.con=y-a.*x==b; [sol,fval,flag,out]=solve(prob),sol.y 求得的最优解为z=1089，即大筐中鸡蛋个数最小为1089个。 3.3 求解下列非线性整数规划问题： max s.t. 解：用蒙特卡洛法计算的Matlab程序如下： clc,clear %rng(shuffle) rng(0) p0=0;n=10^6;tic for i =1:n x=randi([0,99],1,5); [f,g]=mengte(x); if all(g=0) if p0 f x0=x;p0=f end end end x0,p0,toc function[f,g]=mengte(x); f=x(1)^2+x(2)^2+3*x(3)^2+4*x(4)^2+2*x(5)^2-8*x(1)-2*x(2)-3*x(3)-x(4)-2*x(5); g=[sum(x)-400 x(1)+2*x(2)+2*x(3)+x(4)+6*x(5)-800 2*x(1)+x(2)+6*x(3)-200 x(3)+x(4)+5*x(5)-200]; end 求得的解为=46，=98，=1，x4=99，x5=3，对应的最大值为50273。 3.5 求解下列非线性规划问题： max x 解：基于问题求解的Matlab程序如下： clc,clear,prob=optimproblem(ObjectiveSense,max); x=optimvar(x,3); prob.Objective=2*x(1)+3*x(1)^2+3*x(2)+x(2)^2+x(3); con1=[x(1)+2*x(1)^2+x(2)+2*x(2)^2+x(3)=10;x(1)+x(1)^2+x(2)+x(2)^2-x(3)=50;2*x(1)+x(1)^2+2*x(2)+x(3)=40;-x(1)-2*x(2)=1;-x(1)=0] con2=[x(1)^2+x(3)==2] prob.Constraints.con1=con1; prob.Constraints.con2=con2; x0.x=rand(3,1); [sol,fval,flag,out]=solve(p