實(shí)驗(yàn)報(bào)告模板-實(shí)驗(yàn)一

1、實(shí)驗(yàn)一離散時(shí)間信號(hào)的MATLAB實(shí)現(xiàn)課程名稱：數(shù)字信號(hào)處理實(shí)驗(yàn)日期：2016.4.21班級(jí)：數(shù)媒1403姓名：李欣學(xué)號(hào)：1030514327一、實(shí)驗(yàn)?zāi)康?2) 離散序列的基本運(yùn)算與波形變換的MATLAB實(shí)現(xiàn);(3) 線性時(shí)不變系統(tǒng)求解的MATLAB實(shí)現(xiàn); (4) 序列的傅里葉變換與信號(hào)重構(gòu)MATLAB實(shí)現(xiàn)二、實(shí)驗(yàn)要求按照實(shí)驗(yàn)材料內(nèi)容，利用 MATLAB 編程實(shí)現(xiàn)相關(guān)任務(wù)。三、實(shí)驗(yàn)步驟（以實(shí)驗(yàn)材料第四部分“序列的傅里葉變換與信號(hào)重構(gòu)MATLAB實(shí)現(xiàn)”為例， 給出實(shí)驗(yàn)步驟描述、相關(guān)代碼等。）程序如下：Dt=0.00005;t=-0.005: Dt:0.005;xa=exp(-1000*abs(t)

2、;%模擬信號(hào)Wmax=2*pi*2000; K=500;k=0:1:K;W=k*Wmax/K;%連續(xù)時(shí)間傅立葉變換Xa=xa*exp(-1i*t*W)*Dt; Xa=real(Xa);W=-fliplr(W),W(2:501);Xa=fliplr(Xa),Xa(2:501); figure(1)subplot(2,1,1); plot(t*1000,xa,.);xlabel(t in msec); ylabel(xa(t);gtext(模擬信號(hào));subplot(2,1,2); plot(W/(2*pi*1000),Xa*1000,.); xlabel(Frequence in KHz); y

3、label(Xa(jw)*1000);gtext(連續(xù)時(shí)間傅立葉變換); 程序運(yùn)行結(jié)果：程序如下：Dt=0.00005;t=-0.005:Dt:0.005;%模擬信號(hào)xa=exp(-1000*abs(t);Ts=0.0002;n=-25:1:25;%x=exp(-1000*abs(n*Ts);K=500; k=0:1:K;w=pi*k/K;%X=x*exp(-1i*n*w);X=real(X);w=-fliplr(w),w(2:501);X=fliplr(X),X(2:501);figure subplot(2,1,1); plot(t*1000,xa,.);ylabel(xa1(t););

4、hold on ; hold off ; subplot(2,1,2); plot(w/pi,X,.);ylabel(X1(jw);gtext(離散時(shí)間傅立葉變換); gtext(Ts=0.2msec)程序運(yùn)行結(jié)果：程序如下：Ts1=0.0002; Fs1=1/Ts1; n1=-25:1:25;nTs1=n1*Ts1;%x1=exp(-1000*abs(nTs1);Ts2=0.001; Fs2=1/Ts2;n2=-5:1:5;nTs2=n2*Ts2;x2=exp(-1000*abs(nTs2);Dt=0.00005;t=-0.005:Dt:0.005;%模擬信號(hào)重構(gòu)subplot(2,1,1); plot(t*1000,xa1,.);ylabel(xa1(t);title(x1(n)xa(t); hold on;hold off;subplot(2,1,2);plot(t*1000,xa2,.); ylabel(xa2(t);title(從 x2(n)重構(gòu)模擬信號(hào) xa(t); hold on;stem(n2*Ts2*1000,x2); hold off程序運(yùn)行結(jié)果：四、實(shí)驗(yàn)體會(huì)掌握了MATLAB的基本操作以及學(xué)會(huì)了序列的生成;