内科大通信数字信号处理课设.docxVIP

数字信号处理课程设计学院：信息工程学院专业：通信工程班级：12-2班姓名：*** 学号：**********实验一 信号的时域分1.1 连续信号的时域分析一、设计题目用Matlab产生下列信号并与人工分析结果进行比较：(1) x(t)=(5e - t - 5e -3 t )u (t) -1t5(2) x(t)=cos(2πt)cos(20πt) 0t5(3) x(t)=sin(t)/t-10t10二、代码分析及图像（1）代码分析t=-1:0.001:5;%定义域取值m=(t=0); %定义了一个阶跃函数，t=0时m=1，否则为0n=5*exp(-t)-5*exp(-3*t);%定义了一个指数复合函数x=n.*m; %将阶跃函数与指数函数相乘plot(t,x);%输入t值，输出x值（2）代码分析：w0=2*pi;%计算w0值w1=20*pi;%计算w1值t=0:0.001:5;%定义域取值x=cos(w0*t).*cos(w1*t);%将两个cos函数相乘plot(t,x);%输入t值，输出x值 ylabel(x(t));xlabel(t);%将x、y轴定为t和x(t)轴（3）代码分析：t=-10:0.01:10;%定义域取值m=sin(t);%定义了一个sin函数x=m./t;%将两函数相除，形成一个Sa（）函数plot(t,x);%输入t值，输出x值ylabel(x(t));xlabel(t);%将x、y轴定为t和x(t)轴。1.2 离散时间序列的时域分析及信号的运算一、设计题目使用Matlab产生下列序列、作图并与理论值进行比较：（1） x(n)=(0.9) n[sin(0.25πn)+cos(0.25πn)]（2）已知LTI离散系统，x(n)=[1 1 1],h(n)=[0 1 2 3],求y(n)二、代码分析及图像(1)代码分析n=-4:4; %定义域取值x=(0.9).^n;%定义了一个幂函数m=[sin(0.25*pi*n)+cos(0.25*pi*n)];%定义了一个复合正余弦函数y=x.*m;%两函数相乘为一个函数stem(n,y); %输入n值，输出y值ylabel(y);xlabel(n); %将x、y轴定为n和y轴。(2)代码分析x=[1,1,1];%x序列h=[0,1,2,3];%h序列y=conv(x,h);%定义了y 为x和h的卷积，计算离散卷积和subplot(3,1,1);stem([0:length(x)-1],x);%产生三行一列的图形，%并且对第一个图形进行处理ylabel(x);xlabel(Time index n);%将x、y轴定为%Timeindex n和x轴。subplot(3,1,2);stem([0:length(h)-1],h);%对第二个图形进行%理ylabel(h);xlabel(Time index n); %将x、y轴定为%Timeindex n和h轴。subplot(3,1,3);stem([0:length(y)-1],y); %对第三个图形进行%理ylabel(y=x*h);xlabel(Time index n)；%将x、y轴定为%Timeindex n和y=x*h轴。实验二 信号的频域分析一、设计题目利用MATLAB实现由计算有限序列线性卷积，并和线性卷积直接计算的结果进行比较。二、代码分析及图像x=[1 2 0 1];%x的序列值h=[2 2 1 1];%h的序列值L=length(x)+length(h)-1;%其中length(x)表示向量x坐标分量的%个数，即向量的维数XE=fft(x,L);%是补零或截短的n点傅立叶变换，当x(n)的长度小%于L时，在x(n)的尾部补零使x(n)的长度达到n点，当x(n)的长%度大于n时，将x(n)截短使x(n)的长度成n 点，然后对补或者截%短的数据进行快速的傅立叶变换HE=fft(h,L);%进行傅立叶变换y1=ifft(XE.*HE);%进行傅立叶逆变换k=0:L-1;%在0到L-1之间取值subplot(2,1,1);%生成两行一列的图形stem(k,real(y1));%求y1复数的实部，整体画出二维离散数据的%杆图axis([0 6 0 7]);%根据给出的X的最大值最小值，Y的最大值最小%值选择坐标系title(Result of Linear Convolution);%图形名称xlabel(Time index k);ylabel(Amplitude);%横轴为Time %index k纵轴为Amplitudey2=conv(x,h);%计算卷积和error=y1-y2;%误差subplot(2,1,2);%产生两行