随机频率的连续信号matlab,东南大学 信号系统 MATLAB 实验十一 取样定理

电子科学与工程学院 信号与系统 MATLAB 实验

一．

主程序：

clear;

k=0.2;

fn=10;

fs=fn*k;

T=1/fs;

Nc=100/T;

for n=-Nc:1:Nc;

nT(n+Nc+1)=T*n;

f(n+Nc+1)=(sin(10*pi*nT(n+Nc+1))+eps)/(10*pi*nT(n+Nc+1)+eps); end

step=0.005;

t=-0.5:step:0.5;

for N=1:length(t);

for n=1:length(nT);

ff(n)=f(n)*( sin(10*pi*fs*(t(N)-nT(n))+eps)/(10*pi*fs*(t(N)-nT(n))+eps) ); end

fa(N)=sum(ff);

end

y=(sin(10*pi*t)+eps)./(10*pi*t+eps);

subplot(3,1,1);

plot(t,fa,t,y);grid;

hold on;

plot(nT,f,'*');

axis([-0.5,0.5,-0.4,1]);

error=max(abs(fa-y))

hold off

k=0.5;

fn=10;

fs=fn*k;

T=1/fs;

Nc=100/T;

for n=-Nc:1:Nc;

nT(n+Nc+1)=T*n;

f(n+Nc+1)=(sin(10*pi*nT(n+Nc+1))+eps)/(10*pi*nT(n+Nc+1)+eps); end

step=0.005;

t=-0.5:step:0.5;

for N=1:length(t);

for n=1:length(nT);

ff(n)=f(n)*( sin(10*pi*fs*(t(N)-nT(n))+eps)/(10*pi*fs*(t(N)-nT(n))+eps) );