发布时间 : 星期一 文章西工大信号系统实验八 连续信号的采样与重构更新完毕开始阅读c0b32595f18583d04864590c
2,设计性实验
设计一模拟信号：x（t）=3sin（2π*f*t）采样频率fs为5120hz，取信号频率f=150hz（正常采样）和f=3000hz（欠采样）两种情况进行采样分析。 代码： clf;
t=0:0.0000005:0.02; f1=150; f2=3000;
xa1=3*sin(2*pi*f1*t); xa2=3*sin(2*pi*f2*t); fs=5120; T=1/fs;
nn1=-1:T:1; nn2=-1:T:1;
xs1=3*sin(2*pi*f1*nn1); xs2=3*sin(2*pi*f2*nn2); k1=0:length(nn1)-1; k2=0:length(nn2)-1; subplot(3,2,1) plot(t,xa1);grid;
xlabel('时间，msec');ylabel('幅值'); title('连续时间信号 x_{a1}(t)'); axis([0 0.02 -3 3]);
subplot(3,2,2); plot(t,xa2);grid;
xlabel('时间，msec');ylabel('幅值'); title('连续时间信号 x_{a2}(t)'); axis([0 0.001 -3 3]); subplot(3,2,3); stem(k1,xs1);grid;
xlabel('时间，msec');ylabel('幅值'); title('离散时间信号 x1[n]'); axis([0 100 -3 3]); subplot(3,2,4); stem(k2,xs2);grid;
xlabel('时间，msec');ylabel('幅值'); title('离散时间信号 x2[n]'); axis([0 5 -3 3]); subplot(3,2,5) N1=length(xs1); fx1=fft(xs1); df1=fs/N1; n1=0:N1/2; f1=n1*df1;
plot(f1,abs(fx1(n1+1))*2/N1);grid;
set(gca,'XTickMode','manual','XTick',[0,50,100,150,200,250,300]); set(gca,'YTickMode','manual','YTick',[1,2,3,4]); title('离散时间信号 x1[n]和频谱图'); axis([0 300 0 4]) subplot(3,2,6) N2=length(xs2); fx2=fft(xs2); df2=fs/N2; n2=0:N2/2; f2=n2*df2;
plot(f2,abs(fx2(n2+1))*2/N2);grid;
set(gca,'XTickMode','manual','XTick',[100, 500,1000,1500,2120,2500])
set(gca,'YTickMode','manual','YTick',[1,2,3,4]) title('离散时间信号 x2[n]和频谱图'); 结果：