【文章內(nèi)容簡(jiǎn)介】 set(,39。string39。,39。350039。)。%設(shè)置定容排風(fēng)機(jī)的容量 totalC1=str2double(r3C)dC*str2double(r3C)。%集控室剩余空氣量 cC1 = dC*str2double(r3C)。%集控室的空氣量 totalC2=totalC1totalC1*bC。%工作間剩余空氣量 cC3=totalC1*bC。%工作間的空氣量 totalC3=totalC2totalC2*dCShang。 cC2=totalC2*dCShang。%主機(jī)艙上層的空氣量 totalC4=totalC3totalC3*dCDi。 cC4=totalC3*dCDi。%通艙底的空氣量 totalC5=str2double(r3C)dC*str2double(r3C)totalC1*bCtotalC2*dCShangtotalC3*dCDi。 %主機(jī)艙底層的空氣量 set(,39。string39。,num2str(totalC5))。 set(,39。string39。,num2str(cC2))。 set(,39。string39。,num2str(cC4))。 %分別顯示出各個(gè)數(shù)值 Rv=str2double(r3C)*str2double(rC)。 r11C=Rv/(str2double(r4)*str2double(r2))。%每秒鐘的進(jìn)氣量時(shí)間 r16=str2double(t5)*str2double(t6)。 r17=str2double(T0)+str2double(t5)。 r18=r16/r17。 R1=num2str(r18)。%得到外界密度（可以不變） R2C=r11C*r18。%每秒鐘的總質(zhì)量此代碼是依據(jù)排風(fēng)機(jī)的型號(hào)，我們先設(shè)置排風(fēng)機(jī)的容量為‘3500’，并在e4控件中顯示，同時(shí)設(shè)置r3C變量為風(fēng)機(jī)的總?cè)萘?，dC為集控室開關(guān)的開度，從而得到cC1，即空氣的通氣量。同理可以分別得到工作間、主機(jī)艙上層、通艙底以及主機(jī)艙底層的空氣量并用text控件分別顯示出來。然后進(jìn)行固定數(shù)值的計(jì)算，如每秒鐘的進(jìn)氣量以及外界的溫度等。需要注意的是，這里涉及到了供風(fēng)機(jī)的開關(guān)操作的功能設(shè)計(jì)，為了能使風(fēng)機(jī)進(jìn)行開關(guān)操作，本課題設(shè)計(jì)了一個(gè)全局變量global anjian1C，此變量與按鍵TB7共同管理此供風(fēng)機(jī)的開關(guān)操作。下面介紹此風(fēng)機(jī)模塊的兩個(gè)密閉的工作空間，集控室與工作間的功能實(shí)現(xiàn)。集控室的代碼實(shí)現(xiàn)如下：if(k5C==1) %判斷排風(fēng)機(jī)開關(guān)是否開啟 r20C=str2double(T1C)+str2double(t5)。%室內(nèi)溫度34機(jī)房（變） r19C=R2C*dC。%計(jì)算進(jìn)入室內(nèi)的氣體的質(zhì)量 R2mC=num2str(r19C)。%得到每秒鐘進(jìn)入的氣體質(zhì)量 MC=str2double(39。10100039。)*str2double(VC)*str2double(t8)/(str2double(t7)*r20C)。%34機(jī)房的在一大氣壓下的空氣質(zhì)量（初始的室內(nèi)空氣質(zhì)量） m2C=str2double(m34C)。%m34很重要?。。。。?！ m2C=m2C+r19Cstr2double(MzC)。 m34C=num2str(m2C)。m1CP=(MC+str2double(m34C))*str2double(t7)*r20C/(str2double(t8)*str2double(VC))。%得到室內(nèi)的空氣壓力 m22C=m1CP/str2double(r4)。 m22Cc=num2str(m22C)。 set(,39。string39。,num2str(m22C))。%得到相對(duì)壓力end此代碼是該系統(tǒng)模塊的集控室的仿真數(shù)學(xué)模型，首先先獲得室內(nèi)的溫度T1C，然后在獲得每秒鐘進(jìn)入室內(nèi)的空氣的氣體質(zhì)量，再加上原室內(nèi)的空氣的質(zhì)量，得到此時(shí)室內(nèi)總的質(zhì)量。這里應(yīng)該注意一點(diǎn)，室內(nèi)總質(zhì)量是實(shí)時(shí)更新的，由于有供風(fēng)機(jī)與排風(fēng)機(jī)的影響，所以質(zhì)量有增有減，所以室內(nèi)的總質(zhì)量要設(shè)置一個(gè)全局變量，即global m34來記憶室內(nèi)的氣體質(zhì)量。每一個(gè)密閉空間都有這個(gè)變量來記憶室內(nèi)的氣體質(zhì)量。通過進(jìn)而得到室內(nèi)的氣體壓力，然后得到氣壓的相對(duì)壓力。同理，工作間的仿真模型也是一樣的，這里就不一一介紹了。排風(fēng)機(jī)的打開操作仿真完成，同理排風(fēng)機(jī)的關(guān)閉仿真也是相同的實(shí)現(xiàn)，只是對(duì)打開的代碼進(jìn)行了少量的修改，代碼實(shí)現(xiàn)如下：if(xor(get(,39。value39。) , (anjian1C))==1) if(anjian1C==1)set(,39。string39。,39。039。)。 set(,39。string39。,39。039。)。 set(,39。string39。,39。OFF39。)。 set(,39。string39。,39。039。)。 T0=get(,39。string39。)。 totalC1=str2double(39。039。)。 totalC2=str2double(39。039。)。 totalC3=str2double(39。039。)。 totalC4=str2double(39。039。)。 totalC5=str2double(39。039。)。 Rv=str2double(39。039。)。 R2C=str2double(39。039。)。 r20C=str2double(T1C)+str2double(t5)。%室內(nèi)溫度34機(jī)房（變） r19C=R2C*dC。關(guān)閉排風(fēng)機(jī)只是把相應(yīng)的排風(fēng)機(jī)的排風(fēng)量進(jìn)行清零，使此系統(tǒng)模塊進(jìn)行相應(yīng)的無進(jìn)風(fēng)量的操作。實(shí)現(xiàn)密閉空間零輸入的狀態(tài)。同時(shí)，又設(shè)立了一個(gè)global kC變量，此變量的目的是為了進(jìn)行各個(gè)相關(guān)數(shù)據(jù)的更新以及顯示內(nèi)容的更新操作。更新數(shù)據(jù)的代碼如下：if(kCamp。anjian1C==1) T1C=get(,39。string39。)。%更新集控室的溫度 T1Workspace=get(,39。string39。)。%更新工作間的溫度 T0=get(,39。string39。 )。%更新外界的溫度 r4=get(,39。string39。 )。%更新外界的壓強(qiáng) dC=get(,39。value39。)。%更新集控室的開關(guān)開度 bC=get(,39。value39。)。%更新工作間的開關(guān)開度 dCShang=get(,39。value39。)。%更新主機(jī)艙上層的開關(guān)開度 dCDi=get(,39。value39。)。%更新主機(jī)艙通艙底的開關(guān)開度 num2C=get(,39。value39。)。 set(,39。String39。,num2C)。%更新開關(guān)開度并在text中顯示出來 num1C=get(,39。value39。)。 set(,39。String39。,num1C)。 num3C=get(,39。value39。)。 set(,39。String39。,num3C)。 num4C=get(,39。value39。)。 set(,39。String39。,num4C)。 kC=0。%設(shè)置更新標(biāo)志位歸零使其繼續(xù)進(jìn)行操作 set(,39。string39。,mtt2)。 set(,39。string39。,39。設(shè)置中...39。)。end以上代碼實(shí)現(xiàn)了各個(gè)數(shù)據(jù)的更新操作使得仿真的數(shù)據(jù)更能接近實(shí)際的情況，并且為仿真接近真實(shí)性提供了一個(gè)方便的途徑。相應(yīng)的系統(tǒng)模塊1的完整代碼在附錄中可以找到。這里我們就一個(gè)較為復(fù)雜的變?nèi)菖棚L(fēng)機(jī)系統(tǒng)模塊進(jìn)行建模仿真實(shí)現(xiàn)。由上述提到的公式()()()()我們可以得到變?nèi)莸娘L(fēng)機(jī)公式： （）其中是變?nèi)菖棚L(fēng)機(jī)的當(dāng)前壓力，是當(dāng)前排出的氣體量，這是手動(dòng)進(jìn)行調(diào)節(jié)的。是基準(zhǔn)壓力，是基準(zhǔn)的排氣量。這樣通過公式()我們就可以算出氣體變?nèi)莸呐棚L(fēng)機(jī)每分鐘的排出的氣體質(zhì)量。 來介紹變?nèi)莸娘L(fēng)機(jī)系統(tǒng)。 系統(tǒng)模塊3仿真界面。根據(jù)仿真圖，實(shí)現(xiàn)變?nèi)菖棚L(fēng)機(jī)的仿真。這里我們只介紹變?nèi)菖棚L(fēng)機(jī)與定容排風(fēng)機(jī)的異同之處，對(duì)于相同的地方我們就不詳細(xì)介紹。仿真代碼如下：if(kj==0)if (get(,39。value39。)amp。anjian1j==1) if(anjian1j==1) set(,39。string39。,get(,39。value39。))。%設(shè)置變?nèi)菖棚L(fēng)機(jī)的當(dāng)前容量值 totalj=dj*bj。%算出進(jìn)入1,2號(hào)發(fā)電機(jī)中的氣體量 cj = num2str(totalj)。 e1j=djtotalj。%計(jì)算出進(jìn)入分油機(jī)中的氣體量 dengshu52j=num2str(e1j)。 r11j=str2double(r)/(str2double(r4)*str2double(r2))。%時(shí)間 r12j=str2double(r3)*str2double(r3)。 r13j=r11j/r12j。%公式(10)給出的的計(jì)算公式 r16=str2double(t5)*str2double(t6)。 r17=str2double(T0)+str2double(t5)。%計(jì)算出外界的溫度 r18=r16/r17。 R1=num2str(r18)。 %得到外界密度 r14j=dj*dj*dj。 r15j=r13j*r14j。 Rj=num2str(totalj)。%相對(duì)于外界的體積%這里的全局變量global kj是更新數(shù)據(jù)的標(biāo)志位，global anjian1j是判斷供風(fēng)機(jī)的按鍵是否開啟，這些與定容排風(fēng)機(jī)系統(tǒng)相似。根據(jù)公式(10)的計(jì)算公式，算出變?nèi)菖棚L(fēng)機(jī)每秒鐘排出的氣體質(zhì)量。 系統(tǒng)模塊3仿真界面然后介紹一下分油機(jī)間的代碼實(shí)現(xiàn)，由于分油機(jī)間有一個(gè)分油機(jī)間的排風(fēng)機(jī)，所以與前面介紹的密閉空間有一定的不同，分油機(jī)間的排風(fēng)機(jī)仿真代碼如下：if(get(,39。value39。)==1 ) set(,39。string39。,39。50039。)。 Rvshopj=str2double(39。50039。)*str2double(39。60039。)。 r11Cshopj=Rvshopj/(str2double(r4)*str2double(r2))。%每秒鐘的進(jìn)氣量時(shí)間 r16=str2double(t5)*str2double(t6)。 r17=str2double(T0)+str2double(t5)。 r18=r16/r17。 w22shopj=r11Cshopj*r18。%每秒鐘的總質(zhì)量 end if(get(,39。value39。)==0) set(,39。string39。,39。OFF39。)。 w22shopj=0。 end代碼實(shí)現(xiàn)，當(dāng)分油機(jī)間的排風(fēng)機(jī)打開，進(jìn)入氣體的質(zhì)量是w22shopj，排風(fēng)機(jī)關(guān)閉時(shí)進(jìn)入的氣體量為零。以下是分油機(jī)間的仿真代碼實(shí)現(xiàn)：r20jj=str2double(T1jj)+str2double(t5)。%室內(nèi)溫度分油機(jī)間 r19jj=r15j*r18r15j*r18*bj。 R2jj=num2str(r19jj)。 %得到每秒鐘進(jìn)入的氣體質(zhì)量M34jj=str2double(39。10100039。)*str2double(V34jj)*str2double(t8)/(str2double(t7)*r20jj)。%34機(jī)房的空氣質(zhì)量在一大氣壓 m2jj=str2double(m34jj)。 m2jj=m2jj+r19jjw22shopjstr2double(Mz34jj)。 m34jj=num2str(m2jj)。 m1jj=(M34jj+str2double(m34jj))*str2double(t7)*r20jj/(str2double(t8)*str2double(V34jj))。 %得到室內(nèi)的空氣壓力 m22jj=m1jj/str2double(r4)。 set(,39。string39。,num2str(m22jj))。 %得到相對(duì)壓力分油機(jī)間的代碼除了多了一個(gè)排風(fēng)機(jī)外，其他代碼實(shí)現(xiàn)和其他的密閉空間實(shí)現(xiàn)基本相同。 以上一節(jié)介紹的分油機(jī)間的模型為例介紹一下密閉空間模塊的仿真實(shí)現(xiàn)，密閉空間模塊就是運(yùn)用了MATLAB GUI中不同界面?zhèn)鲗?dǎo)數(shù)據(jù)的技術(shù)實(shí)現(xiàn)的。仿真代碼如下：txt1jj=dengshu52j。%每秒鐘的氣體體積 txt2jj=R1。%外界的氣體密度 txt3jj=T1jj。%室內(nèi)溫度34機(jī)房 txt4jj=R2jj。 %得到每秒鐘進(jìn)入的氣體質(zhì)量 txt5jj=num2str(M34jj)。%一大氣壓愿機(jī)