数学建模论文-出版社的资源配置问题

附录附录1：%lin.m利用MATLAB自带的优化工具箱求解clc%各科目各课程的均价jsj=[25.825.528.026.024.725.627.022.925.924.5];jg=[26.427.324.927.523.523.525.732.931.535.3];sx=[21.020.224.819.618.623.313.118.422.525.7];yy=[34.418.733.020.627.921.411.431.323.532.3];lk=[14.718.826.616.713.414.817.524.2];jn=[22.532.320.821.623.035.4];hg=[20.023.625.628.018.926.7];dl=[21.532.424.023.818.222.7];hj=[37.522.220.722.724.332.2];%各课程的书号个数百分比w1=[0.1590909090.1619318180.0397727270.0482954550.0482954550.1534090910.1221590910.0482954550.1477272730.071022727];w2=[0.1356783920.080402010.0904522610.0703517590.0904522610.0904522610.0954773870.160804020.0954773870.090452261];w3=[0.0443548390.1733870970.2392473120.0362903230.1209677420.1438172040.0591397850.0241935480.1021505380.056451613];w4=[0.3559718970.0444964870.0163934430.189695550.060889930.1217798590.0491803280.0444964870.0725995320.044496487];w5=[0.0657894740.1491228070.1052631580.1578947370.0745614040.1403508770.1578947370.149122807];w6=[0.2964824120.0603015080.1105527640.276381910.1005025130.155778894];w7=[0.0761904760.1523809520.20.1238095240.1238095240.323809524];w8=[0.1750.1833333330.1916666670.2166666670.1166666670.116666667];w9=[0.1746031750.2380952380.2222222220.2142857140.0873015870.063492063];c1=[249.4166666363.8375118.3333333238.4375162.4166667276.5818181503.548611183.08333333387.975414.9375];c2=[1425.51875.75386.66666661106.208333286.4166667544.9166667812.9375575.29166662589.93751356.0625];c3=[746.5416667847.37857148609.9033741297.51462.3578432617.573864628.4214286498.8125465.67604172863.553572];c4=[747.8024039525.45907.9166667268.5805555447.8583334274.112559.91666667632.3875347.4464286301.8333333];c5=[5297.56257532.6607148622.24606.10416716071.854167891.254104.6339291241.875];c6=[681.3461539266.875480.25816.1751204.6833341000.75];c7=[1038.75863.95406.975162.9166667782.0833333427.4583333];c8=[226.6666667644.875600.125810.5295.16666671543.666666];c9=[593.1875476.7916667629.1666666414.25718.6251460.75];m=0.7;ys=-[sum(jsj.*w1.*c1)sum(jg.*w2.*c2)sum(sx.*w3.*c3)sum(yy.*w4.*c4)sum(lk.*w5.*c5)sum(jn.*w6.*c6)sum(hg.*w7.*c7)sum(dl.*w8.*c8)sum(hj.*w9.*c9)];y1=m*ys;my=[3.253.243.243.243.23.23.33.273.2];%满意度pf=0.2*sin(pi/4*(my-2));y2=(1-m)*pf.*ys;y=y1+y2;Aeq=[111111111];beq=500;A=[];b=[];vlb=[7035145784038172121];vub=[7844155905042252727];j=[64697769733739433941147148153150146958777838544484250444042403839242520171922212327272721252726];qn=y*j;%01-05年的经济总值[x,fval]=linprog(y,A,b,Aeq,beq,vlb,vub)%x为06年书号分配方案，fval为相应的经济总值附录2：%gagh.m遗传算法主程序clc%定义遗传算法参数NIND=40;%个体数目（Numberofindividuals）MAXGEN=200;%最大遗传代数(Maximumnumber.ofgenerations)NVAR=9;%变量的维数PRECI=20;%变量的二进制位数(Precisionofvariables)GGAP=0.9;%代沟（Generationgap）trace=zeros(MAXGEN,2);%寻优结果的初始值%建立区域描述器(Buildfielddescriptor)vlb=0.5*[110662221187276304040];vub=[6944155795042172522];ub=[vlb;vub];FieldD=[REP([PRECI],[1,NVAR]);ub;REP([1;0;1;1],[1,NVAR])];Chrom=CRTBP(NIND,NVAR*PRECI);%创建初始种群（2进制）P=BS2RV(Chrom,FieldD);%初始种群的10进制化gen=0;%代计数器ObjV=objfunction(P);%计算初始种群个体的目标函数值whilegen<MAXGENf=RANKING(-ObjV);%分配适应度值（AssignfitnessvaluesFitnV=fit(f,P);SelCh=SELECT('SUS',Chrom,FitnV,GGAP);%选择SelCh=RECOMBIN('XOVSP',SelCh,0.7);%重组SelCh=MUT(SelCh);%变异P=BS2RV(SelCh,FieldD);ObjVSel=objfunction(P);[ChromObjV]=REINS(Chrom,SelCh,1,1,ObjV,ObjVSel);%重插入子代的新种群P=BS2RV(Chrom,FieldD);gen=gen+1;%代计数器增加trace(gen,1)=max(ObjV);%遗传算法性能跟踪trace(gen,2)=sum(ObjV)/length(ObjV);endfori=1:NINDifObjV(i)==max(ObjV)value=P(i,:);breakendendobjvalue=max(ObjV);plot(trace(:,1));holdon;plot(trace(:,2),'-.');grid;legend('解的变化','种群均值的变化')functionz=fit(x,y)%适应度函数%x--40组按目标值的大小预分配的适应度%y--40组自变量向量（y1,y2...y9）fori=1:40ifabs(sum(y(i,:))-500)>5x(i)=0.001*x(i);endendz=x;functiony=objfunction(a)%遗传算法目标函数程序%各科目各课程的均价jsj=[25.825.528.026.024.725.627.022.925.924.5];jg=[26.427.324.927.523.523.525.732.931.535.3];sx=[21.020.224.819.618.623.313.118.422.525.7];yy=[34.418.733.020.627.921.411.431.323.532.3];lk=[14.718.826.616.713.414.817.524.2];jn=[22.532.320.821.623.035.4];hg=[20.023.625.628.018.926.7];dl=[21.532.424.023.818.222.7];hj=[37.522.220.722.724.332.2];%各课程的书号个数百分比w1=[0.1590909090.1619318180.0397727270.0482954550.0482954550.1534090910.1221590910.0482954550.1477272730.071022727];w2=[0.1356783920.080402010.0904522610.0703517590.0904522610.0904522610.0954773870.160804020.0954773870.090452261];w3=[0.0443548390.1733870970.2392473120.0362903230.1209677420.1438172040.0591397850.0241935480.1021505380.056451613];w4=[0.3559718970.0444964870.0163934430.189695550.060889930.1217798590.0491803280.0444964870.0725995320.044496487];w5=[0.0657894740.1491228070.1052631580.1578947370.0745614040.1403508770.1578947370.149122807];w6=[0.2964824120.0603015080.1105527640.276381910.1005025130.155778894];w7=[0.0761904760.1523809520.20.1238095240.1238095240.323809524];w8=[0.1750.1833333330.1916666670.2166666670.1166666670.116666667];w9=[0.1746031750.2380952380.2222222220.2142857140.0873015870.063492063];c1=[249.4166666363.8375118.3333333238.4375162.4166667276.5818181503.548611183.08333333387.975414.9375];c2=[1425.51875.75386.66666661106.208333286.4166667544.9166667812.9375575.29166662589.93751356.0625];c3=[746.5416667847.37857148609.9033741297.51462.3578432617.573864628.4214286498.8125465.67604172863.553572];c4=[747.8024039525.45907.9166667268.5805555447.8583334274.112559.91666667632.3875347.4464286301.8333333];c5=[5297.56257532.6607148622.24606.10416716071.854167891.254104.6339291241.875];c6=[681.3461539266.875480.25816.1751204.6833341000.75];c7=[1038.75863.95406.975162.9166667782.0833333427.4583333];c8=[226.6666667644.875600.125810.5295.1666