中南大學(xué) 信號(hào)與系統(tǒng)matlab實(shí)驗(yàn)報(bào)告

《中南大學(xué) 信號(hào)與系統(tǒng)matlab實(shí)驗(yàn)報(bào)告》由會(huì)員分享，可在線閱讀，更多相關(guān)《中南大學(xué) 信號(hào)與系統(tǒng)matlab實(shí)驗(yàn)報(bào)告（23頁(yè)珍藏版）》請(qǐng)?jiān)谘b配圖網(wǎng)上搜索。

1、傳播優(yōu)秀Word版文檔 ，希望對(duì)您有幫助，可雙擊去除！ 實(shí)驗(yàn)一 基本信號(hào)的生成 1．實(shí)驗(yàn)?zāi)康? l 學(xué)會(huì)使用MATLAB產(chǎn)生各種常見(jiàn)的連續(xù)時(shí)間信號(hào)與離散時(shí)間信號(hào)； l 通過(guò)MATLAB中的繪圖工具對(duì)產(chǎn)生的信號(hào)進(jìn)行觀察，加深對(duì)常用信號(hào)的理解； l 熟悉MATLAB的基本操作，以及一些基本函數(shù)的使用，為以后的實(shí)驗(yàn)奠定基礎(chǔ)。 2．實(shí)驗(yàn)內(nèi)容 ⑴ 運(yùn)行以上九個(gè)例子程序，掌握一些常用基本信號(hào)的特點(diǎn)及其MATLAB實(shí)現(xiàn)方法；改變有關(guān)參數(shù)，進(jìn)一步觀察信號(hào)波形的變化。 ⑵ 在 范圍內(nèi)產(chǎn)生并畫出以下信號(hào)： a) ； b) ； c) ； d) 。 源程序： k=-10:10; f1k

2、=[zeros(1,10),1,zeros(1,10)]; subplot(2,2,1) stem(k,f1k) title(f1[k]) f2k=[zeros(1,8),1,zeros(1,12)]; subplot(2,2,2) stem(k,f2k) title(f2[k]) f3k=[zeros(1,14),1,zeros(1,6)]; subplot(2,2,3) stem(k,f3k) title(f3[k]) f4k=2*f2k-f3k; subplot(2,2,4) 22 / 23 stem(k,f4k) title(f4[k]) ⑶ 在

3、范圍內(nèi)產(chǎn)生并畫出以下信號(hào)： a) ； b) ； c) 。 請(qǐng)問(wèn)這三個(gè)信號(hào)的基波周期分別是多少？ 源程序： k=0:31; f1k=sin(pi/4*k).*cos(pi/4*k); subplot(3,1,1) stem(k,f1k) title(f1[k]) f2k=(cos(pi/4*k)).^2; subplot(3,1,2) stem(k,f2k) title(f2[k]) f3k=sin(pi/4*k).*cos(pi/8*k); subplot(3,1,3) stem(k,f3k) title(f3[k]) 其中f1[k]的基波周期是4,

4、f2[k]的基波周期是4, f3[k]的基波周期是16。 實(shí)驗(yàn)二 信號(hào)的基本運(yùn)算 1．實(shí)驗(yàn)?zāi)康? l 學(xué)會(huì)使用MATLAB完成信號(hào)的一些基本運(yùn)算； l 了解復(fù)雜信號(hào)由基本信號(hào)通過(guò)尺度變換、翻轉(zhuǎn)、平移、相加、相乘、差分、求和、微分及積分等運(yùn)算來(lái)表達(dá)的方法； l 進(jìn)一步熟悉MATLAB的基本操作與編程，掌握其在信號(hào)分析中的運(yùn)用特點(diǎn)與使用方式。 2．實(shí)驗(yàn)內(nèi)容 ⑴ 運(yùn)行以上三個(gè)例題程序，掌握信號(hào)基本運(yùn)算的MATLAB實(shí)現(xiàn)方法；改變有關(guān)參數(shù)，考察相應(yīng)信號(hào)運(yùn)算結(jié)果的變化特點(diǎn)與規(guī)律。 ⑵ 已知信號(hào)如下圖所示： a) 用MATLAB編程復(fù)現(xiàn)上圖；

5、 %作業(yè)題2 a： t=-6:0.001:6; ft1=tripuls(t,6,0.5); subplot(2,1,1) plot(t,ft1) title(f(t)) b) 畫出的波形； %b t=-6:0.001:6; ft1=tripuls(2*(1-t),6,0.5); %subplot(1,1,1) plot(t,ft1) title(f(2*(1-t)) c) 畫出的波形； %c h=0.001;t=-6:h:6; yt=tripuls(t,6,0.5); y1=diff(yt)*1/h; plot(t(1:length

6、(t)-1),y1) title(df(t)/dt) d) 畫出的波形。 %d t=-6:0.1:6; for x=1:length(t) y2(x)=quad(tripuls(t,6,0.5),-3,t(x)); end plot(t,y2) title(integral of f(t)) 實(shí)驗(yàn)三 系統(tǒng)的時(shí)域分析 1．實(shí)驗(yàn)?zāi)康? l 學(xué)習(xí)并掌握連續(xù)時(shí)間系統(tǒng)的零狀態(tài)響應(yīng)、沖激響應(yīng)和階躍響應(yīng)的MATLAB求解方法； l 學(xué)習(xí)并掌握離散時(shí)間系統(tǒng)的零狀態(tài)響應(yīng)、沖激響應(yīng)和階躍響應(yīng)的MATLAB求解方法； l 進(jìn)一步深刻理解連續(xù)時(shí)間系統(tǒng)和離散時(shí)間系統(tǒng)的系統(tǒng)函

7、數(shù)零極點(diǎn)對(duì)系統(tǒng)特性的影響； 學(xué)習(xí)并掌握卷積的MATLAB計(jì)算方法。 2．實(shí)驗(yàn)內(nèi)容 ⑴ 運(yùn)行以上五個(gè)例題程序，掌握求解系統(tǒng)響應(yīng)的MATLAB分析方法；改變模型參數(shù)，考察系統(tǒng)響應(yīng)的變化特點(diǎn)與規(guī)律。 ⑵ 設(shè)離散系統(tǒng)可由下列差分方程表示： 計(jì)算時(shí)的系統(tǒng)沖激響應(yīng)。 源程序： k=-20:100; a=[1 -1 0.9]; b=[1]; h=impz(b,a,k); stem(k,h); xlabel(Time(sec)) ylabel(y(t)) ⑶ 設(shè)，輸入，求系統(tǒng)輸出。 （?。? 源程序： k=-10:50; uk=[zeros(1,10),ones(1,

8、51)]; u1k=[zeros(1,20),ones(1,41)]; hk=0.9.^k.*uk; fk=uk-u1k; yk=conv(hk,fk); stem(0:length(yk)-1,yk); ⑷ 已知濾波器的傳遞函數(shù)： 輸入信號(hào)為為隨機(jī)信號(hào)。試?yán)L出濾波器的輸出信號(hào)波形。（?。? 源程序： R=101; d=rand(1,R)-0.5; t=0:100; s=2*sin(0.05*pi*t); f=s+d; subplot(2,1,1); plot(t,d,g-.,t,s,b--,t,f,r-); xlabel(Time index t);

9、legend(d[t],s[t],f[t]); title(處理前的波形) b=[0.22 0];a=[1 -0.8]; y=filter(b,a,f); subplot(2,1,2); plot(t,s,b--,t,y,r-); xlabel(Time index t); legend(s[t],y[t]); title(濾波器輸出波形) 實(shí)驗(yàn)四 周期信號(hào)的頻域分析 1．實(shí)驗(yàn)?zāi)康? l 掌握周期信號(hào)傅立葉級(jí)數(shù)分解與合成的計(jì)算公式 l 掌握利用MATLAB實(shí)現(xiàn)周期信號(hào)傅立葉級(jí)數(shù)分解與綜合方法 l 理解并掌握周期信號(hào)頻譜特點(diǎn) 2.實(shí)驗(yàn)內(nèi)容 1、仿照例程，實(shí)現(xiàn)下述周

10、期信號(hào)的傅立葉級(jí)數(shù)分解與合成： 1 -3 -4 5 4 1 O 要求： （a）首先，推導(dǎo)出求解，，的公式，計(jì)算出前10次系數(shù)； （b）利用MATLAB求解，，的值，其中，求解前10次系數(shù)，并給出利用這些系數(shù)合成的信號(hào)波形。 （a）設(shè)周期信號(hào)的周期為，角頻率，且滿足狄里赫利條件，則該周期信號(hào)可以展開(kāi)成傅立葉級(jí)數(shù)。 （1）三角形式傅立葉級(jí)數(shù) （2）指數(shù)形式傅立葉級(jí)數(shù)

11、 （b）求解，，及合成信號(hào)波形所用程序： function [A_sym,B_sym]=CTFShchsym % 采用符號(hào)計(jì)算求一個(gè)周期內(nèi)連續(xù)時(shí)間函數(shù)f的三角級(jí)數(shù)展開(kāi)系數(shù),再用這些 % 展開(kāi)系數(shù)合成連續(xù)時(shí)間函數(shù)f.傅立葉級(jí)數(shù) % 函數(shù)的輸入輸出都是數(shù)值量 % Nf=6 諧波的階數(shù) % Nn 輸出數(shù)據(jù)的準(zhǔn)確位數(shù) % A_sym 第1元素是直流項(xiàng)，其后元素依次是1,2,3...次諧波cos項(xiàng)展開(kāi)系數(shù) % B_sym 第2,3,4,...元素依次是1,2,3...次諧波sin項(xiàng)展開(kāi)系數(shù) % ta

12、o=1 tao/T=0.2 syms t n k x T=4; tao=T/4; a=-1.5; if nargin<4 Nf=10; end if nargin<5 Nn=32; end x=time_fun_x(t); A0=int(x,t,a,T+a)/T; %求出三角函數(shù)展開(kāi)系數(shù)A0 As=2/T*int(x*cos(2*pi*n*t/T),t,a,T+a); %求出三角函數(shù)展開(kāi)系數(shù)As Bs=2/T*int(x*sin(2*pi*n*t/T),t,a,T+a); %求出三角

13、函數(shù)展開(kāi)系數(shù)Bs A_sym(1)=double(vpa(A0,Nn)); %獲取串?dāng)?shù)組A0所對(duì)應(yīng)的ASC2碼數(shù)值數(shù)組 for k=1:Nf A_sym(k+1)=double(vpa(subs(As,n,k),Nn)); %獲取串?dāng)?shù)組A所對(duì)應(yīng)的ASC2碼數(shù)值數(shù)組 B_sym(k+1)=double(vpa(subs(Bs,n,k),Nn)); %獲取串?dāng)?shù)組B所對(duì)應(yīng)的ASC2碼數(shù)值數(shù)組 end ; if nargout==0 c=A_sym; disp(c); %輸出c為三角級(jí)數(shù)展開(kāi)系數(shù):第1元素是直流項(xiàng)，其后元素依次是1,2,3

14、...次諧波cos項(xiàng)展開(kāi)系數(shù) d=B_sym; disp(d); %輸出d為三角級(jí)數(shù)展開(kāi)系數(shù): 第2,3,4,...元素依次是1,2,3...次諧波sin項(xiàng)展開(kāi)系數(shù) t=-3*T:0.01:3*T; f0=c(1); %直流 f1=c(2).*cos(2*pi*1*t/T)+d(2).*sin(2*pi*1*t/T); % 基波 f2=c(3).*cos(2*pi*2*t/T)+d(3).*sin(2*pi*2*t/T);

15、 % 2次諧波 f3=c(4).*cos(2*pi*3*t/T)+d(4).*sin(2*pi*3*t/T); % 3次諧波 f4=c(5).*cos(2*pi*4*t/T)+d(5).*sin(2*pi*4*t/T); % 4次諧波 f5=c(6).*cos(2*pi*5*t/T)+d(6).*sin(2*pi*5*t/T); % 5次諧波 f6=c(7).*cos(2*pi*6*t/T)+d(7).*sin(2*pi*6*t/T);

16、 % 6次諧波 f7=c(8).*cos(2*pi*7*t/T)+d(8).*sin(2*pi*7*t/T); % 7次諧波 f8=c(9).*cos(2*pi*8*t/T)+d(9).*sin(2*pi*8*t/T); % 8次諧波 f9=c(10).*cos(2*pi*9*t/T)+d(10).*sin(2*pi*9*t/T); % 9次諧波 f10=c(11).*cos(2*pi*10*t/T)+d(11).*sin(2*pi*10*t/T);

17、 % 10次諧波 f11=f0+f1+f2; % 直流+基波+2次諧波 f12=f11+f3; % 直流+基波+2次諧波+3次諧波 f13=f12+f4+f5+f6; % 直流+基波+2次諧波+3次諧波+4次諧波+5次諧波+6次諧波 f14=f13+f7+f8+f9+f10; %0~10次 subplot(2,2,1) plot(t,f0+f1),hold on y=time_fun_e(t); %調(diào)用連續(xù)時(shí)間函數(shù)-周期矩形脈沖 plot(t,y,r:)

18、 title(直流+基波) axis([-8,8,-0.5,1.5]) subplot(2,2,2) plot(t,f12),hold on y=time_fun_e(t); plot(t,y,r:) title(1-3次諧波+直流) axis([-8,8,-0.5,1.5]) subplot(2,2,3) plot(t,f13),hold on y=time_fun_e(t); plot(t,y,r:) title(1-6次諧波+直流) axis([-8,8,-0.5,1.5]) subplot(2,2

19、,4) plot(t,f14),hold on y=time_fun_e(t); plot(t,y,r:) title(1-10次諧波+直流) axis([-8,8,-0.5,1.5]) hold off end function y=time_fun_e(t) % 該函數(shù)是CTFShchsym.m的子函它由符號(hào)函數(shù)和表達(dá)式寫成 a=1.5; T=4; h=1; tao=T/4; t=-3*T:0.01:3*T; e1=1/2+1/2.*sign(t-0.5+tao/2); e2=1/2+1/2.*sign(t-0.5-tao

20、/2); y=h.*(e1-e2); %連續(xù)時(shí)間函數(shù)-周期矩形脈沖 function x=time_fun_x(t) % 該函數(shù)是CTFShchsym.m的子函數(shù)。它由符號(hào)變量和表達(dá)式寫成。 h=1; x1=sym(Heaviside(t))*h; x=x1-sym(Heaviside(t-1))*h; 源程序修改： function [A_sym,B_sym]=CTFShchsym % 采用符號(hào)計(jì)算求一個(gè)周期內(nèi)連續(xù)時(shí)間函數(shù)f的三角級(jí)數(shù)展開(kāi)系數(shù),再用這些 % 展開(kāi)系數(shù)合成連續(xù)時(shí)間函數(shù)f.傅立葉級(jí)數(shù) % 函數(shù)的輸入輸出都是數(shù)值量 %

21、Nf=6 諧波的階數(shù) % Nn 輸出數(shù)據(jù)的準(zhǔn)確位數(shù) % A_sym 第1元素是直流項(xiàng)，其后元素依次是1,2,3...次諧波cos項(xiàng)展開(kāi)系數(shù) % B_sym 第2,3,4,...元素依次是1,2,3...次諧波sin項(xiàng)展開(kāi)系數(shù) % tao=1 tao/T=0.2 syms t n k x T=5; tao=0.2*T; a=0.5; if nargin<4 Nf=6; end if nargin<5 Nn=32; end x=time_fun_x(t); A0=int(x,t,-a,T-a)/T;

22、 %求出三角函數(shù)展開(kāi)系數(shù) A0 As=2/T*int(x*cos(2*pi*n*t/T),t,-a,T-a); %求出三角函數(shù)展開(kāi)系數(shù)As Bs=2/T*int(x*sin(2*pi*n*t/T),t,-a,T-a); %求出三角函數(shù)展開(kāi)系數(shù)Bs A_sym(1)=double(vpa(A0,Nn)); %獲取串?dāng)?shù)組A0所對(duì)應(yīng)的ASC2碼數(shù)值數(shù)組 for k=1:Nf A_sym(k+1)=double(vpa(subs(As,n,k),Nn)); %獲取串?dāng)?shù)組A所對(duì)

23、應(yīng)的ASC2碼數(shù)值數(shù)組 B_sym(k+1)=double(vpa(subs(Bs,n,k),Nn)); %獲取串?dāng)?shù)組B所對(duì)應(yīng)的ASC2碼數(shù)值數(shù)組 end if nargout==0 c=A_sym; disp(c) %輸出c為三角級(jí)數(shù)展開(kāi)系數(shù):第1元素是直流項(xiàng)，其后元素依次是1,2,3...次諧波cos項(xiàng)展開(kāi)系數(shù) d=B_sym; disp(d) %輸出d為三角級(jí)數(shù)展開(kāi)系數(shù): 第2,3,4,...元素依次是1,2,3...次諧波sin項(xiàng)展開(kāi)系數(shù) t=-8*a:0.01:T-a; f0=c(1);

24、 %直流 f1=c(2).*cos(2*pi*1*t/5)+d(2).*sin(2*pi*1*t/5); % 基波 f2=c(3).*cos(2*pi*2*t/5)+d(3).*sin(2*pi*2*t/5); % 2次諧波 f3=c(4).*cos(2*pi*3*t/5)+d(4).*sin(2*pi*3*t/5); % 3次諧波 f4=c(5).*cos(2*pi*4*t/5)+d(

25、5).*sin(2*pi*4*t/5); % 4次諧波 f5=c(6).*cos(2*pi*5*t/5)+d(6).*sin(2*pi*5*t/5); % 5次諧波 f6=c(7).*cos(2*pi*6*t/5)+d(7).*sin(2*pi*6*t/5); % 6次諧波 f7=f0+f1+f2; % 直流+基波+2次諧波 f8=f7+f3; % 直流+基波+2次諧波+3次諧波 f9=f8+f4+f6; % 直流+基波+2次諧波+3次

26、諧波+4次諧波+6次諧波 subplot(2,2,1) plot(t,f0+f1),hold on y=time_fun_e(t); %調(diào)用連續(xù)時(shí)間函數(shù)-周期矩形脈沖 plot(t,y,r:) title(周期矩形波的形成—直流+基波) axis([-4,4.5,-0.5,1.5]) subplot(2,2,2) plot(t,f7),hold on y=time_fun_e(t); plot(t,y,r:) title(周期矩形波的形成—直流+基波+2次諧波) axis([-4,4.5,-0.5,

27、1.5]) subplot(2,2,3) plot(t,f8),hold on y=time_fun_e(t); plot(t,y,r:) title(直流+基波+2次諧波+3次諧波) axis([-4,4.5,-0.5,1.5]) subplot(2,2,4) plot(t,f9),hold on y=time_fun_e(t); plot(t,y,r:) title(基波+2次諧波+3次諧波+4次諧波+6次諧波) axis([-4,4.5,-0.5,1.5]) end function y=time_f

28、un_e(t) % 該函數(shù)是CTFShchsym.m的子函它由符號(hào)函數(shù)和表達(dá)式寫成 a=0.5; T=5; h=1; tao=0.2*T; t=-8*a:0.01:T-a; e1=1/2+1/2.*sign(t+tao/2); e2=1/2+1/2.*sign(t-tao/2); y=h.*(e1-e2); %連續(xù)時(shí)間函數(shù)-周期矩形脈沖 function x=time_fun_x(t) % 該函數(shù)是CTFShchsym.m的子函數(shù)。它由符號(hào)變量和表達(dá)式寫成。 h=1; x1=sym(Heaviside(t+0.5))*h; x=x1-sym(Hea

29、viside(t-0.5))*h; 2、已知周期為T=4的三角波，在第一周期（-2

30、 T T=m*tao， 信號(hào)周期 % Nf 諧波的階數(shù) % m (m=T/tao)周期與脈沖寬度之比,如m=4,8,16,100等 % tao 脈寬:tao=T/m syms t n y if nargin<3 Nf=input(please Input 所需展開(kāi)的最高諧波次數(shù):Nf=); end T=input(please Input 信號(hào)的周期T=); if nargin<5 Nn=32; end y=time_fun_s(t); A0=2/T*int(y,t,0,T/2); As=2/T*i

31、nt(y*cos(2*pi*n*t/T),t,0,T/2); Bs=2/T*int(y*sin(2*pi*n*t/T),t,0,T/2); A_sym(1)=double(vpa(A0,Nn)); for k=1:Nf A_sym(k+1)=double(vpa(subs(As,n,k),Nn)); B_sym(k+1)=double(vpa(subs(Bs,n,k),Nn)); end if nargout==0 An=fliplr(A_sym); %對(duì)A_sym陣左右對(duì)稱交換 An(1,k+1)=A_sym(1);

32、%A_sym的1*k陣擴(kuò)展為1*(k+1)陣 An=fliplr(An); %對(duì)擴(kuò)展后的S1陣左右對(duì)稱交換回原位置 Bn=fliplr(B_sym); %對(duì)B_sym陣左右對(duì)稱交換 Bn(1,k+1)=0; %B_sym的1*k陣擴(kuò)展為1*(k+1)陣 Bn=fliplr(Bn); %對(duì)擴(kuò)展后的S3陣左右對(duì)稱交換回原位置 FnR=An/2-i*Bn/2; % 用三角函數(shù)展開(kāi)系數(shù)A、B值合成付里葉指數(shù)系數(shù) FnL=fliplr(FnR); N=Nf*2*pi/T;

33、 k2=-N:2*pi/T:N; Fn=[FnL,FnR(2:end)]; %subplot(3,3,3) %x=time_fun_e(t); % 調(diào)用連續(xù)時(shí)間函數(shù)-周期矩形脈沖 subplot(2,1,1) stem(k2,abs(Fn)); %畫出周期矩形脈沖的頻譜（T=M*tao） title(連續(xù)時(shí)間函數(shù)周期三角波脈沖的雙邊幅度譜) axis([-80,80,0,0.12]) line([-80,80],[0,0],color,r) line([

34、0,0],[0,0.12],color,r) end function x=time_fun_e(t) % 該函數(shù)是CTFSshbpsym.m的子函數(shù)。它由符號(hào)變量和表達(dá)式寫成。 % t 是時(shí)間數(shù)組 % T 是周期 duty=tao/T=0.2 T=5; t=-2*T:0.01:2*T; tao=T/5; x=rectpuls(t,tao); %產(chǎn)生一個(gè)寬度tao=1的矩形脈沖 subplot(2,2,2) plot(t,x) hold on x=rectpuls(t-5,tao); %產(chǎn)生一個(gè)寬度tao=1的矩形脈,中心位置在t=5處

35、plot(t,x) hold on x=rectpuls(t+5,tao); %產(chǎn)生一個(gè)寬度tao=1的矩形脈,中心位置在t=-5處 plot(t,x) title(周期為T=5，脈寬tao=1的矩形脈沖) axis([-10,10,0,1.2]) function y=time_fun_s(t) syms t y=1-abs(t); x1=sym(Heaviside(t+2)); x=x1-sym(Heaviside(t-2)); y=y*x; ezplot(t,y,[-10,10]) grid 源程序修改： function [A_sy

36、m,B_sym]=CTFSshbpsym(T,Nf) % 采用符號(hào)計(jì)算求[0,T]內(nèi)時(shí)間函數(shù)的三角級(jí)數(shù)展開(kāi)系數(shù)。 % 函數(shù)的輸入輸出都是數(shù)值量 % Nn 輸出數(shù)據(jù)的準(zhǔn)確位數(shù) % A_sym 第1元素是直流項(xiàng)，其后元素依次是1,2,3...次諧波cos項(xiàng)展開(kāi)系數(shù) % B_sym 第2,3,4,...元素依次是1,2,3...次諧波sin項(xiàng)展開(kāi)系數(shù) % T T=m*tao， 信號(hào)周期 % Nf 諧波的階數(shù) % m (m=T/tao)周期與脈沖寬度之比,如m=4,8,16,100等 % tao 脈寬:tao=T/m syms t

37、 n y if nargin<3 Nf=input(please Input 所需展開(kāi)的最高諧波次數(shù):Nf=); end T=input(please Input 信號(hào)的周期T=); if nargin<5 Nn=32; end y=time_fun_s(t); A0=2/T*int(y,t,0,T); As=2/T*int(y*cos(2*pi*n*t/T),t,0,T); Bs=2/T*int(y*sin(2*pi*n*t/T),t,0,T); A_sym(1)=double(vpa(A0,Nn)); for k=1:Nf A_sy

38、m(k+1)=double(vpa(subs(As,n,k),Nn)); B_sym(k+1)=double(vpa(subs(Bs,n,k),Nn)); end if nargout==0 An=fliplr(A_sym); %對(duì)A_sym陣左右對(duì)稱交換 An(1,k+1)=A_sym(1); %A_sym的1*k陣擴(kuò)展為1*(k+1)陣 An=fliplr(An); %對(duì)擴(kuò)展后的S1陣左右對(duì)稱交換回原位置 Bn=fliplr(B_sym); %對(duì)B_sy

39、m陣左右對(duì)稱交換 Bn(1,k+1)=0; %B_sym的1*k陣擴(kuò)展為1*(k+1)陣 Bn=fliplr(Bn); %對(duì)擴(kuò)展后的S3陣左右對(duì)稱交換回原位置 FnR=An/2-i*Bn/2; % 用三角函數(shù)展開(kāi)系數(shù)A、B值合成付里葉指數(shù)系數(shù) FnL=fliplr(FnR); N=Nf*2*pi/T; k2=-N:2*pi/T:N; Fn=[FnL,FnR(2:end)]; subplot(3,3,3) x=time_fun_e(t); % 調(diào)用

40、連續(xù)時(shí)間函數(shù)-周期矩形脈沖 subplot(2,1,1) stem(k2,abs(Fn)); %畫出周期矩形脈沖的頻譜（T=M*tao） title(連續(xù)時(shí)間函數(shù)周期矩形脈沖的雙邊幅度譜) axis([-80,80,0,0.12]) line([-80,80],[0,0]) line([0,0],[0,0.12]) end function x=time_fun_e(t) % 該函數(shù)是CTFSshbpsym.m的子函數(shù)。它由符號(hào)變量和表達(dá)式寫成。 % t 是時(shí)間數(shù)組 % T 是周期 duty=ta

41、o/T=0.2 T=5; t=-2*T:0.01:2*T; tao=T/5; x=rectpuls(t,tao); %產(chǎn)生一個(gè)寬度tao=1的矩形脈沖 subplot(2,2,2) plot(t,x) hold on x=rectpuls(t-5,tao); %產(chǎn)生一個(gè)寬度tao=1的矩形脈,中心位置在t=5處 plot(t,x) hold on x=rectpuls(t+5,tao); %產(chǎn)生一個(gè)寬度tao=1的矩形脈,中心位置在t=-5處 plot(t,x) title(周期為T=5，脈寬tao=1的矩形脈沖) axis([-10,10,0,

42、1.2]) function y=time_fun_s(t) % 該函數(shù)是CTFSshbpsym.m的子函數(shù)。它由符號(hào)變量和表達(dá)式寫成。 syms a a1 T=input(please Input 信號(hào)的周期T=); M=input(周期與脈沖寬度之比M=); A=1; tao=T/M; a=tao/2; y1=sym(Heaviside(t+a1))*A; y=y1-sym(Heaviside(t-a1))*A; y=subs(y,a1,a); y=simple(y); 實(shí)驗(yàn)五 非周期信號(hào)的頻域分析 1．實(shí)驗(yàn)?zāi)康? l 理解非周期信號(hào)的頻域分析方法；

43、l 掌握典型信號(hào)的幅度譜和相位譜； l 理解信號(hào)的調(diào)制特性； l 掌握傅里葉變換的性質(zhì)：尺度變換、時(shí)移、頻移、卷積定理、對(duì)稱性、微分特性。 2．實(shí)驗(yàn)內(nèi)容 1.試計(jì)算寬度為2、幅度為1的三角波信號(hào)在0~fmHz范圍內(nèi)信號(hào)的能量。取fm=0.1~10Hz。 syms t; f=linspace(0.1,10,256); N=length(f); w=zeros(1,N); for k=1:N w(k)=quadl(@(t)2*sinc(t).*sinc(t).*sinc(t).*sinc(t),0,f(k)); end y=0; for k=1:N y

44、=y+w(k); end y plot(f,w); xlabel(Hz); ylabel(E); y =168.2699 2.試計(jì)算寬度和幅度均為1的方波信號(hào)在0~fmHz頻譜范圍內(nèi)所包含的信號(hào)能量。 syms t; f=linspace(0,5,256); N=length(f); w=zeros(1,N); for k=1:N w(k)=quadl(@(t)2*sinc(t).*sinc(t),0,f(k)); end plot(f,w); xlabel(Hz); ylabel(E); 溫馨提示：最好仔細(xì)閱讀后才下載使用，萬(wàn)分感謝！