资源描述：
《数字信号处理(胡广书例题作业程序).doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、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-'>。gridon。ylabel('x(n>andy(n
2、>'>xlabel('n'>1、%---filter求卷积，B(Z>/A(Z>=H(Z>,已知B(Z>和A(Z>，求y(n>=x(n>*h(n>-----b5E2RGbCAPclear。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-'>。g
3、ridon。ylabel('x(n>andy(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>'>。gridon。p1EanqFDPw%率减正弦f0=2。A=2。alf=0.5。fs=16。x2=A*exp(-alf*n/fs>.*sin(2*pi*f0*n/fs>。subplot(323
4、>。plot(n,x2>。xlabel('n'>。ylabel('x2(n>'>。gridon。DXDiTa9E3d%谐波信号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>。RTCrpUDGiTsubplot(322>。plot(n,x3>。xlabel('n'>。ylabel('x3(n>'>。gridon。5PCzVD7HxA%哈明窗f0=10。fs=1000
5、。45/45x4=0.54-0.46*cos(2*pi*f0*n/fs>。subplot(324>。plot(n,x4>。xlabel('n'>。ylabel('x4(n>'>。gridon。jLBHrnAILg%采样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>'>。gridon。xHAQX74J0X2、%产生均匀分布的白噪信号，使均值为0，功率为p%----------
6、-------------------------------------------------------LDAYtRyKfEclear。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>/Nsubplot(211>plot(u1(1:200>>。gridon。ylabel('u(n>'>xlabel('n'>3、%产生高斯分布的白噪信号，使功率为p，并观察数据分布的直方图%------------------
7、-----------------------------------------------Zzz6ZB2Ltkclear。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>>。gridon。ylabel('u(n>'>。xlabel('n'>subplot(212>hist(u,50>。gridon。ylabel('histogramofu(n>'>。4、%产生一sinc函数；%----------
8、-------------------------------------------------------dvzfvkwMI1clear。45/45n=200。stept=4*pi/n。t=-2*pi:stept:2*pi。y=sinc(t>。plot(t,y,t,zeros(size(t>>>。ylabel('sinc(t>
