数字信号处理胡广书例题作业程序.doc

1、 %---filter求卷积，B(Z)/A(Z)=H(Z),已知B(Z)和A(Z)，求y(n)=x(n)*h(n)----- clear; x=ones(100); t=1:100; b=[.001836,.007344,.011016,.007374,.001836]; a=[1,-3.0544,3.8291,-2.2925,.55075]; % y=filter(b,a,x); % 求所给系统的输出，本例实际上是求所给系统的阶跃响应； plot(t,x,r.,t,y,k-);grid on; ylabel(x(n) and y(n)) xlabel(n) 1、 %---filter求卷积，B(Z)/A(Z)=H(Z),已知B(Z)和A(Z)，求y(n)=x(n)*h(n)----- clear; x=ones(100); t=1:100; b=[.001836,.007344,.011016,.007374,.001836]; a=[1,-3.0544,3.8291,-2.2925,.55075]; % y=filter(b,a,x); % 求所给系统的输出，本例实际上是求所给系统的阶跃响应； plot(t,x,r.,t,y,k-);grid on; ylabel(x(n) and y(n)) xlabel(n) 第一章 产生信号，求卷积和自相关函数 1、 %信号产生 n=0:100; %工频 f0=50;A=220;fs=400; x1=A*sin(2*pi*f0*n/fs); subplot(321);plot(n,x1);xlabel(n);ylabel(x1(n)) ;grid on; %率减正弦 f0=2;A=2;alf=0.5;fs=16; x2=A*exp(-alf*n/fs).*sin(2*pi*f0*n/fs); subplot(323);plot(n,x2);xlabel(n);ylabel(x2(n)) ;grid on; %谐波信号 f0=5;A1=1.0;A2=0.5;A3=0.2;fs=100; x3=A1*sin(2*pi*f0*n/fs)+A2*sin(2*pi*2*f0*n/fs)+A3*sin(2*pi*3*f0*n/fs); subplot(322);plot(n,x3);xlabel(n);ylabel(x3(n)) ;grid on; %哈明窗 f0=10;fs=1000; x4=0.54-0.46*cos(2*pi*f0*n/fs); subplot(324);plot(n,x4);xlabel(n);ylabel(x4(n)) ;grid on; %采样 n=-50:50; f0=10;fs=400; w=2*pi*f0*n/fs; x5=sinc(w); subplot(325);plot(n,x5);xlabel(n);ylabel(x5(n)) ;grid on; 2、 % 产生均匀分布的白噪信号，使均值为0，功率为p %----------------------------------------------------------------- clear; p=0.01; N=50000; u=rand(1,N); u=u-mean(u); a=sqrt(12*p); u1=u*a; power_u1=dot(u1,u1)/N subplot(211) plot(u1(1:200));grid on; ylabel(u(n)) xlabel(n) 3、 % 产生高斯分布的白噪信号，使功率为p，并观察数据分布的直方图 %----------------------------------------------------------------- clear; p=0.1; N=500000; u=randn(1,N); a=sqrt(p); u=u*a; power_u=var(u); subplot(211) plot(u(1:200));grid on; ylabel(u(n)); xlabel(n) subplot(212) hist(u,50);grid on; ylabel(histogram of u(n)); 4、 % 产生一 sinc 函数； %----------------------------------------------------------------- clear; n=200; stept=4*pi/n; t=-2*pi:stept:2*pi; y=sinc(t); plot(t,y,t,zeros(size(t))); ylabel(sinc(t)); x