数学建模小题目及答案

1.求下列积分的数值解：
⎰
+∞
+-⋅2
3
2
2
3x x x dx
function y = myfun(x)
y = 1./(x.*(x.^2 - 3*x + 2 ).^(1/3)); warning off all
Q = quad(@myfun,2,100000) Q = quad(@myfun,2,10000000)
Q = quad(@myfun,2,1000000000000000) warning on
当上限为100000,10000000,1000000000时， 定积分的值为x=1.4389,1.4396,1.4396。

因此，可以将1.4396作为此定积分的值。

2.已知
)s i n ()()c o s (),(2h t h t h t e h t f h t ++++=+，dt h t f h g ⎰=10
),()(，画出
]10,10[-∈h 时，)(h g 的图形。

syms t,syms h;
f=exp(t+h)*cos(t+h)+(t+h)^2*sin(t+h); int(f,t,0,10) ans =
1/2*exp(10+h)*cos(10+h)+1/2*exp(10+h)*sin(10+h)-98*cos(10+h)-20*cos(10+h)*h-cos(10+h)*h^2+20*sin(10+h)+2*sin(10+h)*h-1/2*exp(h)*cos(h)-1/2*exp(h)*sin(h)+cos(h)*h^2-2*cos(h)-2*sin(h)*h
ezplot('1/2*exp(10+h)*cos(10+h)+1/2*exp(10+h)*sin(10+h)-98*cos(10+h)-20*cos(10+h)*h-cos(10+h)*h^2+20*sin(10+h)+2*sin(10+h)*h-1/2*exp(h)*cos(h)-1/2*exp(h)*sin(h)+cos(h)*h^2-2*co s(h)-2*sin(h)*h',[-10,10])
3.画出16)5(2
2
=-+y x 绕x 轴一周所围成的图形，并求所产生的旋转体的体积。

主程序：
[y,z]=cylinder(1:0.2:9,100);
mesh(sqrt(16-(sqrt(y.^2+z.^2)-5).^2),y,z); hold on;
mesh(-sqrt(16-(sqrt(y.^2+z.^2)-5).^2),y,z);
求体积dv=
⎰⎰
--2)5(16y dydz
先计算在第Ⅰ卦限的体积 1≤y 2+z 2
≤9 计算⎰
--9
12
)5(16y dy ⎰
--22
911y y dz
=
(
29
1
29)5(16y y ---⎰-)2
21)5(16y y ---
clear
syms z m y
m=sqrt(16-(y-5)^2)*(sqrt(9-y^2)-sqrt(1-y^2)); int(m,1,9)
ans =
-118/3*i-19*3^(1/2)*EllipticK(1/2*2^(1/2))-35/3*3^(1/2)*EllipticPi(1/3,1/2*2^(1/2))+50*3^(1/2)*EllipticE(1/2*2^(1/2))-58/3*i*3^(1/2)*EllipticK(1/2*2^(1/2))+35/6*i*3^(1/2)*EllipticPi(3/4,1/2*2^(1/2))+50*i*3^(1/2)*EllipticE(1/2*2^(1/2))+75/2*log(5)-75/2*log(-3+4*i)
V=8(118/3*i-19*3^(1/2)*EllipticK(1/2*2^(1/2))-35/3*3^(1/2)*EllipticPi(1/3,1/2*2^(1/2))+50*3^(1/2)*EllipticE(1/2*2^(1/2))-58/3*i*3^(1/2)*EllipticK(1/2*2^(1/2))+35/6*i*3^(1/2)*EllipticPi(3/4,1/2*2^(1/2))+50*i*3^(1/2)*EllipticE(1/2*2^(1/2))+75/2*log(5)-75/2*log(-3+4*i))
4.画出下列曲面的图形
（1）旋转单叶双曲面
14
92
22=-+z y x ； x=@(s,t)3.*sec(s).*cos(t);
y=@(s,t)3.*sec(s).*sin(t); z=@(s,t)2.*tan(s); ezmesh(x,y,z) 或者
t=-pi/4:0.1:pi/4; r=0:0.1:2*pi;
[r,t]=meshgrid(r,t); x=3*sec(t).*sin(r); y=3*sec(t).*cos(r); z=2*tan(t); surf(x,y,z)
-5
5
-5
5
-2-1
1
2
（2）马鞍面xy z =； x=-2*pi:0.2:2*pi;
[x,y]=meshgrid(x); z=x.*y; surf(x,y,z);
xlabel('x-axis'),ylabel('y-axis'),zlabel('z-axis'); title('surf')； 或者
x=-2:0.1:2; y=-2:0.1:2;
[xx,yy]=meshgrid(x,y); zz=xx.*yy; surf(xx,yy,zz) 或者
ezsurf(@(x,y)x*y)
5.画出隐函数1cos sin =+y x 的图形。

ezplot('cos(y)+sin(x)-1',[-2*pi,2*pi,-2*pi,2*pi])
6.（1）求函数x
x y -+=12
ln
的三阶导数； clear syms x
diff('log((x+2)/(1-x))',x,3) （2）求向量]425.00
[=a 的一阶向前差分。

a=[0,0.5,2,4]; i=1:3;
b=a(i+1)-a(i)
7.求解非线性方程组
（1）⎪⎩⎪⎨⎧=-+=-+0
60622x y y x
[x,y]=solve('x^2+y-6=0','y^2+x-6=0','x','y')
（2）⎩
⎨⎧=+=++5ln 10tan 10cos sin y x y e y x
[x,y]=solve('exp(x+sin(y))+cos(y)=10','tan(x)+10*log(y)=5','x','y')
8.求函数186)(23-++=x x x x f 的极值点，并画出函数的图形。

clear; syms x
y=x^3+6*x^2+8*x-1; dy=diff(y) x=solve(dy) x=double(x)
作图：fplot('x^3+6*x^2+8*x-1',[-5,1])
9.某单位需要加工制作100套钢架，每套用长为2.9m ，2.1m 和1m 的圆钢各一根。

已知原料长6.9m ，问应如何下料，使用的原材料最省。

解：最简单做法是，在每一根原材料上截取2.9m ，2.1m 和1.5m 的元钢各一根组成一套，每根原材料剩下料头0.9m （7.4-2.9-2.1-1.5=0.9)。

为了做100套钢架，需用原材料100根，共有90m 料头。

若改为用套裁，这可以节约原材料。

下面有几种套裁方案，都可以考虑采用。

见表1-5。

结果 方案 2.9m 2.1m 1m 余料 1 1 1 1 0.9 2 1 0 4 0 3 2 0 1 0.1 4 0 0 6 0.9 5 0 1 4 0.8 6 0 2 2 0.7 7
3
0.6
为了得到100套钢架，需要混合使用各种下料方案。

设按1方案下料的原材料根数为x1,2方案为x 2，3方案为x 3，4方案为x 4,5方案为x 5，6方案为x 6, 7方案为x 7。

根据表1-11的方案，可列出以下数学模型：
{z 1234567min 0.900.10.90.80.7+0.6x x x x x x x =+++++
1231567
1234561234567x +x +2x 100
x +x +2x +3x 100x +4x +x +6x +4x +2x 100x ,x ,x ,x ,x ,x ,x 0
>=⎧⎪>=⎪⎨
>=⎪⎪>⎩ LINGO 程序如下：
min=x1+x2+x3+x4+x5+x6+x7; x1+x2+2*x3>=100; x1+x5+2*x6+3*x7>=100;
x1+4*x2+x3+6*x4+4*x5+2*x6>=100; @gin(x1); @gin(x2); @gin(x3); @gin(x4); @gin(x5); @gin(x6); @gin(x7);
Feasible solution found:
Extended solver steps: 0 Total solver iterations: 0
Variable Value MINZ 27.90000 X1 8.000000 X2 46.00000 X3 31.00000 X4 1.000000 X5 5.000000
Row Slack or Surplus 1 0.000000 2 0.000000 3 0.000000 4 0.000000
N=91根
10. 某部门在今后五年内考虑给下列项目投资，已知：
项目A ，从第一年到第四年每年年初需要投资，并于次年末回收本利115％；
项目B ，从第三年初需要投资，到第五年末能回收本利125％，但规定最大投资额不超过4万元；
项目C ，第二年初需要投资，到第五年末能回收本利140％，但规定最大投资额不超过3万元；
项目D ，五年内每年初可购买公债，于当年末归还，并加利息6％。

该部门现有资金10万元，问它应如何确定给这些项目每年的投资额，使到第五年末拥有的资金的本利总额为最大？
第j 年 项目 1 2 3 4 5 A X1A X2A X3A X4A B X3B C X2C D
X1D
X2D
X3D
X4D
X5D
解：
分析：第一年有10万元资金，10000011=+D X A X ；
第二年投资，D X D X C X A X 106.1222=++；
第三年投资：D X A X D X B X A X 206.1115.1333+=++； 第四年投资：A X D X D X A X 215.1306.144+=+； 第五年投资：D X A X D X 406.1315.15+=
其中，400003<=B X ,300002<=C X ,由于“连续投资问题”要求第五年末部门所拥有的资金的本利总额最大，故目标函数为： D X B X C X A X Z 506.1325.1240.1415.1max +++= 建立“连续投资问题”的线性规划模型：
D X B X C X A X Z 506.1325.1240.1415.1max +++=
⎪⎪⎪
⎪⎩⎪
⎪⎪
⎪⎨⎧<=<=+=+=++=++=++=+300002400003406.1315.15215.1306.144115.1206.1333106.1222100000
11C X B X D X A X D X A X D X D X A X A X D X D X B X A X D X D X C X A X D X A X
用Lingo 求解：
max=1.15*X4A+1.25*X3B+1.4*X2C+1.06*X5D; X1A+X1D=100000;
X2A+X2C+X2D-1.06*X1D=0;
X3A+X3B+X3D-1.15*X1A-1.06*X2D=0;
X4A+X4D-1.15*X2A-1.06*X3D=0;
X5D-1.15*X3A-1.06*X4D=0;
X3B<=40000;
X2C<=30000;
Z=143750元。

第一年：X1A=71698.11元X1D=28301.89元
第二年：X2A=0 X2C=30000元X2D=0
第三年：X3A=0 X3B=40000元X3D=42452.83元
第四年：X4A=45000元X4D=0
第五年：X5D=0
11.已知某工厂计划生产Ⅰ，Ⅱ，Ⅲ三种产品，各种产品需要在A，B，C三种设备上加工生产，具体相关数据如表1。

试研究下列问题。

（1）如何充分发挥已有设备的能力，使生产盈利最大？
（2）如果为了增加产量，可租用其他厂家设备B，每月可租用60台时，租金为1.8万元，试问租用设备B是否合算？
（3）如果该工厂拟增加生产两种新产品Ⅳ和Ⅴ，其中产品Ⅳ需用A设备12台时，B设备5台时，C设备10台时，单位产品盈利21000元；产品Ⅴ需用A设备4台时，B设备4台时，C设备12台时，单位产品盈利1870元。

假如A，B，C三种设备台时不增加，试分别考虑这两种新产品的投产在经济上是否合算？
表1 生产计划的相关数据
ⅠⅡⅢ设备有效台时/每月
A 8 2 10 300
B 10 5 8 400
C 2 13 10 420
单位产品利润/元3000 2000 2900
解：
（1）设使用A为x1台时，B使用x2台时，C是使用x3台时：
程序
max=3000*x1+2000*x2+2900*x3;
8*x1+2*x2+10*x3<=300;
10*x1+5*+5*x2+8*x3<=400;
2*x1+13*x2+10*x3<=420;
Global optimal solution found.
Objective value: 134500.0
Extended solver steps: 6
Total solver iterations: 29
Variable Value Reduced Cost
X1 24.00000 -3000.000
X2 24.00000 -2000.000
X3 5.000000 -2900.000
Row Slack or Surplus Dual Price
1 134500.0 1.000000
2 10.00000 0.000000
3 0.000000 0.000000
4 10.00000 0.000000
（2）max=3000*x1+2000*x2+2900*x3;
8*x1+2*x2+10*x3<=300;
10*x1+5*x2+8*x3<=460;
2*x1+13*x2+10*x3<=420;
@gin(x1);
@gin(x2);
@gin(x3);
综Global optimal solution found.
Objective value: 145000.0
Extended solver steps: 0
Total solver iterations: 4
Variable Value Reduced Cost X1 31.00000 -3000.000
X2 26.00000 -2000.000
X3 0.000000 -2900.000
Row Slack or Surplus Dual Price
1 145000.0 1.000000
2 0.000000 0.000000
3 20.00000 0.000000
4 20.00000 0.000000
合考虑方案（1）和方案（2），方案（2）经济效益为145000不合理.
（3）目标函数：
maxz=3000*x1+2000*x2+2900*x3+21000*x4+1870*x5;
8*x1+2*x2+10*x3+12*x4+4*x5<=300;
10*x1+5*x2+8*x3+5*x4+4*x5<=400;
2*x1+13*x2+10*x3+10*x4+12*x5<=420;
Global optimal solution found.
Objective value: 135960.0
Extended solver steps: 25
Total solver iterations: 202
Variable Value Reduced Cost
X1 26.00000 -3000.000
X2 19.00000 -2000.000
X3 1.000000 -2900.000
X4 1.000000 -2100.000
X5 8.000000 -1870.000
Row Slack or Surplus Dual Price 1 135960.0 1.000000 2 0.000000 0.000000 3 0.000000 0.000000 4 5.000000 0.000000
结论：最大经济效益为1356960元时，x1,x2,x3,x4,x5分别为26、19、1、1、8. 12. 某市政府拟投入一笔资金和一定数量的劳动力建设两类公益项目A 和B ，目的是方便市民的生活，提高城市的生活质量。

根据预测投入1万元资金和1百个劳动力·h （即每个劳动力用1h ）,分别可以建成1个项目A 和两个项目B 。

如果投入1个劳动力·h 需要支付10元，市政府为了用有限的资金和劳动力，并用最快的时间建成这批项目，服务于社会，服务于人民。

市政府依次提出下面的四条要求。

（1）至少要建50个项目A ； （2）至多建设60个项目B ；
（3）至少要利用80万元资金和10000个劳动力·h ； （4）总投入资金不超过预算120万元。

试为该市政府制定一个满意的项目建设方案。

解：
第一目标：11
P d +
第二目标：22P d -
第三目标：33P d +
规划函数：
min 1122334455Z P d P d P d P d P d +-++-=++++
115022601.110.55238010.5241001.110.5525120x d x d x x d x x d x x d +-++-⎧-=⎪+=⎪⎪+-=⎨
⎪+-=⎪
⎪++=⎩
应用lingo 软件，得项目A 数目为70，项目B 数目为60. 程序：
min=p1*d11+p2*d22+p3*d31+p4*d41+p5*d52; x1-d11=50; x2+d22=60;
1.1*x1+0.55*x2-d31=80; x1+0.5*x2-d41=100;
1.1*x1+0.55*x2+d52=120;
运行结果为
Local optimal solution found.
Objective value: 0.000000 Total solver iterations: 4
Variable Value Reduced Cost
P1 0.000000 20.00000 D11 20.00000 0.000000 P2 0.000000 0.000000 D21 0.000000 0.000000 P3 0.000000 30.00000 D31 30.00000 0.000000 P4 0.000000 0.000000 D41 0.000000 0.000000 P5 0.000000 10.00000 D52 10.00000 0.000000 X1 70.00000 0.000000 X2 60.00000 0.000000
Row Slack or Surplus Dual Price 1 0.000000 -1.000000 2 0.000000 0.000000 3 0.000000 0.000000 4 0.000000 0.000000 5 0.000000 0.000000 6 0.000000 0.000000
13.求微分方程组初值问题
⎪⎪⎪
⎩
⎪⎪⎪⎨⎧==+-=-=00)0(,)0(y y x x bxy
sy dt dy
axy rx dt dx
式中，2=r ，1=s ，1=a ，2=b 。

选用ode45函数计算，其相对误差限为5
10-，绝对误差限为6
10-，分别画出初值条件为]3.0,1[],[00=y x ，]5.0,1[，]7.0,1[，]9.0,1[，]1.1,1[解的相平面轨迹图。

程序：新建一个shier.m 文件 function dx=shier(t,x); dx=zeros(2,1);
dx(1)=(2*x(1)-x(1)*x(2));
dx(2)=(-x(2)+2*x(1)*x(2));
options = odeset('RelTol',1e-5,'AbsTol',1e-6); [t,x]=ode45('shier',[0,12],[1,0.3]);
plot(t,x(:,1),'-',t,x(:,2),'*')
[t,x]=ode45('shier',[0,12],[1,0.5]);
plot(t,x(:,1),'-',t,x(:,2),'*')
[t,x]=ode45('shier',[0,12],[1,0.7]);
plot(t,x(:,1),'-',t,x(:,2),'*')
[t,x]=ode45('shier',[0,12],[1,0.9]);
plot(t,x(:,1),'-',t,x(:,2),'*')
[t,x]=ode45('shier',[0,12],[1,1.1]);
plot(t,x(:,1),'-',t,x(:,2),'*')
图一
图二
图三
图四
图五
14．求两个圆10022=+y x ，100)4()3(22=-+-y x 所围公共部分的面积。

程序：
ezplot('x^2+y^2-100',[-20,20]); hold on;
ezplot('(x-3)^2+(y-4)^2-100',[-20,20]); >> syms x;
g=sqrt(100-(x+2.5)^2); int(g,x,0,7.5)
ans =25*pi-25/8*15^(1/2)-50*asin(1/4) double(ans) ans = 53.8027
>> s=4*double(ans) s = 215.2109
15. 已知平面区域56000≤≤x ，48000≤≤y 的高程数据见表3（单位：m ）。

表3
4800 4400 4000 3600 3200 2800 2400 2000 1600 1200 800 400 0
1350 1370 1390 1400 1410 960 940 880 800 690 570 430 290 210 150 **** **** 1410 1430 1440 1140 1110 1050 950 820 690 540 380 300 210 1380 1410 1430 1450 1470 1320 1280 1200 1080 940 780 620 460 370 350 1420 1430 1450 1480 1500 1550 1510 1430 1300 1200 980 850 750 550 500 1430 1450 1460 1500 1550 1600 1550 1600 1600 1600 1550 1500 1500 1550 1550 950 1190 1370 1500 1200 1100 1550 1600 1550 1380 1070 900 1050 1150 1200 910 1090 1270 1500 1200 1100 1350 1450 1200 1150 1010 880 1000 1050 1100 880 1060 1230 1390 1500 1500 1400 900 1100 1060 950 870 900 936 950 830 980 1180 1320 1450 1420 400 1300 700 900 850 810 380 780 750 740 880 1080 1130 1250 1280 1230 1040 900 500 700 780 750 650 550 650 760 880 970 1020 1050 1020 830 800 700 300 500 550 480 350 510 620 730 800 850 870 850 780 720 650 500 200 300 350 320 370 470 550 600 670 690 670 620 580 450 400 300 100 150 250 X Y /
0 400 800 1200 1600 2000 2400 2800 3200 3600 4000 4400 4800 5200 5600
试用二维插值求y x ,方向间隔都为10的高程，画出该区域的等高线和三维视图，并求该区域的表面积。

建一个moutain.m 文件 x=0:400:5600; y=0:400;4800;
z=[1350 1370 1390 1400 1410 960 940 880 800 690 570 430 290 210 150
1370 1390 1410 1430 1440 1140 1110 1050 950 820 690 540 380 300 210
1380 1410 1430 1450 1470 1320 1280 1200 1080 940 780 620 460 370 350
1420 1430 1450 1480 1500 1550 1510 1430 1300 1200 980 850 750 550 500
1430 1450 1460 1500 1550 1600 1550 1600 1600 1600 1550 1500 1500 1550 1550
950 1190 1370 1500 1200 1100 1550 1600 1550 1380 1070 900 1050 1150 1200
910 1090 1270 1500 1200 1100 1350 1450 1200 1150 1010 880 1000 1050 1100
880 1060 1230 1390 1500 1500 1400 900 1100 1060 950 870 900 936 950
830 980 1180 1320 1450 1420 400 1300 700 900 850 810 380 780 750
740 880 1080 1130 1250 1280 1230 1040 900 500 700 780 750 650 550
650 760 880 970 1020 1050 1020 830 800 700 300 500 550 480 350
510 620 730 800 850 870 850 780 720 650 500 200 300 350 320
370 470 550 600 670 690 670 620 580 450 400 300 100 150 250]
figure(1);
meshz(x,y,z)
xlabel('X'),ylabel('Y'),zlabel('Z')
xi=0:10:5600;
yi=0:10:4800;
>> figure(2)
z1i=interp2(x,y,z,xi,yi','nearest');
surfc(xi,yi,z1i)
xlabel('X'),ylabel('Y'),zlabel('Z')
>> figure(3)
z2i=interp2(x,y,z,xi,yi');
surfc(xi,yi,z2i)
xlabel('X'),ylabel('Y'),zlabel('Z')
>> figure(4)
z3i=interp2(x,y,z,xi,yi','cubic');
surfc(xi,yi,z3i)
xlabel('X'),ylabel('Y'),zlabel('Z')
figure(5)
subplot(1,3,1),contour(xi,yi,z1i,10,'r');
subplot(1,3,2),contour(xi,yi,z2i,10,'r');
subplot(1,3,3),contour(xi,yi,z3i,10,'r');
16. 已知矩阵⎥⎥
⎥⎦
⎤
⎢⎢⎢⎣⎡=365624541A ，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=321x x x x ，求二次型Ax x x x x f T =),,(321在单位球面12
32221=++x x x 上的最小值。

MATLAB 程序：
syms x1 x2 x3;
A=[1,4,5;4,2,6;5,6,3]; x=[x1;x2;x3]; t=[x1,x2,x3]; f=t*A*x f =
(x1+4*x2+5*x3)*x1+(4*x1+2*x2+6*x3)*x2+(5*x1+6*x2+3*x3)*x3
LINGO 程序：
min=(x1+4*x2+5*x3)*x1+(4*x1+2*x2+6*x3)*x2+(5*x1+6*x2+3*x3)*x3; x1^2+x2^2+x3^2=1;
过MATLAB 计算得出
f= (x1+4*x2+5*x3)*x1+(4*x1+2*x2+6*x3)*x2+(5*x1+6*x2+3*x3)*x3 代入LINGO 的出 X1=x2=0 x3=1 f=3；
17. 求解线性规划问题：
5432113070809020 max x x x x x Z ++++=
⎪⎪⎪⎩⎪
⎪⎪⎨⎧=≥≤++≤+≥+≥++5
,,1 ,048
23120
2330305
423143521 i x x x x x x x x x x x i 程序：max=20*x1+90*x2+80*x3+70*x4+30*x1; s.t. x1+x2+x5>=30; x3+x4>=30;
3*x1+2*x3<=120; 3*x2+2*x4+x5<=48;
18.18.求解数学规划问题
∑=1000
1
max i i
x
s.t.⎩
⎨⎧≥≤⨯011500x Ax 其中
1000500)(⨯=ij a A ，这里ij a 是服从均值为5，标准差为2的正态分布的随
机数；15001⨯表示500个元素全部为1的列向量。

：
标注c 为110001⨯表示1000个元素全部为1的行向量；
A 为服从均值为5，标准差为2的正态分布的随机数组成的矩阵； D 为15001⨯表示500个元素全部为1的列向量。

主程序：
经计算得出max=0.2022
f=-1*[1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;
1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1 ;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1; 1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1 ;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1; 1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1 ;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1; 1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1 ;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1; 1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1 ;1;1];
R=random('Normal',5,2,500,1000);
A=R;
b=[1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1, 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1 ,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1, 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1 ,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1, 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1 ,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1, 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1 ,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1, 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1 ,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1, 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1 ,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1, 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1 ,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1]
lb=zeros(1000,1);
[x,fval]=linprog(f,A,b,[],[],lb)
经计算得出max=-0.2022
19.根据表4某猪场24头育肥猪4个胴体性状的数据资料，试进行瘦肉量y对眼肌面积（x1）、腿肉量(x2)、腰肉量(x3)的多元回归分析。

表4 某养猪场数据资料
序号瘦肉量
y(kg)
眼肌面积
x1(cm2)
腿肉量
x2(kg)
腰肉量
x3(kg)
序
号
瘦肉量
y(kg)
眼肌面积
x1(c m2)
腿肉量
x2(kg)
腰肉量
x3(kg)
1 15.0
2 23.7
3 5.49 1.21 1
4 15.94 23.52 5.18 1.98
2 12.62 22.34 4.32 1.35 15 14.3
3 21.86 4.86 1.59
3 14.86 28.8
4 5.04 1.92 16 15.11 28.9
5 5.18 1.37
4 13.98 27.67 4.72 1.49 17 13.81 24.53 4.88 1.39
5 15.91 20.83 5.35 1.5
6 18 15.58 27.65 5.02 1.66
6 12.4
7 22.27 4.27 1.50 19 15.85 27.29 5.55 1.70 7 15.80 27.57 5.25 1.85 20 15.2
8 29.07 5.26 1.82 8 14.32 28.01 4.62 1.51 21 16.40 32.47 5.18 1.75
9 13.76 24.79 4.42 1.46 22 15.02 29.65 5.08 1.70 10 15.18 28.96 5.30 1.66 23 15.73 22.11 4.90 1.81 11 14.20 25.77 4.87 1.64 24 14.75 22.43 4.65 1.82 12 17.07 23.17 5.80 1.90 13
15.40
28.57
5.22
1.66
要求
（1）求y 关于321,,x x x 的线性回归方程
3322110x c x c x c c y +++=，
计算3210,,,c c c c 的估计值；
（2）对上述回归模型和回归系数进行检验（要写出相关的统计量）；
（3）试建立y 关于321,,x x x 的二项式回归模型，并根据适当统计量指标选择一个较好的模型。

（4）利用表4的数据分别利用Matlab 的命令lsqcurvefit 和nlinfit 拟合非线性函数
)tan(cos sin 24321113x a x e a x a y x a ++=
(1)程序：y=[15.02 12.62 14.86 13.98 15.91 12.47 15.8 14.32 13.76 15.18 14.2 17.07 15.4 15.94 14.33 15.11 13.81 15.58 15.85 15.28 16.4 15.02 15.73 14.75];
x1=[23.73 22.34 28.84 27.67 20.83 22.27 27.57 28.01 24.79 28.96 25.77 23.17 28.57 23.52 21.86 28.95 24.53 27.65 27.29 29.07 32.47 29.65 22.11 22.43];
x2=[5.49 4.32 5.04 4.72 5.35 4.27 5.25 4.62 4.42 5.3 4.87 5.8 5.22 5.18 4.86 5.18 4.88 5.02 5.55 5.26 5.18 5.08 4.9 4.65];
x3=[1.21 1.35 1.92 1.49 1.56 1.5 1.85 1.51 1.46 1.66 1.64 1.9 1.66 1.98 1.59 1.37 1.39 1.66 1.7 1.82 1.75 1.7 1.81 1.82]; X=[ones(24,1), x1',x2',x3'];
[c,bint,r,rint,s]=regress(y',X); c,bint,s,rcoplot(r,rint) 结果 c =
1.0272 0.0074
2.1141
1.8980
所以
y=1.2072+0.0074x1+2.1141x2+1.8980x3
5
101520
-1.5
-1
-0.5
0.5
1
1.5
Residual Case Order Plot
R e s i d u a l s
Case Number。