数字信号处理第一次上机

数字信号处理第一次上机作业1.p131.10代码：dt=0.0001;%连续函数时域步长tf=6;t=0:dt:tf;xa=sqrt(t)+cos(t);T=0.5;%采样间隔n=0:1:tf/T;x=sqrt(n*T)+cos(n*T);deltax=0.5;xq=deltax*round(x/deltax);subplot(1,2,1);plot(t,xa,':'),holdon,gridon;%连续时间信号,绘制点线plot(n*T,x,'o');%离散时间信号,绘制圆圈stem(n*T,xq,'*'),gridon;%数字信号,绘制星号legend('连续时间信号xa','离散时间信号x','数字信号xq')subplot(1,2,2)stairs(n*T,xq),gridonlegend('数字信号采样保持')set(gcf,'color','w')运行结果：2.输出以下三种信号：1）x(t)=sin⁡(1.25πt+0.25π)，t∈[0,10];连续信号。2）对x(t)采样，采样频率f_s=0.5Hz；离散时间信号。3）对x(t)采样，采样频率f_s=2.5Hz；离散时间信号。4）对x(t)采样，采样频率f_s=4Hz；离散时间信号。代码：clearn1=0:2:10;x1=sin(1.25*pi*n1+0.25*pi);subplot(4,1,1)stem(n1,x1);xlabel('n1');ylabel('x1');n2=0:0.4:10;x2=sin(1.25*pi*n2+0.25*pi);subplot(4,1,2)stem(n2,x2);xlabel('n2');ylabel('x2');n3=0:0.25:10;x3=sin(1.25*pi*n3+0.25*pi);subplot(4,1,3)stem(n3,x3);xlabel('n3');ylabel('x3');t=0:0.01:10;x=sin(1.25*pi*t+0.25*pi);subplot(4,1,4)plot(t,x);xlabel('t');ylabel('x');结果3.计算习题2.13（b）中两个离散信号的线性卷积（书面作业题中的一道）。输出卷积前两个序列x_1(n),x_2(n)及卷积后序列x_3(n)=x_1(n)*x_2(n)的图形。注：要求不使用matlab中的卷积函数，自己编制卷积程序，并和卷积函数计算的结果进行比较！注意程序适用性（即：输入任意两个有限长序列，程序都能计算其线性卷积结果）！代码：x=[1111];y=[654321];N1=length(x);N2=length(y);N=N1+N2-1;xl=0:1:N1-1;yl=0:1:N2-1;zl=0:1:N-1;[A]=juanji(x,y);figure(1)subplot(3,1,1)scatter(xl,x,'o');subplot(3,1,2)scatter(yl,y,'o');subplot(3,1,3)scatter(zl,A,'o');function[A]=juanji(xn,hn)N1=length(xn);N2=length(hn);N=N1+N2-1;for(n=1:N)y(n)=0;for(m=1:N2)k=n-m+1;if(k=1&k=N1)y(n)=y(n)+hn(m)*xn(k);endendendA=y;end运行结果