資源描述:
《書中matlab源程序》由會(huì)員上傳分享,免費(fèi)在線閱讀,更多相關(guān)內(nèi)容在教育資源-天天文庫。
1、書中Matlab源程序第1章緒論【例1-1】有一名學(xué)生,期末有5門功課要考試,可用的復(fù)習(xí)時(shí)間有18小時(shí)。假定這五門課程分別是數(shù)學(xué)、英語、計(jì)算機(jī)基礎(chǔ)、畫法幾何和專業(yè)概論。如果不復(fù)習(xí)直接參加考試,這五門功課預(yù)期的考試成績分別為65分、60分、70分、60分和65分。復(fù)習(xí)以1小時(shí)為一單元,每增加1小時(shí)復(fù)習(xí)時(shí)間,各門功課考試成績就有可能提高,每復(fù)習(xí)1小時(shí)各門功課考試成績提高的分?jǐn)?shù)分別為3分、4分、5分、4分和6分。問如何安排各門功課的復(fù)習(xí)時(shí)間可使平均成績不低于80分,并且數(shù)學(xué)和英語成績分別不低于70分和75分。解:設(shè)分配在數(shù)學(xué)、英語、計(jì)算機(jī)基礎(chǔ)、畫法幾何和專業(yè)概論這五門功課
2、的復(fù)習(xí)時(shí)間分別為,則可列出如下的目標(biāo)函數(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)計(jì)算結(jié)果為:x=1.66673.75005.00000.00005.8333fval=16.250050【例1-2】某工廠要生產(chǎn)兩種規(guī)格的電冰箱,分別用Ⅰ和Ⅱ表示。生產(chǎn)電冰箱需要兩
3、種原材料A和B,另外需設(shè)備C。生產(chǎn)兩種電冰箱所需原材料、設(shè)備臺時(shí)、資源供給量及兩種產(chǎn)品可獲得的利潤如表1-1所示。問工廠應(yīng)分別生產(chǎn)Ⅰ、Ⅱ型電冰箱多臺,才能使工廠獲利最多?表1.1資源需求與限制資源ⅠⅡ資源限制設(shè)備111200臺時(shí)原料A211800千克原料B011000千克單位產(chǎn)品獲利220元250元求最大收益產(chǎn)品Ⅰ用原料限制800千克解:設(shè)生產(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ù)的曲線。解:應(yīng)用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)化模型,并求解。解:該繪制目標(biāo)函數(shù)曲面、約束函數(shù)曲線及求解程序如下:(1)繪制目標(biāo)函數(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