計算機(jī)仿真試驗(yàn)報告采樣系統(tǒng)仿真

1、博/ 大挈課程名稱：計算機(jī)仿真實(shí)驗(yàn)項目：采樣控制系統(tǒng)數(shù)字仿真專業(yè)班級：姓 名：學(xué) 號：實(shí)驗(yàn)室號：實(shí)驗(yàn)組號：實(shí)驗(yàn)時間：2016.04.27批閱時間： 指導(dǎo)教師：成 績：沈陽工業(yè)大學(xué)實(shí)驗(yàn)報告（適用計算機(jī)程序設(shè)計類）專業(yè)班級：_ 學(xué)號： 姓名： 實(shí)驗(yàn)名稱：采樣控制系統(tǒng)數(shù)字仿真.實(shí)驗(yàn)?zāi)康模菏煜げ蓸涌刂葡到y(tǒng)數(shù)值積分法仿真流程，程序編制及仿真結(jié)果分析。.實(shí)驗(yàn)內(nèi)容：（1）學(xué)習(xí)MATLAB語言實(shí)現(xiàn)采樣控制系統(tǒng)仿真建模。（2）學(xué)習(xí)使用仿真程序進(jìn)行系統(tǒng)分析。.實(shí)驗(yàn)方案（程序設(shè)計說明）建立一個采樣系統(tǒng)，被控對象的狀態(tài)空間模型，在階躍輸入信號 下進(jìn)行采樣處理，并且觀察采樣后的曲線和采樣前的曲線，進(jìn)行對比, 4.實(shí)驗(yàn)步

2、驟或程序（經(jīng)調(diào)試后正確的源程序）（填寫主要步驟與程序代碼等。見附件 A）.程序運(yùn)行結(jié)果（結(jié)果見附件A）.出現(xiàn)的問題及解決方法沒有理解怎么用離散相似法先對式子進(jìn)行離散化處理，在程序的 設(shè)計和運(yùn)用上遇到困難。通過查閱網(wǎng)上資料和相關(guān)書籍，了解并學(xué)會相應(yīng)的知識，最后成 功運(yùn)用到這次采樣中。附件A沈陽工業(yè)大學(xué)實(shí)驗(yàn)報告（適用計算機(jī)程序設(shè)計類）專業(yè)班級： 學(xué)號： 姓名： 實(shí)驗(yàn)步驟或程序：例1：被控對象的狀態(tài)空間模型為x1 =ux2 =x1-x2y=x2卜的狀態(tài)響應(yīng)。采樣周期為T=0.1so初始狀態(tài)x1（0）=x2（0）=0 ,在階躍輸入信號 4）=1作用（1）程序如下:G=1 0;0.01 0.99;H=0

3、.01;0;c=0 1;h=0.01;T=0.1;N=round(T/h);TF=15;%在一個采樣周期內(nèi)，連續(xù)對被控對象的離散化狀態(tài)方程遞推N次（內(nèi)循環(huán)）M=round(TF/T);x=0;0;u=1;t=0;%數(shù)字控制器的差分方程遞推 （外循環(huán)）xt=x;for i=1:My=c*x;r(i)=1;e(i)=r(i)-y;u=0.5*u+e(i);for j=1:N x=G*x+H*u;xt=xt,x;%t內(nèi)循環(huán)內(nèi)循環(huán)中遞推求受控對象狀態(tài)響應(yīng)（步長為為記載各采樣（KT）時刻的行向量n)t=t,j*h+(i-1)*T;end endx1t=xt(1,:);x2t=xt(2,:);plot(t

4、,x1t, :k,t,x2t,k);legend(x1, x2);grid;xlabel(time(s);(2)得到的結(jié)果：num=0.005,0.005,0;%den=1,-2.4,1.863,-0.453; % yk,x,n=dstep(num,den); %yk%T=0.1;%for k=1:n如圖為采樣后的結(jié)果，我們看到采樣后的曲線比較真實(shí)曲線存在誤差。在10s范圍內(nèi)誤差誤差較大，振幅比真實(shí)曲線大。隨著時間的推移，誤差逐漸減小，最后趨于穩(wěn)定。例2. 給出返回系統(tǒng)輸出階躍響應(yīng)的離散數(shù)據(jù)，在單位階躍響應(yīng)下的采樣。（1）程序如下： clear all脈沖傳遞函數(shù)分子多項式按z的降募系數(shù)排列的行向量脈沖傳遞函數(shù)分母多項式按z的降募系數(shù)排列的行向量為存放在輸出離散序列的數(shù)組，n為dstep函數(shù)自動設(shè)定的采樣點(diǎn)已知系統(tǒng)采樣周期為 0.1splot(k*T,yk(k),*k）; %k為采樣序列號，k*T為第k次采樣對應(yīng)的時刻hold onendxlabel( 時間(s),);grid時間(S)(3)系統(tǒng)