fft1.m
function
result=fft1
(w,hanshu,n)
fori=1:length(w)
m=hanshu.*((exp(-1
i*(i-1)*pi/100)).^n);
a=sum(m);
endfor
i=1:length(w)
result(i)=a;
end
x1=0:pi/10:(8*pi);
w=linspace(0,8*pi,length(x1));
figure
subplot(211)
plot(x1,sin(x1)+2*cos(x1),'k'); %原時域連續訊號y=sin(t)+2cost
line([0
30],[0
0]);
xlabel('t');ylabel('x(t)');
title('原時域連續訊號y=sin(t)+2cost');
grid
sin1=sin(x1)+2*cos(x1);
n=0:(length(x1)-1);
subplot(212)
plot(w,abs(fft1(w,sin1,n))); %其對應頻域訊號y=fft（sin(t)）
xlabel('w');ylabel('x(w)');
title('其對應頻域訊號y=ft（sin(t)+2cos(t)）');
grid
n1=input('請輸入取樣點數n:');
n=0:n1;
zb=size(n);
figure
%sinf=sin(8
*pi*n/zb(2));
sinf=sin(8
*pi*n/zb(2))+2
*cos(8
*pi*n/zb(2));
subplot(211);
stem(n,sinf,'.');
xlabel('n');ylabel('x(n)');
title('取樣後的時域訊號y=x(n)');
w=0:(pi/100):4
*pi;
subplot(212)
plot(w,abs(fft1(w,sinf,n)));
xlabel('w');ylabel('x(w)');
title('取樣後的頻域訊號y=ft（sin(n)）');
wp=0.1;ws=0.7;rp=2;as=30;
[n,wc]=buttord(wp,ws,rp,as);
[b,a]=butter(n,wc);
[h,w]=freqz(b,a,512,1000); %返回在取樣頻率下頻率向量
figure;
plot(w*1000/(2
*pi),abs(h));
xlabel('hz');ylabel('頻率響應幅度');
title('低通濾波器');
grid;
figure
y=filter(b,a,sinf);
subplot(2,1,1);plot(y,'k')%恢復後的連續訊號
xlabel('t');ylabel('x(t)');
title('恢復後的連續訊號y=sin(t)+2cost');
grid;
y=fft(y,512);
w=(0:255)/256
*500;
subplot(2,1,2),
plot(w,abs([y(1:256)]));%繪製頻譜圖
xlabel('hz');ylabel('頻率響應幅度');
title('頻譜圖');
grid;
在matlab中理解取樣定理
奈奎斯特取樣定理 眾所周知一句話，取樣頻率大於等於訊號最高頻率的兩倍才能使採得的訊號不失真，也就輸專業術語，在頻域發生混疊，不便於後續的操作。a 1 f 2 f sample 4 t 0 1 f sample 1 y a sin 2 pi f t subplot 3,2,1 plot t,y tit...
關於取樣定理
取樣定理，又稱夏農取樣定理，奈奎斯特取樣定理，是 資訊理論，特別是通訊與 訊號處理 學科中的乙個重要基本結論 e.t.whittaker 1915 年發表的統計理論 克勞德 夏農 與harry nyquist 都對它作出了重要貢獻。另外 v.a.kotelnikov 也對這個定理做了重要貢獻。取樣是...
夏農取樣定理
編輯夏農取樣定理，又稱奈奎斯特取樣定理，是資訊理論，特別是通訊與訊號處理學科中的乙個重要基本結論。1924年奈奎斯特 nyquist 就推導出在理想低通訊道的最高大碼元傳輸速率的公式 理想低通訊道的最高大碼元傳輸速率 2w log 2n 其中w是理想低通訊道的頻寬，n是電平強度 中文名 夏農取樣定理...