【文章內(nèi)容簡(jiǎn)介】 ot(t,h3t),title(39。\sigma=,\omega=039。),xlabel(39。\sigma39。),ylabel(39。j\ome 20 ga39。) （ 2）由系統(tǒng)函數(shù)零、極點(diǎn)分布決定頻率特性 給出 零極點(diǎn)分布圖，討論他們是何種類型的濾波器，畫(huà)出其幅頻特性曲線。 方法：首先根據(jù) 零極點(diǎn)分布圖 寫(xiě)出零極點(diǎn)，因 H(s)的分母多項(xiàng)式之根構(gòu)成極點(diǎn)，所以可寫(xiě)出 H(s)的表達(dá)式，由此寫(xiě)出 )(?jH 的表達(dá)式，應(yīng)用 MATLAB 軟件，用 ilaplace 拉普拉斯逆變換 函數(shù)，把 )(?jH 代換為 )( ?jH ， 再利用 plot 函數(shù)，畫(huà)出 )( ?jH 的圖 形。最后根據(jù)圖形來(lái)判斷屬于 何種類型的濾波器 。 算法（舉圖 a 為例）： s=0::5。 %s 的范圍為 [0， 5]，步長(zhǎng)為 H=1./(s.*i+1)+1./(s.*i+2)。 %列出 H的表達(dá)式 plot(s,H) %用 plot函數(shù)畫(huà)出圖形 Warning: Imaginary parts of plex X and/or Y arguments ignored. xlabel(39。\omega39。) %標(biāo)注 x軸坐標(biāo) ylabel(39。|H(j\omega)|39。) %標(biāo)注 x軸坐標(biāo) title(39。a 圖幅頻特性曲線39。) %標(biāo)注圖名 主要技術(shù)問(wèn)題 問(wèn)題 一： 在分析 由系統(tǒng)函數(shù)零、極點(diǎn)分布決定時(shí)域特性 題目中，并不能很好的理解題目的含義，欠缺理論知識(shí)。 問(wèn)題二： MATLAB 創(chuàng)建子圖 subplot 函數(shù)的正確應(yīng)用。 問(wèn)題三：波形分析。在分析 由系統(tǒng)函數(shù)零、極點(diǎn)分布決定時(shí)域特性 題目中， 21 共畫(huà)有 19 張小圖（有重復(fù)），如何在這么多的圖中找出極點(diǎn)分布與原函數(shù)波形的對(duì)應(yīng)關(guān)系。 問(wèn)題四：在 由系統(tǒng)函數(shù)零、極點(diǎn)分布決定頻率特性 的題目中，如何寫(xiě)出 )(?jH的表達(dá)式。 難題的解決方法 鞏固信號(hào)與系統(tǒng)以及 MATLAB 的相關(guān)知識(shí)，把題目的含義理順，整體把握設(shè)計(jì)的要求，運(yùn)用控制變量法進(jìn)類比，找出波形的規(guī)律。 )(?jH 的表達(dá)式可由H(s)得出。 系統(tǒng)設(shè)計(jì)的亮點(diǎn)或創(chuàng)新點(diǎn) 亮點(diǎn)一： 拉普拉斯變換與 MATLAB 軟件的結(jié)合使用，是問(wèn)題得以簡(jiǎn)化。 拉普拉斯變換 是 為簡(jiǎn)化計(jì)算而建立的實(shí)變量函數(shù)和復(fù)變量函數(shù)間的一種函數(shù)變換。 而 MATLAB 軟件特點(diǎn)是簡(jiǎn)潔和智能化。兩個(gè)工具結(jié)合使用，避免了復(fù)雜的計(jì)算以及不夠精確的手工畫(huà)圖。 亮點(diǎn)二： 創(chuàng)建子 圖， 分析波形時(shí) 進(jìn)行橫向與縱向的類比。 在設(shè)計(jì) 由系統(tǒng)函數(shù)零、極點(diǎn)分布決定時(shí)域特性 題目中，運(yùn)用了創(chuàng)建子圖的方法，并且控制變量的保持與改變，使所畫(huà)出來(lái)的圖對(duì)比更加鮮明，為分析波形的規(guī)律做好鋪墊。 第四章 、總結(jié)與體會(huì) 這次兩個(gè)星期的 MATLAB 軟件的學(xué)習(xí)及應(yīng)用 ，轉(zhuǎn)眼就過(guò)去了，在思考與探索中，我發(fā)覺(jué)每一天都是那么的短暫。在設(shè)計(jì)過(guò)程中，我們遇到了些許難題，克服難題的過(guò)程是一個(gè)痛苦與歡樂(lè)并存的過(guò)程。 首先，通過(guò)設(shè)計(jì)研究，我對(duì)這些理論知識(shí)有了全新的認(rèn)識(shí)，因?yàn)橥ㄟ^(guò)自學(xué)，可以找出自己想要理解的內(nèi)容。我也發(fā)現(xiàn)，在設(shè)計(jì)之前，我 們必須把理論知識(shí)打 22 扎實(shí)，否則很難對(duì)題目進(jìn)行整體的把握。這次的課程設(shè)計(jì)也提醒著我要在理論知識(shí)上花點(diǎn)功夫，為以后 畢業(yè)論文撰寫(xiě) 打下基礎(chǔ)，我也希望自己以后把這些知識(shí)運(yùn)用到實(shí)踐中去。 在設(shè)計(jì)之前，我閱讀了一些書(shū)籍來(lái)獲得一定的感性知識(shí)，然后才有自己的想法和觀點(diǎn)， 這次的題目并不難， 最大的困難是步驟多， 考驗(yàn)的是耐心。 通過(guò)設(shè)計(jì)研究， 我 還 了解了 運(yùn)用 MATLAB 設(shè)計(jì)的基本步驟，也了解了MATLAB 方便，簡(jiǎn)潔，高效的優(yōu)點(diǎn) ，掌握此軟件定能為以后的學(xué)習(xí)打下了基礎(chǔ) 。 在研究 由系統(tǒng)函數(shù)零、極點(diǎn)分布決定時(shí)域特性 的題目時(shí)，由于自己的粗心，導(dǎo) 致坐標(biāo)標(biāo)錯(cuò)，辛苦畫(huà)出的 5張圖全部得重新畫(huà)，浪費(fèi)大量的時(shí)間和精力。粗心的結(jié)果很殘酷，在以后的學(xué)習(xí)生活中，我會(huì)記住這次的教訓(xùn)。 在設(shè)計(jì)中，我也遇到了些許困難，例如，剛開(kāi)始設(shè)計(jì)題目的意思理解錯(cuò)誤，走入了一個(gè)很復(fù)雜的歪路，導(dǎo)致的結(jié)果是不知道如何分析波形，停滯在那不知如何，經(jīng)過(guò)同學(xué)的指點(diǎn)，才正確理解了題目的意思。這個(gè)教訓(xùn)也在告訴這我，不能鉆牛角尖，抓緊某一點(diǎn)不放手，應(yīng)該開(kāi)拓思路，可以與同學(xué)討論，也可以查閱相關(guān)書(shū)籍，從而得到靈感。 此外，本次課程設(shè)計(jì)報(bào)告是按照畢業(yè)設(shè)計(jì)報(bào)告的要求來(lái)的，因此我抓住此次機(jī)會(huì)，對(duì) word 的使用 作了進(jìn)一步的熟悉 ， 學(xué)會(huì)了 公式 產(chǎn)生，目錄的制作以及文檔的排版等，如果在畢業(yè)設(shè)計(jì)時(shí)再去熟悉這些小細(xì)節(jié)，勢(shì)必會(huì)花費(fèi)較多的時(shí)間，所以在平時(shí)應(yīng)該注意這些細(xì)節(jié)。 總之，這次實(shí)訓(xùn)受益頗多，在以后的學(xué)習(xí)中，我會(huì)好好利用這次的學(xué)習(xí)經(jīng)驗(yàn)，結(jié)合實(shí)際在專業(yè)領(lǐng)域內(nèi)進(jìn)行更深入的學(xué)習(xí)。 參考文獻(xiàn) [1] 鄭君里 (楊為理 ). 信號(hào)與系統(tǒng) [M]（第二版）上冊(cè) 北京： 高等教育出版社 ， 23 2021． [2]張志涌（楊祖櫻） . MATLAB 教程 [M] 北京： 北京航空航天大學(xué)出版社 ， 2021． [3]徐金明 （ 張孟喜、丁濤 ） . MATLAB 實(shí)用教程 [M] 北京： 清華大學(xué)出版社 ， 2021． [4]黨宏社 . 信號(hào)與系統(tǒng)實(shí)驗(yàn)（ MATLAB 版） [M] 西安 ： 西安電子科技大學(xué)出版社 ， 2021． [5] 劉會(huì)燈 （ 朱飛 ） . 信號(hào)與系統(tǒng) :MATLAB 綜合實(shí)驗(yàn) [M] 北京： 高等教育出版社 ，2021． 24 附錄 分析 由系統(tǒng)函數(shù)零、極點(diǎn)分布決定時(shí)域特性 的代碼如下： ? =， 0:: ??? syms t s。 h1=1/(s+*i)。 h2=1/(s+*i)。 h3=1/(s+)。 h1t=ilaplace(h1) h1t = exp(1/25*t)*(cos(3/5*t)+i*sin(3/5*t)) h2t=ilaplace(h2) h2t = exp(1/25*t)*(cos(3/10*t)+i*sin(3/10*t)) h3t=ilaplace(h3) h3t = exp(1/25*t) t=50::50。 h1t =exp(1/25.*t).*(cos(3/5.*t)+i.*sin(3/5.*t))。 h2t =exp(1/25.*t).*(cos(3/10.*t)+i.*sin(3/10.*t))。 h3t =exp(1/25.*t)。 subplot(3,1,1),plot(t,h1t),title(39。\sigma=,\omega=39。),xlabel(39。\sigma39。),ylabel(39。j\omega39。) Warning: Imaginary parts of plex X and/or Y arguments ignored. subplot(3,1,2),plot(t,h2t),title(39。\sigma=,\omega=39。),xlabel(39。\sigma39。),ylabel(39。j\omega39。) Warning: Imaginary parts of plex X and/or Y arguments ignored. 25 subplot(3,1,3),plot(t,h3t),title(39。\sigma=,\omega=039。),xlabel(39。\sigma39。),ylabel(39。j\omega39。) ? =0， 0:: ??? syms t s。 h1=1/(*i)。 h2=1/(*i)。 h3=1/s。 h1t=ilaplace(h1) h1t = cos(3/5*t)+i*sin(3/5*t) h2t=ilaplace(h2) h2t = cos(3/10*t)+i*sin(3/10*t) h3t=ilaplace(h3) h3t = 1 t=50::50。 h1t =cos(3/5.*t)+i.*sin(3/5.*t)。 h2t =cos(3/10.*t)+i.*sin(3/10.*t)。 h3t =1。 subplot(3,1,1),plot(t,h1t),title(39。\sigma=0,\omega=39。),xlabel(39。\sigma39。),ylabel(39。j\omega39。) Warning: Imaginary parts of plex X and/or Y arguments ignored. subplot(3,1,2),plot(t,h2t),title(39。\sigma=0,\omega=39。),xlabel(39。\sigma39。),ylabel(39。j\omega39。) Warning: Imaginary parts of plex X and/or Y arguments ignored. 26 subplot(3,1,3),plot(t,h3t),title(39。\sigma=,\omega=039。),xlabel(39。\sigma39。),ylabel(39。j\omega39。) ? =， 0:: ??? syms t s。 h1=1/(*i)。 h2=1/(*i)。 h3=1/()。 h1t=ilaplace(h1) h1t = exp(1/25*t)*(cos(3/5*t)+i*sin(3/5*t)) h2t=ilaplace(h2) h2t = exp(1/25*t)*(cos(3/10*t)+i*sin(3/10*t)) h3t=ilaplace(h3) h3t = exp(1/25*t) t=50::50。 h1t =exp(1/25.*t).*(cos(3/5.*t)+i.*sin(3/5.*t))。 h2t =exp(1/25.*t).*(cos(3/10.*t)+i.*sin(3/10.*t))。 h3t =exp(1/25*t)。 subplot(3,1,1),plot(t,h1t),title(39。\sigma=,\omega=39。),xlabel(39。\sigma39。),ylabel(39。j\omega39。) Warning: Imaginary parts of plex X and/or Y arguments ignored. subplot(3,1,2),plot(t,h2t),title(39。\sigma=,\omega=39。),xlabel(39。\sigma39。),ylabel(39。j\omega39。) Warning: Imaginary parts of plex X and/or Y arguments ignored. 27 subplot(3,1,3),plot(t,h3t),title(39。\sigma=,\omega=039。),xlabel(39。\sigma39。),ylabel(39。j\omega39。) ? =::,? = syms t s。 h1=1/(s+*i)。 h2=1/(s+*i)。 h3=1/(*i)。 h4=1/(*i)。 h5=1/(*i)。 h1t=ilaplace(h1) h1t = exp(1/25*t)*(cos(3/5*t)+i*sin(3/5*t)) h2t=ilaplace(h2) h2t = exp(1/50*t)*(cos(3/5*t)+i*sin(3/5*t)) h3t=ilaplace(h3) h3t = cos(3/5*t)+i*sin(3/5*t) h4t=ilaplace(h4) h4t = exp(1/50*t)*(cos(3/5*t)+i*sin(3/5*t)) h5t=ilaplace(h5) h5t = exp(1/25*t)*(cos(3/5*t)+i*sin(3/5*t)) t=50::50。 h1t =exp(1/25.*t).*(cos(3/5.*t)+i.*sin(3/5.*t))。 h2t =exp(1/50.*t).*(cos(3/5.*t)+i*sin(3/5.*t))。 h3t =cos(3/5.*t)+i.*sin(3/5*t)。 28 h4t =exp(1/50.*t).*(cos(3/5.*t)+i.*sin(3/5