信号与系统仿真实验报告
- 格式:docx
- 大小:580.96 KB
- 文档页数:29
subplot(2,2,3),stem(n,yy,'r'),title('N=16的FFT频谱');
subplot(2,2,2),stem(n,abs(y)) ,title('N=16的DFT幅频特性');
>> axis([k1,k2,0,1.5])
>> title('单位冲激序列')
图像如下:
2.单位阶跃序列 =
程序如下:
>> %lab2.2
>> k1=-5;k2=10;k0=0;
>> k=k1:-k0-1;
>> kk=-k0:k2;
>> n=length(k);
>> nn=length(kk);
>> u=zeros(1,n);
>> k=k1:k2;
>> f=exp(a*k);
>> stem(k,f,'filled')
>> title('负指数序列')
图像如下:
5.单位矩形序列 =
程序如下:
>> %lab2.4
>> k1=-5;k2=25;k0=0;N=20;
>> k=k1:k2;
>> n=length(k);
>> u=zeros(1,n);
>> figure(1)
>> impz(b,a,-10:50)
>> ylabel('h(n)');
>> title('单位冲激响应')
>> figure(2)
>> w=[0:1:500]*pi/500;
>> freqz(b,a,w)
>> title('H(w)频率特性')
单位冲击响应图形如下:
H(w)频率特性图如下:
>> uu=[(k>=k0)&(k<=k0+N-1)];
>> stem(k,u,'filled');
>> hold on
>> stem(k,uu,'filled')t;> axis([k1,k2,0,1]);
>> title('单位矩形序列')
图像如下:
6.取样函数序列
程序如下:
subplot(2,2,2),stem(n,abs(y)) ,title('N=10的DFT幅频特性');
subplot(2,2,4),stem(n,abs(fft(x)),'r'),title('N=10的FFT幅频特性');
N=16;
n=0:N-1;
x=cos((pi./8)*n);
y=dft(x,N);yy=fft(x);
subplot(2,2,3),stem(n,yy,'r'),title('N=22的FFT频谱');
subplot(2,2,2),stem(n,abs(y)) ,title('N=22的DFT幅频特性');
>> %stairs(k,f)
>> plot(k,f)
>> axis([k1,k2,-0.4,1.1])
>> title('取样函数')
取样函数图像如下:
7.正弦函数 )
周期T=20时,程序如下:
>> %lab1.6.1
>> k1=0;k2=60;w=pi/10;
>> k=k1:0.01:k2;
>> f=sin(w*k);
当w=pi/15时,程序如下:
>> %lab2.6.1
>> k1=0;k2=60;w=pi/15;
>> k=k1:k2;
>> f=sin(w*k);
>> stem(k,f,'filed')
>> %axis([k1,k2,-1.1,1.1])
>> title('正弦函数序列')
图像如下:
当w=pi/10时,程序如下:
>> plot(k,f)
>> title('正弦函数')
图像如下:
周期T=30时,程序如下:
>> %lab1.6.2
>> k1=0;k2=60;w=pi/15;
>> k=k1:0.01:k2;
>> f=sin(w*k);
>> plot(k,f)
>> title('正弦函数')
图像如下:
实验二、离散时间信号的表示及可视化
= ; = ; = (分别取a>0及a<0);
= ; ; );(分别取不同的 值);
1.单位冲击序列 =
程序如下:
>> %lab2.1
>> k1=-5;k2=5;k0=0;
>> k=k1:k2;
>> n=length(k);
>> f=zeros(1,n);
>> f(1,-k0-k1+1)=1;
>> stem(k,f,'filled')
function[x,n]=stepseq(n0,ns,nf)
n=[ns:nf];
x=[(n-n0)>=0];
**************************************
function[y,ny]=conv_m(x,nx,h,nh)
ny1=nx(1)+nh(1);
ny2=nx(length(x))+nh(length(h));
>> stairs(k,u)
>> plot([-k0,-k0],[0,1])
>> hold off
>> axis([k1,k2,0,1.1])
>> title('单位阶跃函数')
单位阶跃函数图像如下:
结果分析:由图像可知,在t=0时,函数值跳变到1.
3.正指数函数 = (a>0)
程序如下:
>> %lab1.3.1
单位矩形函数图像如下:
结果分析:上图显示[-5,25]区间的图像。在0<t<10时出现矩形。
6.取样函数
程序如下:
>> %lab1.5
>> k1=-60;k2=60;w=pi/20;
>> k=k1:0.01:k2;
>> %k=k1:k2;
>> f=sin(w*k)./(w*k);
>> %stem(k,f,'filled')
程序如下:
>> %lab2.3.1
>> k1=-5;k2=20;a=0.2;
>> k=k1:k2;
>> f=exp(a*k);
>> stem(k,f,'filled')
>> title('正指数序列')
图像如下:
4.负指数序列 = (a<0)
程序如下:
>> %lab2.3.2
>> k1=-5;k2=20;a=-0.2;
>> k=k1:k2;
>> n=length(k);
>> u=zeros(1,n);
>> uu=[(k>=k0)&(k<=k0+N-1)];
>> stairs(k,u)
>> hold on
>> stairs(k,uu)
>> hold off
>> axis([k1,k2,0,1.2])
>> title('单位矩形函数')
>> uu=ones(1,nn);
>> %f(1,-k0-k1+1)=1;
>> stem(kk,uu,'filled')
>> hold on
>> stem(k,u,'filled')
>> hold off
>> axis([k1,k2,0,1.5])
>> title('单位阶跃序列')
图像如下:
3.正指数序列 = (a>0)
>> k1=-5;k2=10;a=0.2;
>> k=k1:k2;
>> f=exp(a*k);
>> plot(k,f)
>> %axis([k1,k2,0,1.5])
>> title('正指数函数')