matlab数字信号处理作业.doc

Matlab数字信号处理实验第1次作业 院系：理学院2010级电子信息科学与技术 学号：2010142112 姓名：李家宁 阅读例子程序，观察输出波形，理解每条语句的含义。 程序1： clear all;close all;clc;%清除原所有变量，关闭窗口，对命令窗口请屏 n=0:9;x=(0.8).^n;%设定函数 X=fft(x,4096);%将信号做FFT运算 Xk1=DFT(n,x,10);Xk2=DFT(n,x,20);%分别将学x(n)做 DTT运算 subplot(211),plot((0:4095)/4095*10,abs(X));%画图排列 hold on;%前轴及图形保持而不被刷新stem(0:9,abs(Xk1),LineWidth,2);%画火柴柱形 legend(x(n)的幅频特性,x(n)的10点DFT的幅度);%画Xk1的频谱、幅度谱 xlabel(k);axis([-1 11 -1 5]);grid;%横坐标的位置 subplot(212),plot((0:4095)/4095*20,abs(X));hold on; 画图排列 stem(0:19,abs(Xk2),LineWidth,2); %画火柴柱形 legend(x(n)的幅频特性,x(n)的20点DFT的幅度); %画Xk1的频谱、幅度谱 xlabel(k);axis([-1 21 -1 5]);grid; %横坐标的位置 set(gcf,color,w);%颜色设定 程序2： xn=[0,1,2,3,4,5,6,7];%设定n的取值 N=length(xn);%将N设定为n取值总长度 n=0:N-1;k=0:N-1;%设定n的变化范围 Xk=xn*exp(-j*2*pi/N).^(n*k); %离散傅里叶变换 x=(Xk*exp(j*2*pi/N).^(n*k))/N; %离散傅里叶逆变换 subplot(2,2,1),stem(n,xn,k);%画x(n)的频谱图形 title(x(n));axis([-1,N,1.1*min(xn),1.1*max(xn)]);%设定标题和横坐标的取值 subplot(2,2,2),stem(n,abs(x),k); %显示逆变换结果 title(IDFT|X(k)|);%将幅度谱命名 axis([-1,N,1.1*min(x),1.1*max(x)]);%设定横坐标的范围 subplot(2,2,3),stem(k,abs(Xk),k); %显示序列的幅度谱 title(|X(k)|); %将幅度谱命名 axis([-1,N,1.1*min(abs(Xk)),1.1*max(abs(Xk))]); 设定纵坐标的范围 subplot(2,2,4),stem(k,angle(Xk),k);%显示序列的相位谱 title(arg|X(k)|); 将相位谱命名 axis([-1,N,1.1*min(angle(Xk)),1.1*max(angle(Xk))]); 将相位谱横坐标范围 程序3： xn=[0,1,2,3,4,5,6,7]; %设定n的取值 N=length(xn); %将N设定为n取值总长度 n=0:4*N-1;k=0:4*N-1; %设定n的变化范围 xn1=xn(mod(n,N)+1);%xn1=[xn,xn,xn,xn]; Xk=xn1*exp(-j*2*pi/N).^(n*k); %离散傅里叶变换 subplot(2,2,1),stem(xn,k); %显示序列主值 title(原主值信号x(n)); %将原信号命名 subplot(2,2,2),stem(n,xn1,k); %显示周期序列 title(周期序列信号); %将周期号命名 axis([-1,4*N,1.1*min(xn1),1.1*max(xn1)]); 设定横坐标的范围 subplot(2,2,3),stem(k,abs(Xk),k); %显示序列的幅度谱 title(|X(k)|); %将幅度谱命名 axis([-1,4*N,1.1*min(abs(Xk)),1.1*max(abs(Xk))]); 设定纵坐标的范围 subplot(2,2,4),stem(k,angle(Xk),k);%显示序列的相位谱 title(arg|X(k)|); %将相位谱命名 axis([-1,4*N,1.1*min(angle(Xk)),1.1*max(angle(Xk))]); 设定纵坐标的范围 程序4： xn=[0,1,2,3,4,5,6,7]; %设定n的取值 N=length(xn); %将N设定为n取值总长度 n=0:N-1; %设定n的变化范围 w=linspac