电子科学与工程学院 信号与系统 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) );