【正文】 ,10,10])。axis([30,10,10,10])。axis([20,20,10,10])。axis([20,20,10,10])。axis([20,20,10,10])。 %創(chuàng)建繪圖窗口1 % 子圖1：作出原函數(shù)圖f（t）subplot(2,3,1)。grid on3（2）執(zhí)行結(jié)果編寫如圖3的函數(shù)并用MATLAB繪出滿足下面要求的圖形。 %（自變量t，脈寬w）plot(t,f)。)。抽樣信號Sa(t)39。) 。omega0=pi。xlabel(39。grid on%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% f=A*cos(omega0*t+phi) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%t=10::10。f(t)=A*exp(alpha*t)39。用MATLAB表示抽樣信號(sinc(t))、矩形脈沖信號(rectpuls(t, width))及三角脈沖信號(tripuls(t, width, skew))。t=10::10。ylabel(39。plot(t,f)title(39。) 。)。（1）源程序%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% f=sinc(t) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%t=20::20。 %（自變量t，脈寬w）plot(t,f)。grid on%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 三角脈沖信號 y=tripuls(t,width, skew) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% t=10::10。t39。 %使用單位矩形脈沖和單位三角脈沖來對圖象進(jìn)行描述f2=3*tripuls(t6,w2)。xlabel(39。xlabel(39。xlabel(39。xlabel(39。xlabel(39。5（1）源程序%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % （1）f= a.^k.*stepfun(k,0)。 axis([1,10, 0, 2])。f39。filled39。k39。 %設(shè)定離散自變量u=stepfun(k,0)。 title(39。f39。f=f[k2]*u[k]。) 。)。stem(k+2,f,39。k39。)。ylabel(39。 %起點(diǎn)的和是起點(diǎn)，終點(diǎn)的和是終點(diǎn)[0,5]+[0,5]=[0,10]stem(k,f,39。k39。二、實(shí)驗(yàn)內(nèi)容1. 分別用函數(shù)lsim和卷積積分兩種方法求如圖8所示系統(tǒng)的零狀態(tài)響應(yīng)。b=[1]。Yzs(t) 39。 %用卷積的定義式計算系統(tǒng)的零狀態(tài)響應(yīng)ezplot(i,[0 5])。Yzs(t) 39。b=[1]。系統(tǒng)的沖激響應(yīng)h(t)39。step(b,a,5)xlabel(39。)。h=().^k.*stepfun(k,0)。零狀態(tài)響應(yīng)39。x=ones(1,length(n))。 + b1*x(iN:iN+M)39。ylabel(39。%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%clear。*t2)。)。 subplot(2,1,2)。title(39。輸入疊加的諧波次數(shù):39。f(t)39。w(rad)39。 %|Fw|ezplot(FFw)。調(diào)制信號為一取樣信號，利用MATLAB分析幅度調(diào)制(AM)產(chǎn)生的信號頻譜，比較信號調(diào)制前后的頻譜并解調(diào)已調(diào)信號。xlabel(39。)。Fx39。 % 調(diào)制后信號的數(shù)值表示subplot(4,2,3)。title(39。 %作圖，調(diào)制后的信號頻譜xlabel(39。)。)。subplot(4,2,6)。)。 ff2=subs(f2,t,t1)。title(39。ylabel(39。在接收端，將收到的信號乘以解調(diào)信號，同樣是cos(w0 t)，頻譜搬移，注意在零頻有頻譜疊加 ，因此幅度仍為原信號的1/2（子圖6）。（1）源程序：18%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 畫出系統(tǒng)的幅頻曲線 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%系統(tǒng)的頻率響應(yīng)函數(shù)為 H(jw)=1/(jw*RC+1) clear。 H=1/(i*w*(1e3)+1)。t/s39。 subplot(2,1,2)。title(39。（1）源程序：20%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 線性不失真條件%分析以下幾種情況下信號的頻譜和波形變化 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %取信號exp(2*abs(t)) %(1)滿足線性無失真條件時clear。 subplot(2,3,2)。 i=i+1。 %輸出信號r1=ifourier(R1,t)。幅度譜39。 subplot(2,3,6)。e=exp(2*abs(t))。)。 plot(b,C)title(39。 subplot(2,3,4)。)title(39。 end a=3::3。e=exp(2*abs(t))。)。 plot(b,C)title(39。 ezplot(r1,[3,3]) title(39。)。plot(b,C) axis([3,3,3,3])。 figure(4)。 axis([3,3,0,2])。相位譜39。兩個條件均不滿足時39。)。相位譜39。注意：當(dāng)不滿足線性不失真條件時（情況2,3,4），可能會產(chǎn)生新的頻率分量，波形會失真，不能用y(t)= |H(jw)|*cos(wt+θ+∠|H(jw)|) 來算。syms t w。 ezplot(ff,[5,5]) title(39。采樣信號頻譜，39。)。 fs=5。 由1中采樣得到的離散信號重建對應(yīng)的連續(xù)時間信號：情況一；情況二、。ezplot(f,[5,5])title(39。 % 情況一：fs2*fm % fs=3。 %低通濾波器Fc=fc*G。重建時域信號，fs=2fm(fs=3)39。)。 subplot(4,2,8)。二、實(shí)驗(yàn)原理及方法本實(shí)驗(yàn)所研究的系統(tǒng)函數(shù)H(s)是有理函數(shù)形式，分子、分母分別是m、n階多項式。bode(b,a)。 % Butterworth 低通濾波器，階數(shù)[b,a]=zp2tf(z,p,k)。 % 三階Butterworth濾波器 b=[0,1]。 % 高通濾波器% Chebyshev Highpass% 截止頻率 f=10 kHz，即 Wp= 2 * pi * 10k krad/s，% 阻帶為0~ 8 kHz ，通帶波紋 Rp 3 Rs 60 dB clear。 %模擬濾波器定階 [z,p,k]=cheb1ap(n,Rp)。grid。grid。grid。grid。 magH1=20*log10(magH1)。Frequency w(rad/s)39。 ，設(shè)計指標(biāo)為：截止頻率 f=10 kHz，即 Wp= 2 * pi * 10k krad/s，阻帶為0~ 8 kHz ，通帶波紋 Rp 3 Rs 60 dB。階數(shù)越高，元件越復(fù)雜，成本越高，應(yīng)根據(jù)系統(tǒng)對濾波器截止特性的需求（最大通帶起伏，最小阻帶衰減等）選擇合適的階數(shù)。s39。ylabel(39。 plot(w1,magH1)。n=0。a=[1, 1]。a=[1, 4 , 1]。a=[1,3,2]。 %零極點(diǎn)式轉(zhuǎn)化為多項式構(gòu)成的傳遞函數(shù)形式 [bt,at]= lp2hp(b,a,W0)。Ws=2*pi*8e3。figure(3)。 % 將歸一化頻率對應(yīng)的多項式系數(shù)轉(zhuǎn)成截止頻率W0所對應(yīng)的多項式系數(shù)figure(2)。二階Butterworth低通濾波器39。繪制波特圖，輻頻曲線和相頻曲線。重建時域信號,fs2fm(fs=)39。 Fc=fc*G。 % 情況二：fs2*fm % fs=。 %原信號的頻譜Fourier變換subplot(4,2,5)。 for n=3:3 %采樣后的信號的頻域表示 fc=fc+(pi/2)*fs*(heaviside(w2n*2*pi*fs)+heaviside(w+2n*2*pi*fs))。)。（1）源程序%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 由1中采樣得到的離散信號重建% 對應(yīng)的連續(xù)時間信號：情況一、 ；情況二、% 提示 根據(jù)內(nèi)插公式重建連續(xù)時間信號，或者% 根據(jù)采樣后頻譜通過低通濾波器的形式% 重建連續(xù)時間信號。 for n=3:3 fc=fc+(pi/2)*fs*(heaviside(w2n*2*pi*fs)+heaviside(w+2n*2*pi*fs))。fc=0。 fs=1。)。 subplot(3,2,1)。掌握信號采樣后的頻譜。 axis([3,3,3,3])。 subplot(2,3,6)i=1。axis([3,3,7])。 H1=v^2*(1j*v)/(1+j*v)。for a=3::3 R11=subs(Fe,v,a)。 ezplot(e,[3,3])axis([3,3,2])。相位譜39。 i=1。)。)。i=1。 subplot(2,3,1)。相位譜39。)。只滿足恒定幅值條件時39。)。 i=1。subplot(2,3,1)。相位譜39。 axis([3,3,0,])。滿足線性無失真條件時39。 subplot(2,3,3)plot(b,C) title(39。幅度譜39。syms t v。)。xlabel(39。ylabel(39。fy=simplify(ifourier(Fy,t))。 a