第一篇：中南大學(xué)化工原理仿真實(shí)驗(yàn)報(bào)告

化工原理計(jì)算機(jī)仿真實(shí)驗(yàn)

班級(jí)：化學(xué)工程與工藝1102班

姓名：王翔

學(xué)號(hào)：1505110321

日期：2014年1月1日

本套軟件系統(tǒng)包括8個(gè)單元仿真實(shí)驗(yàn)：

實(shí)驗(yàn)一 離心泵性能的測(cè)試

實(shí)驗(yàn)二 管道阻力實(shí)驗(yàn)

實(shí)驗(yàn)三 傳熱實(shí)驗(yàn)

實(shí)驗(yàn)四 吸收實(shí)驗(yàn)

實(shí)驗(yàn)五 流體流動(dòng)形態(tài)的觀測(cè)

實(shí)驗(yàn)六 柏努利方程實(shí)驗(yàn)

實(shí)驗(yàn)七 干燥實(shí)驗(yàn)

實(shí)驗(yàn)八 精餾實(shí)驗(yàn)

以下是實(shí)驗(yàn)?zāi)M觀測(cè)過程和計(jì)算機(jī)生成的實(shí)驗(yàn)報(bào)告。

圖1 離心泵性能的測(cè)試 觀察氣蝕現(xiàn)象（1）

圖2 離心泵性能的測(cè)試 觀察氣蝕現(xiàn)象（2）

圖3 離心泵性能的測(cè)試 離心泵特性曲線測(cè)定實(shí)驗(yàn)報(bào)告（1）

圖4 離心泵性能的測(cè)試 離心泵特性曲線測(cè)定實(shí)驗(yàn)報(bào)告（2）

圖5 離心泵性能的測(cè)試 離心泵特性曲線測(cè)定實(shí)驗(yàn)報(bào)告（3）

圖6 離心泵性能的測(cè)試 離心泵特性曲線測(cè)定實(shí)驗(yàn)報(bào)告（4）

圖7 管道阻力的測(cè)定實(shí)驗(yàn)報(bào)告（1）

圖8 管道阻力的測(cè)定實(shí)驗(yàn)報(bào)告（2）

圖9 管道阻力的測(cè)定實(shí)驗(yàn)報(bào)告（3）

圖10 傳熱實(shí)驗(yàn)

圖11 傳熱實(shí)驗(yàn)報(bào)告（1）

圖12 傳熱實(shí)驗(yàn)報(bào)告（2）

圖13 傳熱實(shí)驗(yàn)報(bào)告（3）

圖14 傳熱實(shí)驗(yàn)報(bào)告（4）

圖15 吸收實(shí)驗(yàn) 觀察液泛現(xiàn)象

圖16 吸收實(shí)驗(yàn)報(bào)告

圖17 液體流動(dòng)形態(tài)的觀測(cè) 觀察滯留形態(tài)

圖18 液體流動(dòng)形態(tài)的觀測(cè)實(shí)驗(yàn)報(bào)告

圖19 柏努利方程實(shí)驗(yàn) 觀察測(cè)壓孔與水流方向方位角與水位變化（1）

圖20 柏努利方程實(shí)驗(yàn) 觀察測(cè)壓孔與水流方向方位角與水位變化（2）

圖21 干燥實(shí)驗(yàn)報(bào)告（1）

圖22 干燥實(shí)驗(yàn)報(bào)告（2）

圖23 干燥實(shí)驗(yàn)報(bào)告（3）

圖24 干燥實(shí)驗(yàn)報(bào)告（4）

圖25 精餾實(shí)驗(yàn) 動(dòng)態(tài)平衡調(diào)整

圖26 精餾實(shí)驗(yàn)報(bào)告（1）

圖27 精餾實(shí)驗(yàn)報(bào)告（2）

第二篇：中南大學(xué)系統(tǒng)仿真實(shí)驗(yàn)報(bào)告

ans =

實(shí)驗(yàn)一 MATLAB 中矩陣與多項(xiàng)式的基本運(yùn)算 實(shí)驗(yàn)任務(wù) 1.了解 MATLAB 命令窗口和程序文件的調(diào)用。?熟悉如下 MATLAB 的基本運(yùn)算：

① 矩陣的產(chǎn)生、數(shù)據(jù)的輸入、相關(guān)元素的顯示； ② 矩陣的加法、乘法、左除、右除； ③ 特殊矩陣：?jiǎn)挝痪仃?、?1 ”矩陣、“ 0 ”矩陣、對(duì)角陣、隨機(jī)矩陣的產(chǎn)生和 運(yùn)算； ④ 多項(xiàng)式的運(yùn)算：多項(xiàng)式求根、多項(xiàng)式之間的乘除。

基本命令訓(xùn)練 1、>> eye(2)ans = 1 0 0 1 >> eye(4)ans = 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 2、>> ones(2)1 1

ans =1 >> ones(4)ans = 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 >> ones(2,2)ans = 1 1 1 1 >> ones(2,3)ans = 1 1 1 1 1 1 >> ones(4,3)ans = 1 1 1 1 1 1 1 1 1 1 1 1 3、>> zeros(2)

0 0 0 0 >> zeros(4)ans = 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 >> zeros(2,2)ans = 0 0 0 0 >> zeros(2,3)ans = 0 0 0 0 0 0 >> zeros(3,2)ans = 0 0 0 0 00 4、隨機(jī)陣 >> rand(2,3)ans = 0.2785 0.9575 0.1576 0.5469 0.9649 0.9706 >> rand(2,3)

ans = 0.9572 0.8003 0.4218 0.4854 0.1419 0.9157 5、>> diag(5)ans = 5 >> diag(5,5)ans = 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 5 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 >> diag(2,3)ans = 0 0 0 2 0 0 0 0 0 0 0 0 0 0 0 0 6、(inv(A)為求 A 的逆矩陣)>> B=[5 3 1;2 3 8;1 1 1],inv(B)3 1

1

ans =

0.6250 0.2500-2.6250-0.7500-0.5000 4.7500 0.1250 0.2500-1.1250

>> A=[2 3;4 4],B=[5 3;3 8],inv(A),inv(B);AB,A/B,inv(A)*B,B*inv(A)A = 2 3 4 4 B = 5 3 3 8 ans =-1.0000 0.7500 1.0000-0.5000 ans =-2.7500 3.0000 3.5000-1.0000 ans = 0.2258 0.2903 0.6452 0.2581 ans =-2.7500 3.0000

3.5000 -1.0000 ans =-2.0000 2.2500 5.0000-1.7500 7、>> p =[1,-6,-72,-27], roots(p)p = 1-6-72-27 ans = 12.1229-5.7345-0.3884 >> p=[2,3,6],roots(p)p = 2 3 6 ans =-0.7500 + 1.5612i-0.7500-1.5612i 8、(A 為 n*n 的方陣)>> A=[0 1 0;-4 4 0;-2 1 2],poly(A),B=sym(A),poly(B)A = 0 1 0-4 4 0-2 1 2 ans =-6 12-8

B = [ 0, 1, 0] [-4, 4, 0] [-2, 1, 2]

ans = x A 3-6*x A 2+12*x-8 9,、(conv 是多項(xiàng)式相乘，deconv 是多項(xiàng)式相除)>> u=[1 2 4 6 ],v=[5 0 0-6 7],conv(u,v)u = 1 2 4 6 v =0 0-6 7

ans =10 20 24-5-10-8 42

>> v=[1 2 4 6 ],u=[5 0 0-6 7],deconv(u,v)v = 1 2 4 6 u = 5 0 0-6 7 ans = 5-10 10、（點(diǎn)乘是數(shù)組的運(yùn)算，沒有點(diǎn)的乘是矩陣運(yùn)算）

>> a = [2 5;3 4], b =[3 1;4 7],a.*b,a*b a = 2 5 3 4 b = 3 1 4 7 ans = 6 5 12 28 ans = 26 37 25 31 >> a = [2 3];b = [4 7];a.*b = [8 21];a*b %錯(cuò)誤 a*b“ = 29;11、（who 可以看到你用過的一些變量，來了）

>> who Your variables are: A B a ans b p u >> whos Name Size Bytes whos 是把該變量及所存儲(chǔ)的大小等信息都顯示出 Class Attributes 2x2 32 double

B 2x2 32 double a 1x2 16 double ans 1x2 16 double b 1x2 16 double p

1x3 24 double u 1x5 40 double v 1x4 32 double12、>> A=[2 5 3;6 5 4],disp(A),size(A),length(A)A = 2 5 3 6 5 4 2 5 3 6 5 4 ans = 2 3 ans = 3 實(shí)驗(yàn)二 MATLAB 繪圖命令 實(shí)驗(yàn)任務(wù) 熟悉 MATLAB 基本繪圖命令，掌握如下繪圖方法：．坐標(biāo)系的選擇、圖形的繪制；?圖形注解（題目、標(biāo)號(hào)、說明、分格線）的加入;3 ?圖形線型、符號(hào)、顏色的選取 基本命令訓(xùn)練 1、>>t=[0:pi/360:2*pi];x=cos(t)+ cos(t*4);y=si n(t)+ sin(t*4);xlabel(”x 軸“);ylabel(”y 軸“);plot(y,x),grid;2、>> t=0:0.1:100;x=3*t;y=4*t;z=si n(2*t);plot3(x,y,z, ”g:“)

■15 i 0 5 0 05 1 1 5 2

3、>>x = linspace(-2*pi,2*pi,40);y=si n(x);stairs(x,y)4、>> t=[0:pi/360:2*pi];x=cos(t)+ cos(t*4)+ sin(t*4);y=si n(t)+ si n(t*4);plot(y,x, ”r:“);

xlabel(”x 軸“);ylabel(”y 軸“);

6、>>th=[0:pi/20:2*pi];x=exp(j*th);plot(real(x),imag(x),”r-.“);grid;text(0,0,”中心“);

5、>> th=[0:pi/1000:2*pi]”;r=cos(2*th);polar(th,r);title(“四葉草圖”)270 四葉草圖

7、>>x=-2:0.01:2;8、y=-2:0.01:2;9、[X,Y] = meshgrid(x,y);Z = Y.*exp(-X.A 2-Y.A 2);[C,h] = con tour(X, YZ);set(h,“ShowText”,“o n”,“TextStep”,get(h,“LevelStep”)*2)_ 1OS ?I ,5 2 ■n.s o o.s 8、>>x = 0:0.2:10;y = 2*x+3;subplot(411);plot(x,y);grid;title(“y 的原函數(shù)”);subplot(412);semilogy(x,y);grid;title(“對(duì) y 取對(duì)數(shù)”);丫的原畫數(shù) 40----------1------------------------------------------1--------------------1--------------------1-----------------------------------------1--------------------1-------------------1--------

l| I p | il ■ | i| I 九 _____ 1-___ — ___ I ____ L ___ :…… ：

___ J _

_______ L ___ u i| I |l I , il _-■」 “ j I ■I __ h-_____________ I I ■ Q 」 【 I 1

F I I I I II I I IT 10 1

□ 1 2 3 4 5 6 r 6 9 10 , 對(duì) y 取對(duì)數(shù) 對(duì)弋觀對(duì)數(shù) subplot(413);semilogx(x,y);10 10 id 10 1lZ 10 w 10 10 40 20 0 對(duì)好對(duì)數(shù) ?掌握循環(huán)、分支語句的編寫，學(xué)會(huì)使用 look for、help 命令

grid;title(”對(duì) x 取對(duì)數(shù)“);subplot(414);loglog(x,y);grid;title(”對(duì) xy 均取對(duì)數(shù)“);9、>>x =-3:0.3:3;bar(x,exp(-x.*x),”g“)實(shí)驗(yàn)三 MATLAB 程序設(shè)計(jì) 實(shí)驗(yàn)任務(wù) 1 ?熟悉 MATLAB 程序設(shè)計(jì)的方法和思路；

程序舉例 1、>> f=[1,1];i=1;while f(i)+f(i+1)<1000 f(i+2)=f(i)+f(i+1);i=i+1;end f,i f = Columns 1 through 14 1 1 2 3 5 8 13 21 34 55 89 144 233 377 Columns 15 through 16 610 987 i = 15 2、>> m=3;n=4;for i=1:m for j=1:n a(i,j)=1/(i+j-1);end end format rat1/2 1/3 1/2 1/3 1/4 1/3 1/4 1/5(分?jǐn)?shù)格式形式。用有理數(shù)逼近顯示數(shù)據(jù))>>m=5;n=4;for i=1:m for j=1:n a(i,j)=1/(i+j-1);end end format rat a a = 1 1/2 1/3 1/4 1/2 1/3 1/4 1/5 1/3 1/4 1/5 1/6 1/4 1/5 1/6 1/7 1/5 1/6 1/7 1/8 3、程序中沒有 format rat 命令時(shí)，如果上次運(yùn)行結(jié)果沒有清除，輸出的結(jié)果就是上次運(yùn)行的 結(jié)果！但是運(yùn)用 clear 命令清楚之前的運(yùn)行結(jié)果之后就會(huì)正常運(yùn)行。

4、>> x=input（”請(qǐng)輸入 x 的值:’）;if x==10 y=cos(x+1)+sqrt(x*x+1);else y=x*sqrt(x+sqrt(x));1/4 1/5 1/6

end y

請(qǐng)輸入 x 的值 :2 y = 2391/647 x=input(“請(qǐng)輸入 x 的值:’);if x==10 y= fprintf(” 不在定義域內(nèi)，請(qǐng)重新輸入：

“);return else y=1/(x-10);end y 請(qǐng)輸入 x 的值 :2-1/8 5、>> p=[0 0 0 1 3 0 2 0 0 9];for i=1:length(p), if p(1)==0,p=p(2:length(p));end;end;p

p = Columns 1 through 5 3

0

0 Columns 6 through 7

0 9

>>p=[0 0 0 1 3 0 2 0 0 9];p(p==0)=[];p

p =3

6、>> e2(500)

ans =1 2 3 5 8 13 21 34 55 89 144 233

377 >> lookfor ffibno e2-ffibno 計(jì)算斐波那契亞數(shù)列的函數(shù)文件 >> help e2 ffibno 計(jì)算斐波那契亞數(shù)列的函數(shù)文件 n 可取任意自然數(shù) 程序如下（用法：

lookfor 關(guān)鍵詞 在所有 M 文件中找 “關(guān)鍵詞 ”，比如：

lookfor max（即尋找關(guān)鍵詞 “ max”）

其實(shí)就和我們平時(shí)用 CTRL+F 來查找關(guān)鍵詞”是一樣的 而 help 是顯示 matlab 內(nèi)置的幫助信息 用法：

help 命令，比如 help inv，作用就是調(diào)用 inv 這個(gè)命令的幫助）

程序設(shè)計(jì)題 用一個(gè) MATLAB 語言編寫一個(gè)程序：輸入一個(gè)自然數(shù)，判斷它是否是素?cái)?shù)，如果是，輸出“ It is one prime ”，如果不是，輸出“ It is n ot o ne prime.”。要求通 過調(diào)用子函數(shù)實(shí)現(xiàn)。

最好能具有如下功能：

①設(shè)計(jì)較好的人機(jī)對(duì)話界面，程序中 含有提示

性的輸入輸出語句。

②能實(shí)現(xiàn)循環(huán)操作，由操作者輸入相關(guān)命令來控制 是否繼續(xù)進(jìn)行素?cái)?shù)的判斷。

如果操作者希望停止這種判斷，則可以退出程序。

③ 如果所輸入的自然數(shù)是一個(gè)合數(shù)，除了給出其不是素?cái)?shù)的結(jié)論外，還應(yīng)給出至少 一種其因數(shù)分解形式。例：輸入 6，因?yàn)?6 不是素?cái)?shù)。則程序中除了有“ It is not one prime ”的結(jié)論外，還應(yīng)有：“ 6=2*3 ”的說明。

function sushu while 1 x=input（” 請(qǐng)輸入一個(gè)自然數(shù) “）;if x<2 disp（” 既不是質(zhì)數(shù)又不是合數(shù) “）;Else if isprime(x)==1 disp（” 這是一個(gè)素?cái)?shù) “）;Else disp(” 這是一個(gè)合數(shù)，可以因式分解為：

“)for n=2:sqrt(x)if rem(x,n)==0 num3=x;num1=n;num2=x/n;disp([num2str(num3),”=“,num2str(num1),”x“,num2str(num2)])End End End

end y=input(” 是否繼續(xù)判斷？繼續(xù)請(qǐng)按 1 ,按任意鍵退出 ：

“)； if y~=1 break end end 實(shí)驗(yàn)四 MATLAB 的符號(hào)計(jì)算與 SIMULINK 的使用 實(shí)驗(yàn)任務(wù) 1.掌握 MATLAB 符號(hào)計(jì)算的特點(diǎn)和常用基本命令;2.掌握 SIMULINK 的使用。

程序舉例 1.求矩陣對(duì)應(yīng)的行列式和特征根 >> a=sym(”[a11 a12;a21 a22]“);da=det(a)ea=eig(a)da = a11*a22-a12*a21 ea = 1/2*a11+1/2*a22+1/2*(a11 A 2-2*a11*a22+a22 A 2+4*a12*a21)A(1 /2)1/2*a11+1/2*a22-1/2*(a11 A 2-2*a11*a22+a22 A 2+4*a12*a21)A(1 /2)a=sym(”[2 3;1 5]“);

da=det(a)ea=eig(a)da = 7 ea = 7/2+1/2*21A(1/2)7/2-1/2*21A(1/2)2.求方程的解(包括精確解和一定精度的解)>> r1=solve(”xA2+x-1“)rv=vpa(r1)rv4=vpa(r1,4)

rv20=vpa(r1,20)r1 = 1/2*5 八(1/2)-1/2-1/2*5 八(1/2)-1/2 rv =.6894848260-1.689484826 rv4 =.6180-1.618 rv20 =.689484820-1.68948482 3 ． a=sym(”a“);b=sym(”b“);c=sym(w=10;x=5;y=-8;z=11;A=[a,b;c,d] B=[w,x;y,z] det(A)det(B)A = [ a, b] ”c“);d=sym(”d“);%定義 4 個(gè)符號(hào)變量 %定義 4 個(gè)數(shù)值變量 %建立符號(hào)矩陣 A %建立數(shù)值矩陣 B %計(jì)算符號(hào)矩陣 A 的行列式 %計(jì)算數(shù)值矩陣 B 的行列式

[ c, d]

X

10-8 11 ans = a*d-b*c ans = 150 4.? syms x y;s=(-7*x A 2-8*y A 2)*(-x A 2+3*y A 2);expand(s)% 對(duì) s 展開 collect(s,x)% 對(duì) s 按變量 x 合并同類項(xiàng)（無同類項(xiàng)）

factor(ans)% 對(duì) ans 分解因式

ans = 7*xA4-13*xA2*yA2-24*yA4 ans = 7*xA4-13*xA2*yA2-24*yA4 ans =(8*yA2+7*xA2)*(xA2-3*yA2)5.對(duì)方程 AX=b 求解 >> A=[34,8,4;3,34,3;3,6,8];b=[4;6;2];X=linsolve(A,b)% 調(diào)用 linsolve 函數(shù)求解 Ab % 用另一種方法求解

0.0675 0.1614

diff(f)%未指定求導(dǎo)變量和階數(shù)，按缺省規(guī)則處理

0.1037 ans = 0.0675 0.1614 0.1037 6 ． 對(duì)方程組求解 a11*x1+a12*x2+a13*x3=b1 a21*x1+a22*x2+a23*x3=b2 a31*x1+a32*x2+a33*x3=b3 >>syms a11 a12 a13 a21 a22 a23 a31 a32 a33 b1 b2 b3;A=[a11,a12,a13;a21,a22,a23;a31,a32,a33];b=[b1;b2;b3];XX=Ab %用左除運(yùn)算求解(X=linsolve(A,b)%調(diào)用 linsolve 函數(shù)求的解)XX =(a12*a23*b3-a12*b2*a33+a13*a32*b2-a13*a22*b3+b1*a22*a33-b1*a32*a23)/(a11*a22*a33-a11*a32*a23-a12*a21*a33+a32*a21*a13-a22*a31*a13+a31*a12*a 23)-(a11*a23*b3-a11*b2*a33-a21*a13*b3-a23*a31*b1+b2*a31*a13+a21*b1*a33)/(a11*a22*a33-a11*a32*a23-a12*a21*a33+a32*a21*a13-a22*a31*a13+a31*a12* a23)(a32*a21*b1-a11*a32*b2+a11*a22*b3-a22*a31*b1-a12*a21*b3+a31*a12*b2)/(a11*a22*a33-a11*a32*a23-a12*a21*a33+a32*a21*a13-a22*a31*a13+a31*a12*a 23)7 ． syms a b t x y z;f=sqrt(1+exp(x));

%求 f 對(duì) x 的二階導(dǎo)數(shù) %求 f 對(duì) x 的三階導(dǎo)數(shù) %按參數(shù)方程求導(dǎo)公式求 y 對(duì) x 的導(dǎo)數(shù) ans = 1/2/(1+exp(x))A(1/2)*exp(x)ans =-2*si n(x)-x*cos(x)ans =-3*cos(x)+x*si n(x)ans =-b*cos(t)/a/si n(t)三、SIMULINK 的使用

f=x*cos(x);diff(f,x,2)diff(f,x,3)f1=a*cos(t);f2=b*si n(t);diff(f2)/diff(f1)

diff(f)%未指定求導(dǎo)變量和階數(shù)，按缺省規(guī)則處理

仿真圖: 波形圖: 其中: G i(R

0.14

實(shí)驗(yàn)五 MATLAB 在控制系統(tǒng)分析中的應(yīng)用 實(shí)驗(yàn)任務(wù) 1.掌握 MATLAB 在控制系統(tǒng)時(shí)間響應(yīng)分析中的應(yīng)用;2.掌握 MATLAB 在系統(tǒng)根軌跡分析中的應(yīng)用； 3.掌握 MATLAB 控制系統(tǒng)頻率分析中的應(yīng)用； 4.掌握 MATLAB 在控制系統(tǒng)穩(wěn)定性分析中的應(yīng)用 基本命令 1.step 2.impulse 3.in itial 4.Isim 5.rlocfi nd 6.bode 7.margin 8.nyquist 9.Nichols 10.cloop 程序舉例 1.求下面系統(tǒng)的單位階躍響應(yīng) 1.5 G(s)

t >> num=[4];den=[1 , 1 , 4];step(num , den)[y , x , t]=step(num , den);tp=spli ne(y , t , max(y))% 計(jì)算峰值時(shí)間 max(y)% 計(jì)算峰值 tp = Step Response 0.5 Time(sec)p m 1.6062 ans = 1.4441 0.18 0.16

2.求如下系統(tǒng)的單位階躍響應(yīng) X i 0 1 X i 0 u X 2 6 5 X 2 1

>> a=[0,1;-6,-5];b=[0;1];c=[1,0];d=0;[y,x]=step(a,b,c,d);Plot(y)

1,0 X i X 2 3.求下面系統(tǒng)的單位脈沖響應(yīng): G(s)兒 >> num=[4];den=[1 , 1 ,4];impulse(nu m,de n)Time(sec)Response to Initial Conditions 4.已知二階系統(tǒng)的狀態(tài)方程為: e u p m A 0.25 0.2 0.15 0.1 0.05 0-0.05

L I 1-

r

-廠、、/ Impulse Response

c=[1 , 0] ；

d=[0];xO=[1 ,0] ； subplot(1 , 2,1);in itial(a , b , c ,d,x0)subplot(1 , 2,2);impulse(a , b , c , d):系統(tǒng)傳遞函數(shù)為: G(s)占 輸入正弦信號(hào)時(shí)，觀察輸出 信號(hào)的相位差。

>> num=[1];den=[1 ,1];t=0 : 0.01 : 10;u=s in(2*t);hold on plot(t,u, ”r“)lsim(nu m,de n,u,t)

Real Axis

6.有一二階系統(tǒng)，求出周期為 4 秒的方波的輸出響應(yīng) 2s 2

5s 1 ~2 s 2s 3 >>num=[2 5 1];den=[1 2 3];t=(0:.1:10);period=4;u=(rem(t,period)>=period./2);% 看 rem 函數(shù)功能 lsim(nu m,de n,u,t);7.已知開環(huán)系統(tǒng)傳遞函數(shù)，繪制系統(tǒng)的根軌跡，并分析其穩(wěn)定性

G(s)k(s 2)(s 2

4s 3)2

>>num=[1 2];den 仁 [1 4 3];den=conv(de n1,de n1);figure(1)-6 rlocus(num,den)[k,p]= rlocfi nd(nu m,de n)l L l ■ 匚一 L — C--------------------

---------

-r t r r

-10-8-6-4-2 4 6 Root Locus 0 2 2 2

0 0 2 2--CPXA vyanma n-8 2.5 G(s)Lin ear Simulati on Results1.5 e 0.5-0.5-1-1.5-2 Time(sec)p m

impulse response(k=55)

5

figure(2)k=55;num 仁 k*[1 2];den=[1 4 3];den 1=c onv(de n,den);[nu m,de n]=cloop(nu m1,de n1,-1);impulse(nu m,de n)title(”impulse resp onse(k=55)“)-1.5-1 0 200 400 600 Time(sec)800 1000 1200 5 5

0 05

0 0

o

--x 10 8!-------

impulse resp on se(k=56)figure(3)k=56;num 仁 k*[1 2];den=[1 4 3];den 1=c onv(de n,den);[nu m,de n]=cloop(nu m1,de n1,-1);impulse(nu m,de n)-2-4-6-8 L

0 500 1000 Time(sec)title(”impulse resp on se(k=56)“)1500 2000 2500

Select a point in the graphics win dow selected_po int =-2.5924-0.0248i 0.7133-3.4160-2.5918-0.9961 + 0.4306i-0.9961-0.4306i Bode Diagram 8.作如下系統(tǒng)的 bode 圖 G(s)>> n=[1 , 1];d=[1 , 4,11 , 7];bode(n , d),grid on Frequency(rad/sec)s 1 s 3

4s 2

11s 7 9.系統(tǒng)傳函如下 G(s)s 1 0.5s e(s 2)3 求有理傳函的頻率響應(yīng)，然后在同一張圖上繪出以四階伯德近似表示的系統(tǒng)頻 率響應(yīng) >> num=[1];de n=co nv([1 2],co nv([1 2],[1 2]));

w=logspace(-1,2);t=0.5;

ylabel(”gai n“);subplot(2,1,2);semilogx(w,p1,w,p2,”g--“);grid on;xlabel(”freque ncy“);ylabel(”phase“);10.已知系統(tǒng)模型為

求它的幅值裕度和相角裕度 >> n=[3.5];d=[1 2 3 2];[Gm,Pm,Wcg,Wcp]=margi n(n,d)

G(s)3.5 s 3

2s 2

3s 2 bode plot [m1,p1]=bode(nu m,de n, 2);p1=p1-t*w”*180/pi;[n 2,d2]=pade(t,4);nu mt=c onv(n2,nu m);den t=(c onv(de n, d2));2 freque xlabel(“freque

Gm =

1.1433 Pm =

7.1688 Wcg =

1.7323 Wcp =

1.6541

nyq uist(n, d1);hold on nyq uist(n, d2);nyq uist(n, d3);nyq uist(n, d4);Nyquist Diagram 80 |--------------------

--------------------[--------------------[---------

11.二階系統(tǒng)為: G(s)n s 2n S S 令 wn=1, 分別作出 E =2,1 , 0.707 , C A 5 y

a 時(shí)的 nyquist 曲線。

m >> n=[1];d 仁 [1 , 4,1];d2=[1 , 2 , 1] ； d3=[1 , 1.414,1];d4=[1,1,1];Nyquist Diagram Real Axis 12.已知系統(tǒng)的開環(huán)傳遞函數(shù)為

_一-”—一 一-一_一_ ~-60

005

S S

2)

S S32

S S

v vr ra an nk ky ya a

14.一多環(huán)系統(tǒng)，其結(jié)構(gòu)圖如下，使用 Nyquist 頻率曲線判斷系統(tǒng)的穩(wěn)定性。

16.7s(0.85s 1)(0.25s 1)(0.0625s 1)繪制系統(tǒng)的 Nyqusit 圖, 并討論系統(tǒng)的穩(wěn)定性.>> G=tf(1000,co nv([1,3,2],[1,5]));nyquist(G);axis(”square“)13.分別由 w 的自動(dòng)變量和人工變量作下列系統(tǒng)的 nyquisF 曲線：

m

G(s)1 s(s 1)>> n=[1];d=[1 , 1 ,0];nyquist(n ,d);% 自動(dòng)變量 n=[1];d=[1 , 1 ,0];w=[0.5 : 0.1 : 3];nyq uist(n , d , w);% 人工變量-2-1-1.5-0.9-0.8-0.7-0.6-0.5 Real Axis-0.4-0.3-0.2-0.1 5 5

--0 0

S SI IX XA A

u u^ ^a an nT T9 9a a

卩

G(s)

30 20 10 0

figure(2)[nu m2,de n2]=cloop(nu m,de n);impulse(nu m2,de n2);

>> k1=16.70.0125;z1=[0];p1=[-1.25-4-16];[nu m1,de n1]=zp2tf(z1,p1,k1);[nu m,de n]=cloop(nu m1,de n1);[z,p,k]=tf2zp(nu m,de n);p figure(1)nyq uist(nu m,de n)-2-1 Nyquist Diagram-0.5 0 0.5 Real Axis 1 1.5 1 0.5 s a 0 n g m-0.5-1-1.5 1.5 20 Impulse Resp onse-10.5969 +36.2148i-10.5969-36.2148i-0.0562 15.已知系統(tǒng)為: eanML—pm

-5-10-15 0 0.1 0.2 0.3 Time(sec)Nichols Chart 0.4 0.5 0.6-2 10 0 10 2

Open-Loop Phase(deg)Frequency(rad/sec)

G(s)s(s 1)作該系統(tǒng)的 nichols 曲線。

>> n=[1];d=[1 , 1 , 0];ni chols(n , d);16.已知系統(tǒng)的開環(huán)傳遞函數(shù)為: G(s)k s(s 1)(s 2)當(dāng) k=2 時(shí)，分別作 nichols 曲線和波特圖 >> num=1;den=conv(co nv([1 0],[1 1]),[0.5 1]);subplot(1,2,1);ni chols(nu m,de n);grid;B

% n ichols 曲線 G subplot(1,2,2);2 g=tf(nu m,de n);bode(feedback(g,1,-1));grid;Nichols Chart koflr ea M^hnaaM

msuvesaB—

50

--

00

-270 Bode Diagram

90--

% 波特圖 17.系統(tǒng)的開環(huán)傳遞函數(shù)為:

分別確定 k=2 和 k=10 時(shí)閉環(huán)系統(tǒng)的穩(wěn)定性 >> d 仁 [1 , 3,2,0];n 仁 [2];[nc1 , dc1]=cloop(n1 , d1 ,-1);roots(dc1)d2=d1;n2=[10];[nc2 , dc2]=cloop(n2 , d2,-1);roots(dc2)ans =-2.5214-0.2393 + 0.8579i-0.2393-0.8579i ans =-3.3089 0.1545 + 1.7316i 0.1545-1.7316i 18.系統(tǒng)的狀態(tài)方程為：

X 1 4 3 0 X 1 1 X 2 1 0 0 X 2 0 u X 3 0 1 0 X 3 0

X 1

y 0 1 2 x 2

X 3

G(s)k s(s 1)(s 2)

試確定系統(tǒng)的穩(wěn)定性。

>> a=[-4,-3,0;1,0,0;0,1,0];b=[1;0;0];c=[0,1,2];d=0;eig(a)% 求特征根 ran k(ctrb(a,b))ans = 0-1-3 ans = 3 實(shí)驗(yàn)六連續(xù)系統(tǒng)數(shù)字仿真的基本算法 實(shí)驗(yàn)任務(wù) 1.理解歐拉法和龍格-庫塔法的基本思想； 2 ?理解數(shù)值積分算法的計(jì)算精度、速度、穩(wěn)定性與步長(zhǎng)的關(guān)系;程序舉例 1.取 h=0.2 ,試分別用歐拉法、RK2 法和 RK4 法求解微分方程的數(shù)值解，并 比較計(jì)算精度。

注：解析解：y(t).1 2t >> clear t(1)=0 ； y(1)=1;y_euler(1)=1;y_rk2(1)=1;y_rk4(1)=1;h=0.001;% 步長(zhǎng)修改為 0.001 for k=1:5 y(t)y(t)y(0)2t y(t)1

t(k+1)=t(k)+h;y(k+1)=sqrt(1+2*t(k+1));end for k=1:5 y_euler(k+1)=y_euler(k)+h*(y_euler(k)-2*t(k)/y_euler(k));end for k=1:5 k1= y_rk2(k)-2*t(k)/y_rk2(k);k2=(y_rk2(k)+h*k1)-2*(t(k)+h)/(y_rk2(k)+h*k1);y_rk2(k+1)=y_rk2(k)+h*(k1+k2)/2;end for k=1:5 k1= y_rk4(k)-2*t(k)/y_rk4(k);k2=(y_rk4(k)+h*k1/2)-2*(t(k)+h/2)/(y_rk4(k)+h*k1/2);k3=(y_rk4(k)+h*k22)-2*(t(k)+h/2)/(y_rk4(k)+h*k2/2);k4=(y_rk4(k)+h*k3)-2*(t(k)+h)/(y_rk4(k)+h*k3);

y_rk4(k+1)=y_rk4(k)+h*(k1+2*k2+2*k3+k4)/6;end disp(” 時(shí)間 解析解 yt=[t“, y”, y_euler“, y_rk2”, y_rk4“];disp(yt)歐拉法 RK2 法 RK4 法 ”)時(shí)間 解析解 歐拉法 RK2 法 RK4法0 1.0000 1.0000 1.0000 1.0000 0.0010 1.0010 1.0010 1.0010 1.0010 0.0020 1.0020 1.0020 1.0020 1.0020 0.0030 1.0030 1.0030 1.0030 1.0030 0.0040 1.0040 1.0040 1.0040 1.0040 0.0050 1.0050 1.0050 1.0050 1.0050

y(t)2Ry(t)y(t)0 在 0 t 10 上的數(shù)字仿真解(已知：

y(0)y(0)0)，并將不同步長(zhǎng)下的仿真結(jié)果與解析解進(jìn)行精度比較。

說明：

已知該微分方程的解析解分別為：

100，yt y t 100cost(當(dāng) R 0)100e 2t c°s 仝 t 10 ^ e*s in 仝 t 2 3 2(當(dāng) R 0.5)采用 RK4 法進(jìn)行計(jì)算，選擇狀態(tài)變量: 2.考慮如下二階系統(tǒng)：

x 1 y x 2 y 則有如下狀態(tài)空間模型及初值條件 x 1 x 2 x 1(0)100 x 2 x 1 2Rx 2 x 2(0)0 y x 1 采用 RK4 法進(jìn)行計(jì)算。

>> clear h=input(“ 請(qǐng)輸入步長(zhǎng) h=”);M=round(10/h);t(1)=0;y_0(1)=100;y_05(1)=100;和 y_05 分別對(duì)應(yīng)于為 R=0 和 R=0.5)x1(1)=100;x2(1)=0;y_rk4_0(1)=x1(1);y_rk4_05(1)=x1(1);% 求解析解 for k=1:M t(k+1)=t(k)+h;y_0(k+1)=100*cos(t(k+1));y_05(k+1)=100*exp(-t(k+1)/2).*cos(sqrt(3)/2*t(k+1))+100*sqrt(3)/3*exp(-t(k+1)/2).*si n(sqrt(3)/2*t(k+1));end% 輸入步長(zhǎng) % 置總計(jì)算步數(shù) % 置自變量初值 % 置解析解的初始值(y_0 % 置狀態(tài)向量初值 % 置數(shù)值解的初值

% 利用 RK4 法求解 % R=0 for k=1:M k11=x2(k);k12=-x1(k);k21=x2(k)+h*k12/2;k22=-(x1(k)+h*k11/2);k31=x2(k)+h*k22/2;k32=-(x1(k)+h*k21/2);k41=x2(k)+h*k32;k42=-(x1(k)+h*k31);x1(k+1)=x1(k)+h*(k11+2*k21+2*k31+k41)/6;x2(k+1)=x2(k)+h*(k12+2*k22+2*k32+k42)/6;y_rk4_0(k+1)=x1(k+1);end % R=0.5 for k=1:M k11=x2(k);k12=-x1(k)-x2(k);k21=x2(k)+h*k12/2;k22=-(x1(k)+h*k11/2)-(x2(k)+h*k12/2);k31=x2(k)+h*k22/2;k32=-(x1(k)+h*k21/2)-(x2(k)+h*k22/2);k41=x2(k)+h*k32;k42=-(x1(k)+h*k31)-(x2(k)+h*k32);x1(k+1)=x1(k)+h*(k11+2*k21+2*k31+k41)/6;x2(k+1)=x2(k)+h*(k12+2*k22+2*k32+k42)/6;y_rk4_05(k+1)=x1(k+1);end % 求出誤差最大值 err_0=max(abs(y_0-y_rk4_0));err_05=max(abs(y_05-y_rk4_05));% 輸出結(jié)果 disp（“ 最大誤差（R=0）

最大誤差（R=0.5）”）

err_max=[err_0,err_05];disp（err_max）

請(qǐng)輸入步長(zhǎng) h=0.5 最大誤差（R=0）

最大誤差（R=0.5）

0.4299 0.0460 沒運(yùn)行一次程序，輸入一個(gè)步長(zhǎng)，記得可到相應(yīng)的最大誤差，將結(jié)果可列表如 下：

步長(zhǎng) h 0.0001 0.0005 0.001 0.005 0.01 0.05 0.1 0.5

F =0 5.4330 1.6969 1.0574 4.1107 6.6029 4.1439 6.6602 4.2988 差 最大誤 X 10-10 X 10-10 X 10-10 X 10-9 X 10-8 X 10-5 X 10-4 X 10-1

R=0.5 2.7649 6.8123 5.3753 4.0902 6.5425 4.1365 6.7152 4.5976 差 最大誤 X 10-11 X 10-12 X 10-12 X 10-10 X 10-9 X 10-6 X 10-5 X 10-2 從上表中可以看出，當(dāng)步長(zhǎng) h=0.001 時(shí)，總誤差最??；當(dāng)步長(zhǎng) h 小于 0.001 時(shí)，由于舍入誤差變大而使總誤差增加；當(dāng)步長(zhǎng) h 大于 0.001 時(shí)，則由于截?cái)嗾` 差的增加也使得總誤差加大。另外，當(dāng)系統(tǒng)的解變化激烈時(shí)（如 R=0），誤差對(duì) 步長(zhǎng)的變化較為敏感；當(dāng)系統(tǒng)的解變化平穩(wěn)時(shí)，步長(zhǎng)的變化對(duì)誤差的影響就要緩 和得多。數(shù)值積分算法確定以后，在選擇步長(zhǎng)時(shí)，需要綜合考慮。

第三篇：化工原理實(shí)驗(yàn)報(bào)告封面格式

課程名稱： 化工原理實(shí)驗(yàn) 課題

系院：

專業(yè)：

實(shí)驗(yàn)人員：

學(xué)號(hào)：

同組人員：

實(shí)驗(yàn)時(shí)間：

指導(dǎo)教師：流體流動(dòng)阻力的測(cè)定生物與化學(xué)工程系11級(jí)應(yīng)用化學(xué)朱穎妍1***

47、戴軍凡

49、陳謝萍

51、徐曉麗

52、丁忱2014年2 月 25 日楊金杯

第四篇：化工原理實(shí)驗(yàn)報(bào)告評(píng)分標(biāo)準(zhǔn)

化工原理實(shí)驗(yàn)報(bào)告評(píng)分標(biāo)準(zhǔn)

化工原理實(shí)驗(yàn)報(bào)告評(píng)分分為五部分：

預(yù)習(xí)20%，考勤10%，實(shí)驗(yàn)操作10%，實(shí)驗(yàn)數(shù)據(jù)處理40%，思考題20%。

1.預(yù)習(xí)評(píng)分標(biāo)準(zhǔn)

進(jìn)實(shí)驗(yàn)室之前，仔細(xì)閱讀實(shí)驗(yàn)講義相關(guān)內(nèi)容，認(rèn)真寫好實(shí)驗(yàn)預(yù)習(xí)報(bào)告。完整的預(yù)習(xí)報(bào)告應(yīng)包括如下內(nèi)容：實(shí)驗(yàn)名稱、實(shí)驗(yàn)班級(jí)、姓名、學(xué)號(hào)以及實(shí)驗(yàn)?zāi)康?、?shí)驗(yàn)內(nèi)容、基本原理、設(shè)備流程(設(shè)備工藝流程圖)、實(shí)驗(yàn)步驟及注意事項(xiàng)等。缺內(nèi)容者酌情扣分。

2.考勤評(píng)分標(biāo)準(zhǔn)

按時(shí)到實(shí)驗(yàn)室簽到，著裝整齊者為10分。每遲到一分鐘扣一分。

3.實(shí)驗(yàn)操作評(píng)分標(biāo)準(zhǔn)

正確操作實(shí)驗(yàn)裝置者得10分。

4.實(shí)驗(yàn)數(shù)據(jù)處理評(píng)分標(biāo)準(zhǔn)

實(shí)驗(yàn)數(shù)據(jù)及中間處理過程表格化，每一計(jì)算步驟需要有詳細(xì)的計(jì)算舉例，計(jì)算公式正確，結(jié)果合理40分。

5.思考題評(píng)分標(biāo)準(zhǔn)

結(jié)合實(shí)驗(yàn)，認(rèn)真完整地回答實(shí)驗(yàn)講義后面的思考題，且回答正確20分。

第五篇：化工原理仿真系統(tǒng)研究論文

一、化工原理仿真系統(tǒng)的制作

化工原理實(shí)驗(yàn)包括流體流動(dòng)阻力測(cè)定、離心泵性能測(cè)定、傳熱、精餾、吸收與解吸、干燥、萃取等基本單元操作，分別由不同的儀器儀表和管道組合而成。在仿真軟件中，把各種設(shè)備和管道用flash畫出，再根據(jù)每一套裝置流程圖的要求，以真實(shí)、立體的效果來實(shí)現(xiàn)。

1.整體結(jié)構(gòu)。實(shí)驗(yàn)仿真系統(tǒng)的開發(fā)過程分為三個(gè)階段：實(shí)驗(yàn)前的準(zhǔn)備、實(shí)驗(yàn)過程及數(shù)據(jù)記錄和數(shù)據(jù)處理。前兩個(gè)階段在Flash動(dòng)畫制作軟件上完成，第三階段在VisualStudio2005軟件開發(fā)工具上完成，并且使用Ac-cess數(shù)據(jù)庫進(jìn)行數(shù)據(jù)的存儲(chǔ)與交換。

2.仿真系統(tǒng)的實(shí)現(xiàn)。在計(jì)算機(jī)模擬化工原理實(shí)驗(yàn)時(shí)，需要通過動(dòng)態(tài)數(shù)學(xué)模型來模擬真實(shí)的實(shí)驗(yàn)操作，該模型主要包括實(shí)驗(yàn)指導(dǎo)、素材演示、仿真操作、數(shù)據(jù)處理、考題測(cè)試、幫助功能等內(nèi)容。下面以離心泵性能測(cè)定為例詳細(xì)說明仿真系統(tǒng)的制作過程。在實(shí)驗(yàn)準(zhǔn)備階段與實(shí)驗(yàn)開始階段的Flash動(dòng)畫的制作過程中，考慮到實(shí)驗(yàn)步驟有先后，以及更好地做到人機(jī)交互，必須使用專門為Flash腳本開發(fā)的ActionScript語言。如點(diǎn)擊水泵開啟按鈕必須在閥門開啟以后才能啟動(dòng)，直至水灌滿后，才可以點(diǎn)擊關(guān)閉水泵按鈕。為了使實(shí)驗(yàn)更具有真實(shí)性，需設(shè)置閥門的流量控制，分為10個(gè)級(jí)別，可以逐漸增大或減小。運(yùn)用VisualStudio.Net開發(fā)環(huán)境編寫C#程序，可以通過拖動(dòng)添加組件，并自動(dòng)生成組件需要的代碼。在制作化工原理實(shí)驗(yàn)?zāi)M課件時(shí)，可通過VisualStudio屬性窗口設(shè)置各種開發(fā)元素屬性如外觀、名稱等，且屬性窗口中顯示的內(nèi)容，隨著選擇開發(fā)元素的不同而動(dòng)態(tài)改變。利用VisualStudi“o工具箱”，可以向應(yīng)用程序添加標(biāo)準(zhǔn)控件。在設(shè)置好窗體和控件后，利用Vi-sualStudio的代碼編輯器編寫程序代碼。在命令窗口中，可以直接輸入并執(zhí)行各種命令，調(diào)試應(yīng)用程序，并通過在即時(shí)窗口的命令行中輸入表達(dá)式或變量名，可以得到它們的值。編寫程序過程中，難免會(huì)遇到一些錯(cuò)誤，開發(fā)人員需要對(duì)應(yīng)用程序進(jìn)行調(diào)試，查找錯(cuò)誤的根源，以期達(dá)到設(shè)計(jì)要求。離心泵性能測(cè)定實(shí)驗(yàn)涉及到流體流動(dòng)、水泵運(yùn)轉(zhuǎn)、儀表變化、閥門打開或關(guān)閉等動(dòng)作，在仿真系統(tǒng)中通過Flash動(dòng)畫來實(shí)現(xiàn)這些動(dòng)作的動(dòng)態(tài)效果，使整個(gè)實(shí)驗(yàn)過程表現(xiàn)得更加真實(shí)。用Flash中的按鈕實(shí)現(xiàn)動(dòng)畫交互效果，控制整個(gè)實(shí)驗(yàn)的操作并對(duì)數(shù)據(jù)進(jìn)行采集，同時(shí)將數(shù)據(jù)傳入C#，由C#對(duì)數(shù)據(jù)庫進(jìn)行讀寫操作，然后作出離心泵特性曲線圖。

3.實(shí)驗(yàn)數(shù)據(jù)產(chǎn)生及處理。化工原理實(shí)驗(yàn)過程中往往要測(cè)定溫度、壓強(qiáng)、濃度、流速等數(shù)據(jù)，同時(shí)必須對(duì)這些參數(shù)進(jìn)行整理和分析，并運(yùn)用相關(guān)的理論公式進(jìn)行計(jì)算，才能達(dá)到實(shí)驗(yàn)預(yù)期目的?；ぴ韺?shí)驗(yàn)實(shí)測(cè)數(shù)據(jù)多，繪圖耗時(shí)費(fèi)力，計(jì)算公式復(fù)雜，有時(shí)甚至需要進(jìn)行迭代計(jì)算，借助計(jì)算機(jī)輔助程序可圓滿解決這些問題。在仿真軟件中，通過C#語言設(shè)計(jì)數(shù)據(jù)處理程序。根據(jù)各化工單元操作理論建立數(shù)學(xué)模型，使仿真數(shù)據(jù)在實(shí)際操作的數(shù)據(jù)范圍內(nèi)隨機(jī)產(chǎn)生，以保證每個(gè)學(xué)生在進(jìn)行仿真實(shí)驗(yàn)時(shí)即使初始條件相同，也不會(huì)得到完全相同的實(shí)驗(yàn)結(jié)果，更接近真實(shí)操作狀況。試驗(yàn)完成后，點(diǎn)擊“記錄數(shù)據(jù)”按鈕，計(jì)算機(jī)會(huì)自動(dòng)記錄數(shù)據(jù)，并在后臺(tái)進(jìn)行數(shù)據(jù)傳遞，然后根據(jù)預(yù)先輸入的計(jì)算公式進(jìn)行數(shù)據(jù)處理。數(shù)據(jù)處理后被保存到Access數(shù)據(jù)庫中，再通過調(diào)用，將數(shù)據(jù)輸出在DataGridView進(jìn)行顯示，或據(jù)此數(shù)據(jù)繪制實(shí)驗(yàn)曲線。

二、操作過程及功能概述

主界面使用VisualStudioC#中的窗體，通過添加菜單欄來控制試驗(yàn)的選擇。其特點(diǎn)是方便、簡(jiǎn)單易用，更重要的是為今后仿真系統(tǒng)的逐步完善提供了空間。首先水泵的開關(guān)按鈕是不可用的，必須在打開閥門以后，才能啟動(dòng)水泵。當(dāng)水泵與閥門同時(shí)啟動(dòng)后，便開始灌水，在這期間禁用系統(tǒng)中所有的按鈕。待灌水過程結(jié)束，先關(guān)閥門再關(guān)水泵。點(diǎn)擊“開始實(shí)驗(yàn)”按鈕，可以開啟下一個(gè)界面繼續(xù)實(shí)驗(yàn)。先打開水泵，然后打開閥門，通過閥門調(diào)節(jié)流量級(jí)別，儀表數(shù)值會(huì)隨之變化。點(diǎn)擊“記錄數(shù)據(jù)”按鈕，將儀表的數(shù)值記錄在數(shù)據(jù)庫中。當(dāng)數(shù)據(jù)記錄完畢，點(diǎn)擊“查看數(shù)據(jù)”按鈕，屏幕上顯示10組數(shù)據(jù)以及由公式計(jì)算得出的“揚(yáng)程”、“有效功率”、“效率”數(shù)值。點(diǎn)擊“繪圖”按鈕，可直接繪制出H-Q、P-Q及η-Q三條特性曲線。無論是實(shí)驗(yàn)結(jié)束還是中途關(guān)閉實(shí)驗(yàn)窗體，都將出現(xiàn)一個(gè)對(duì)話框以提示實(shí)驗(yàn)者“是否保存當(dāng)前數(shù)據(jù)？”操作者可根據(jù)提示對(duì)實(shí)驗(yàn)數(shù)據(jù)進(jìn)行取舍。集合Flash動(dòng)畫和C#語言優(yōu)點(diǎn)開發(fā)的化工原理實(shí)驗(yàn)仿真系統(tǒng)，具有界面直觀、操作簡(jiǎn)單、支持人機(jī)交互、占用空間小等特點(diǎn)，能顯著提高化工原理實(shí)驗(yàn)教學(xué)的效果，減少實(shí)驗(yàn)設(shè)備投資和損壞，降低實(shí)驗(yàn)投入成本，避免實(shí)驗(yàn)事故的發(fā)生。計(jì)算機(jī)輔助教學(xué)，特別是計(jì)算機(jī)仿真實(shí)驗(yàn)在化工教學(xué)過程中的應(yīng)用，使學(xué)生接觸了一種全新的實(shí)驗(yàn)手段，激發(fā)了學(xué)生學(xué)習(xí)的積極性和主動(dòng)性，使學(xué)生創(chuàng)新意識(shí)得到培養(yǎng)，從而提高了整體教學(xué)質(zhì)量。