数字信号处理(胡广书例题作业程序)

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)');

xlabel('t=-2*pi~2*pi');

grid on;

5、

% 产生一chirp 信号；

% chirp(T0,F0,T1,F1):

% T0: 信号的开始时间；F0:信号在T0时的瞬时频率，单位为Hz；% T1: 信号的结束时间；F1:信号在T1时的瞬时频率，单位为Hz；%-----------------------------------------------------------------

clear;

t=0:0.001:1;

x=chirp(t,0,1,125);

plot(t,x);

ylabel('x(t)')

xlabel('t')

6、

% 计算两个序列的线性卷积；

%-----------------------------------------------------------------

clear;

N=5;%第一个序列的长度

M=6;%第二个序列的长度

L=N+M-1;