鐘曉蕓
(廣西師范大學(xué)電子與信息工程學(xué)院,廣西 桂林 541004)
信號(hào)與系統(tǒng)課程是電子信息類專業(yè)必修的重要專業(yè)課程之一,高等數(shù)學(xué)、電路和復(fù)變函數(shù)等是它的先修課程,而其后續(xù)課程則有數(shù)字信號(hào)處理和自動(dòng)控制原理等[1-2]。該課程既是信號(hào)處理領(lǐng)域的入門(mén)課程,也包含后續(xù)課程學(xué)習(xí)及學(xué)生考取研究生的重要基礎(chǔ)理論,其特點(diǎn)是概念抽象較難理解,數(shù)學(xué)推導(dǎo)復(fù)雜較難掌握,結(jié)果靈活較難運(yùn)用。該課程的理論性和實(shí)踐性都較強(qiáng),其課程目標(biāo)主要是使學(xué)生獲得連續(xù)系統(tǒng)的時(shí)域分析和s域分析、離散系統(tǒng)的時(shí)域分析、傅里葉變換和系統(tǒng)的頻域分析、拉普拉斯變換、系統(tǒng)狀態(tài)動(dòng)力學(xué)分析等,同時(shí)提高學(xué)生的實(shí)踐能力??紤]到該課程涉及較多的現(xiàn)代數(shù)學(xué)知識(shí),概念抽象,內(nèi)容靈活,沒(méi)有一定的養(yǎng)成訓(xùn)練和實(shí)驗(yàn)試探,相當(dāng)多的同學(xué)會(huì)對(duì)其原理和方法感到難以正確把握,特別是對(duì)于課程中諸如時(shí)域分析和頻域分析以及系統(tǒng)動(dòng)力學(xué)特征等重要的知識(shí)點(diǎn)和較難掌握的內(nèi)容,僅通過(guò)教師傳統(tǒng)課堂的教學(xué),容易造成學(xué)生們不愿意手工求解,或者對(duì)原理的運(yùn)用不得要領(lǐng)。
為此,文章在課程實(shí)驗(yàn)的教學(xué)環(huán)節(jié)中有針對(duì)性地結(jié)合若干經(jīng)過(guò)簡(jiǎn)化的科研前沿問(wèn)題,引入若干Matlab仿真實(shí)驗(yàn),引導(dǎo)學(xué)生通過(guò)計(jì)算機(jī)編程完成數(shù)值仿真,以便可視化地演示課程的一些基本原理及方法,達(dá)到科研實(shí)踐與課程教學(xué)相互融合,增強(qiáng)學(xué)生在計(jì)算機(jī)軟件的配合下正確把握課本知識(shí),適當(dāng)開(kāi)展一定的科學(xué)探索的目的。實(shí)踐表明,由于結(jié)合一些適當(dāng)簡(jiǎn)化的科研前沿問(wèn)題,加上Matlab程序設(shè)計(jì)語(yǔ)言非常簡(jiǎn)潔而高效,可以很快掌握,極大調(diào)動(dòng)了學(xué)生自己進(jìn)行課程實(shí)驗(yàn)的積極性。大家主動(dòng)參與實(shí)驗(yàn)、自己試著編寫(xiě)M程序?qū)崿F(xiàn)對(duì)所學(xué)基本原理和方法進(jìn)行實(shí)驗(yàn)驗(yàn)證,包括對(duì)于文獻(xiàn)上的一些可類比的新結(jié)果進(jìn)行檢驗(yàn),有效地促進(jìn)了學(xué)生創(chuàng)新能力的提升??陀^上,由于Matlab軟件不僅運(yùn)算能力強(qiáng)大,而且計(jì)算精度特別高,其數(shù)值計(jì)算誤差可以調(diào)節(jié)。這些特點(diǎn)使得Matlab 仿真結(jié)果能夠真實(shí)地反映與理論分析的一致性,相比傳統(tǒng)的課程實(shí)驗(yàn)教學(xué)具有明顯的創(chuàng)新優(yōu)勢(shì)。通過(guò)這些實(shí)驗(yàn),將有助于學(xué)生理解和運(yùn)用現(xiàn)代信號(hào)與系統(tǒng)課程的基本原理和方法,獲得利用Matlab仿真解決現(xiàn)代分?jǐn)?shù)階系統(tǒng)的初步科研訓(xùn)練。
基于科研實(shí)踐與課程教學(xué)相互融合理念[3],結(jié)合Matlab仿真開(kāi)展信號(hào)與系統(tǒng)課程實(shí)驗(yàn)的設(shè)計(jì),有效打破了以在實(shí)驗(yàn)中大量使用硬件和常規(guī)問(wèn)題為特征的傳統(tǒng)實(shí)驗(yàn)教學(xué)的局限性,實(shí)驗(yàn)項(xiàng)目設(shè)計(jì)可以很容易地觸及科研前沿內(nèi)容,同時(shí)聯(lián)系到連續(xù)或離散系統(tǒng)的時(shí)域仿真和頻域分析等理論課程中的難點(diǎn)和重點(diǎn),學(xué)生們通過(guò)編寫(xiě)M程序并自己上機(jī)進(jìn)行測(cè)試分析,不僅增加了實(shí)驗(yàn)設(shè)計(jì)的前沿性和實(shí)驗(yàn)項(xiàng)目選擇的靈活性,而且使得信號(hào)與系統(tǒng)的計(jì)算可以直觀演示。這樣,在開(kāi)展信號(hào)與系統(tǒng)課堂理論學(xué)習(xí)的同時(shí),通過(guò)計(jì)算機(jī)的配合,不僅學(xué)會(huì)Matlab 基本的編程技巧和實(shí)現(xiàn)方法,而且嘗試將體現(xiàn)課程基本原理的重難點(diǎn)習(xí)題及通過(guò)將課本中的一些連續(xù)信號(hào)系統(tǒng)的狀態(tài)方程進(jìn)行分?jǐn)?shù)階推廣。具體來(lái)說(shuō),我們結(jié)合若干課程實(shí)驗(yàn)內(nèi)容的分?jǐn)?shù)階推廣并且借用Matlab 程序和仿真演示,給學(xué)生布置一些應(yīng)用Matlab語(yǔ)言的課程習(xí)題作為實(shí)驗(yàn)作業(yè),不僅有利于提高學(xué)生的學(xué)習(xí)興趣,培養(yǎng)學(xué)生綜合運(yùn)用計(jì)算機(jī)輔助分析實(shí)際控制系統(tǒng)的科研創(chuàng)新能力,而且有利于加深對(duì)信號(hào)與系統(tǒng)課程的基本原理、基本方法的理解。
在課程實(shí)驗(yàn)中應(yīng)用Matlab 進(jìn)行分?jǐn)?shù)階的信號(hào)和系統(tǒng)的可視化繪圖,主要使用到plot,stairs,stem,jieyue,即階躍函數(shù)等,以及一些簡(jiǎn)單的拉普拉斯變換及其逆變換的Matlab 求解命令laplace,ilaplace 的分?jǐn)?shù)階實(shí)現(xiàn)。
與經(jīng)典的整數(shù)階系統(tǒng)類似,簡(jiǎn)單分?jǐn)?shù)階非線性系統(tǒng)的穩(wěn)定性分析方法最主要的也是根軌跡法和李雅普諾夫方法。根軌跡法主要通過(guò)雅可比迭代法分解法計(jì)算系統(tǒng)雅可比矩陣A的特征值λ(A),并且根據(jù)其分布情況判定系統(tǒng)的穩(wěn)定性,如果涉及導(dǎo)數(shù)的階為分?jǐn)?shù)階α∈(0,1],且arg(λ(A)) >πα/2 時(shí),可判定系統(tǒng)是穩(wěn)定的。李雅普諾夫方法主要通過(guò)判別滿足PA+ATP=-Q的矩陣P的性質(zhì)判斷系統(tǒng)是否穩(wěn)定。這些分析方法涉及Matlab 函數(shù)主要有:時(shí)域分析函數(shù),比如,impulse.m、step.m 和lsim.m。頻域分析分析函數(shù),比如bode.m、nyquist.m和freqs.m等。其中,impulse 函數(shù)和step函數(shù)來(lái)求解LTI系統(tǒng)的沖激響應(yīng)和階躍響應(yīng)數(shù)值解,而lsim 函數(shù)來(lái)對(duì)LTI 系統(tǒng)的零狀態(tài)響應(yīng)進(jìn)行數(shù)值仿真。step(sys)函數(shù)既可以求解系統(tǒng)的單位階躍響應(yīng),也可以繪制響應(yīng)波形圖,impulse(sys)函數(shù)既可以求解系統(tǒng)的單位沖激響應(yīng),也可以繪制響應(yīng)的圖形,其他M 函數(shù)還有mtimes.m、plus.m 和feedback.m等,還可以利用Simulink,通過(guò)仿真求解。
分?jǐn)?shù)階微積分誕生于三百年前,主要研究任意階次的微分和積分的特性及其應(yīng)用問(wèn)題。分?jǐn)?shù)階控制系統(tǒng)已成為現(xiàn)代非線性控制理論的一個(gè)不可或缺的部分[4-5]。隨著分?jǐn)?shù)階系統(tǒng)在現(xiàn)代實(shí)際工程領(lǐng)域的作用日益凸顯,目前,分?jǐn)?shù)階非線性系統(tǒng)的研究已經(jīng)引起了越來(lái)越多的學(xué)者的興趣。其主要原因在于分?jǐn)?shù)階系統(tǒng)自身具有的特性,即分?jǐn)?shù)階系統(tǒng)具有記憶特性,這是整數(shù)階系統(tǒng)所沒(méi)有的,所以研究含有分?jǐn)?shù)階微分或積分環(huán)節(jié)的控制系統(tǒng)在工程應(yīng)用中愈發(fā)重要。
本文以一個(gè)經(jīng)簡(jiǎn)化后的工程問(wèn)題中的分?jǐn)?shù)階的均勻RC線連接的電壓源模擬電路[5]的波動(dòng)為例,說(shuō)明基于Matlab的教學(xué)實(shí)驗(yàn)的具體設(shè)計(jì)內(nèi)容。圖1顯示的電路系統(tǒng)表示了電壓源的內(nèi)部傳統(tǒng)電感L= 1H與一個(gè)非常長(zhǎng)的(l≈∞)均勻RC傳輸線的連接。該電路的波動(dòng)特征可由以下簡(jiǎn)單的微分方程表示-σ0u(t) +σ0e(t),其中,σ0=L-1C-1R, 0 <α≤2。如果α= 1.5,相應(yīng)的系統(tǒng)的動(dòng)力學(xué)特征在文獻(xiàn)[5]中進(jìn)行了詳細(xì)研究,如果系統(tǒng)的階α= 1.0,相應(yīng)的整數(shù)階系統(tǒng)就是熟悉的經(jīng)典RC電路模型。在σ0= 0.1,e(t) = 10e-t的情形,由此實(shí)現(xiàn)傳統(tǒng)的教材內(nèi)容與前沿的分?jǐn)?shù)階系統(tǒng)研究適當(dāng)結(jié)合起來(lái)。
圖1 仿真系統(tǒng)電路圖
通過(guò)如下的Matlab 語(yǔ)言:
可以得到均勻RC線連接的電壓源模擬電路的仿真結(jié)果如圖2(a) 所示。改變系統(tǒng)的設(shè)置,當(dāng)σ0= 0.1,e(t) =cos(0.4t)時(shí),該電壓源模擬電路的時(shí)間波形如圖2(b)所示。
圖2 電路系統(tǒng)的動(dòng)力學(xué)特征仿真波形圖
本課程實(shí)驗(yàn)中電路系統(tǒng)的動(dòng)力學(xué)特征仿真在計(jì)算機(jī)上模擬進(jìn)行,同時(shí)將表征系統(tǒng)特征的微分方程的階處理為一個(gè)在0到1之間取值的分?jǐn)?shù)參數(shù)進(jìn)行實(shí)驗(yàn)結(jié)果分析。參加實(shí)驗(yàn)的學(xué)生在教師的指導(dǎo)下,針對(duì)不同的分?jǐn)?shù)階取值選擇,通過(guò)仿真波形圖直觀地發(fā)現(xiàn)分?jǐn)?shù)階的改變對(duì)于電路特征的影響。實(shí)驗(yàn)過(guò)程中,通過(guò)采用圖1所示的仿真系統(tǒng)電路盡可能地降低實(shí)驗(yàn)復(fù)雜度,忽略了實(shí)驗(yàn)環(huán)境噪聲以及隨機(jī)性因素的影響,體現(xiàn)在系統(tǒng)的微分方程模型中,方程的右邊關(guān)于電路狀態(tài)變量的函數(shù)是簡(jiǎn)單的一次函數(shù)。鼓勵(lì)學(xué)生在完成實(shí)驗(yàn)后,探索其他的方法進(jìn)一步提高仿真精度,比如考慮對(duì)應(yīng)的隨機(jī)系統(tǒng)或者將高斯白噪聲考慮進(jìn)去,還可以進(jìn)一步考慮模糊控制。結(jié)合實(shí)驗(yàn)教學(xué)讓學(xué)生注意,在復(fù)雜而艱辛的科學(xué)研究過(guò)程中,對(duì)科學(xué)問(wèn)題進(jìn)行適當(dāng)?shù)暮?jiǎn)化無(wú)論是對(duì)于成功獲得理論研究成果還是有效獲得實(shí)驗(yàn)現(xiàn)象都是非常必要的。學(xué)生可以像科學(xué)家那樣,結(jié)合自己的學(xué)習(xí)專業(yè)及自己感興趣的方向,選擇自己比較熟悉的課題,在一定的理論分析和結(jié)果預(yù)測(cè)的基礎(chǔ)上,通過(guò)編寫(xiě)計(jì)算機(jī)程序,比如M 程序,進(jìn)行計(jì)算機(jī)仿真研究。嘗試按照科研論文的寫(xiě)作要求和格式,把課程的實(shí)驗(yàn)報(bào)告作為一篇基本符合投稿要求的教學(xué)科研論文進(jìn)行撰寫(xiě),不僅給出課程實(shí)驗(yàn)的結(jié)果,而且分析產(chǎn)生這些結(jié)果的內(nèi)在機(jī)理及其表現(xiàn)形式,探討這些實(shí)驗(yàn)結(jié)果與文獻(xiàn)上已經(jīng)發(fā)表的科研結(jié)果的聯(lián)系與區(qū)別。同時(shí),訓(xùn)練學(xué)生的科研語(yǔ)言表達(dá)能力,讓他們通過(guò)圖書(shū)館期刊及CNKI 中國(guó)知網(wǎng)了解科技論文的寫(xiě)作格式和要求,激發(fā)他們初步了解和思考自己所學(xué)專業(yè)的相關(guān)科研前沿問(wèn)題,有效提高科研意識(shí),培養(yǎng)他們從課程內(nèi)容出發(fā)通過(guò)引申、變通等方式發(fā)現(xiàn)科學(xué)問(wèn)題、通過(guò)實(shí)驗(yàn)探索科學(xué)問(wèn)題、通過(guò)理論分析解決科學(xué)問(wèn)題的科研創(chuàng)新能力。在這個(gè)過(guò)程中,教師要及時(shí)收集學(xué)生對(duì)于課程實(shí)驗(yàn)的建議和實(shí)驗(yàn)反饋意見(jiàn),進(jìn)一步完善實(shí)驗(yàn)設(shè)計(jì),把課程實(shí)驗(yàn)的教學(xué)質(zhì)量提升到一個(gè)新的臺(tái)階。
科研實(shí)踐與課程教學(xué)相互融合的實(shí)驗(yàn)教學(xué)設(shè)計(jì)是信號(hào)與系統(tǒng)課程實(shí)驗(yàn)教學(xué)改革,實(shí)現(xiàn)“專創(chuàng)融合”和“思創(chuàng)融合”的重要舉措之一。本文提出基于科研實(shí)踐與課程教學(xué)相互融合理念結(jié)合Matlab 的信號(hào)與系統(tǒng)課程實(shí)驗(yàn)教學(xué)設(shè)計(jì),并給出一個(gè)典型實(shí)驗(yàn)說(shuō)明科研前沿內(nèi)容與課程基本理論內(nèi)容如何結(jié)合有效進(jìn)行Matlab 仿真的程序設(shè)計(jì)與用法。本文所述的分?jǐn)?shù)階電路系統(tǒng)的動(dòng)力學(xué)特征仿真實(shí)驗(yàn)設(shè)計(jì)是“科研實(shí)踐與課程教學(xué)相互融合”課程實(shí)驗(yàn)的典型案例之一。它使得同學(xué)們?cè)诖髮W(xué)本科階段就接觸到與自己專業(yè)相同或者相近的學(xué)科前沿的若干熱點(diǎn)課題,比如,通過(guò)編寫(xiě)計(jì)算機(jī)M 程序和對(duì)比該實(shí)驗(yàn)電路的整階及分?jǐn)?shù)階動(dòng)力學(xué)特征,獲得機(jī)會(huì)通過(guò)本科課程實(shí)驗(yàn)教學(xué)親身體驗(yàn)一個(gè)相對(duì)完整的科研過(guò)程。這既有助于對(duì)《信號(hào)與系統(tǒng)》課程內(nèi)容,特別是系統(tǒng)函數(shù)及系統(tǒng)的狀態(tài)變量分析的深入理解和掌握,也對(duì)學(xué)生在計(jì)算機(jī)軟件的配合下通過(guò)實(shí)驗(yàn)探索的體驗(yàn)促進(jìn)初步的科研訓(xùn)練有重要作用。