MATLAB实验报告7.docVIP

1．Carpenter Problem. A carpenter makes tables and bookcases. He is trying to determine how many each type of furniture he should make each week The carpenter realizes a net unit profit of $25 per table and $30 per bookcase. He has up to(多达) 600 board feet of lumber to devote weekly to the project and up to 120 hours of labor. He estimates that it requires 20 board feet of lumber and 5 hours of labor to complete a table and 30 board feet of lumber and 4 hour of labor for a bookcase. Moreover, he signed a contract to deliver four tables and two bookcases each week. The carpenter wishes to determine a weekly production schedule for tables and bookcases that maximize his profits. 设桌子的个数为x1，书架的个数为x2 f=[-25,-30]; A=[20 30;5 4]; b=[600;120]; Aeq=[];beq=[]; lb=[4;2]; ub=[inf,inf]; [x,feval]=linprog(f,A,b,Aeq,beq,lb,ub) Optimization terminated. x = 17.1429 8.5714 feval = -685.7143 则最优解为桌子17.1429 书架8.5714 最大利润为685.7143 2. Solve the Ladder Problem by using three methods. (Animation is optional) Suppose there is a big house with two levels. A flower house made of glass is next to it. On the side where the two houses are next to each other, there is a window. Assume we put a ladder onto the wall over the glass house, so that we can clean the window by standing on top of the ladder. What is the shortest length of the ladder? 设房屋的高度为H=10米，花房高度为h=5米，宽度s=3米。 （1）第一种解法： 设梯子总长为z，上半段长为x，下半段长为y，则有 条件： x/(x+y)=s/(s+(y^2-h^2)^(1/2))即为x*(y^2-25)(1/2)-3*y=0 x^2-s^2=(H-h)^2,y0即为x^2=34 求解： Min z=x+y 代码： function f=fun(x) f=x(1)+x(2) function [c,ceq]=confun(x) c=[] ceq=[x(1)*(x(2)^2-25)^(1/2)-3*x(2)] [x,fval]=fmincon(@fun,[4,6],[],[],[],[],[0,0],[34^(1/2),inf],@confun) 结果 No active inequalities. x = 4.6531 6.5410 fval = 11.1941 （2）第二种解法： 设梯上小三角形高度为x，下小三角形的底为y，则 条件： x/(x+h)=s/(s+y)即为xy=sh 0xH-h,y0 求解： Min z=(x^2+s^2)^(1/2)+(y^2+h^2)^(1/2) 代码： function f=fun(x) f=(x(1)+9)^(1/2)+(x(2)+25)^(1/2) function [c,ceq]=confun(x) c=[] ceq=[x(1)*x(2)-15] [x,fval]=fmincon(@fun,