西南交大信号与系统实验报告..docxVIP

信号与系统实验报告XXX201411XXXX微电一班实验一 连续时间信号的采样 一、实验目的: 用不同的采样频率采样并进行比较，验证采样定理。 对不同信号采样验证时域与频域的性质。 二、误差分析:题目采用数值方法估计，是把用一个栅格序列来近似，由于，题中用 之间的有限长度信号来近似表示。此处存在误差。 三、实验结论：严格地说，除了用符号处理工具箱(Symbolics)外，不可能用MATLAB来分析模拟信号。然而如果用时间增量足够小的很密的网格对采样，就可得到一根平滑的曲线和足够长的最大时间来显示所有的模态。这样就可以进行近似分析。令是栅网的间隔且，则 （2）可以用一个数组来仿真一个模拟信号。不要混淆采样周期和栅网间隔，因为后者是MATLAB中严格地用来表示模拟信号的。类似地，付利叶变换关系也可根据（2）近似为： （3）现在，如果（也就是）是有限长度的。则公式（3）与离散付利叶变换关系相似，因而可以用同样的方式以MATLAB来实现，以便分析采样现象。采样频率必须大于2即奈奎斯特采样频率（为有限带宽信号的带宽），该信号才可以由它的采样值重构。 若小于2则会在x(n)中产生混叠，为欠采样。产生的傅立叶变换图形是把互相交叠的的复制品叠加的结果，所以产生了失真。 四、思考题： 1．通过实验说明信号的时域与频域成反比的关系。 从实验一做出的图形看出： 的时域是时域的1/2倍，的频域是的频域的2倍。 2．分别求出奈奎斯特采样间隔，并与例一的信号的奈奎斯特采样间隔比较。 的奈奎斯特采样频率为8000样本每秒；则奈奎斯特采样间隔为 1/8000; 是例一的信号的1/2;的奈奎斯特采样频率为2000样本/秒；则奈奎斯特采样间隔为1/2000, 是例一的信号的2倍。 五. 完整书写你的全部MATLAB程序，并给出程序运行的结果图。1.Dt=0.00005; t=-0.005:Dt:0.005; xa=exp(-1000*abs(t));%连续时间傅立叶变换Wmax=2*pi*2000;K=500;k=0:1:K;W=k*Wmax/K; Xa=xa*exp(-j*t*W)*Dt; Xa=real(Xa);W=[-fliplr(W),W(2:501)]; %频率从-Wmax to WmaxXa=[fliplr(Xa),Xa(2:501)]; %Xa 介于 -Wmax和 Wmax之间% subplot(1,1,1)subplot(2,1,1);plot(t*1000,xa);xlabel(t 毫秒); ylabel(xa(t)); title(模拟信号)subplot(2,1,2);plot(W/(2*pi*1000),Xa*1000);xlabel(频率（单位：KHz)); ylabel(Xa(jW)*1000)title(连续时间傅立叶变换)2.% 模拟信号Dt=0.00005;t=-0.005:Dt:0.005;xa=exp(-1000*abs(t));%离散时间信号Ts=0.0002;n=-25:1:25;x=exp(-1000*abs(n*Ts));%离散时间傅立叶变换K=500;k=0:1:K;w=pi*k/K;X=x*exp(-j*n*w);X=real(X);w=[-fliplr(w),w(2:K+1)];X=[fliplr(X),X(2:K+1)];subplot(1,1,1)subplot(2,1,1);plot(t*1000,xa);xlabel(t 毫秒);ylabel(x1(n));title(离散信号);hold on stem(n*Ts*1000,x);% gtext(Ts=0.2毫秒);hold offsubplot(2,1,2); plot(w/pi,X);xlabel(以pi为单位的频率);ylabel(X1(w));title(离散时间傅立叶变换);clearclc[y,fs,bits]=wavread(P_n8); n=length(y);y_zp=fft(y,n); %对n点进行傅里叶变换到频域f=fs*(0:n/2-1)/n; % 对应点的频率fp=1500;fc=1700;As=100;Ap=1;%通带边缘频率fp,阻带边缘频率fc wc=2*pi*fc/fs; wp=2*pi*fp/fs; %数字频率 = 模拟频率/采样频率wdel=wc-wp; beta=0.112*(As-8.7); N=ceil((As-8)/2.285/wdel); wn= kaiser(N+1,beta); ws=(wp+wc)/2/pi; b=fir1(N,ws,wn); figure(2); freqz(b,1);x=fftfilt(b,y); X=fft(x,n); figur