這是我做的數(shù)字信號處理實驗的“抽樣”部分
用了一個星期時間獨立完成，保證原創(chuàng)
內(nèi)含實驗結(jié)果圖fig格式
源代碼 m格式doc格式
以及實驗報告
因為不支持此類fig，m格式文件，需要評論私我
是用matlab做的
實驗三 3.3.1抽樣引起的混疊 一實驗內(nèi)容 f0=5000Hz,fs分別取100Hz,1kHz,10kHz,繪出x[n]及其DTFT 二實驗程序 n=0:20; N=1000; f0=500; fai=pi/4; fs1=100;fs2=1000;fs3=10000; x1=sin(2*pi*f0*n/fs1+fai); x2=sin(2*pi*f0*n/fs2+fai); x3=sin(2*pi*f0*n/fs3+fai); [X1,W1]=dtft(x1,N); [X2,W2]=dtft(x2,N); [X3,W3]=dtft(x3,N); subplot(331),stem(n,x1) grid,title('fs=100Hz的x[n]圖形') xlabel('n'),ylabel('x1[n]') subplot(332),plot(W1,abs(X1)) grid,title('fs=100Hz的x[n]的DTFT幅頻特性') xlabel('W1'),ylabel('|X1|') subplot(333),plot(W1,angle(X1)) grid,title('fs=100Hz的x[n]的DTFT相頻特性') xlabel('W1'),ylabel('fai(X1)') subplot(334),stem(n,x2) grid,title('fs=1000Hz的x[n]圖形') xlabel('n'),ylabel('x2[n]') subplot(335),plot(W2,abs(X2)) grid,title('fs=1000Hz的x[n]的DTFT幅頻特性') xlabel('W2'),ylabel('|X2|') subplot(336),plot(W2,angle(X2)) grid,title('fs=1000Hz的x[n]的DTFT相頻特性') xlabel('W2'),ylabel('fai(X2)') subplot(337),stem(n,x3) grid,title('fs=10000Hz的x[n]圖形') xlabel('n'),ylabel('x3[n]') subplot(338),plot(W3,abs(X3)) grid,title('fs=100Hz的x[n]的DTFT幅頻特性') xlabel('W3'),ylabel('|X3|') subplot(339),plot(W3,angle(X3)) grid,title('fs=100Hz的x[n]的DTFT相頻特性') xlabel('W3'),ylabel('fai（X3）') 三實驗結(jié)果 3.3.2抽樣的頻域視圖 一實驗內(nèi)容 a.已知 Xa（t）=e-1000|t| 求出并繪制其傅里葉變換xa(jΩ) b.以5000HZ和1000Hz分別對其進(jìn)行采樣得到x1(n),x2(n)；畫出他們的DTFT并與xa(jΩ)。 二實驗程序 %a dt=0.00004; t=-0.04:0.00004:0.04; xa=exp(-1000*abs(t)); fmagplot(xa,dt)
%b n=0:20; N=1000; t1=n/5000; t2=n/1000; x1=exp(-1000*abs(t1)); x2=exp(-1000*abs(t2)); [X1,W1]=dtft(x1,N); [X2,W2]=dtft(x2,N); subplot(211),plot(W1,real(X1)) grid,title('5000Hz采樣') xlabel('w1') subplot(212),plot(W2,real(X2)) grid,title('5000Hz采樣') xlabel('w2') 三實驗結(jié)果
3.3.3從樣本重建信號 一實驗內(nèi)容 1、 擬合正弦波 【實驗內(nèi)容】 假設(shè)三個樣本符合正弦波，能不能確定其各個參數(shù)，w需要什么條件。不能的話，給出理由。選取幾個樣本繪制圖像。 2、 線性與多項式插值 【實驗內(nèi)容】 a. 使用matlab用直線連接樣本，解釋plot如何繪制圖像。 b. 將三角形沖擊與樣本卷積。證明假設(shè)t=-1和t=3的樣本是0，上面結(jié)果與線性插值 相同 二實驗程序 %3.3.3.1 t=[0:0.001:6] w1=pi/3; w2=5*pi/3; x1=2*cos(w1*n*dt); x2=2*cos(w2*n*dt); subplot(211),plot(t,x1) grid on,set(gca,'Xtick',t) title('w1=pi/3'); xlabel('t') subplot(212),plot(t,x2) grid on,set(gca,'xtick',t) title('w2=5*pi/3'); xlabel('t')
%3.3.3.2a x=[0,2]; y=[1,1];z=[2,-1] plot([x(1) y(1) z(1)],[x(2) y(2) z(2)]); grid,set(gca,'xtick',0:0.01:2) %3.3.3.2b
h=[0.2 0.4 0.6 0.8 1.0 0.8 0.6 0.4 0.2]; t=-4:14; x=[2 0 0 0 0 1 0 0 0 0 -1]; hx=conv(h,x); %卷積 n=[0 2 1 -1 0]; %插值點 x=-1:3; xt=-1:0.1:3; %被插值點 y=interp1(x,n,xt,'linear') %線性插值結(jié)果 subplot(211),plot(t,hx) grid title('卷積結(jié)果') xlabel('t') subplot(212),plot(xt,y) grid title('線性插值結(jié)果') xlabel('xt') 三實驗結(jié)果 3.3.3.1 3.3.3.2a
3.3.3.2b 四結(jié)果分析 不能根據(jù)三點得到正弦信號各個量。
完整的Word格式文檔51黑下載地址：
實驗三.docx
(1.1 MB, 下載次數(shù): 15)
2018-12-11 13:39 上傳
點擊文件名下載附件
實驗報告
| 
