數(shù)學建模實驗答案_建立數(shù)學模型

1、實驗01 建立數(shù)學模型（4學時）（第1章 建立數(shù)學模型）教材中給出原始數(shù)據(jù)，結合模型，得到結果。但如何求得結果這一過程沒有給出，實際上要用MATLAB軟件編寫程序來求得，這應該交給實驗課來完成?？紤]到上學期同學們剛學習MATLAB語言，編程能力不強，所以有關的程序給出來供同學們進行驗證。要求同學們要讀懂程序。1.（求解，編程）如何施救藥物中毒p1011人體胃腸道和血液系統(tǒng)中的藥量隨時間變化的規(guī)律（模型）：其中，x(t)為t時刻胃腸道中的藥量，y(t)為t時刻血液系統(tǒng)中的藥量，t=0為服藥時刻。1.1（求解）模型求解p1011要求： 用MATLAB求解微分方程函數(shù)dsolve求解該微分方程（符號

2、運算）。 用MATLAB的化簡函數(shù)simplify化簡所得結果。 結果與教材P11上的內(nèi)容比較。提示：dsolve和simplify的用法可用help查詢。建議在命令窗口中操作。 求解的語句及運行結果：>> x,y= dsolve('Dx=-a*x','Dy=a*x-b*y','x(0)=1100','y(0)=0');>> disp(x,y) 1100*exp(-a*t), exp(-a*t)*exp(-b*t)*(1100*a*exp(a*t)/(a - b) - (1100*a*exp(b*t)/(a

3、 - b)>> disp(simplify(x,y); 1100*exp(-a*t), (1100*a*exp(-t*(a + b)*(exp(a*t) - exp(b*t)/(a - b)1.2（編程）結果分析p11已知=0.1386, =0.1155，將上題中得到x(t)和y(t)兩條曲線畫在同一個圖形窗口內(nèi)。參考圖形如下。提示：MATLAB命令plot, fplot, hold on/off, grid on/off, xlabel, ylabel, text。 編寫的程序和運行結果：程序1：用plotclc;clear;a=0.1386; b=0.1155;t=0:0.01

4、:25;x=1100./exp(a*t);y=-(1100*a*(1./exp(a*t) -1./exp(b*t)/(a - b);plot(t,x,t,y);grid on;xlabel('itt /h'); ylabel('itx,ity /mg');text(2,1100/exp(a*2),' itx(itt)');text(3,-(1100*a*(1/exp(a*3) - 1/exp(b*3)/(a - b),' ity(itt)');程序2：用fplot和匿名函數(shù)clc;clear;a=0.1386; b=0.1155;

5、fplot(t)1100/exp(a*t),-(1100*a*(1/exp(a*t) - 1/exp(b*t)/(a - b),0 25);grid on;xlabel('itt /h'); ylabel('itx,ity /mg');text(2,1100/exp(a*2),' itx(itt)');text(3,-(1100*a*(1/exp(a*3) - 1/exp(b*3)/(a - b),' ity(itt)');2.（編程，驗證）商人們怎樣安全過河p89三名商人各帶一個隨從乘船渡河，一只小船只能容納二人，由他們自己劃行

6、。隨從們密約，在河的任一岸，一旦隨從的人數(shù)比商人多，就殺人越貨。但是如何乘船的大權掌握在商人們手中。商人們怎樣才能安全渡河呢？模型構成決策： 每一步（此岸到彼岸或彼岸到此岸）船上的人員。要求：在安全的前提下（兩岸的隨從數(shù)不比商人多），經(jīng)有限步使全體人員過河。xk第k次渡河前此岸的商人數(shù)yk第k次渡河前此岸的隨從數(shù) xk , yk=0,1,2,3; k=1,2,¼過程的狀態(tài)sk=(xk , yk)允許狀態(tài)集合S=(x, y)| x=0, y=0,1,2,3; x=3, y=0,1,2,3; x=y=1,2uk第k次渡船上的商人數(shù)vk第k次渡船上的隨從數(shù) uk , vk=0,1,2; k

7、=1,2,¼決策dk=(uk , vk)允許決策集合D=(u , v)| u+v =1, 2狀態(tài)轉移律sk+1=sk+(-1)kdk多步?jīng)Q策問題求dkÎD(k=1, 2, ¼, n), 使skÎS, 并按轉移律由 s1=(3,3) 到達sn+1=(0,0)。2.1（編程）求允許決策集合D和允許狀態(tài)集合SD是2行多列矩陣，每一列是一個決策。S是2行多列矩陣，每一列是一種狀態(tài)。要求： 編寫程序求D和S，并輸出。 S的第一列是3,3'，最后一列是0,0 '。 編寫的程序和運行結果：程序：clear; clc;%求允許決策集合D（2×n

8、1，n1種決策）D=;for u=0:2 for v=0:2 if u+v=1|u+v=2 D=D,u;v; end endend%求允許狀態(tài)集合S（2×n2，n2種狀態(tài)）S=;for x=3:-1:0 for y=3:-1:0 if x=0|x=3|x=y S=S,x;y; end endend%首列狀態(tài) ( 商人數(shù), 仆從數(shù))' = ( 3, 3 )'，末列為( 0, 0 )'D, S運行結果：2.2（驗證）求動態(tài)允許狀態(tài)集合SS和狀態(tài)轉移矩陣A上面允許狀態(tài)集合S沒有指明當時船是在此岸還是在彼岸，應該將S中的每一種狀態(tài)再分為兩種狀態(tài)，需增加一個元素（值為-

9、1或1）放在第三行。定義動態(tài)允許狀態(tài)集合SS=(x, y, z)'| x=0, y=0,1,2,3; x=3, y=0,1,2,3; x=y=1,2;z=-1,1(x, y, -1)，表示從此岸渡河前此岸的允許狀態(tài)(x, y)。(x, y, 1)，表示從彼岸渡河前此岸的允許狀態(tài)(x, y)。SS是三行多列矩陣，每一列表示一種狀態(tài)，列下標為其編號。定義狀態(tài)轉移矩陣A，其中，A(i, j)=1表示D中存在決策使狀態(tài)i轉到j，否則A(i, j)=0。程序如下（輸入時，不必把注釋也輸入）：程序運行結果（參考）：要求： 將程序接在上題的程序之后（去掉最后多余的輸出語句），程序最后給出顯示SS和A

10、的語句。 運行程序，輸出SS和A。對照參考答案，如數(shù)值不一致，請檢查程序。MATLAB函數(shù)all的用法見提示。 運行的完整程序和運行結果：程序：clear; clc;%求允許決策集合D（2×n1，n1種決策）D=;for u=0:2 for v=0:2 if u+v=1|u+v=2 D=D,u;v; end endend%求允許狀態(tài)集合S（2×n2，n2種狀態(tài)）S=;for x=3:-1:0 for y=3:-1:0 if x=0|x=3|x=y S=S,x;y; end endend%首列狀態(tài) ( 商人數(shù), 仆從數(shù))' = ( 3, 3 )'，末列為( 0

11、, 0 )'%動態(tài)允許狀態(tài)集合SS（3×n3，n3種狀態(tài)）%-1，從此岸渡河前此岸的允許狀態(tài)%1，從彼岸渡河前此岸的允許狀態(tài)SS=S;-ones(1,size(S,2),S;ones(1,size(S,2);SSnum=size(SS,2);%狀態(tài)總數(shù)，SS中的列下標對應狀態(tài)的編號%SS(:,1)=3,3,-1'（起點），SS(:,end)=0,0,1'（終點）%狀態(tài)轉移矩陣A，A(i,j)=1表示存在決策使狀態(tài)i轉到j，其它為0A=zeros(SSnum);for i=1:SSnum for j=1:SSnum for d=D%順序取D的每一列給d s=SS

12、(1:2,i)+SS(3,i)*d;-SS(3,i); if all(s=SS(:,j)%所有元素不為0時為真 A(i,j)=1; break; end end endendSS, A運行結果：2.3（驗證）給出一個商人們安全過河的方案程序如下（輸入時，不必把注釋也輸入）：程序運行結果：要求： 讀懂以上程序，掌握all和any函數(shù)的應用。 將程序接在之前的程序后（去掉最后多余的輸出語句），運行程序。對照答案，如不一致，請檢查程序。MATLAB函數(shù)any的用法見提示。 給出運行的完整程序和運行結果：程序：clear; clc;%求允許決策集合D（2×n1，n1種決策）D=;for u=

13、0:2 for v=0:2 if u+v=1|u+v=2 D=D,u;v; end endend%求允許狀態(tài)集合S（2×n2，n2種狀態(tài)）S=;for x=3:-1:0 for y=3:-1:0 if x=0|x=3|x=y S=S,x;y; end endend%首列狀態(tài) ( 商人數(shù), 仆從數(shù))' = ( 3, 3 )'，末列為( 0, 0 )'%動態(tài)允許狀態(tài)集合SS（3×n3，n3種狀態(tài)）%-1，從此岸渡河前此岸的允許狀態(tài)%1，從彼岸渡河前此岸的允許狀態(tài)SS=S;-ones(1,size(S,2),S;ones(1,size(S,2);SSnum

14、=size(SS,2);%狀態(tài)總數(shù)，SS中的列下標對應狀態(tài)的編號%SS(:,1)=3,3,-1'（起點），SS(:,end)=0,0,1'（終點）%狀態(tài)轉移矩陣A，A(i,j)=1表示存在決策使狀態(tài)i轉到j，其它為0A=zeros(SSnum);for i=1:SSnum for j=1:SSnum for d=D%順序取D的每一列給d s=SS(1:2,i)+SS(3,i)*d;-SS(3,i); if all(s=SS(:,j)%所有元素不為0時為真 A(i,j)=1; break; end end endend%SK為多行SSnum列的矩陣，第i列對應SS的第i列的狀態(tài)。

15、%第i行表示第i次渡河前有那些狀態(tài)，1表示有對應下標列號的狀態(tài)，0則無。%第i+1行是第i行狀態(tài)轉移過來的狀態(tài)%當?shù)赟Snum狀態(tài)出現(xiàn)時停止，表示有解。%或者新行的元素全0時停止，表示無解。%新行不允許出現(xiàn)之前出現(xiàn)過的狀態(tài)。s=1,zeros(1,SSnum-1);%從狀態(tài)3,3,-1'開始SK=s;while any(s) && s(SSnum)=1 k=find(s);%求最后一行為1的列下標 s=any(A(k,:),1);%得到新行，元素值為1或0 s=s&any(SK,1);%去掉新行中之前出現(xiàn)的狀態(tài) SK=SK;s;%新行加入，作為SK的末行end%

16、多步?jīng)Q策，直到找到目標。有難度！if any(s)=0 %新行全0disp('無解！'); return;end%最后一行的第SSnum狀態(tài)(j)開始往回找轉換到它的狀態(tài)(i)，直到狀態(tài)1。%得到一個狀態(tài)號的順序sk（行向量），為一個解。sk=zeros(1,size(SK,1);sk(1)=1; sk(1,end)=SSnum; j=SSnum;for k=size(SK,1)-1:-1:2 for i=find(SK(k,:) if A(i,j)=1 sk(k)=i; j=i; break; end endenddisp(1:length(sk);SS(1:2,sk);%輸

17、出結果給出程序的運行結果：3.（求解）商人們怎樣安全過河（修改）p9對第2題的問題改動，用類似的方法求解。3.1 在第2題中修改商人數(shù)和隨從數(shù)有四名商人各帶一個隨從，其它同第2題。修改第2題中的程序求解。 修改的程序部分和完整程序的運行結果（安全過河的方案）：%求允許狀態(tài)集合S（2×n2，n2種狀態(tài)）S=;for x=4:-1:0 for y=4:-1:0 if x=0|x=4|x=y S=S,x;y; end endend%首列狀態(tài) (商人數(shù), 仆從數(shù))' = ( 4, 4 )'，末列為( 0, 0 )'3.2 在3.1題中修改船容納的人數(shù)船能容納3人。修改

18、3.1題中的程序求解。 修改的程序部分和完整程序的運行結果（安全過河的方案）：%求允許決策集合D（2×n1，n1種決策）D=;for u=0:3 for v=0:3 if u+v>=1 && u+v<=3 D=D,u;v; end endend4.（編程）安全過河問題（人、貓、雞、米）p21習題5人帶著貓、雞、米過河，除需要人劃船之外，船至多能載貓、雞、米三者之一，而當人不在場時貓要吃雞、雞要吃米。試設計一個安全過河方案。模仿商人過河問題的程序編寫程序解決本問題。 程序：clear; clc;%求允許決策集合D（2×n1，n1種決策）D=1,0,

19、0,0',1,1,0,0',1,0,1,0',1,0,0,1'%求允許狀態(tài)集合S（2×n2，n2種狀態(tài)）S=;j=1;for a=23+22+21+1:-1:0 s=fix(a/23);mod(fix(a/22),2);mod(fix(a/2),2);mod(a,2);%取二進制數(shù)字 if s(1)&&(s(2)&&s(4)|s(3)|s(1)&&(s(2)&&s(4)|s(3) S(:,j)=s; j=j+1; end%s=人，貓，雞，米end%首列狀態(tài)=(1,1,1,1)'，末

20、列為(0,0,0,0)'%動態(tài)允許狀態(tài)集合SS（3×n3，n3種狀態(tài)）%-1，從此岸渡河前此岸的允許狀態(tài)%1，從彼岸渡河前此岸的允許狀態(tài)SS=S;-ones(1,size(S,2),S;ones(1,size(S,2);SSnum=size(SS,2);%狀態(tài)總數(shù)，SS中的列下標對應狀態(tài)的編號%SS(:,1)=1,1,1,1,-1'（起點）,SS(:,end)=0,0,0,0,1'（終點）%狀態(tài)轉移矩陣A，A(i,j)=1表示存在決策使狀態(tài)i轉到j，其它為0A=zeros(SSnum);for i=1:SSnum for j=1:SSnum for d=D%順

21、序取D的每一列給d s=SS(1:4,i)+SS(5,i)*d;-SS(5,i); if all(s=SS(:,j)%所有元素不為0時為真 A(i,j)=1; break; end end endend%SK為多行SSnum列的矩陣，第i列對應SS的第i列的狀態(tài)。%第i行表示第i次渡河前有那些狀態(tài)，1表示有對應下標列號的狀態(tài)，0則無。%第i+1行是第i行狀態(tài)轉移過來的狀態(tài)%當?shù)赟Snum狀態(tài)出現(xiàn)時停止，表示有解。%或者新行的元素全0時停止，表示無解。%新行不允許出現(xiàn)之前出現(xiàn)過的狀態(tài)。s=1,zeros(1,SSnum-1);%從狀態(tài)1,1,1,1,-1'開始SK=s;while any

22、(s) && all(SSnum-find(s) k=find(s);%求最后一行為1的列下標 s=any(A(k,:),1);%得到新行，元素值為1或0 s=s&any(SK,1);%去掉新行中之前出現(xiàn)的狀態(tài) SK=SK;s;%新行加入，作為SK的末行endif any(s) %新行全0disp('無解！'); return;end%最后一行的第SSnum狀態(tài)(j)開始往回找轉換到它的狀態(tài)(i)，直到狀態(tài)1。%得到一個狀態(tài)號的順序sk（行向量），為一個解。sk=zeros(1,size(SK,1);sk(1)=1; sk(1,end)=SSnum; j

23、=SSnum;for k=size(SK,1)-1:-1:2 for i=find(SK(k,:) if A(i,j)=1 sk(k)=i; j=i; break; end endenddisp(1:length(sk);SS(1:4,sk);%輸出結果 程序的運行結果5.（選做，編程）商人們怎樣安全過河（自編）2.3中的程序有一定難度，所以給出了參考程序，初學MATLAB的同學可能想不到用其中的一些語句。第2題給出的是求解商人們安全過河問題的一種算法程序，還有其它的算法程序，比如回溯法、遞歸算法等。要求： 按照自己的理解和想法，改寫2.3中的程序。 或用別的算法編寫解決商人們安全過河的程序。

24、 程序和運行結果：解1（遞歸）程序function river()global mark D SHOW;%全局變量mark=2*ones(4);mark(1,:)=0; mark(4,:)=0;%mark(i,j)=0為允許狀態(tài)mark(2,2)=0; mark(3,3)=0;D=0 1; 0 2; 1 0; 1 1; 2 0;%允許決策集合mark(4,4)=2; SHOW=;fun(3,3,-1);disp(3,3;SHOW);returnfunction f=fun(s,t)global mark D SHOW;%全局變量if s=0f=1; return;endfor k=1:5s1=

25、s+t*D(k,:);i=s1(1)+1; j=s1(2)+1;if all(s1>=0) && all(s1<=3) && (mark(i,j)=0 | mark(i,j)=t)if mark(i,j)=0%標記狀態(tài)(i-1,j-1)的使用情況mark(i,j)=-t;elsemark(i,j)=2;end if fun(s1,-t)SHOW=s1;SHOW; f=1; return; endendendif s=3 disp('無解');endf=0;運行結果3 33 13 23 0 %船?？勘税叮劝阉须S從渡過彼岸3 11 1

26、2 20 2 %船?？勘税?，再把所有商人渡過彼岸0 3 %船停靠此岸，此岸只有隨從，彼岸只有商人0 10 20 0 %船?？勘税?，最后把所有隨從渡過彼岸三個階段：(3,3,船停此岸)(3,0,船停彼岸)(0,3,船停此岸)(0,0,船停彼岸)第一階段：所有隨從渡過彼岸（商人不動）；第二階段：使所有商人渡過彼岸，并將隨從渡回此岸；第三階段：所有隨從渡過彼岸（商人不動）。解2（借用棧）思路定義狀態(tài)：此岸的商人數(shù)，隨從數(shù)，小船停靠在此岸或彼岸。從初始狀態(tài)開始，找到一個決策得到下一個狀態(tài)，再從新狀態(tài)開始，找到一個決策得到下一個狀態(tài)，.。在這個過程中，出現(xiàn)過的狀態(tài)將不可再用。若從某個狀態(tài)開始，所有決策都

27、得不到新狀態(tài)，則返回上一狀態(tài)，并從下一個決策開始，試探是否能得到新的狀態(tài)。直到到達狀態(tài)(0,0)，或無解。算法1. 給出允許狀態(tài)集合；%用4階方陣mark標記，其中mark(i,j)=0表示(i-1,j-1)為允許狀態(tài)，mark(i,j)=2表示(i-1,j-1)為不可達狀態(tài)2. 給出允許決策集合；% D=0 1; 0 2; 1 0; 1 1; 2 03. 初始狀態(tài)進棧；%包括：當前狀態(tài)，出發(fā)標記（-1/1為從此岸/彼岸出發(fā)），選擇第幾個決策4. 初始狀態(tài)標記；%mark(i,j)=-1/1表示從此岸/彼岸出發(fā)轉移到狀態(tài)(i-1,j-1)已走過，mark(i,j)=2表示上述兩種情況都出現(xiàn)或不

28、可達狀態(tài)5. 循環(huán)直到?？栈虻竭_狀態(tài)(0,0) 5.1 出棧，將出棧的狀態(tài)作為當前狀態(tài)，并選擇下一個決策；%原決策失敗，回溯 5.2 循環(huán)直到?jīng)Q策都用過或到達最終狀態(tài) 5.2.1 由當前狀態(tài)、出發(fā)標記和決策得到轉移到的新狀態(tài)； 5.2.2 若新狀態(tài)是允許狀態(tài)且沒走過，則執(zhí)行5.2.2.1；否則選擇下一個決策 5.2.2.1 標記新狀態(tài)的訪問信息；%mark(i,j)=-1/1/2 5.2.2.2 當前狀態(tài)重新進棧（改變了決策）；%存儲決策過程 5.2.2.4 將新狀態(tài)作為當前狀態(tài)，置出發(fā)標記，選擇第1個決策；6. 若有解，則輸出決策過程；否則輸出“無解”；程序clear; clc;mark=2

29、*ones(4);mark(1,:)=0; mark(4,:)=0;%mark(i,j)=0為允許狀態(tài)mark(2,2)=0; mark(3,3)=0;D=0 1; 0 2; 1 0; 1 1; 2 0;%允許決策集合sp=1; st(sp,:)=3,3,-1,0; %進棧，(3,3)為當前狀態(tài)，-1表示小船?？看税?，決策0（未做決策）mark(4,4)=2;%標記為不可達狀態(tài)s0=3,3;while sp>0&& any(s0>0)%棧未空且s0中至少有一個元素>0s0=st(sp,1,2); ss=st(sp,3);%當前狀態(tài)，ss=-1/1小船?？看税?

30、彼岸d=st(sp,4);%之前用過的決策，該決策失敗sp=sp-1;%出棧 d=d+1;%試探下一個決策 while d<6 && any(s0>0)s1=s0+ss*D(d,:);%下一個狀態(tài)，不包括小船?？繝顟B(tài)i=s1(1)+1;j=s1(2)+1;%狀態(tài)轉換為下標if all(s1>=0) && all(s1<=3) && (mark(i,j)=0 | mark(i,j)=-ss)if mark(i,j)=0%標記狀態(tài)(i-1,j-1)的使用情況mark(i,j)=ss;elsemark(i,j)=2;endsp=

31、sp+1;%進棧st(sp,:)=s0,ss,d;%當前狀態(tài)重新進棧，改變了決策s0=s1; ss=-ss; d=1;%下一個狀態(tài)作當前狀態(tài)，選擇第1個決策else%決策d失敗d=d+1;endendendif s0=0disp(1:sp+1',st(1:sp,1:2);s0);%有解elsedisp('無解'); end運行結果 1 3 3 2 3 1 3 3 2 4 3 0 5 3 1 6 1 1 7 2 2 8 0 2 9 0 3 10 0 1 11 0 2 12 0 0解3clear; clc;%求允許狀態(tài)集合SS=;for x=0:3 for y=0:3 if

32、 (x=0&&(y>=0&&y<=3). |(x=3&&(y>=0&&y<=3). |(x=1&&y=1|x=2&&y=2) S=S,x;y; end endendSA=S;%從此岸(A)到彼岸(B)渡河前此岸的允許狀態(tài)集合SB=S;%從彼岸(B)到此岸(A)渡河前此岸的允許狀態(tài)集合%求允許決策集合DD=;for u=0:2 for v=0:2 if u+v=1|u+v=2 D=D,u;v; end endend k=1;%第1次渡河s=3;3;%第1次渡河前的狀態(tài)Snum(

33、1)=1;%第1次渡河前的狀態(tài)總數(shù)SA(:,end)=; SB=SA;%去掉狀態(tài)(3,3)while Snum(k)>0 Snum(k+1)=0;%第k+1次渡河前的狀態(tài)總數(shù)初始化 k1=size(s,2)-Snum(k)+1;k2=size(s,2); for i=k1:k2 xk=s(1,i); yk=s(2,i); for d=D xk1=xk+(-1)k*d(1); yk1=yk+(-1)k*d(2); if (-1)k=1 for j=1:size(SB,2) if xk1=SB(1,j)&&yk1=SB(2,j) s=s,SB(:,j); SB(:,j)=;

34、Snum(k+1)=Snum(k+1)+1; break; end end else for j=1:size(SA,2) if xk1=SA(1,j)&&yk1=SA(2,j) s=s,SA(:,j); SA(:,j)=; Snum(k+1)=Snum(k+1)+1; break; end end end end end k=k+1;endSnums(:,1:Snum(1)m=Snum(1);for k=1:11 disp(s(:,m+1:m+Snum(k+1); m=m+Snum(k+1);end解4i=1;S=;D=;S(:,i)=3;3;%開始S(:,1)=3,3

35、9;while S(1,i)=0|S(2,i)=0%當s(:,i)=0,0'完成 D(:,i)=round(2*rand(2,1);%隨機取0,1,2給決策D(:,i) d=D(:,i); if d(1)+d(2)=1&&d(1)+d(2)=2%成立則不是允許決策 continue;%重新選取決策 end S(:,i+1)=S(:,i)+(-1)i*D(:,i);%得到下一個狀態(tài)S(:,i+1) s=S(:,i+1); if (s(1)=1&&s(2)=1)|(s(1)=2&&s(2)=2)|(s(1)>3|s(1)<0)|(

36、s(2)>3|s(2)<0)%成立則不是允許狀態(tài) continue;%重新選取決策 end i=i+1; if i=105%等于此循環(huán)次數(shù)判無解 disp('無解');return; endendfor i=1:size(S,2)%添加狀態(tài)的第3行元素，-1為船在此岸，1時在彼岸 SS(:,i)=S(:,i);(-1)i;end%區(qū)分過程中的不同狀態(tài)for i=1:size(SS,2)-1%去掉狀態(tài)轉移過程中多余的狀態(tài) a=0; for j=i+1:size(SS,2)%找出與第i狀態(tài)相同的最后一個狀態(tài) if any(SS(:,i)-SS(:,j)=0 a=j; e

37、nd end if i<a%刪除第i到a-1列，多余的狀態(tài)轉移 SS(:,i:a-1)=; endenddisp(1:size(SS,2);SS)%該程序運行效率含隨機性，可能誤判為無解！解5syms k;%定義一個關于k的，儲存各個狀態(tài)之間轉化的符號矩陣as=3,2,1,0;ss=;%能夠存在的狀態(tài)集合for i=1:length(s)%計算ss for j=1:length(s) if i=length(s)|i=1|s(i)=s(j) sss=s(i),s(j); ss=ss;sss; end endendignore,p=sort(sum(ss,2),'descend&#

38、39;);%按總數(shù)排序,便于計算狀態(tài)矩陣ss=ss(p,:);a=sym(zeros(length(ss),length(ss);for i=1:length(ss) for j=i+1:length(ss) m=abs(sum(ss(i,:)-sum(ss(j,:); if m=1|m=2%兩個狀態(tài)的總人數(shù)之差只能是1或者2 if ss(i,:)>=ss(j,:)%保證總人數(shù)多的狀態(tài)每個分量都要大 a(i,j)=2+(-1)k; a(j,i)=(-1)k; end end endenda(1,2)=0;%當?shù)谝淮斡?3,3)變到(3,2)時，下一次必為(3,2) 變到(3,3)，避免無

39、用功k=1;dd=1;s0=ss(1,:);%初始狀態(tài)S=ss(1,:);%儲存總的狀態(tài)轉移方案D=;%每一步的決策集合while any(s0)&k<=100 b=eval(a);%確定此時的a矩陣 for i=1:length(ss)%判斷s0處于s中的第幾行 if s0(1,:)=ss(i,:) break end end m=find(b(i,:)=1); %確定此時的方案集合 if k>1&length(m)>=2 m(find(m=dd(k-1)=;%避免走回頭路 end ignore,p=sort(rand(1,length(m);%在剩下的狀態(tài)

40、中隨機找一個可行狀態(tài) s0=ss(m(p(1),:); k=k+1; dd=dd;m(p(1);S=S;s0;D=D;abs(S(k,:)-S(k-1,:);endSD40附1：實驗提示第2.2題MATLAB函數(shù)all的用法：對于向量v，all(v)當v的所有元素都非0時，返回true（值為1）；否則返回false（值為0）。（類似v的所有元素“與”）對于矩陣X，all(X)返回一個行向量，每一個元素值為all作用于X對應列的結果。all(X, 1)作用于列（同all(X)），all(X, 2)作用于行。第2.3題MATLAB函數(shù)any的用法：對于向量v，any(v)當v中有元素非0時，返回t

41、rue（值為1）；否則返回false（值為0）。（類似v的所有元素“或”）對于矩陣X，any(X)返回一個行向量，每一個元素值為any作用于X對應列的結果。any(X, 1)作用于列（同any(X)），any(X, 2)作用于行。第3題要給出安全過河的方案。第4題在設計允許狀態(tài)集合和允許決策集合時，要把“人”考慮進去附2：第1章 建立數(shù)學模型數(shù)學模型出現(xiàn)在生產(chǎn)、工作和社會活動中：電氣工程師必須建立所要控制的生產(chǎn)過程的數(shù)學模型，用這個模型對控制裝置作出相應的設計和計算，才能實現(xiàn)有效的過程控制。氣象工作者為了得到準確的天氣預報，一刻也離不開根據(jù)氣象站、氣象衛(wèi)星匯集的氣壓、雨量、風速等資料建立的數(shù)學

42、模型。生理醫(yī)學專家有了藥物濃度注人體內(nèi)隨時間和空間變化的數(shù)學模型，就可以分析藥物的療效，有效地指導臨床用藥。城市規(guī)劃工作者需要建立一個包括人口、經(jīng)濟、交通、環(huán)境等大系統(tǒng)的數(shù)學模型，為領導層對城市發(fā)展規(guī)劃的決策提供科學根據(jù)。廠長經(jīng)理們要是能夠根據(jù)產(chǎn)品的需求狀況、生產(chǎn)條件和成本、貯存費用等信息，籌劃出一個合理安排生產(chǎn)和銷售的數(shù)學模型，一定可以獲得更大的經(jīng)濟效益。對于科學技術人員和應用數(shù)學工作者來說，建立數(shù)學模型是溝通實際問題與數(shù)學工具之間聯(lián)系的一座橋梁。1.1 從現(xiàn)實對象到數(shù)學模型原型和模型原型 指人們在現(xiàn)實世界里關心、研究或者從事生產(chǎn)、管理的實際對象。在科技領域通常使用系統(tǒng)、過程等詞匯：如機械系

43、統(tǒng)、電力系統(tǒng)、生態(tài)系統(tǒng)、生命系統(tǒng)、社會經(jīng)濟系統(tǒng)；又如鋼鐵冶煉過程、導彈飛行過程、化學反應過程、污染擴散過程、生產(chǎn)銷售過程、計劃決策過程等。本書所述的現(xiàn)實對象、研究對象、實際問題等均指原型。模型 則指為了某個特定目的將原型的某一部分信息簡縮、提煉而構造的原型替代物。模型的基本特征是由構造模型的目的決定的：模型不是原型原封不動的復制品，原型有各個方面和各種層次的特征，而模型只要求反映與某種目的有關的那些方面和層次。一個原型，為了不同的目的可以有許多不同的模型：放在展廳里的飛機模型應該在外形上逼真，但是不一定會飛。參加航模競賽的模型飛機要具有良好的飛行性能，在外觀上不必苛求。在飛機設計、試制過程中用

44、到的數(shù)學模型和計算機模擬，只要求在數(shù)量規(guī)律上真實反映飛機的飛行動態(tài)特性，毫不涉及飛機的實體。用模型替代原型的方式分類，模型可分為物質模型（形象模型）：直觀模型、物理模型等。理想模型（抽象模型）：思維模型、符號模型、數(shù)學模型等。直觀模型：那些供展覽用的實物模型，以及玩具、照片等，通常把原型的尺寸按比例縮小或放大，主要追求外觀上的逼真。物理模型：科技工作者為一定目的根據(jù)相似原理構造的模型，它不僅可以顯示原型的外形或某些特征，而且可以用來進行模擬實驗，間接地研究原型的某些規(guī)律。如波浪水箱中的艦艇模型用來模擬波浪沖擊下艦艇的航行性能，風洞中的飛機模型用來試驗飛機在氣流中的空氣動力學特性。有些現(xiàn)象直接用

45、原型研究非常困難，更可借助于這類模型，如地震模擬裝置、核爆炸反應模擬設備等。思維模型：指通過人們對原型的反復認識，將獲取的知識以經(jīng)驗形式直接貯存于人腦中，從而可以根據(jù)思維或直覺作出相應的決策。如汽車司機對方向盤的操縱、一些技藝性較強的工種（如鉗工）的操作。通常說的某些領導者憑經(jīng)驗作決策也是如此。Ø 優(yōu)缺點：思維模型便于接受，也可以在一定條件下獲得滿意的結果，但是它往往帶有模糊性、片面性、主觀性、偶然性等缺點，難以對它的假設條件進行檢驗，且不便于人們的相互溝通。符號模型：是在一些約定或假設下借助于專門的符號、線條等，按一定形式組合起來描述原型。如地圖、電路圖、化學結構式等。具有簡明、方

46、便、目的性強及非量化等特點。數(shù)學模型：是由數(shù)字、字母或其它數(shù)學符號組成的，描述現(xiàn)實對象數(shù)量規(guī)律的數(shù)學公式、圖形或算法。什么是數(shù)學模型？航行問題：甲乙兩地相距750km，船從甲到乙順水航行需30h，從乙到甲逆水航行需50h，問船速、水速各若干？用x，y分別代表船速和水速，可以列出方程( x + y ) · 30 = 750，( x y ) · 50 = 750這組方程就是航行問題的數(shù)學模型。方程的解x = 20 km/h，y = 5 kn/h，最終給出了航行問題的答案。建立數(shù)學模型的基本內(nèi)容：根據(jù)建立數(shù)學模型的目的和問題的背景作出必要的簡化假設（航行中設船速和水速為常數(shù)）；用

47、字母表示待求的未知量（x，y代表船速和水速）；利用相應的物理或其它規(guī)律（勻速運動的距離等于速度乘以時間)，列出數(shù)學式子（二元一次方程）；求出數(shù)學上的解答（x = 20，y = 5）；用這個答案解釋原問題（船速和水遞分別為20 km/h和5 km/h）；最后用實際現(xiàn)象驗證上述結果。數(shù)學模型可以描述為，對于現(xiàn)實世界的一個特定對象，為了一個特定目的，根據(jù)特有的內(nèi)在規(guī)律，做出一些必要的簡化假設，運用適當?shù)臄?shù)學工具，得到的一個數(shù)學結構。本書的重點不在于介紹現(xiàn)實對象的數(shù)學模型是什么樣子，而是要討論建立數(shù)學模型的全過程。建立數(shù)學模型簡稱數(shù)學建?；蚪！?shù)學模擬：運用數(shù)字式計算機的計算機模擬。它根據(jù)實際系統(tǒng)或

48、過程的特性，按照一定的數(shù)學規(guī)律用計算機程序語言模擬實際運行狀況，并依據(jù)大量模擬結果對系統(tǒng)或過程進行定量分析。例如，通過各種工件在不同機器上按一定工藝順序加工的模擬，能夠識別生產(chǎn)過程中的瓶頸環(huán)節(jié)；通過高速公路上交通流的模擬，可以分析車輛在路段上的分布特別是堵塞的狀況。與用物理模型的模擬實驗相比，計算機模擬有明顯的優(yōu)點：成本低、時間短、重復性高、靈活性強。數(shù)學建模與計算機模擬有人把計算機模擬作為建立數(shù)學模型的手段之一。但是數(shù)學模型在某種意義下描述了對象內(nèi)在特性的數(shù)量關系，其結果容易推廣，特別是得到了解析形式答案時，更易推廣。而計算機模擬則完全模仿對象的實際演變過程，難以從得到的數(shù)字結果分析對象的內(nèi)

49、在規(guī)律。當然，對于那些因內(nèi)部機理過于復雜、目前尚難以建立數(shù)學模型的實際對象，用計算機模擬獲得一些定量結果，可稱是解決問題的有效手段。1.2 數(shù)學建模的重要意義P4數(shù)學，作為一門研究現(xiàn)實世界數(shù)量關系和空間形式的科學，在它產(chǎn)生和發(fā)展的歷史長河中，一直是和人們生活的實際需要密切相關的。作為用數(shù)學方法解決實際問題的第一步，數(shù)學建模自然有著與數(shù)學同樣悠久的歷史。兩千多年以前創(chuàng)立的歐幾里德幾何，17世紀發(fā)現(xiàn)的牛頓萬有引力定律，都是科學發(fā)展史上數(shù)學建模的成功范例。進入20世紀以來，隨著數(shù)學以空前的廣度和深度向一切領域的滲透，和電子計算機的出現(xiàn)與飛速發(fā)展，數(shù)學建模越來越受到人們的重視。1) 在一般工程技術領域

50、，數(shù)學建模仍然大有用武之地。在以聲、光、熱、力、電這些物理學科為基礎的諸如機械、電機、土木、水利等工程技術領域中，數(shù)學建模的普遍性和重要性不言而喻。雖然這里的基本模型是已有的，但是由于新技術、新工藝的不斷涌現(xiàn)，提出了許多需要用數(shù)學方法解決的新問題；高速、大型計算機的飛速發(fā)展，使得過去即便有了數(shù)學模型也無法求解的課題（如大型水壩的應力計算、中長期天氣預報等）迎刃而解；建立在數(shù)學模型和計算機模擬基礎上的CAD技術，以其快速、經(jīng)濟、方便等優(yōu)勢，大量地替代了傳統(tǒng)工程設計中的現(xiàn)場實驗、物理模擬等手段。2) 在高新技術領域，數(shù)學建模幾乎是必不可少的工具。無論是發(fā)展通訊、航天、微電子、自動化等高新技術本身，還是將高新技術用于傳統(tǒng)工業(yè)去創(chuàng)造新工藝、開發(fā)新產(chǎn)品，計算機技術支持下的建模和模擬都是經(jīng)常使用的有效手段。數(shù)學建模、數(shù)值計算和計算機圖形學等相結合形成的計算機軟件，已經(jīng)被固化于產(chǎn)品中，在許多高新技術領域起著核心作用，被認為是高新技術的特征之一。在這個意義上，數(shù)學不再僅僅作為一門科學，是許多技術的基礎，而且直接定向了技術的前臺。高技術本質上是一種數(shù)學技術。3) 數(shù)學迅速進入一些新領域，為數(shù)學建模開拓了許多新的處女她。隨著數(shù)學向諸如經(jīng)濟、人口、生態(tài)、地質等所謂非物理領域的滲透，一些交叉學科如計量經(jīng)濟學、人口控制論、