信号与系统实习报告

《信号与系统实习报告》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、实验三3•clearall;dt=0.01;咎采样时间间隔T=2;tao=l;N=100;%T脉冲周期，tao脉宽，N谐波次数omega0=2*pi/T;t=-6:dt:6;d=-6:T:6;fxx=pulstran(t,d,1rectpuls',tao);figure(1);plot(t，fxx);n=-N:N;F=fxx*exp(-j*t1*n*omega0)*dt/T;figure(2);Fl=abs(F);subplot(121)plot(n,Fl);phaF=angle(F);subplot(122)stem(n,phaF,1fi

2、ll1)实验四］对连续信号xa(t)=Ae~asin(Q0/)w(^),A=444.128,a=5oV2^,Qo=50血龙进行理想采样，可得采样序列x(n)=X^nT)=Ae(X，，T05^550。下图给出了连续信号七⑴的幅频特性曲线，由此图可以确定对七⑴采用的采样频率。要求：分别取采样频率几为IKHz、300Hz和200Hz,用matlab编程求采样序列兀⑺)的幅X(eja))频特性'7O观察所得采样序列的幅频特性，说明冇无频谱混叠现象。参考程序代码:%产生序列x(n)n=0:50;A二444.128;a=50*sqrt(2.0)*pi;

3、T=0.001;%T分别取1/1000、1/300、1/200w0=50*sqrt(2.0)*pi;x二A*exp(-a*n*T).*sin(wO*n*T);%函数f的表达式subplot(1,2,1),stem(n,x)titleC理想采样序列fs=1000Hz，)%绘制x(n)的幅度谱k二-250:250;W=pi/125*k;X=x*(exp(-j*pi/125))「(n'*k);%由公式计算DTFTmagX=abs(X);subplot(1,2,2),plot(W/pi,magX)titleC理想釆样序列的幅度谱')2.设x(n)=

4、cos(0.48加)+cos(0.52加)(1)取兀(力(0N,则将x截短为N点序列，再作N点DFT(2)若x的长度M

5、：⑴和⑵的代码：n=0:10;M=length(n);xl=cos(0.48*pi*n)+cos(0.52*pi*n);subplot(2,2,l)stem(n,xl)xlabel('n')title(*x(n)0<=n<=10,)k=0:250;N=length(k);w=2*pi/N*k;WN=exp(-j*2*pi/N);kn=n'*k;WNkn=WN.Akn;X=xl*WNkn;subplot(2,2,2)plot(w/pi,abs(X))xlabelCw/pi1)titleCx(n)ffi里叶变换的近似幅度)k=0:10;N=le

6、ngth(k);Xl=fft(xl);w=2*pi/N*k;subplot(2,2,3)stem(v/pi,abs(Xl))holdonplot(w/pi,abs(Xl),'r:f)xlabelCw/pi*)titleCX(k)的幅度(变换区间长度N=ll)*)k=0:20;N=length(k);X2=fft(xl,N);w=2*pi/N*k;subplot(2,2,4)stem(w/pi,abs(X2))holdonplot(w/pi,abs(X2)/r:r)xlabel('w/pijtitleCX(k)的幅度(变换区间长度N=21)

7、*)x(n)傅里叶变换的近似幅度15x(n)0<=n<=10w/pi(1)和⑵的结果:••/■：0O"p・I1鮎•爲竹勰•f崭I0.511.52w/pi可见，通过补零加长序列长度，使得DFT变换后的谱线更密，但只是观察(计算)分辨率提高了，从图中很难看出信号的频谱部分。⑶的代码：n=0:100;M=length(n);x3=cos(0.48*pi*n)+cos(0.52*pi*n);subplot(2J,l)stem(n,x3)xlabelCn')titleCx(n)0<=n<=100,)k=0:100;N=length(k);X3=fft

8、(x3);w=2*pi/N*k;subplot(2,l,2)plot(w/pi,abs(X3))xlabel('w/pi‘)titleCX(k)的幅度)(3)的结果：x(n)0<