時域信號的matlab仿真系統(tǒng)(參考版)
2025-05-30 22:51本頁面
【正文】 參考文獻(xiàn)[1]管致中、夏恭恪、孟橋：信號與線性系統(tǒng)，高等教育出版社，,P63P81[2]鄭君里、應(yīng)啟衍、楊為理：信號與系統(tǒng)，高等教育出版社，P93P99[3]丁玉美、高西全：數(shù)字信號處理，西安電子科技大學(xué)出版社，P45P66[4]程佩青：數(shù)字信號處理教程，清華大學(xué)出版社，P92P113[5]王沐然：MATLAB與科學(xué)計算，電子工業(yè)出版社，P281P316[6]劉保柱、蘇彥華、張宏林：，人民郵電出版社，P267P288[7]張玉華、馬旭波：MATLAB在信號分析與處理中的應(yīng)用，山西大學(xué)工程學(xué)院，P3P4成績評定表學(xué)生姓名x學(xué)生學(xué)號x學(xué)生專業(yè)通信工程廣播電視方向?qū)W生班級廣播電視班指導(dǎo)教師x職 稱講師設(shè)計題目時域信號的MATLAB仿真系統(tǒng)設(shè)計評語評定成績等級分?jǐn)?shù)指導(dǎo)教師簽字： 年 月 日說明：分?jǐn)?shù)采用百分制，相應(yīng)的等級為：優(yōu)：90～100：良：80～89；中：70～79；及格：60～69；不及格：60分以下。經(jīng)過老師的點(diǎn)評，對論文的結(jié)構(gòu)和格式做了修改，更好的論述了本文的題目，增加了一些內(nèi)容，充實(shí)了文本。在論文完成過程中，在圖書館查閱了很多資料，并且在網(wǎng)上搜集了很多資料，最終完成了這篇論文。由于現(xiàn)代計算機(jī)技術(shù)的發(fā)展，數(shù)字信號處理越來越占有重要地位，所以在有限的篇幅內(nèi)，簡要介紹了信號的抽樣，它是鏈接模擬信號處理和數(shù)字信號處理的橋梁。對于連續(xù)系統(tǒng)，簡要介紹了信號描述及運(yùn)算，連續(xù)時域系統(tǒng)分析和連續(xù)時域系統(tǒng)的仿真。set(gcf,’color’,’w’)。title(’impz函數(shù)求出的單位沖激響應(yīng)’)。ylabel(’h(n)’)。subplot(2,1,1),stem(n1,h1,’.’)。h1=filter(b,a,x)。n1=0:89。b=[1 ]。close all。 %直接型沖激響應(yīng)圖39 沖激響應(yīng)波形圖（2）系統(tǒng)由差分方程描述 (28)分別用filter函數(shù)和impz函數(shù)繪出其單位沖激響應(yīng)h(n)的線圖。stem(n,yout2,’.’)。圖38 離散系統(tǒng)的仿真圖MATLAB程序如下：n=0:20。title(’抽樣信號的頻譜’)。grid on。subplot(2,2,4)。hold off。ylabel(’fs(t)’)。axis([4 4 ])。stem(t2,fst)。plot(t1,ft,’:’)。fst=((1+cos(t2))/2).*(uCT(t2+pi)uCT(t2pi))。title(’升余弦脈沖信號的頻譜’)。grid on。subplot(2,2,2)。W=pi*k/(N*dt)。N=500。ylabel(’f(t)’)。axis([4 4 ])。plot(t1,ft)。ft=((1+cos(t1))/2).*(uCT(t1+pi)uCT(t1pi))。dt=。當(dāng)采樣間隔時，其實(shí)現(xiàn)的MATLAB程序代碼如下：clear all。升余弦脈沖信號為 (27)下面用MATLAB編程實(shí)現(xiàn)該信號經(jīng)沖激脈沖抽樣后得到的抽樣信號及其頻譜。從數(shù)學(xué)上講，抽樣過程就是抽樣脈沖p(t)和原連續(xù)信號f(t)相乘的過程，即： (26)因此，可以用傅里葉變換的頻域卷積性質(zhì)來求抽樣信號的頻譜。其仿真波形如下圖：圖36 單位取樣響應(yīng) 信號抽樣信號抽樣是連續(xù)時間信號分析向離散時間信號分析、連續(xù)信號處理向數(shù)字信號處理的第一步廣泛應(yīng)用于實(shí)際的各類系統(tǒng)中。grid on。n=0:30。a=[3 4 2]。有一離散線性系統(tǒng)的差分方程如下： (25)下面利用MATLAB的impz函數(shù)繪出該系統(tǒng)的單位取樣響應(yīng)。其仿真波形如下圖：圖35 零狀態(tài)響應(yīng) 單位取樣響應(yīng)系統(tǒng)的單位取樣響應(yīng)定義為系統(tǒng)在激勵下系統(tǒng)的零狀態(tài)響應(yīng)，用h(n)表示。xlabel(’n’)。stem(n,y,’fill’)。x=(1/2).^n。b=[1 2]。其實(shí)現(xiàn)的MATLAB程序代碼如下：clear all。函數(shù)filter的調(diào)用格式如下： (23)其中x為輸入的離散序列；y為輸出的離散序列；y的長度和x的長度一樣；b和a分別為差分方程右端與左端的系數(shù)向量。兩個序列的加減乘除是對應(yīng)離散樣點(diǎn)值的加減乘除，因此，可通過MATLAB點(diǎn)乘和點(diǎn)除、序列移位和反折來實(shí)現(xiàn)，與連續(xù)時間信號處理方法基本一樣。axis([0 2 100 200])。plot(w/pi,HP*180/pi)。axis([0 2 ])。xlabel(‘\omega 單位：pi’)。HM=abs(H)HP=angle(H)subplot(211)。w=0::2*pi。%應(yīng)用freqs函數(shù)求頻率響應(yīng)并作圖b=[1,0]。 %定義頻率采樣間隔subplot(224)pz_plxy(f1,f2,k,p,z)。f2=80。 %定義系統(tǒng)（d）零點(diǎn)行向量p=[100]。k= %定義頻率采樣間隔subplot(223)pz_plxy(f1,f2,k,p,z)。f1=0。 %繪制系統(tǒng)幅頻響應(yīng)曲線%繪制系統(tǒng)（c）的幅頻響應(yīng)曲線z=[20+i*40 20i*40]。k=。f1=0。 %繪制系統(tǒng)幅頻響應(yīng)曲線%繪制系統(tǒng)（b）的幅頻響應(yīng)曲線Z=[i*150 –i*150]。k=。f1=0。%xh07%繪制系統(tǒng)（a）的幅頻響應(yīng)曲線z=[0]。subplot(212)。[p,z]=pzmap(H)subplot(2111)。a=[1 2 5]。其波形如圖所示:圖25 零輸入響應(yīng)，零狀態(tài)響應(yīng)及完全響應(yīng)3 時域信號分析及MATLAB仿真 連續(xù)線性系統(tǒng)的時域分析及仿真（1）已知連續(xù)系統(tǒng)的系統(tǒng)函數(shù),繪制系統(tǒng)函數(shù)的零極點(diǎn)分布圖和計算系統(tǒng)函數(shù)的零極點(diǎn)位置和沖激響應(yīng)時域波形。grid on。subplot(3,1,3)。grid on。subplot(3,1,2)。grid on。yzs=simplify()yzs=1/4*exp(1/2*t)*Heaviside(t)yt=simplify(yzi+yzs)
公司管理相關(guān)推薦
鄂ICP備17016276號-1