【正文】 階數(shù)越高，元件越復雜，成本越高，應根據(jù)系統(tǒng)對濾波器截止特性的需求（最大通帶起伏，最小阻帶衰減等）選擇合適的階數(shù)。 具有陷波作用，相位也在陷波頻率上發(fā)生反轉(zhuǎn)。 %模擬濾波器定階 得到n=11。s39。 ，設計指標為：截止頻率 f=10 kHz，即 Wp= 2 * pi * 10k krad/s，阻帶為0~ 8 kHz ，通帶波紋 Rp 3 Rs 60 dB。grid。)。ylabel(39。Frequency w(rad/s)39。)。 % 畫出幅頻曲線end title(39。 plot(w1,magH1)。 magH1=20*log10(magH1)。 % 將傳遞函數(shù)由零極點式轉(zhuǎn)化成多項式形式 [H1,w1]=freqs(b,a,w)。 for n=1:1:8 [z,p,k]=buttap(n)。n=0。grid。一階全通濾波器39。bode(b,a)。a=[1, 1]。grid。帶阻濾波器39。bode(b,a)。a=[1, 4 , 1]。grid。二階帯通濾波器39。bode(b,a)。a=[1,3,2]。grid。Chebyshev Highpass39。 bode(bt,at)。 %零極點式轉(zhuǎn)化為多項式構(gòu)成的傳遞函數(shù)形式 [bt,at]= lp2hp(b,a,W0)。 %模擬濾波器定階 [z,p,k]=cheb1ap(n,Rp)。s39。Rs=60。Ws=2*pi*8e3。 % 高通濾波器% Chebyshev Highpass% 截止頻率 f=10 kHz，即 Wp= 2 * pi * 10k krad/s，% 阻帶為0~ 8 kHz ，通帶波紋 Rp 3 Rs 60 dB clear。)。title(39。figure(3)。 % 三階Butterworth濾波器 b=[0,1]。)。title(39。 % 將歸一化頻率對應的多項式系數(shù)轉(zhuǎn)成截止頻率W0所對應的多項式系數(shù)figure(2)。 % Butterworth 低通濾波器，階數(shù)[b,a]=zp2tf(z,p,k)。W0=1。grid。二階Butterworth低通濾波器39。bode(b,a)。a=[1,2^(1/2),1]。繪制波特圖。繪制波特圖，輻頻曲線和相頻曲線。二、實驗原理及方法本實驗所研究的系統(tǒng)函數(shù)H(s)是有理函數(shù)形式，分子、分母分別是m、n階多項式。27實驗八 傳輸函數(shù)與系統(tǒng)特性一、實驗目的學習由系統(tǒng)函數(shù)確定系統(tǒng)頻率特性的方法。26（2）執(zhí)行結(jié)果（3）分析： ，原信號最高頻率wm=2，,即fm=，當采樣頻率fs2fm時（），信號頻譜結(jié)構(gòu)與原信號頻譜結(jié)構(gòu)不同，信號發(fā)生混疊，波形失真，采樣不理想；只有采樣頻率fs2fm時(圖中采樣頻率3)，采樣后的信號頻譜結(jié)構(gòu)與原信號相同，波形不失真，采樣成功。重建時域信號,fs2fm(fs=)39。 subplot(4,2,8)。重建信號頻譜,fs2fm(fs=)39。 subplot(4,2,6)。 Fc=fc*G。)。ezplot(fc,[40,40]) title(39。for n=3:3 fc=fc+(pi/2)*fs*(heaviside(w2n*2*pi*fs)+heaviside(w+2n*2*pi*fs))。 % 情況二：fs2*fm % fs=。重建時域信號，fs=2fm(fs=3)39。 subplot(4,2,7)。重建信號頻譜,fs=2fm(fs=3)39。 %原信號的頻譜Fourier變換subplot(4,2,5)。 %低通濾波器Fc=fc*G。)。 ezplot(fc,[40,40])title(39。 for n=3:3 %采樣后的信號的頻域表示 fc=fc+(pi/2)*fs*(heaviside(w2n*2*pi*fs)+heaviside(w+2n*2*pi*fs))。 % 情況一：fs2*fm % fs=3。頻域信號39。 subplot(4,2,2)。)。ezplot(f,[5,5])title(39。 f=sin(2*t)/(2*t)。clear。（1）源程序%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 由1中采樣得到的離散信號重建% 對應的連續(xù)時間信號：情況一、 ；情況二、% 提示 根據(jù)內(nèi)插公式重建連續(xù)時間信號，或者% 根據(jù)采樣后頻譜通過低通濾波器的形式% 重建連續(xù)時間信號。 由1中采樣得到的離散信號重建對應的連續(xù)時間信號：情況一；情況二、。)。ezplot(fc,[40,40])title(39。 for n=3:3 fc=fc+(pi/2)*fs*(heaviside(w2n*2*pi*fs)+heaviside(w+2n*2*pi*fs))。 fs=5。采樣信號頻譜，采樣頻率339。end subplot(3,2,5)。fc=0。)。ezplot(fc,[40,40])title(39。for n=6:6 fc=fc+(pi/2)*fs*(heaviside(w2n*2*pi*fs)+heaviside(w+2n*2*pi*fs))。 fs=1。采樣信號頻譜，39。 end subplot(3,2,3)。 %采樣頻率fsfc=0。)。 ezplot(ff,[5,5]) title(39。 ff=fourier(f)。時域信號39。 subplot(3,2,1)。syms t w。%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% clear。三、實驗內(nèi)容：已知，以為采樣頻率，對進行采樣得到，觀察隨著由小變大，頻譜的變化，最后得出與兩者頻譜之間的關系。掌握信號采樣后的頻譜。注意：當不滿足線性不失真條件時（情況2,3,4），可能會產(chǎn)生新的頻率分量，波形會失真，不能用y(t)= |H(jw)|*cos(wt+θ+∠|H(jw)|) 來算。.，信號在通過兩個條件均不滿足的系統(tǒng)后，波形失真，幅度譜和相位譜均改變。 ，信號通過只滿足幅值條件 的系統(tǒng)，波形失真，幅度譜不變，相位譜變化。 axis([3,3,3,3])。相位譜39。enda=3::3。 C(i)=angle(R11)。 subplot(2,3,6)i=1。)。)title(39。 ezplot(39。axis([3,3,7])。兩個條件均不滿足時39。 subplot(2,3,4)。 D=abs(R1)。 H1=v^2*(1j*v)/(1+j*v)。相位譜39。 subplot(2,3,3)。 i=i+1。for a=3::3 R11=subs(Fe,v,a)。 axis([3,3,0,2])。幅度譜39。 subplot(2,3,2)。 ezplot(e,[3,3])axis([3,3,2])。 figure(4)。syms t v。 %（4）兩個條件均不滿足時clear。相位譜39。plot(b,C) axis([3,3,3,3])。 end b=3::3。 C(i)=angle(R11)。 i=1。)。)title(39。 ezplot(39。)。 ezplot(r1,[3,3]) title(39。 r1=ifourier(R1,t)。 H1=v^2*exp(j*v*1)。)。 plot(b,C)title(39。 endb=3::3。 C(i)=angle(R11)。i=1。)。 ezplot(Fe,[3,3]) title(39。Fe=fourier(e,v)。 subplot(2,3,1)。e=exp(2*abs(t))。close。axis([3,3,3,3])。相位譜39。 end a=3::3。 C(i)=angle(R11)。 subplot(2,3,6) i=1。)。)title(39。 ezplot(39。 axis([3,3,1,2])。只滿足恒定幅值條件時39。 subplot(2,3,4)。D=abs(R1)。 H1=(1j*v)/(1+j*v)。)。 plot(b,C)title(39。 end b=3::3。 C(i)=angle(R11)。 i=1。)。 ezplot(Fe,[3,3])title(39。 Fe=fourier(e,v)。subplot(2,3,1)。e=exp(2*abs(t))。close。axis([3,3,3,3])。相位譜39。 subplot(2,3,6)。 i=i+1。 %畫輸出信號相位譜for a=3::3 R11=subs(R1,v,a)。 axis([3,3,0,])。幅度譜39。abs(8/(4+v^2)*exp(i*v))39。 axis([3,3,2])。滿足線性無失真條件時39。 %輸出信號r1=ifourier(R1,t)。 H1=2*exp(j*v*1)。)。 subplot(2,3,3)plot(b,C) title(39。 i=i+1。 %數(shù)值法畫相位譜for a=3::3 R11=subs(Fe,v,a)。 axis([3,3,0,2])。幅度譜39。 subplot(2,3,2)。 %第一幅圖像ezplot(e,[3,3]) %ezplot畫符號函數(shù)的幅度譜axis([3,3,2])。 %輸入信號subplot(2,3,1)。syms t v。（1）源程序：20%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 線性不失真條件%分析以下幾種情況下信號的頻譜和波形變化 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %取信號exp(2*abs(t)) %(1)滿足線性無失真條件時clear。信號任選，分析以下幾種情況下信號的頻譜和波形變化：(1)系統(tǒng)滿足線性不失真條件時；(