%时域采样与频域采样实验figure(1)A=444.18;a=50*pi*sqrt(2);omega=50*pi*sqrt(2);%模拟信号Tp=64/1000;Fs=1000;%抽样频率n=0:Tp*Fs-1;y=A*exp(-a*n*1/Fs).*sin(omega*n*1/Fs);subplot(2,2,1);stem(n,y,.);xlabel(抽样频率为1000HZ时信号Y的抽样图形);y1=fft(y,64);%频域抽样用fft计算DFT。subplot(2,2,2);n=0:length(y1)-1;stem(n,y1,.);xlabel(抽样频率为1000HZ时信号y的64点得FFT);subplot(2,2,3);n1=0:length(angle(y)-1;plot(n1,angle(y);xlabel(抽样频率为1000HZ时信号Y的相频特性图);figure(2)A=444.18;a=50*pi*sqrt(2);omega=50*pi*sqrt(2);%模拟信号Tp=64/1000;Fs=300;%抽样频率n=0:Tp*Fs-1;y=A*exp(-a*n*1/Fs).*sin(omega*n*1/Fs);/两向量相乘用点乘y2=y,zeros(1,46);subplot(2,2,1);n2=0:length(y2)-1;stem(n2,y2,.);xlabel(抽样频率为300HZ时信号Y的抽样图形);y1=fft(y2,64);%频域抽样subplot(2,2,2);n3=0:length(y1)-1;stem(n3,y1,.);xlabel(抽样频率为300HZ时信号y的64点得FFT);subplot(2,2,3);n1=0:length(angle(y2)-1;plot(n1,angle(y2);%相位xlabel(抽样频率为300HZ时信号Y的相频特性图);figure(3)A=444.18;a=50*pi*sqrt(2);omega=50*pi*sqrt(2);Tp=64/1000;抽样的实现Fs=200;%抽样频率n=0:Tp*Fs-1;%信号长度y=A*exp(-a*n*1/Fs).*sin(omega*n*1/Fs);/t换成nTy2=y,zeros(1,52);subplot(2,2,1);n2=0:length(y2)-1;stem(n2,y2,.);xlabel(抽样频率为200HZ时信号Y的抽样图形);y1=fft(y2,64);%频域抽样subplot(2,2,2);n3=0:length(y1)-1;stem(n3,y1,.);xlabel(抽样频率为200HZ时信号y的64点得FFT);subplot(2,2,3);n1=0:length(angle(y2)-1;plot(n1,angle(y2);%相位xlabel(抽样频率为200HZ时信号Y的相频特性图);figure(4)n1=0:13;n2=14:26;x1=n1+1;x2=27-n2;x0=x1,x2,zeros(1,5);%信号xnx00=x1,x2;xk=fft(x0,32);%信号xn的32点频域抽样n3=0:length(xk)-1;subplot(2,2,1);stem(n3,abs(xk),.);%信号幅频特性xlabel(信号XN的32点FFT);xk1=fft(x00,16);%16点得频域抽样n4=0:length(xk1)-1;subplot(2,2,2);stem(n4,abs(xk1),.);xlabel(信号XN的16点FFT);n5=0:length(x0)-1;subplot(2,2,3)stem(n5,x0,.);%原始信号xnxlabel(信号X0的波形);figure(5)y5=ifft(xk,32);n5=0:length(xk)-1;subplot(2,2,1);stem(n5,y5,.);xlabel(32点IFFT);y51=ifft(xk1,16);n51=0:length(y51)-1;subplot(2,2,2);stem(n51,y51,.);xlabel(16点IFFT);%实验二figure(6)x1n=ones(1,4);%信号X1Nn1=0:3;n2=4:7;x1=n1+1;x2=8-n2;x2n=x1,x2;%信号X2Nn3=0:3;n4=4:7;x3=4-n;x4=n4-3;x3n=x3,x4;%信号X3N即分段信号的表示方法，用向量表示的。y1n=fft(x1n,8);%信号X1N的8点ffty11n=fft(x1n,16);/离散傅立叶变换求频谱即求信号的频谱。你懂的注意点数要满足的条件,不然就有失真啦就不接近DTFT的波形啦。subplot(2,2,1);n1n=0:length(y1n)-1;/长度要匹配，才能显示出来！这个必须懂。stem(n1n,abs(y1n),.);%幅频图xlabel(信号X1N的8点幅频图);subplot(2,2,2);n11n=0:length(y11n)-1;stem(n11n,abs(y11n),.);xlabel(信号x1n的16点幅频图);y2n=fft(x2n,8);%信号x2n的频域8点抽样subplot(2,2,3);n2n=0:length(y2n)-1;stem(n2n,abs(y2n),.);xlabel(信号x2n的8点幅频图);y2n=fft(x2n,16);%信号x2n的频域16点抽样subplot(2,2,4);n2n=0:length(y2n)-1;stem(n2n,abs(y2n),.);%信号x2n的16点幅频图xlabel(信号x2n的16点幅频图);figure(7)y3n=fft(x3n,8);subplot(2,2,1);n3n=0:length(y3n)-1;stem(n3n,abs(y3n),.);%信号x3n的8点幅频图xlabel(信号x3n的8点幅频图);y3n=fft(x3n,16);subplot(2,2,2);n3n=0:length(y3n)-1;stem(n3n,abs(y3n),.);%信号x3n的16点幅频图xlabel(信号x3n的16点幅频图);figure(8)n1=0:7;n2=0:15;y1=cos(1/4*pi*n1);y2=cos(1/4*pi*n2);y11=fft(y1);%信号x4n的8点FFTy22=fft(y2);%信号x4n的16点FFTn11=0:length(y11)-1;n22=0:length(y22)-1;subplot(2,2,1);stem(n11,abs(y11),.);xlabel(信号x4n的8点幅频图);subplot(2,2,2);stem(n22,abs(y22),.);xlabel(信号x4n的16点得幅频图);n3=0:7;n4=0:15;x11=cos(1/4*pi*n3)+cos(1/8*pi*n3);%信号x5n长度为8x22=cos(1/4*pi*n4)+cos(1/8*pi*n4);%信号x5n长度为16y1=fft(x11);y2=fft(x22);subplot(2,2,3);n33=0:length(y1)-1;n44=0:length(y2)-1;stem(n33,abs(y1),.);xlabel(信号x5n的8点幅频图);subplot(2,2,4);stem(n44,abs(y2),.);xlabel(信号x5n的16点幅频图);figure(9)Fs=64;T=1/Fs;n1=0:15;t=n1*T;x1=cos(8*pi*t)+cos(16*pi*t)+cos(20*pi*t);%模拟信号的16点抽样y1=fft(x1);n11=0:length(y1)-1;subplot(2,2,1);stem(n11,abs(y1),.);xlabel(抽样信号的16点幅频图);Fs=64;T=1/Fs;n2=0:31;t1=n2*T;x2=cos(8*pi*t1)+cos(16*pi*t1)+cos(20*pi*t1);%模拟信号的32点抽样y2=fft(x2);n22=0:length(y2)-1;subplot(2,2,2);stem(n22,abs(y2),.);xlabel(抽样信号的32点幅频图);Fs=64;T=1/Fs;n3=0:63;t2=n3*T;x3=cos(8*pi*t2)+cos(16*pi*t2)+cos(20*pi*t2);%模拟信号的64点抽样y3=fft(x3,64);n33=0:length(y3)-1;subplot(2,2,3);stem(n33,abs(y3),.);xlabel(抽样信号的64点幅频图);滤波器的设计T=1;Fs=1/T;wp=0.25*pi,0.45*pi;%转换成模拟角频率ws=0.15*pi,0.55*pi;rp=3;rs=40;n,wc=buttord(wp,ws,rp,rs,s);b,a=butter(n,wc,s);bz,az=impinvar(b,a,Fs);%脉冲响应法h,w=freqz(bz,az,whole);/数字系统的频谱求法用freqz(),模拟系统的频谱用freqs();在数字频谱中，whole将显示0到2*pi内的频谱图，没写则是0到pi。实验频谱图可以说明。subplot(2,2,1);plot(w/pi,abs(h);xlabel(w/pi);ylabel(|H(w)|);T=1;Fs=1/T;wpz=0.25,0.45;wsz=0.15,0.55;wp=2*tan(wpz*pi/2);%双线性法中的预畸变ws=2*tan(wsz*pi/2);rp=3;rs=40;n,wc=buttord(wp,ws,rp,rs,s);/间接法，里面的角频率都是模拟量b,a=butter(n,wc,s);bz,az=bilinear(b,a,Fs);%双线性法h,w=freqz(bz,az,whole);subplot(2,2,2);plot(w/pi,abs(h);xlabel(w/pi);ylabel(|H(w)|);