科技大学数字信号处理课程设计.doc

文档介绍：数字信号处理课程设计
第1章信号的时域分析
连续信号的时域分析
用Matlab产生下列信号并与人工分析结果进行比较:
(1) r (t)= t u (t) -1<t<10
(2) x (t)=1+cos10 t -1<t<1
(3) x (t)=(5e - t - 5e -3 t ) u (t) -1<t<5
(4) x (t)=cos(2πt) cos(20πt) 0<t<5
(5) x (t)=sin(t)/t -10<t<10
(1)t=-1::5;
x=(t>=0);
plot(t,x);
axis([-2,6,-,]);
(2)t=-1::1;
x=1+cos(10*t);
plot(t,x);
ylabel('x(t)');xlabel('t');
(3)
t=0::5;
x=t(t>=0);
plot(t,x);
axis([-2,6,-,]);
t=0::10;
m=(t>=0);
n=5*exp(-t)-5*exp(-3*t);
x=n.*m;
plot(t,x);
(4)
w0=2*pi;
w1=20*pi;
t=0::5;
x=cos(w0*t).*cos(w1*t);
plot(t,x);
ylabel('x(t)');xlabel('t');
(5)
t=-10::10;
m=sin(t);
x=m./t;
plot(t,x);
ylabel('x(t)');xlabel('t');
离散时间序列的时域分析及信号的运算
、作图并与理论值进行比较:
(1)x(n)=2δ(n+n0)
(2) x(n)=() n [sin()+cos()]
n=-4:4;
x=().^n;
y=[sin(*pi*n)+cos(*pi*n)];
m=x.*y;
stem(n,m);
(3)已知LTI离散系统,x(n)=[1 1 1],h(n)=[0 1 2 3],求y(n)
x=[1,1,1];
h=[0,1,2,3];
y=conv(x,h);
subplot(3,1,1);stem([0:length(x)-1],x);
ylabel('x');xlabel('Time index n');
subplot(3,1,2);stem([0:length(h)-1],h);
ylabel('h');xlabel('Time index n');
subplot(3,1,3);stem([0:length(y)-1],y);
ylabel('y=x*h');xlabel('Time index n');
(4)已知x(t)=e –2 tu (t), y(t)=e - tu (t) , 求:x(t) * y(t)
t=-10:10;
u=(t>=0);
m=exp(-2*t);
n=exp(-1*t);
x=m.*u;
y=n.*u;
h=conv(x,y);
subplot(3,1,1);stem([0:length(x)-1],x);
ylabel('x(n))');xlabel('Time index n');
subplot(3,1,2);stem([0:length(y)-1],y);
ylabel('y(n)');xlabel('Time index n');
subplot(3,1,3);stem([0:length(h)-1],h);
ylabel('h(n)=x(n)*y(n)');xlabel('Time index n');
(5)已知信号x(t)=(1+t/2)[u(t+2)-u(t-2)], 求x(t+2),x(t-2),x(-t),x(2t),-x(t)
t=-10:10;
m=(t>=2);
n=(t>=-2);
x=(1+(t./2)).*(n-m);
plot(t+2,x);
t=-10:10;
m=(t>=2);
n=(t>=-2);
x=(1+(t./2)).*(n-m);
plot(t-2,x);
t=-10:10;
m=(t>=2);
n=(t>=-2);
x=(1+(t./2)).*(n-m);
plot(-t,x);
t=-10:10;
m=(t>=2);
n=(t>=-2);
x=(1+(t./2)).*(n-m);
plot(2*t,x);
t=-10:10;
m=(t>=2);
n=(t>=-2