大連理工大學(xué)實(shí)驗(yàn)預(yù)習(xí)報(bào)告信號(hào)實(shí)驗(yàn)

大連理工大學(xué)實(shí)驗(yàn)預(yù)習(xí)報(bào)告學(xué)院（系）：電信專業(yè)：電子信息工程（英語(yǔ)強(qiáng)化）班級(jí)：電英1001班姓名：學(xué)號(hào)：組：___實(shí)驗(yàn)時(shí)間：實(shí)驗(yàn)室：實(shí)驗(yàn)臺(tái)：指導(dǎo)教師簽字：成績(jī)：實(shí)驗(yàn)三信號(hào)抽樣一、實(shí)驗(yàn)?zāi)康膶W(xué)會(huì)運(yùn)用MATLAB完成信號(hào)抽樣及對(duì)抽樣信號(hào)的頻譜進(jìn)行分析；學(xué)會(huì)運(yùn)用MATLAB改變抽樣間隔，觀察抽樣后信號(hào)的頻譜變化；學(xué)會(huì)運(yùn)用MATLAB對(duì)抽樣后的信號(hào)進(jìn)行重建。二、習(xí)題設(shè)有三個(gè)不同頻率的正弦信號(hào)，頻率分別為?，F(xiàn)在用抽樣頻率對(duì)這三個(gè)正弦信號(hào)進(jìn)行抽樣，用MATLAB命令畫(huà)出各抽樣信號(hào)的波形及頻譜，并分析頻率混疊現(xiàn)象。1)clcclearTs=1/3800;dt=0.0001;t1=-0.01:dt:0.01;ft=sin(2*pi*100*t1);subplot(221)plot(t1,ft),gridonxlabel('Time'),ylabel('f(t)')title('f1')N=500;k=-N:N;W=2*pi*k/((2*N+1)*dt);Fw=dt*ft*exp(-j*t1'*W);subplot(222)plot(W,abs(Fw)),gridonxlabel('\omega'),ylabel('F(w)')title('f1頻譜')t2=-0.01:Ts:0.01;fst=sin(2*pi*100*t2);subplot(223)plot(t1,ft,':'),holdonxlabel('Time'),ylabel('fs(t)')title('抽樣信號(hào)')stem(t2,fst),gridonFsw=Ts*fst*exp(-j*t2'*W);subplot(224)plot(W,abs(Fsw)),gridonxlabel('\omega'),ylabel('Fs(w)')title('抽樣信號(hào)頻譜')2)clcclearTs=1/3800;dt=0.0001;t1=-0.01:dt:0.01;ft=sin(2*pi*200*t1);subplot(221)plot(t1,ft),gridonxlabel('Time'),ylabel('f(t)')title('f1')N=500;k=-N:N;W=2*pi*k/((2*N+1)*dt);Fw=dt*ft*exp(-j*t1'*W);subplot(222)plot(W,abs(Fw)),gridonxlabel('\omega'),ylabel('F(w)')title('f1頻譜')t2=-0.01:Ts:0.01;fst=sin(2*pi*200*t2);subplot(223)plot(t1,ft,':'),holdonxlabel('Time'),ylabel('fs(t)')title('抽樣信號(hào)')stem(t2,fst),gridonFsw=Ts*fst*exp(-j*t2'*W);subplot(224)plot(W,abs(Fsw)),gridonxlabel('\omega'),ylabel('Fs(w)')title('抽樣信號(hào)頻譜')3)clcclearTs=1/3800;dt=0.0001;t1=-0.01:dt:0.01;ft=sin(2*pi*3800*t1);subplot(221)plot(t1,ft),gridonxlabel('Time'),ylabel('f(t)')title('f1')N=500;k=-N:N;W=2*pi*k/((2*N+1)*dt);Fw=dt*ft*exp(-j*t1'*W);subplot(222)plot(W,abs(Fw)),gridonxlabel('\omega'),ylabel('F(w)')title('f1頻譜')t2=-0.01:Ts:0.01;fst=sin(2*pi*3800*t2);subplot(223)plot(t1,ft,':'),holdonxlabel('Time'),ylabel('fs(t)')title('抽樣信號(hào)')stem(t2,fst),gridonFsw=Ts*fst*exp(-j*t2'*W);subplot(224)plot(W,abs(Fsw)),gridonxlabel('\omega'),ylabel('Fs(w)')title('抽樣信號(hào)頻譜')結(jié)合抽樣定理，用MATLAB編程實(shí)現(xiàn)信號(hào)經(jīng)沖激脈沖抽樣后得到的抽樣信號(hào)及其頻譜，并利用重構(gòu)信號(hào)。clcclearwm=2;wc=1.2*wm;Ts=0.1;n=-100:100;nTs=n*Ts;fs=sinc(nTs/pi);t=-4:0.1:4;ft=fs*Ts*wc/pi*sinc((wc/pi)*(ones(length(nTs),1)*t-nTs'*ones(1,length(t))));t1=-4:0.1:4;f1=sinc(t/pi);subplot(311)plot(t1,f1,':'),holdonstem(nTs,fs),gridonxlabel('nTs'),ylabel('f(nTs)');title('Ts=0.1fs（t）')axis([-88-1.11.1])holdoffsubplot(312)plot(t,ft),gridonxlabel('t'),ylabel('f(t)');title('fs（t）重建Sa（t）')axis([-88-1.11.1])error=abs(ft-f1);subplot(313)plot(t,error),gridonxlabel('t'),ylabel('error(t)');title('重建信號(hào)與Sa（t）信號(hào)絕對(duì)誤差')axis([-88-1.11.1])【實(shí)驗(yàn)體會(huì)】通過(guò)這次實(shí)