【正文】 ，時域信號 fy1=subs(fy,t,t1)。)。 %輸入信號 subplot(2,3,1)。)。 %畫輸出信號相位譜 for a=3::3 R11=subs(R1,v,a)。 Fe=fourier(e,v)。 H1=(1j*v)/(1+j*v)。 subplot(2,3,6) i=1。 Fe=fourier(e,v)。 H1=v^2*exp(j*v*1)。 C(i)=angle(R11)。 subplot(2,3,2)。 D=abs(R1)。 C(i)=angle(R11)。 三、實驗內(nèi)容： 已知 ，以 為采樣頻率，對 進行采樣得到 ，觀察隨著由小變大， 頻譜的變化，最后得出 與 兩者頻譜之間的關(guān)系。 %采樣頻率 fs fc=0。 end subplot(3,2,5)。clear。 ezplot(fc,[40,40]) title(39。重建時域信號， fs=2fm(fs=3)39。 % 情況一： fs2*fm % fs=3。 由 1 中采樣得到的離散信號重建對應(yīng)的連續(xù)時間信號：情況一、 ms ff 2? ；情況二、 ms ff 2? 。)。 ezplot(ff,[5,5]) title(39。 注意：當(dāng)不滿足線性不失真條件時（情況 2,3,4），可能會產(chǎn)生新的頻率分量，波形會失真，不能用 y(t)= |H(jw)|*cos(wt+θ +∠ |H(jw)|) 來算。)。相位譜 39。 figure(4)。)。 plot(b,C) title(39。 e=exp(2*abs(t))。) title(39。 plot(b,C) title(39。 e=exp(2*abs(t))。幅度譜 39。 i=i+1。 101100 1 0 0 8 0 6 0 4 0 2 0F r e q u e n c y ( r a d / s )Phase (degrees)101100102101100F r e q u e n c y ( r a d / s )Magnitude0 0 . 0 2 0 . 0 4 0 . 0 6 0 . 0 8 0 . 1 0 . 1 2 0 . 1 4 0 . 1 6 0 . 1 8 0 . 221012t/sfx1 輸入信號 f( t ) 0 0 . 0 2 0 . 0 4 0 . 0 6 0 . 0 8 0 . 1 0 . 1 2 0 . 1 4 0 . 1 6 0 . 1 8 0 . 221012t/sfy1 輸出信號 y ( t ) )()()( ???? ejHjH ? 20 （ 1）源程序： %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % 線性不失真條件 %分析以下幾種情況下信號的頻譜和波形變化 % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %取信號 exp(2*abs(t)) %(1)滿足線性無失真條件時 clear。 subplot(2,1,2)。 f ( t )+RC y ( t )( ) c os( 10 0 ) c os( 20 00 )f t t t??)(tf )(ty ~0?t)( ?jH 22 c?)( ?jH c? c? 19 H=1/(i*w*(1e3)+1)。 在接收端，將收到的信號乘以解調(diào)信號，同樣是 cos(w0 t)，頻譜搬移，注意在零頻有頻譜疊加 ，因此幅度仍為原信號的 1/2（子圖 6）。 title(39。)。)。 %作圖，調(diào)制后的信號頻譜 xlabel(39。 % 調(diào)制后信號的數(shù)值表示 subplot(4,2,3)。)。 4 3 2 1 0 1 2 3 40 . 0 30 . 0 40 . 0 50 . 0 60 . 0 70 . 0 80 . 0 90 . 10 . 1 1w ( r a d )1 / abs ( 3 + w i )2|F(w)|3 2 1 0 1 2 3123456789w ( r a d )2 / abs ( w )|F(w)| 16 調(diào)制信號為一取樣信號，利用 MATLAB 分析幅度調(diào)制 (AM)產(chǎn)生的信號頻譜，比較信號調(diào)制前后的頻譜并解調(diào)已調(diào)信號。w(rad)39。輸入疊加的諧波次數(shù) :39。 subplot(2,1,2)。*t2)。 ylabel(39。 x=ones(1,length(n))。 h=().^k.*stepfun(k,0)。 step(b,a,5) xlabel(39。b=[1]。 %用卷積的定義式計算系統(tǒng)的零狀態(tài)響應(yīng) ezplot(i,[0 5])。b=[1]。k39。) 。 ylabel(39。k39。)。)。f39。 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% clear。 %設(shè)定自變量 u=stepfun(k,0)。 axis([30,10,10,10])。 axis([20,20,10,10])。 %創(chuàng)建繪圖窗 2 0 1 5 1 0 5 0 5 10 15 20 0 . 500 . 51抽樣信號 S a ( t )2 1 . 5 1 0 . 5 0 0 . 5 1 1 . 5 2202矩形脈沖信號 f ( t ) = r e c t p u l s ( t , w )t2 1 . 5 1 0 . 5 0 0 . 5 1 1 . 5 2202三角脈沖信號 f ( t ) = t r i p u l s ( t , w )t 4 口 1 % 子圖 1：作出原函數(shù)圖 f（ t） subplot(2,3,1)。)。 w=1。 plot(t,f) title(39。phi=pi/3。 (1)源程序 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % f=A*exp(alpha*t) % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% A=1。alpha=。 f=A*cos(omega0*t+phi)。f=A*sin(omega0*t+phi)39。 %脈寬為 1 f=rectpuls(t,w)。 xlabel(39。 plot(t,f)。 grid on % 子圖 3： 作出 f（ t2） subplot(2,3,3)。 grid on % 子圖 5： 作出 f（ *t+1） subplot(2,3,5)。 %設(shè)定單位階躍序列 f= a.^k.*stepfun(k,0)。 A=。) 。 xlabel(39。 xlabel(39。) 。f39。 grid on （ 2）執(zhí)行結(jié)果 如上圖，繪出了指定的序列 波形 若 [ ] {1,1,1,1, 0, 0}xn ? ，計算 [ ] [ ]xn xn? 。) 。 y=lsim(b,a,f,t)。 %符號函數(shù)作圖 xlabel(39。 impulse(b,a,5) xlabel(39。Time(sec)39。 e=[1 1 1 1]。 % recur function begin N=length(a)。y[n]39。 Fn=*trapz(gt1,2)*t1(2)*trapz(gt2,2)*t2(2) %求出傅里葉級數(shù)的系數(shù)（復(fù)指數(shù)形式） %梯形積分法，積分運算近似成求和計算 ,trapz(被積函數(shù)，維數(shù) ) %代表每一點 t的函數(shù)值乘以 delta t的疊加，這里 gt1為 10行 100列的二維數(shù)組， %deltat用 t1(2)表示 ， t1數(shù)組的第二個元素，即為(10)/100= sig=abs(Fn)。 stem(n,angle(Fn))。) %輸入疊加的諧波次數(shù) for k=1:Mark y=y+2/(k*pi)*(1cos(k*pi))*sin(k*pi*t)。)。設(shè)載波信號的頻率為 100Hz。 Fx=fourier(x)。 % 畫調(diào)制后信號的時域圖像 plot(t1,y1)。w (rad/s)39。 ylabel(39。 title(39。低通取出的信號波形 39。然后 通過一個幅值為 2 截止頻率為 ws（ wm ws w0wm）的低通濾波器后取出原信號頻譜，再轉(zhuǎn)化成時域信號，即完成了信號重建。 %低通濾波 Fy=Fx*H。 plot(t1,fy1)。 close 。 end b=3::3。)。 figure(2)。相位譜 39。幅度譜 39。 figure(3)。相位譜 39。 axis([3,3,3,3])。 subplot(2,3,1)。) axis([3,3,1,1])。 axis([3,3,0,2])。 2 0 200 . 511 . 52te x p ( 2 abs ( t ))2 0 200 . 511 . 52v幅度譜2 0 21 0 . 500 . 51相位譜2 0 243210t滿足相位條件2 0 2202R1幅度譜2 0 2202相位譜2 0 200 . 511 . 52te x p ( 2 abs ( t ))2 0 200 . 511 . 52v幅度譜2 0 21 0 . 500 . 51相位譜2 0 20246t兩個條件均不滿足時2 0 200 . 511 . 52v幅度譜2 0 2202相位譜 24 實驗七 信號采樣與重建 一、實驗?zāi)康? 深刻理解采樣定理的內(nèi)容。頻域信號 39。 fs=3。提示 根據(jù)內(nèi)插公式重建連續(xù)時間信號，或者根據(jù)采樣后頻譜通過低通濾波器的形式重建連續(xù)時間信號。fc=0。)。 end subplot(4,2,3)。 % 5 0 500 . 51t時域信號5 0 500 . 511 . 5w頻域信號 4 0 2 0 0 20 4000 . 511 . 5w采樣信號頻譜，采樣頻率 0 . 5 4 0 2 0 0 20 4000 . 511 . 5w采樣信號頻譜，采樣頻率 1 4 0 2 0 0 20 40024w采樣信號頻譜，采樣頻率 3 4 0 2 0 0 20 4002468w采樣信號頻譜，采樣頻率 5 26 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % 原信號 % close。 for n=3:3 fc=fc+(pi/2)*fs*(heaviside(w2n*2*pi*fs)+heaviside(w+2n*2*pi*fs))。 fs=。 二、實驗設(shè)備 安裝有 以上版本的 PC 機一臺。 for a=3::3 R11=subs(R1,v,a)。 R1=Fe*H1。 Fe=fourier(e,v)。 for a=3::3 R11=subs(R1,v,a)。 axis([3,3,1,1])。 ezplot(e,[3,3]) axis([3,3,2])。 axis([3,