数字信号处理课程设计

XXXXX姓名:XXXXXXXXXXXXX1一、课程设计题目代码:W=10;f=1000;n=-10:W-1;t=n/f;Xexp00*abs(t));subplotplottX模拟信号曲线xlabeltsylabelxan)');title%标题模拟信号tfNdtNt(1:N)*dt;wf25;Nf=50;w=linspace(0,wf,Nf);dwwf/(Nf-1);W1=-50:50;Xatexp(-1000*abs(t));F=Xat*exp(-1i*t'*w1)*dt;w[-fliplr(w1),w1(2:Nf)];式叶变换率的频谱Yexp)-1)./(exp(2)-exp(1-1i*W1)-exp(1+1i*W1)+1);FfliplrF,F1(2:Nf)];t=[-fliplr(t),t];plotwFlinewidth);%画傅立叶变换曲线xlabelwpiylabelXajΩ)');title%标题傅立叶变换10.80.60.40.201-0.00500.0050.014321000-2002002(2)Dt00005;tDt5;xaexp*abs(t));%取时间从-0.005s到0.005s这段模拟信号Ts01;Ts2=0.0002;n:25;xexp0*abs(n*Ts1));xexp0*abs(n*Ts2));Kk0:1:K;w=pi*k/K;拟角频率Xxexp-j*n'*w);傅立叶变换Xxexp-j*n'*w);傅立叶变换XrealXX12=real(X2);w[-fliplr(w),w(2:101)];频率范围扩展为从-到+XfliplrX),X11(2:101)];X12=[fliplr(X12),X12(2:101)];xlabelwpiylabelX(jw)');xlabelwpiylabelX(jw)');fs=1000Hz信DTFT3210-1-0.8-0.6-0.4-0.200.20.40.60.81fs=5000Hz信DTFT50-1-0.8-0.6-0.4-0.200.20.40.60.813n);ycospin+cos(0.52*pi*n);n);x=y(1:1:10);titlexnnxlabel'n');ylabel('x(n)');axis[0,10,-2.5,2.5]);wlinspace(0,2*pi,length(x));xwx*exp(-j*[1:length(x)]'*w);magxabs(xw);otwmagxxlabelwylabelxjw');axispi]);bplottxmagxabsx1);xaxis[xminxmaxyminymax])axis[xminxmaxyminymax])4w20-21050xlabelkylabelxk');axis0]);w20-21050axis[xminxmaxyminymax])x(n)(0<=n<=9)01093861475201093861475nDTFT10500316254031625wDFT01638210475901638210475kn);ycospin+cos(0.52*pi*n);n9);xy),zeros(1,90)];subplottitlexnnzeros');xlabel'n');ylabel('x(n)');axis[0,100,-2.5,2.5]);wlinspace(0,2*pi,length(x));xwx*exp(-j*[1:length(x)]'*w);magxabs(xw);otwmagxaxis[xminxmaxyminymax])5jwxlabelwylabelxjw');axispi]);bplotjwtxmagxabsx1);xxlabelkylabelxk');axis10]);axis[xminxmaxyminymax])axis[xminxmaxyminymax])x(n)(0<=n<=9+90zeros)20-2090308060705010200903080607050102040nDTFT10500316254031625wDFT1050506070308020900105060703080209001040kn9);ycospin+cos(0.52*pi*n);n9);x=y(1:1:100);titlexnn;xlabelnylabelxn');6axis5,2.5]);wlinspace(0,2*pi,length(x));xwx*exp(-j*[1:length(x)]'*w);magxabs(xw);otwmagxxlabelwylabelxjw');axispi]);bplottxmagxabsx1);xxlabelkylabelxk');axis54]);axis[xminxmaxyminymax])axis[xminxmaxyminymax])axis[xminxmaxyminymax])x(n)(0<=n<=99)20-2010308060705020103080607050204090nDTFT5000316254031625wDFT500010506070308020105060703080209040k7n27);ycospin+cos(0.52*pi*n);n27);x=y(1:1:128);titlexnn);xlabel'n');ylabel('x(n)');axis[0,128,-2.5,2.5]);wlinspace(0,2*pi,length(x));xwx*exp(-j*[1:length(x)]'*w);magxabs(xw);otwmagxxlabelwylabelxjw');axispi]);bplottxmagxabsx1);xxlabelkylabelxk');axis65]);axis[xminxmaxyminymax])axis[xminxmaxyminymax])axis[xminxmaxyminymax])x(n)(0<=n<=127)20-2020406080100120nDTFT5000123456wDFT500020406080100120k8203020100203020100n9);ycospin+cos(0.52*pi*n);n9);x=y(1:1:50);subplot3,1,1);stem(n1,x);title'x(n)(0<=n<=49)');xlabel'n');ylabel('x(n)');axis[0,50,-2.5,2.5]);wlinspace(0,2*pi,length(x));xwx*exp(-j*[1:length(x)]'*w);magxabs(xw);lotwmagxxlabelwylabelxjw');axispi]);bplottxmagxabsx1);xxlabelkylabelxk');axis0]);axis[xminxmaxyminymax])axis[xminxmaxyminymax])axis[xminxmaxyminymax])x(n)(0<=n<=49)0102520501540010252050154035455nDTFT30201001036254103625wDFT1510200253550151020025355040455k9ws=[0.3,0.75];wp=[0.45,0.65];RpRsplotwpiabsHabelwpielplotwpiangleHabelwpiel巴特沃斯滤带通波器响频曲线频曲线N=7wc=0.42840.6693B0.0003003A=1.0000019005700950IIR幅度和相位频响曲线：1.510.5000.40.610.50.900.40.610.50.90.80.10.70.24200.50.40.60.20.70.10.50.40.60.20.70.10.800.91ws1=0.3*pi;ws2=0.75*pi;wp1=0.45*pi;wp2=0.65*pi;wb(wp1-ws1+ws2-wp2)/2;过渡带宽wp/pi+(ws2-wp2)/(2*pi)];理想带通截止频率Nceil(3.3*2*pi/wb);窗口长度bfirNwc滤波器系数hnbn激响应率响应plotwpiabsH频曲线Rxlabelwpiylabel度');plot(w/pi,angle(H));%画相频曲线Rxlabelwpiylabel位');0.00110.00050.00010.0006-0.0031-0.00080.005102160998-0.001845400.025097970.11260.15914545-0.0119091602-0.0002-0.0008-0.00310.00060.00010.00050.00111.510.5000.10.20.40.5000.10.20.40.50.60.70.80.91