信號與線性系統課程設計

上傳人：奔*** IP屬地：河北更新時間：2024-03-01 格式：doc 頁數：48 大?。?14.50KB 人氣指數：12 舉報 版權申訴

已閱讀1頁，還剩47頁未讀，繼續(xù)免費閱讀

版權說明：本文檔由用戶提供并上傳，收益歸屬內容提供方，若內容存在侵權，請進行舉報或認領

文檔簡介

1、　　信號與線性系統　　課程設計　　班級：　　學號：　　姓名：&l

2、t;/p>　　目錄　　實驗一連續(xù)信號的時域分析 …………………………… 3　　實驗二連續(xù)時間系統的時域分析 ……………………… 8　　實驗三連續(xù)信號的頻域分析 …………………………… 14　　實驗四

3、連續(xù)系統的頻域分析 …………………………… 22　　實驗五信號采樣與重建 ………………………………… 33　　實驗六離散時間信號和系統分析 ……………………… 39　　實驗總結 ………………………………………………… 46　　實驗一連續(xù)信號的時域分析</p

4、>　　一、實驗目的　　1、熟悉MATLAB 軟件。　　2、掌握常用連續(xù)信號與離散信號的MATLAB 表示方法。　　二、實驗設備　　安裝有matlab6.5 以上版本的P

5、C 機一臺。　　三、實驗原理　　四、實驗內容　　1、用MATLAB表示連續(xù)信號：，Acos(ω0 t +? )，Asin(ω0 t +? )。　　源程序：&l

6、t;/p>　　clc　　clear　　close　　syms t;　　f1=2*exp(t)

7、;　　f2=2*cos(3*t+4);　　f3=2*sin(3*t+4);　　subplot(2,2,1);　　ezplot(f1,[-10,2]);　　xlabel('t');<p&

8、gt;　　title('f(t)=2e^t');　　grid on;　　subplot(2,2,2);　　ezplot(f2,[-5,5]);　　xlabel('t');<

9、p>　　title('f(t)=2cos(3t+4)');　　grid on;　　subplot(2,2,3);　　ezplot(f3,[-5,5]);　　xlabel('t');</p&g

10、t;　　title('f(t)=2sin(3t+4)');　　grid on　　2、用MATLAB表示抽樣信號(sinc(t))、矩形脈沖信號(rectpuls(t, width))　　及三角脈沖信號(tripuls(t, width,

11、skew))。　　源程序：　　clc　　clear　　close　　t=-5:

12、0.01:5;　　f1=sinc(t);　　f2=3*rectpuls(t,4);　　f3=3*tripuls(t,4,0);　　subplot(2,2,1);　　plot(t,f1);

13、　xlabel('t');　　title('f(t)=sinc(t)');　　grid on;　　subplot(2,2,2)　　plot(t,f2);　　xl

14、abel('t');　　title('f(t)=3rectpuls(t,4)');　　grid on;　　axis([-5,5,-1,4]);　　subplot(2,2,3);&

15、lt;p>　　plot(t,f3);　　xlabel('t');　　title('f(t)=3rectpuls(t,4,0)');　　grid on;　　axis([-5,5,-1,4]);<

16、;/p>　　3、編寫如圖3 的函數并用MATLAB 繪出滿足下面要求的圖形。　　(1) f (?t); (2) f (t ? 2); (3) f (1? 2t); (4) f (0.5t +1).　　源程序：　　clc<

17、;/b>　　clear　　close　　t1=-14:0.05:2;　　t2=0:0.05:16;　　t=-6.5:0.05:1.5;</p&

18、gt;　　t4=-6:0.05:26;　　f1=4*rectpuls(-t1-6,12)+3*tripuls(-t1-6,4,0);　　f2=4*rectpuls(t2-8,12)+3*tripuls(t2-8,4,0);　　f3=4*rectpuls(1-2*t-6,12)+3*tripu

19、ls(1-2*t-6,4,0);　　f4=4*rectpuls(0.5*t4+1-6,12)+3*tripuls(0.5*t4+1-6,4,0);　　subplot(2,2,1);　　plot(t1,f1);　　xlabel('t');

20、　　title('f(-t)');　　subplot(2,2,2);　　plot(t2,f2);　　xlabel('t');　　title('f(t-2)');

21、　　subplot(2,2,3);　　plot(t,f3);　　xlabel('t');　　title('f(1-2t)');　　subplot(2,2,4);　　plot(t4,f4);</

22、p>　　xlabel('t');　　title('f(0.5t+1)');　　實驗二連續(xù)時間系統的時域分析　　一、實驗目的　　1、掌握卷積計算方法。

23、　　2、掌握函數lsim，impulse，step 的用法，lsim 為求取零狀態(tài)響應，　　impulse 為求取單位脈沖響應，step 為求取單位階躍響應。　　3、運用課堂上學到的理論知識，從RC、RL 一階電路的響應中正確　　區(qū)分零輸入響應、零狀態(tài)響應、自由響應與受迫響應。

24、;　　二、實驗設備　　安裝有matlab6.5 以上版本的PC 機一臺。　　三、實驗原理　　四、實驗內容　　1. 分別用函數lsim 和

25、卷積積分兩種方法求如圖7 所示系統的零狀態(tài)　　響應。其中L=1，R=2，e(t) = ε(t)，i(0? ) = 2。　　源程序：　　方法一：　　clear&

26、lt;/p>　　close　　clc　　t=0:0.01:10;　　f=exp(-t);　　a=[1 2];

27、　　b=[1];　　y=lsim(b,a,f,t);　　plot(t,y);　　xlabel('Time(sec)');　　ylabel('i(t)');&l

28、t;p>　　axis([0,10,-0.025,0.275]);　　grid on;　　方法二：　　易求得系統的沖激響應為ε(t)　　clear;<

29、;/p>　　clc;　　close;　　syms t x;　　e=exp(-x);　　h=exp(-2.*(t-x));</

30、p>　　i=int(e.*h,x,0,t);　　ezplot(i,[0,10]);　　xlabel('Time(sec)');　　ylabel('i(t)');　　title('exp(-t)*exp(-

31、2t)');　　grid on;　　2. 求上述系統的沖激響應與階躍響應。　　沖激響應源程序：　　clear;

32、;　　close;　　clc;　　a=[1 2];　　b=[1];　　impulse(b,a,10);

33、;　　xlabel('Time(sec)');　　ylabel('i(t)');　　axis([0,10,-0.1,1]);　　grid on;　　階躍響應源程序

34、：　　clear;　　close;　　clc;　　a=[1 2];<p

35、>　　b=[1];　　step(b,a,10);　　xlabel('Time(sec)');　　ylabel('i(t)');　　axis([0,10,0,0.55]);&

36、lt;p>　　grid on;　　五、思考題　　1. 為什么連續(xù)時間系統的零狀態(tài)響應為激勵與沖擊響應的卷積？　　答：根據卷積的定義，函數e(t)與函數h(t）相卷積后，就是在變量由負無窮到正無窮范圍內，對于某一t值時乘積e(τ)h(

37、t-τ)曲線下的面積，也就是：r(t)=e(t)*h(t），又零狀態(tài)響應與系統的特性和外加激勵有關，所以如問題。　　2. 利用卷積積分法計算系統響應應從幾個方面進行?　　答：利用卷積積分法先要將系統的沖擊響應求出，之后再將其與激勵卷積即可　　實驗三連續(xù)信號的頻域分析&l

38、t;b>　　一、實驗目的　　1. 掌握周期信號的頻譜—— Fourier 級數的分析方法。　　2. 深入理解信號頻譜的概念，掌握典型信號的頻譜以及Fourier 變換　　的主要性質。　　3. 掌握調制與解調的

39、基本原理及濾波器的使用。　　二、實驗設備　　安裝有matlab6.5 以上版本的PC 機一臺。　　三、實驗原理　　四、實驗內容<

40、p>　　1. 求如圖所示周期矩形脈沖信號的Fourier級數表達式，畫出頻譜圖，并用前N次諧波合成的信號近似。　　源程序：　　close　　clear&

41、lt;b>　　clc　　syms t n;　　T=2*pi; %設T為2*pi　　N=9; %設N為9　　f=heaviside(t)-2*heaviside(t-T/2)+heaviside(t-T); s

42、ubplot(2,2,1);　　ezplot(f,[0,2.1*pi]);　　title('原函數');　　h=exp(-j*n*2*pi/T*t); %用指數傅里葉級數表示，角頻率2*pi/T　　A1=int(f.*h,t,0,T);

43、　　A=2/T*A1;　　for n1=-N:-1　　C(n1+10)=subs(A,n,n1);　　end　　for n1=1:N

44、　　C(n1+10)=subs(A,n,n1);　　end　　C(10)=0;%N=0時系數　　subplot(2,2,3);　　k=-N:N;　　ste

45、m(k,abs(C)); %作出幅度譜　　ylabel('Cn 的幅度');　　xlabel('\Omega');　　subplot(2,2,4);　　stem(k,angle(C)); %作出相位譜

46、;　　ylabel('Cn的相位');　　xlabel('\Omega');　　f1=0;　　for m=-N:N　　f1=f1+1/2*C(m+10)*exp(j*m*t); %前N次諧波合成的信號近

47、似　　end　　subplot(2,2,2);　　ezplot(f1,[0,2.01*pi]);　　title('前N次諧波合成的信號近似');　　2、試用fourier()函

48、數求下列信號的傅里葉變換F( jω) ,并畫出F( jω)　　(1) f (t) = te?3tε (t)　　(2) f (t) = sgn(t)　　(1)源程序：　　clear

49、;　　close　　clc　　syms t x ;　　x=fourier(t*exp(-3*t)*heaviside(t));　　x&l

50、t;/p>　　z=abs(x);　　ezplot(z);%符號函數作圖函數　　xlabel('Time(sec)');　　ylabel('|F(jω)|');

51、;　　grid on　　得F(jω)= 1/(3+ j*ω)^2　　(2) 源程序：　　clear　　close&

52、lt;p>　　clc　　syms t x ;　　x=fourier(2*heaviside(t)-1);%2*heaviside(t)-1即為sgn(t)　　x　　z=ab

53、s(x);　　ezplot(z);%符號函數作圖函數　　xlabel('Time(sec)');　　ylabel('|F(jω)|');　　grid on&l

54、t;p>　　得F(jω)= -2*j/ω　　3、調制信號為一取樣信號，利用MATLAB 分析幅度調制(AM)產生　　的信號頻譜，比較信號調制前后的頻譜并解調已調信號。設載波信號　　的頻率為100Hz。　　源程序：</p&g

55、t;　　clc;　　clear;　　close;　　Fm=10;　　t1=0:0.00002:0.

56、2;　　syms t v;　　x=sin(2.0*pi*Fm*t)/(2.0*pi*Fm*t);　　subplot(3,2,1);　　ezplot(x,[0,0.2]);　　title(

57、9;原函數');　　Fx=fourier(x,v);　　subplot(3,2,2);　　ezplot(Fx,[-50*pi,50*pi]);　　axis([-50*pi,50*pi,-0.05,0.1]);　　title

58、('頻譜');　　y=x*cos(200*pi*t);　　subplot(3,2,3);　　b=subs(y,t,t1);　　plot(t1,b);　　title('調制后');

59、　　axis([0,0.2,-1,1]);　　Fy=fourier(y,v);　　subplot(3,2,4);　　ezplot(Fy,[-250*pi,250*pi]);　　axis([-250*pi,250*pi,-0.05,0.1]);<

60、/p>　　title('頻譜');　　z=y*cos(200*pi*t);　　Fz=fourier(z,v);　　G=-heaviside(v-20*pi)+heaviside(v+20*pi);%門函數　　F

61、x1=Fz*G;　　x1=2*ifourier(Fx1,v);%濾波過程中幅度減半且反向。　　subplot(3,2,5);　　ezplot(x1,[0,0.2]);　　title('解調后');

62、;　　subplot(3,2,6);　　ezplot(2*Fx1,[-50*pi,50*pi]);　　axis([-50*pi,50*pi,-0.05,0.1]);　　title('頻譜');　　五、思考題</p

63、>　　1、根據試驗1 的結果，解釋Gibbs 現象。　　答：因為對于具有不連續(xù)點的函數，即使級數的項無限增大，在不連續(xù)處，級數之和不收斂于函數f(t)；在躍變點附近的波形，總是不可避免的存在有起伏震蕩，從而使躍變點的值超過一形成過沖，造成吉布斯現象。　　2、比較周期信號與非周期信號的頻譜。<p&g

64、t;　　區(qū)別：　　1．周期信號的頻譜為離散頻譜，非周期信號的頻譜為連續(xù)頻譜。　　2.周期信號的頻譜為Fn的分布，表示每個諧波分量的復振幅；而周期信號的頻譜為F(jω)的分布，(F(jω)/2)ω表示合成諧波分量的復振幅，所以也將稱為頻譜密度函數。　　聯系：&

65、lt;/b>　　1.都是反映將時域信號表示為正弦類信號時各諧波分量的分布特性。　　2.若周期信號是連續(xù)非周期信號的周期延拓，則兩者的關系為　　F(jω)= ； =　　3、調制與解調的基本原理是什么？為什么要進行調制？&l

66、t;p>　　調制：調制就是使信號f(t)控制載波的某一個或某些參數（如振幅、頻率、相位等），是這些參數按照信號f(t)的規(guī)律變化的過程。載波可以是正弦波或脈沖序列。以正弦型信號作載波的調制叫做連續(xù)波調制。調制后的載波就載有調制信號所包含的信息，稱為已調波。對于連續(xù)波調制，已調信號可以表示為：　　它有振幅頻率和相位三個參數構成。改變三個參數中的任何一個都可以攜帶同樣的信息。因此連續(xù)波的調制可

67、分為調幅、調相、和調頻。　　解調：解調是調制的逆過程，它的作用是從已調波信號中取出原來的調制信號。調制過程是一個頻譜搬移的過程，它將低頻信號的頻譜搬移到載頻位置。如果要接收端回復信號，就要從已調信號的頻譜中，將位于載頻的信號頻譜再搬回來。　　之所以進行解調，是因為無線電通信是通過空間輻射方式傳輸信號的，調制過程可以將信號的頻譜搬移到容易以電磁波形勢輻射的

68、較高頻率范圍；此外，調制過程可以將不同的信號通過頻譜搬移托付至不同頻率的載波上，實現多路復用，不至于互相干擾。　　實驗四連續(xù)系統的頻域分析　　一、實驗目的：　　掌握連續(xù)時間系統變換域分析的基本方法。 　　二、實驗

69、設備：　　安裝有matlab6.5以上版本的PC機一臺。　　三、實驗原理　　四、實驗內容：　　如圖所示系統：<

70、;p>　　對不同的RC值，用freqs函數畫出系統的幅頻曲線。　　易求得H(jω)=1/（1+jωRC），RC的取值依次為100、10、……、0.00001時的幅頻曲線。　　源程序：　　close<p&

71、gt;　　clear　　clc　　b=[0 1];　　for c=-5:2　　RC=10^c;&

72、lt;p>　　a=[RC 1];　　freqs(b,a);　　axis([10^(-2),10^(5),0.1,1]);　　hold on　　end&l

73、t;/p>　　得到一系列幅頻曲線，從左到右依次為RC的取值依次為100、10、……、0.00001時的幅頻曲線。圖中褐色虛線表示縱坐標取值為0.707，紅線表示橫坐標為100，綠線橫坐標為2000。　　(b) 信號包含了一個低頻分量和一個高頻分量。確定適當的RC值，濾除信號中的高頻分量并畫出信號和在s范圍內的波形。提示： |H( jω)|為最大值的/2處對應的頻率為通

74、帶截止頻率ωc,首先求取|H( jω)|并找到ωc和RC關系，然后根據題意選定ωc即可確定RC值。　　由（a）中的圖可知，當RC=-2時符合題意。　　源程序：　　close　　c

75、lear　　clc　　t=0:0.001:0.2;　　f=cos(100*t)+cos(2000*t);　　subplot(2,1,1);　　plot(t,f);</

76、p>　　y1=cos(100*t)/(1+j*100*10^(-2))+cos(2000*t)/(1+j*2000*10^(-2));　　subplot(2,1,2);　　plot(t,y1)　　2、信號任選，分析以下幾種情況下信號的頻譜和波形變化：<p&

77、gt;　　（1）系統滿足線性不失真條件時；　?。?）系統只滿足恒定幅值條件時；　　（3）系統只滿足相位條件時；　?。?）系統兩個條件均不滿足時。　　提示：利用fourier求取信號的傅立葉變換E(jω)，然后設計　　H( j

78、ω) = H( jω) eφ (ω )使之滿足不同條件，計算R( jω)= E(jω) H(jω)并畫頻譜圖。　?。?）源程序：　　clc　　clear<

79、;b>　　close　　syms t v;　　e=exp(-2*abs(t));　　subplot(2,3,1);　　ezplot(e,[-3,3]);　　axis(

80、[-3,3,-0.2,2]);　　Fe=fourier(e,v);　　subplot(2,3,2);　　ezplot(Fe,[-3,3]);　　title('幅度譜');　　axis([-3,3,0,2]);<

81、;/p>　　i=1;　　for a=-3:0.02:3　　R11=subs(Fe,v,a);　　C(i)=angle(R11);　　i=i+1;

82、　　end　　b=-3:0.02:3;　　subplot(2,3,3);　　plot(b,C);　　title('相位譜');　　axis([-3,3,-1,

83、1]);　　H1=2*exp(-j*v*1);　　R1=Fe*H1;　　r1=ifourier(R1,t);　　subplot(2,3,4);　　ezplot(r1,[-3,3]);</

84、p>　　title('滿足線性不失真條件');　　axis([-3,3,-0.2,2])　　subplot(2,3,5);　　ezplot('abs(8/(4+v^2)*exp(-i*v))');　　title(

85、'幅度譜');　　axis([-3,3,0,2.2]);　　i=1;　　for a=-3:0.02:3　　R11=subs(R1,v,a);　　C(i)=angle(R11);

86、　　i=i+1;　　end　　b=-3:0.02:3;　　subplot(2,3,6);　　plot(b,C);

87、　　title('相位譜');　　axis([-3,3,-3,3]);　?。?）源程序：　　clc　　clear<

88、;p>　　close　　syms t v;　　e=exp(-2*abs(t));　　subplot(2,3,1);　　ezplot(e,[-3,3]);<p&g

89、t;　　axis([-3,3,-0.2,2]);　　Fe=fourier(e,v);　　subplot(2,3,2);　　ezplot(Fe,[-3,3]);　　title('幅度譜');　　axis([-3,3,

90、0,2]);　　i=1;　　for a=-3:0.02:3　　R11=subs(Fe,v,a);　　C(i)=angle(R11);　　i=i+1;<

91、;/p>　　end　　b=-3:0.02:3;　　subplot(2,3,3);　　plot(b,C);　　title('相位譜');　　axis(

92、[-3,3,-1,1]);　　H1=(1-j*v)/(1+j*v);　　R1=Fe*H1;　　D=abs(R1);　　r1=ifourier(R1,t);　　subplot(2,3,4);&l

93、t;/p>　　ezplot(r1,[-3,3]);　　title('只滿足恒定幅值');　　axis([-3,3,-1,2]);　　subplot(2,3,5);　　ezplot('4*abs(1/(4+v^2)*(1-

94、i*v)/(1+i*v))');　　title('幅度譜');　　axis([-3,3,0,2]);　　subplot(2,3,6)　　i=1;　　for a=-3:0

95、.02:3　　R11=subs(R1,v,a);　　C(i)=angle(R11);　　i=i+1;　　end　　a=-3:0.02:3;</p&g

96、t;　　plot(a,C);　　title('相位譜');　　axis([-3,3,-3,3]);　?。?）源程序：　　clc<

97、p>　　clear　　close　　syms t v;　　e=exp(-2*abs(t));　　subplot(2,3,1);&

98、lt;p>　　ezplot(e,[-3,3]);　　axis([-3,3,-0.2,2]);　　Fe=fourier(e,v);　　subplot(2,3,2);　　ezplot(Fe,[-3,3]);　　title(

99、9;幅度譜');　　axis([-3,3,0,2]);　　i=1;　　for a=-3:0.02:3　　R11=subs(Fe,v,a);　　C(i)=angle(R11);</

100、p>　　i=i+1;　　end　　b=-3:0.02:3;　　subplot(2,3,3);　　plot(b,C);　　tit

101、le('相位譜');　　axis([-3,3,-1,1]);　　H1=v^2*exp(-j*v*1);　　R1=Fe*H1;　　r1=ifourier(R1,t);　　subp

102、lot(2,3,4);　　ezplot(r1,[-3,3]);　　title('滿足相位條件');　　axis([-3,3,-4,0.2])　　subplot(2,3,5);　　ezplot('R1'

103、;);　　title('幅度譜');　　axis([-3,3,-3,3]);　　i=1;　　for a=-3:0.02:3　　R11=subs(R1,v,a);</p&g

104、t;　　C(i)=angle(R11);　　i=i+1;　　end　　b=-3:0.02:3;　　subplot(2,3,6);

105、plot(b,C);　　axis([-3,3,-3,3]);　　title('相位譜');　?。?）源程序：　　clc

106、　clear　　close　　syms t v;　　e=exp(-2*abs(t));　　subplot(2,3,1);　　ezplo

107、t(e,[-3,3]);　　axis([-3,3,-0.2,2]);　　Fe=fourier(e,v);　　subplot(2,3,2);　　ezplot(Fe,[-3,3]);　　title('幅度譜');<

108、;/p>　　axis([-3,3,0,2]);　　i=1;　　for a=-3:0.02:3　　R11=subs(Fe,v,a);　　C(i)=angle(R11);

109、;　　i=i+1;　　end　　b=-3:0.02:3;　　subplot(2,3,3);　　plot(b,C);　　title('相位譜'

110、;);　　axis([-3,3,-1,1]);　　H1=v^2*(1-j*v)/(1+j*v);　　R1=Fe*H1;　　D=abs(R1);　　r1=ifourier(R1,t);</

111、p>　　subplot(2,3,4);　　ezplot(r1,[-3,3]);　　title('兩個條件均不滿足');　　axis([-3,3,-0.5,7]);　　subplot(2,3,5);<

112、p>　　ezplot('4*abs(1/(4+v^2)*v^2*(1-i*v)/(1+i*v))');　　title('幅度譜');　　axis([-3,3,0,2]);　　subplot(2,3,6)　　i=1;&

113、lt;/b>　　for a=-3:0.02:3　　R11=subs(R1,v,a);　　C(i)=angle(R11);　　i=i+1;　　end&

114、lt;/p>　　a=-3:0.02:3;　　plot(a,C);　　title('相位譜');　　axis([-3,3,-3,3]);　　五、思考題

115、　　1、連續(xù)系統頻域與復頻域分析的基本方法是什么？　　頻域分析方法：將激勵信號分解為正弦分量，找出聯系響應與激勵的系統函數H(jw)，求每一頻率分量的響應，從響應的頻譜函數R(jw)求傅里葉反變換從而求得響應r(t)。　　復頻域分析方法：通過拉普拉斯變換將時域中的積分微分方程變成復頻域中的代數方程，在復頻域中進行代數運算后則可得到系統響應的復頻域解，將此

116、解在經反變換則得到最終的時域解。　　2、若信號經過系統不發(fā)生失真，則對系統頻響有何要求？　　系統頻響的幅頻特性在整個頻率范圍中為一常數，即系統具有無限頻寬的響應均勻的通頻帶，系統的相頻特性是經過原點的直線。　　實驗五信號采樣與重建　　一、實驗目的<

117、;/b>　　1、深刻理解采樣定理的內容。　　2、掌握信號采樣后的頻譜。　　二、實驗設備　　安裝有matlab6.5 以上版本的PC 機一臺。　　三、實驗原理</

118、b>　　四、實驗內容：　　1、已知f (t) = Sa(2t)，以fs為采樣頻率，對f (t)進行采樣得到fa(t) ，觀　　察隨著fs 由小變大， fa(t)頻譜的變化，最后得出fa(t)與f (t)兩者頻譜

119、　　之間的關系。　　源程序：　　clc　　close　　clear<p

120、>　　syms t w;　　f=sin(2*t)/(2*t);　　subplot(3,2,1);　　ezplot(f,[-5,5]);　　title('時域信號');

121、ff=fourier(f);　　subplot(3,2,2);　　ezplot(ff,[-5,5]);　　title('頻譜信號');　　fs=1;fc=0;　　for n=-6:6<

122、p>　　fc=fc+(pi/2)*fs*(-heaviside(w-2-n*2*pi*fs)+heaviside(w+2-n*2*pi*fs));　　end　　subplot(3,2,4);　　ezplot(fc,[-40,40]);<p&

123、gt;　　title('采樣信號頻譜，采樣頻率1');　　fs=3;fc=0;　　for n=-3:3　　fc=fc+(pi/2)*fs*(-heaviside(w-2-n*2*pi*fs)+heaviside(w+2-n*2*pi*fs));<b&

124、gt;　　end　　subplot(3,2,5);　　ezplot(fc,[-40,40]);　　title('采樣信號頻譜，采樣頻率3');　　fs=5;fc=0;　　for n=-3:3&l

125、t;/p>　　fc=fc+(pi/2)*fs*(-heaviside(w-2-n*2*pi*fs)+heaviside(w+2-n*2*pi*fs));　　end　　subplot(3,2,6);　　ezplot(fc,[-40,40]);<

126、/p>　　title('采樣信號頻譜，采樣頻率5');　　fs=0.5;fc=0;　　for n=-3:3　　fc=fc+(pi/2)*fs*(-heaviside(w-2-n*2*pi*fs)+heaviside(w+2-n*2*pi*fs));

127、　　end　　subplot(3,2,3);　　ezplot(fc,[-40,40]);　　title('采樣信號頻譜，采樣頻率0.5');　　2、由實驗1 中采樣得到的離散信號重建對應的連續(xù)時間

128、信號：情況　　一、f s ≥ 2f m；情況二、f s < 2f m。提示根據內插公式重建連續(xù)時間信　　號，或者根據采樣后頻譜通過低通濾波器的形式重建連續(xù)時間信號。　　源程序：　　clc&l

129、t;/p>　　close　　clear　　syms t w;　　f=sin(2*t)/(2*t);　　subplot(4,2,1)

130、;　　ezplot(f,[-5,5]);　　title('時域信號');　　ff=fourier(f);　　subplot(4,2,2);　　ezplot(ff,[-5,5]);<p&

131、gt;　　title('頻譜信號');　　fs=3;fc=0;　　for n=-3:3　　fc=fc+(pi/2)*fs*(-heaviside(w-2-n*2*pi*fs)+heaviside(w+2-n*2*pi*fs));　　end

132、　　subplot(4,2,3);　　ezplot(fc,[-40,40]);　　title('采樣信號頻譜，fs≥2fm(fs=3)');　　G=-heaviside(w-4)+heaviside(w+4);&

133、lt;p>　　Fc=fc*G;　　f1=ifourier(Fc)/fs;　　subplot(4,2,5);　　ezplot(Fc,[-40,40]);　　title('重建信號頻譜,fs≥2fm(fs=3)');

134、　　subplot(4,2,7);　　ezplot(f1,[-5,5]);　　title('重建時域信號,fs≥2fm(fs=3)');　　fs=0.5;fc=0;　　for n=-3:3&

135、lt;p>　　fc=fc+(pi/2)*fs*(-heaviside(w-2-n*2*pi*fs)+heaviside(w+2-n*2*pi*fs));　　end　　subplot(4,2,4);　　ezplot(fc,[-40,40]);<

136、;p>　　title('采樣信號頻譜，fs<2fm(fs=0.5)');　　G=-heaviside(w-4)+heaviside(w+4);　　Fc=fc*G;　　f1=ifourier(Fc)/fs;

137、　subplot(4,2,6);　　ezplot(Fc,[-40,40]);　　title('重建信號頻譜,fs<2fm(fs=0.5)');　　subplot(4,2,8);　　ezplot(f1,[-5,5]);<p

138、>　　title('重建時域信號,fs<2fm(fs=0.5)');　　五、思考題　　1、隨著采樣頻率 fs 從小到大變化，fa(t) 的頻譜發(fā)生怎樣的變化，與　　f (t)頻譜之間的關系如何？與理論計算結果之間是否完全一致？如

139、　　果不一致，請分析可能是什么原因導致的？　　隨著采樣頻率 fs 從小到大變化，fa(t) 的頻譜剛開始時出現混疊現象，之后會出現不相互疊加的頻譜，并且相鄰頻譜之間的間隔越來越大，幅度會增大一些。　　2、采樣頻率f s 分別滿足情況一與情況二時，由fa (t)重建的(t)與原信號f (t)是否完全相同？如果不相同，試分析原因。&l

140、t;/p>　　情況一：重建的圖像衰減一定幅度得到的圖像和原圖像完全相同。　　情況二：不一樣。因為fa (t) 出現混疊現象，使其解調后的圖像與原圖像不同。　　實驗六離散時間信號和系統分析　　一、實驗目的

141、　　1. 掌握常用離散信號的MATLAB 表示方法。　　2、掌握用MATLAB 計算卷積和及零狀態(tài)響應的方法。　　二、實驗設備　　安裝有matlab6.5 以上版本的PC 機一臺。　　三、實驗原理 </b&g

142、t;　　四、實驗內容　　1、用MATLAB表示離散信號：，Asin(k)。　　源程序：　　clc

143、　close　　clear　　k=-3:3;　　fk=2.^(k);　　subplot(1,2,1);　　stem(k,fk,'

144、;filled');　　axis([-4,4,0,9]);　　h=-7:7;　　fh=2*sin(h);　　subplot(1,2,2);　　stem(h,fh,'filled

眾賞文庫> 全部分類> 畢業(yè)設計

溫馨提示

1. 本站所有資源如無特殊說明，都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
2. 本站的文檔不包含任何第三方提供的附件圖紙等，如果需要附件，請聯系上傳者。文件的所有權益歸上傳用戶所有。
3. 本站RAR壓縮包中若帶圖紙，網頁內容里面會有圖紙預覽，若沒有圖紙預覽就沒有圖紙。
4. 未經權益所有人同意不得將文件中的內容挪作商業(yè)或盈利用途。
5. 眾賞文庫僅提供信息存儲空間，僅對用戶上傳內容的表現方式做保護處理，對用戶上傳分享的文檔內容本身不做任何修改或編輯，并不能對任何下載內容負責。
6. 下載文件中如有侵權或不適當內容，請與我們聯系，我們立即糾正。
7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

信號與線性系統課程設計

文檔簡介

溫馨提示

最新文檔

評論

信號與線性系統課程設計

文檔簡介

溫馨提示

最新文檔

評論

免費下載