計算機仿真試驗報告采樣系統(tǒng)仿真

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

2、驟或程序（經(jīng)調(diào)試后正確的源程序）（填寫主要步驟與程序代碼等。見附件 A）.程序運行結果（結果見附件A）.出現(xiàn)的問題及解決方法沒有理解怎么用離散相似法先對式子進行離散化處理，在程序的 設計和運用上遇到困難。通過查閱網(wǎng)上資料和相關書籍，了解并學會相應的知識，最后成 功運用到這次采樣中。附件A沈陽工業(yè)大學實驗報告（適用計算機程序設計類）專業(yè)班級： 學號： 姓名： 實驗步驟或程序：例1：被控對象的狀態(tài)空間模型為x1 =ux2 =x1-x2y=x2卜的狀態(tài)響應。采樣周期為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)響應（步長為為記載各采樣（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)得到的結果：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如圖為采樣后的結果，我們看到采樣后的曲線比較真實曲線存在誤差。在10s范圍內(nèi)誤差誤差較大，振幅比真實曲線大。隨著時間的推移，誤差逐漸減小，最后趨于穩(wěn)定。例2. 給出返回系統(tǒng)輸出階躍響應的離散數(shù)據(jù)，在單位階躍響應下的采樣。（1）程序如下： clear all脈沖傳遞函數(shù)分子多項式按z的降募系數(shù)排列的行向量脈沖傳遞函數(shù)分母多項式按z的降募系數(shù)排列的行向量為存放在輸出離散序列的數(shù)組，n為dstep函數(shù)自動設定的采樣點已知系統(tǒng)采樣周期為 0.1splot(k*T,yk(k),*k）; %k為采樣序列號，k*T為第k次采樣對應的時刻hold onendxlabel( 時間(s),);grid時間(S)(3)系統(tǒng)