nyquist.m

来自「Matlab基础程序联系：成绩统计,练习switch和if else组合；抽样定」· M 代码 · 共 22 行

function Nyquist(ori_freq,sample_rate)
sample_signal=zeros(1,64);
ori_signal= 100*sin(2*pi*1000*[0:.00001:0.001]);
for k=1:64     
    sample_signal(k) = 100*sin(2*pi*k*ori_freq/sample_rate); 
end
figure,
title('time');
subplot(211)
plot(ori_signal );
subplot(212)
plot(sample_signal );
f_ori = abs(fft(ori_signal));
f_sample = abs(fft(sample_signal));
figure,
title('frequency');
subplot(211)
plot(f_ori );
subplot(212)
plot(f_sample );