1、tite（） 系统单位脉冲响应h1（n）sublot（2,2,2）;y=y1（）te（y1,y）;（e）1（n）与R（）的卷积y21（n）subot（2, ，3）; yh2（n）;tstem（2n,y）; %调用函数se绘图tte（） 系统单位脉冲响应h2（n）suplo（2，2, 4）; yy22（n）e（y22n, ）;（g）2（n）与R（n）的卷积22（n）=%内容3： 谐振器分析u=ones（1, 256）; %产生信号n=0：255;snsn（0.04*n）+si（0.4*n）； %产生正弦信号A=,-123,0.801；=1/100,0,-1/100.9; %系统差分方程系数向量B
2、和A31n=filtr（B，A,n）; %谐振器对u的响应y31n2n=filr（,，xi）;%谐振器对正弦信号的响应y3figur（3）subplo（,1,1）;y3（）；tstem（y31n,y）titl（h） 谐振器对（n）的响应1（n）sult（2,1,）;y3（n）;s（,y）;tl（i） 谐振器对正弦信号的响应y32（n）function tstem（xn,） = 0：leng（xn）;st（n,xn,.）；xlabe（n）;yll（yn%xlabel（ ）;ylabe（yn）;is（0,n（ed）,min（x）,1.2*x（xn）；实验二时域采样理论验证程序exp2a.Tp641
3、00;观察时间Tp=64微秒%产生长采样序列x（）% Fs=00; T=1Fs;Fs=1000; =1s;M=Tp*Fs； n=0：M-1;A=444.18;alph=pi0*20.5；omega=i*205;xn=Axp（-alph*T）*in（omen*T）;k=T*fft（xnt,M）;M点FFTxnt）ynxa（T）sbplt（,2,）;tsem（xnt,n）; %调用自编绘图函数tsem绘制序列图bo on； titl（a） F=100Hzk0:M-;fk=k/Tp;subpo（3,2,）；plt（fk,abs（Xk）；（a）TFTxa（nT），Fs=1000Hlabel（f（z）
4、yabe（幅度xis（0,Fs,0，1.*ma（abs（X）= Fs=300z和 Fs=200Hz的程序与上面s=100z的程序完全相同。% s30Hzp64/1000;%观察时间p=64微秒%产生M长采样序列x（n） s=10; T=1/Fs;s=30； T1/Fs;M=Tp*Fs; n=0:A=.128; lph=pi*0*205; mega=pi*5020.5；xn=*x（alph*T）sin（omega*n*T）;Xk=T*fft（nt,M）;M点FFTxt）y=a（nT）;fgre（）;sbpot（,2,1）;tstem（xt,yn）;%调用自编绘图函数tse绘制序列图ox on;
5、tte（b） F=300Hk=0:-1；fkkTp;sublo（3,2,2）; pot（fk,abs（X）;tte（b） T*FTa（）,Fs=1000Hxlab（f（H） b（axis（0,Fs,0,1.2*ax（as（X）% fs=200HzTp=64/000;%观察时间p=4微秒%产生M长采样序列（n）% s=100； =1/Fs;s=00;=1Fs；MT*F； =:M1；A=448; ah=pi0*20.5; omeg=pi*50205；xnA*exp（-alph*nT）.*sn（*n*）;Xk=Tfft（xnt,）;%M点FTxn）n=x（n）iure（3）;splot（3,2,1）
6、;tstem（n，y）; %调用自编绘图函数em绘制序列图box n; tite（c） F=200Hz0:-； k=k/Tp;spot（3,2,2）； plot（f,bs（Xk）;title（） TFTxa（），Fs100Hzxlae（f（） ylabel（幅度ais（0,Fs,0,12max（abs（k）频域采样理论验证程序exp2.mM=7;N=32;=0:M;产生M长三角波序列x（n）xa=0:flo（M/2）; x= il（M/2）-:-:; n=a,b;Xkft（x，1024）;%1024点FTx（n）, 用于近似序列x（n）的TFX2kff（n,32）;32点FFTx（n）3n=f
7、t（X32k）;%32点FTX2（k）得到x3（n）X16k=X32k（1:：N）;%隔点抽取X2k得到1（K）16n=ifft（X16k,/）; 1点ITX16（）得到x16（n）ublt（3，2,2）;te（n,xn,.b on（b）三角波序列x（n）;xlael（nyael（n）xi（,0）1023;wk=2*k1024；%subpt（3，,1）;lot（w,abs（Xk）；（a）Tx（n）;lab（ogpi）;ylbl（|X（ejomeg）|ax（,1,0，20）/2-1；ubpot（3,2,3）;ste（k,bs（X16），）;box otitl（） 16点频域采样xla（k）;yl
8、al（|X_1_6（k）|）;xis（0,，0,200）n=0:N/2-1;sblot（3,2,4）;stm（1,16n,.）;bo ontile（）16点IDTX_6（k）;xel（ylabel（x_1_（n）ax（0,32,0,0）N-1;sb（3,2,）；sem（k,bs（32k）,boontie（e） 32点频域采样xlabl（kael（|_（k）|）;axis（,6,0）n1=0:-1;sulo（3,2,）;stm（n1，x32n,box on（f）32点IDT_2（k）xlabe（yabel（x_3_（）;axs（,2,0,2）unctiottem（xn,n） = 0:lengt（
9、xn）1;stem（n,x，.）;lel（yabl（n%xla（ ）;ylael（y）;axis（0,（ed）,min（xn）,.2x（xn）;实验三%实验三程序exp3.%用FFT对信号作频谱分析clara;cse all实验内容（）=x1=ones（1,4）;产生序列向量x（n）=R（n）M=8;xa1:（/2）;xb=（/2）：1:;x2nxa，xb; %产生长度为8的三角波序列x2（）x=xb,xa；1k8=fft（x1n,）; %计算x1n的8点DFTX1k6ft（n,16）; 计算n的点DT28=ft（x2n,8）; 计算x的8点DFTX216ft（x2n，16）; 计算x1n的6
10、点DF3k8fft（x3n，）； %计算x1n的点DFTk6=ft（3n，1）; 计算x1n的点DFT以下绘制幅频特性曲线subpot（2,2，）;mtem（Xk8）;%绘制8点DFT的幅频特性图tle（1a） 点DFTx_1（n）xlabel（/）;ais（,0,12*max（s（X）subpot（2,3）；mstem（1k16）;绘制6点FT的幅频特性图ttle（1b）1点DFTx1（n）;xlael（/axs（,2,0,.*m（ab（X1k1）figure（2）o（2，2，1）;mste（X2k）;%绘制8点FT的幅频特性图titl（2a） 点DTx_2（）；xlbe（ ab（幅度）；a
11、is（,2,1.*max（a（X28）suplo（，2,）;stem（X2k1）;绘制16点DFT的幅频特性图（b）16点DFTx_2（n）xlabel（/ ylabl（幅度ais（0,2，0,12*max（as（X2k6）ubpot（，2，3）;mst（X3k8）;（） 8点F_3（）;labe（ ylabel（axis（，0,2*x（as（X3k8）subplot（2，）;mstem（X3k1）;%绘制16点DT的幅频特性图（3）1点DTx_（n）xlal（ ylbe（幅度）;xi（,0,12*（abs（31）实验内容（2） 周期序列谱分析=N=8;n=0: %F的变换区间N8x=cos（
12、n/4）;xos（i*4）+cos（pi8）;X4k8=fft（x4n）; %计算x4n的8点DTX5fft（）; %计算x5n的8点FTN=;n0: %FT的变换区间6x4=s（*n/4）;x5n=cos（i*n/4）+cos（pi*/8）;Xk16=t（x）； %计算x的1点X56=fft（x5n）; %计算x5的1点DFTfue（3）sublo（2，2,1）;msem（X4k）;%绘制8点DT的幅频特性图itle（4a） 8点DFx_（n） xlal（axis（0，2，0，1.ma（abs（X4k8）subplot（2,2,3）；stm（X4k1）;（4b）16点DTx4（） lbe（/
13、ylal（as（0,2，0,2*ma（abs（X46）plot（2,2）;mstem（X5k）；%绘制8点D的幅频特性图（5）8点Fx_5（n）label（/）; yabel（幅度axi（0,2,0,1.2*max（abs（5k8）subpot（2，2,4）;mtem（516）；（b）16点DFT（） ylael（ais（,2,0,1.2ax（b（X）%实验内容（3） 模拟周期信号谱分析=figur（4）Fs4;T=1Fs；N=16; %FFT的变换区间N=16xno（8*pi*n*）+cos（6*pi*）+cos（20p*n）;%对x6（t）16点采样X66fft（xnT）； %计算xnT的
14、6点DFT616=fftift（X6k16）; %将零频率移到频谱中心 TN*T；F=/Tp； 频率分辨率Fk=-/2：N/2-1；fk=k*; 产生1点DF对应的采样点频率（以零频率为中心）sblot（3,1）；ste（fk,as（X6k6）,box on %绘制8点DFT的幅频特性图tt（a） 16点|DTx_6（nT）|f（Hz）;%ylab（ais（N*/2-,N*F2,0,1.2max（ab（Xk1）N=3;-; %FFT的变换区间N=1xnT=os（8pi*）cs（16pin*）+os（20pi*T）；对x6（t）2点采样63fft（xnT）； %计算xnT的3点DFTX632ff
15、tshif（62）; %将零频率移到频谱中心 Tp=N*;=/p; %频率分辨率FkN2:N/21;fkk*F;产生16点F对应的采样点频率（以零频率为中心）subplot（3, 1, 2）;se（k， abs（X6k2）,oxo%绘制8点DF的幅频特性图（6b） 32点DFTx_（）|）；abe（axs（N*F/21, */-1，, 1.2*ma（ab（6k2）N=64;n=0： F的变换区间N=16x6nTcos（8*n*T）+cos（16*pi*n*）+os（2*pi*n*）;对（t）4点采样X6k64=fft（x6nT）; %计算x6nT的4点DFTX6k4=ffthift（X6k6）
16、; %将零频率移到频谱中心 =N*T；1Tp; %频率分辨率k=-N2：N2-1;f=kF; %产生6点DFT对应的采样点频率（以零频率为中心）sul（3,1,3）；stm（,ab（X664）,bo on绘制点DFT的幅频特性图（a） 64点|DFTx_6（）label（f（Hz）ai（-N*F/2-,N*F/2-1，, .2*max（bs（X6k6）fncto mtem（Xk）%sem（Xk）绘制频域采样序列向量Xk的幅频特性图lgt（Xk）;-1;wk=2k/M;%产生M点DFT对应的采样点频率（关于i归一化值）stem（w,abs（Xk）,%绘制点FT的幅频特性图?wi）；xi（0,，0
17、,1.2mx（ab（Xk）;实验四%实验四程序ep4% I数字滤波器设计及软件实现%lear all; close allF=10000; =1/Fs； %采样频率%调用信号产生函数mstg产生由三路抑制载波调幅信号相加构成的复合信号st tmsg;%低通滤波器设计与实现=f=20; fs=45;wp=2*f/Fs;ws2*fs/Fs; rp=.1; s=6;DF指标（低通滤波器的通、 阻带边界频率）N,wp=ellpord（wp，ws,p,rs）； %调用lipord计算椭圆DF阶数和通带截止频率wpB,Alp（，rp,s,wp）;%调用ellp计算椭圆带通D系统函数系数向量B和y1t=fi
18、ter（B,st）; 滤波器软件实现%低通滤波器设计与实现绘图部分fgure（）; sublot（3，,1）；mylo（,A）;调用绘图函数mlt绘制损耗函数曲线t=y_（t）sbpt（3,1，2）; lt（yt,T,yt）; %调用绘图函数tplo绘制滤波器输出波形%=带通滤波器设计与实现=l=40； fu=60; fsl=75; fu=90;w=2*fpl/Fs,2*fu/Fs；=*flFs,*su/Fs;=0.1； rs60；，wpellipord（p，ws,rp，rs）;调用ellipor计算椭圆D阶数N和通带截止频率w,=elli（N,rp,rs,w）；调用l计算椭圆带通DF系统函数
19、系数向量B和Ay2t=filt（B,A,st）; %滤波器软件实现%带通滤波器设计与实现绘图部分（省略）补充figure（）；sbplot（3，1，3）;pot（,A）；%-gr（3）;yt=y2（t）；sublot（3,，1）； tplot（y2t，T，yt）;%=%高通滤波器设计与实现=fp=890; =0;p=2f/F;s=*fs/F;rp=.1； rs； ,w=liord（wp,ws,rp,rs）;%调用elliord计算椭圆F阶数N和通带 %截止频率wpB,A=lip（N,p,r,wp,gh%调用ei计算椭圆带通DF系统函数系数向量和A3t=fr（B,s）; %滤波器软件实现% 高低通滤波器设计与实现绘图部分（省略）%补充2fe（3）; subplo（，,）;myot（,A）;-fiue（3）;3（t）suplt（3，,3）; tplo（3t,T,t）;%
