資源描述：
《書中matlab源程序》由會員上傳分享，免費在線閱讀，更多相關內(nèi)容在教育資源-天天文庫。
1、書中Matlab源程序第1章緒論【例1-1】有一名學生，期末有5門功課要考試，可用的復習時間有18小時。假定這五門課程分別是數(shù)學、英語、計算機基礎、畫法幾何和專業(yè)概論。如果不復習直接參加考試，這五門功課預期的考試成績分別為65分、60分、70分、60分和65分。復習以1小時為一單元，每增加1小時復習時間，各門功課考試成績就有可能提高，每復習1小時各門功課考試成績提高的分數(shù)分別為3分、4分、5分、4分和6分。問如何安排各門功課的復習時間可使平均成績不低于80分，并且數(shù)學和英語成績分別不低于70分和75分。解：設分配在數(shù)學、英語、計算機基礎、畫法幾何和專業(yè)概論這五門功課
2、的復習時間分別為，則可列出如下的目標函數(shù)和限制條件為：本例具體程序如下：%li_1_1f=[11111];A=[11111;-3-4-5-4-6;-30000;0-4000;30000;04000;00500;00040;00006];b=[18;-80;-5;-15;35;40;30;40;35];lb=zeros(6,1)[x,fval]=linprog(f,A,b,[],[],lb)計算結果為：x=1.66673.75005.00000.00005.8333fval=16.250050【例1-2】某工廠要生產(chǎn)兩種規(guī)格的電冰箱，分別用Ⅰ和Ⅱ表示。生產(chǎn)電冰箱需要兩
3、種原材料A和B，另外需設備C。生產(chǎn)兩種電冰箱所需原材料、設備臺時、資源供給量及兩種產(chǎn)品可獲得的利潤如表1-1所示。問工廠應分別生產(chǎn)Ⅰ、Ⅱ型電冰箱多臺，才能使工廠獲利最多？表1.1資源需求與限制資源ⅠⅡ資源限制設備111200臺時原料A211800千克原料B011000千克單位產(chǎn)品獲利220元250元求最大收益產(chǎn)品Ⅰ用原料限制800千克解：設生產(chǎn)Ⅰ、Ⅱ兩種產(chǎn)品的數(shù)量分別為。則可獲得的最大收益為Matlab求解程序如下：%li_1_2clc;closeall;f=-[220250];A=[11;21;10;01];b=[1200;1800;800;1000];xl=[0
4、0];[x,fval]=linprog(f,A,b,[],[],xl)x1=[0:1800];x2=[0:2000];[xm1,xm2]=meshgrid(x1,x2);x21=1200-x1;x22=1800-2*x1;x23=(-fval-220*x1)/250;50plot(x1,x21,x1,x22,[0:1:1000],1000,800,[0:1:1500],x1,x23,'r')axis([0,1400,0,2000])xlabel('x1');ylabel('x2');holdonz=200*xm1+250*xm2;[C,h]=contour(xm1,
5、xm2,z);text_handle=clabel(C,h);set(text_handle,'BackgroundColor',[11.6],'Edgecolor',[.7.7.7]);holdoff【例1-3】繪制下面函數(shù)的曲線。解：應用plot()函數(shù)繪制該函數(shù)曲線的程序如下：%li_1_3f=inline('2*sin(x)+log(x)','x')x=linspace(0.1,2*pi,15);y=feval(f,x);plot(x,y,'-rs','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor',
6、'g','MarkerSize',10)xlabel('0.1leqThetaleq2pi')ylabel('2sin(Theta)+ln(Theta)');title('Plotof2sin(Theta)+ln(Theta)')text(pi/4,sin(-pi/4),'leftarrow2sin(Theta)+ln(Theta)','HorizontalAlignment','left')legend('-')gridon【例1-4】用圖形表示如下優(yōu)化模型，并求解。解：該繪制目標函數(shù)曲面、約束函數(shù)曲線及求解程序如下：（1）繪制目標函數(shù)曲面的
7、程序%li_1_4_1functionli_1_4_1()clc;clearall;closeall;n=20;50x1=linspace(0,2,n);x2=linspace(0,6,n);[xm1,xm2]=meshgrid(x1,x2);fori=1:nforj=1:nxx=[xm1(i,j),xm2(i,j)];zm(i,j)=fun_obj(xx);endendfigure(1)meshc(xm1,xm2,zm)xlabel('x1');ylabel('x2');zlabel('zm')（2）繪制約束函數(shù)曲線及求解的程序%li_1_4_2functio
