李 杰

(承德石油高等?？茖W(xué)校熱能工程系，河北承德 067000)

數(shù)值積分法［1］發(fā)展到今天，已有許多種擬合手段，迄今為止，應(yīng)用最廣泛的當(dāng)屬四階龍格－庫(kù)塔法［2］，并以變步長(zhǎng)的龍格庫(kù)塔法為優(yōu)。但鑒于其計(jì)算量和計(jì)算精度有待進(jìn)一步提高，有必要進(jìn)一步尋找出最佳步長(zhǎng)，以便在不無(wú)限地減小步長(zhǎng)條件下達(dá)到較高的精度，從而得出一個(gè)較可靠的步長(zhǎng)值，供分析實(shí)際工程系統(tǒng)的穩(wěn)定性作參考使用，對(duì)此本文做了大量的實(shí)驗(yàn)數(shù)據(jù)及實(shí)驗(yàn)圖像，進(jìn)而分析得出了最佳仿真時(shí)間與仿真步長(zhǎng)［2］。

1 實(shí)驗(yàn)原理及仿真步驟

1.1 實(shí)驗(yàn)原理

運(yùn)用matlab［3］中的命令語(yǔ)句ctrb(a，b)與obsv(a，c)，分別對(duì)其用命令語(yǔ)句rank取秩，以此判斷其可控性與可觀察性。

1.2 仿真步驟

編制RK-4仿真程序，繪制出單輸入單輸出及多輸入多輸出曲線圖［4］。仿真程序流程如圖1所示:

1.3 軟件程序代碼

仿真用原程序代碼清單如下:

2 實(shí)驗(yàn)仿真結(jié)果

在實(shí)驗(yàn)中，選取了大量不同的仿真時(shí)間tf以及步長(zhǎng)h，從中挑選出如下幾種最具代表性的實(shí)驗(yàn)圖像(圖2—圖11)。

由圖2、圖3可以看出當(dāng)仿真時(shí)間tf=2，仿真步長(zhǎng)分別為:h=0.05，h=0.01時(shí)，波形明顯變緩，過渡過程時(shí)間明顯加長(zhǎng)，系統(tǒng)達(dá)到穩(wěn)定時(shí)所用時(shí)間過長(zhǎng)，嚴(yán)重影響控制系統(tǒng)的實(shí)時(shí)性。

由圖4可以看出當(dāng)仿真時(shí)間tf=10，仿真步長(zhǎng)h=0.5時(shí)，波形開始變陡，過渡過程時(shí)間縮短，但波形明顯失真，不利于控制系統(tǒng)的穩(wěn)定性。由圖5可以看出，隨著仿真步長(zhǎng)的減小，當(dāng)h=0.02時(shí)，波形明顯無(wú)失真。

由圖6可以看出當(dāng)仿真時(shí)間tf=20，仿真步長(zhǎng)h=0.005時(shí)，波形開始變緩，過渡過程時(shí)間明顯縮短，波形明顯變好，有利于控制系統(tǒng)的穩(wěn)定性，并可提高系統(tǒng)的實(shí)時(shí)性。由圖7可以看出，隨著仿真步長(zhǎng)的增大，tf=30，步長(zhǎng)h=0.005時(shí)，波形進(jìn)一步變陡，可縮短系統(tǒng)達(dá)到穩(wěn)態(tài)時(shí)所用時(shí)間，但不利于控制系統(tǒng)的穩(wěn)定性。

由圖8可以看出當(dāng)仿真時(shí)間tf=60，仿真步長(zhǎng)h=0.5時(shí)，仿真速度明顯提高，過渡過程時(shí)間明顯縮短，但波形明顯變陡，部分圖形失真，仿真精度降低，不利于控制系統(tǒng)的穩(wěn)定性。由圖9可以看出，隨著仿真步長(zhǎng)的減小，步長(zhǎng)h=0.05時(shí)，波形無(wú)失真，雖然可提高仿真的精度，但波形明顯變陡，同樣不利于控制系統(tǒng)的穩(wěn)定性。

由圖10，圖11明顯看出，隨著仿真步長(zhǎng)的增大，當(dāng)仿真步長(zhǎng)h=1時(shí)，不論仿真時(shí)間增大還是減小，波形嚴(yán)重失真，控制系統(tǒng)變的不穩(wěn)定性，難以控制，更談不上仿真的速度與精度。

3 結(jié)論

本文通過大量的仿真圖形與數(shù)據(jù)分析得出結(jié)論:當(dāng)仿真時(shí)間tf大于20秒時(shí)波形明顯變陡，但當(dāng)仿真時(shí)間tf小于20秒時(shí)波形明顯變緩，但過渡過程時(shí)間明顯加長(zhǎng)，綜合分析最佳仿真時(shí)間為tf=20，當(dāng)步長(zhǎng)為h等于0.2時(shí)波形開始失真，等于1時(shí)算法不穩(wěn)定，當(dāng)?shù)扔?.001時(shí)運(yùn)行時(shí)間明顯加長(zhǎng)，當(dāng)步長(zhǎng)等于0.005時(shí)為最佳步長(zhǎng)，同時(shí)把本實(shí)驗(yàn)得出的結(jié)果應(yīng)用于實(shí)際工程，結(jié)果較為理想。所以得出:最佳仿真時(shí)間tf=20，最佳仿真步長(zhǎng)h=0.005，當(dāng)h=0.2時(shí)波形開始失真h=1算法不穩(wěn)定。

［1］劉廷建.數(shù)值計(jì)算方法［M］.成都:成都科技大學(xué)出版社，1997.

［2］沈艷霞，周平.MATLAB在電力電子學(xué)仿真中的應(yīng)用［J］.無(wú)錫輕工大學(xué)學(xué)報(bào)，2001，20(1):92～94.

［3］劉衛(wèi)國(guó)，陳昭平，張穎.MATLAB程序設(shè)計(jì)與應(yīng)用［M］.北京:高等教育出版社，2002.

［4］鄧國(guó)揚(yáng)，盛義發(fā).基于Matlab Simulink的電力電子系統(tǒng)的建模與仿真［J］.南華大學(xué)學(xué)報(bào)(理工版)，2003，17(1):1～6.