這是我做的數字信號處理實驗的“抽樣”部分
用了一個星期時間獨立完成，保證原創
內含實驗結果圖fig格式
源代碼 m格式doc格式
以及實驗報告
因為不支持此類fig，m格式文件，需要評論私我
是用matlab做的
實驗三 3.3.1抽樣引起的混疊 一實驗內容 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）') 三實驗結果 3.3.2抽樣的頻域視圖 一實驗內容 a.已知 Xa（t）=e-1000|t| 求出并繪制其傅里葉變換xa(jΩ) b.以5000HZ和1000Hz分別對其進行采樣得到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') 三實驗結果
3.3.3從樣本重建信號 一實驗內容 1、 擬合正弦波 【實驗內容】 假設三個樣本符合正弦波，能不能確定其各個參數，w需要什么條件。不能的話，給出理由。選取幾個樣本繪制圖像。 2、 線性與多項式插值 【實驗內容】 a. 使用matlab用直線連接樣本，解釋plot如何繪制圖像。 b. 將三角形沖擊與樣本卷積。證明假設t=-1和t=3的樣本是0，上面結果與線性插值 相同 二實驗程序 %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') %線性插值結果 subplot(211),plot(t,hx) grid title('卷積結果') xlabel('t') subplot(212),plot(xt,y) grid title('線性插值結果') xlabel('xt') 三實驗結果 3.3.3.1 3.3.3.2a
3.3.3.2b 四結果分析 不能根據三點得到正弦信號各個量。
完整的Word格式文檔51黑下載地址：
實驗三.docx
(1.1 MB, 下載次數: 15)
2018-12-11 13:39 上傳
點擊文件名下載附件
實驗報告 下載積分: 黑幣 -5
| 
QQ好友和群
QQ空間
騰訊微博
騰訊朋友
收藏
淘帖
頂
踩
