下载文档原格式
实验题目：用Matlab语言编程实现学生成绩管理系统一、实验目的1.掌握建立和执行M文件的方法，学会M函数的编写及其调用，掌握多种函数类型；2.学习Matlab中的关系运算和逻辑运算，掌握它们的表达形式和用法；3.掌握Matlab中的选择结构和循环结构；4.学习使用帮助文档查找更多的Matlab功能；5.掌握Matlab的基本绘图命令，学习运用Matlab绘制一维、二维、三维图形的方法。
二、实验要求1.使用函数（多个函数M文件）进行分析和设计；2.要使用到以下函数类型：主函数、子函数、私有函数。
三、实验内容功能内容：1.学生信息为：学号、姓名、三门课成绩（英语、数学、物理）；2.从键盘输入学生成绩，建立信息库；3.从文件读入学生信息；4.显示信息库中的所有学生的信息；5.将学生成绩信息库保存为文件；6.根据学号查询某一个学生的成绩统计信息：三门课成绩、总分、平均分；7.根据学号将某一个学生的信息从信息库删除；8.新增一个学生的信息；9.按总分成绩从高到低排序输出学生信息；10.图形输出所有学生成绩统计信息（饼图pie、柱状图bar）。
四、实验设计及分析1 主函数的设计及分析主函数主要用于显示学生管理系统的主界面，以及通过对其他几个私有函数的调用来实现相应的功能。
主函数名为zonghe()，首先利用disp()屏幕输出函数来提示用户如何使用该管理系统。
然后通过input()输入函数输入字符，通过判断输入的字符与哪个私有函数对应的字符匹配来调用相应的私有函数实现相应的功能。
其中将条件语句放入while循环中，使系统能循环运行直到用户输入‘H‘才跳出循环。
主函数zonghe()的源代码：function zonghe()disp('Welcome to the students'' information library!');disp('Input A to input students'' information');disp('Input B to display all students'' information');disp('Input C to search a stuedent''s information');disp('Input D to delete a stuedent''s information');disp('Input E to add a new stuedent''s information');disp('Input F to display stuedents'' information from high to low of total score');disp('Input G to graphical display all stuedents'' information'); disp('Input H to exit');x=input('Please input the FUNCTION NUMBER:','s');while x~='H'if x=='A'disp('Input students'' information.');ztest();elseif x=='B'disp('Display all students'' information.');displaytest();elseif x=='C'disp('Search a stuedent''s information.');search();elseif x=='D'disp('Delete a stuedent''s information.');shanchu();elseif x=='E'disp('Add a new stuedent''s information.');add();elseif x=='F'disp('Display stuedents'' information from high to low of total score.');paixu();elseif x=='G'disp('Graphical display all stuedents'' information.');tuxing();elseif x=='H'disp('Thank you for using the students'' information library!');break;elsedisp('Invalid function number. Program is terminated!');endx=input('\nPlease input the FUNCTION NUMBER:','s');endend2 七个私有函数的设计与分析2.1 function ztest()的设计与分析ztest()函数实现从键盘输入学生成绩建立信息库的功能。
１. 假定一因果系统为y(n)-0.4y(n-1)+0.75y(n-2)=2.2403x(n)+2.4908x(n-1)+2.2403x(n-2)用MA TLAB 程序仿真该系统，输入三个不同的输入序列：)1.02c o s()(1n n x ⋅=π，)4.02cos()(2n n x ⋅=π，)(3)(221n x n x x -= 计算并并显示相应的输出)(1n y ， )(2n y 和)(n y 。
clf;n=0:40;a=2;b=-3;x1=cos(2*pi*0.1*n);x2=cos(2*pi*0.4*n);x=a*x1+b*x2;num=[2.2403, 2.4908,2.2403];den=[1 -0.4 0.75];ic=[0 0]; %设置零初始条件y1=filter(num,den,x1,ic); %计算输入为x1(n)时的输出y1(n)y2=filter(num,den,x2,ic); %计算输入为x2(n)时的输出y2(n)y=filter(num,den,x,ic); %计算输入为x (n)时的输出y(n)yt= a*y1+b*y2;%画出输出信号subplot(2,1,1)stem(n,y);ylabel('振幅');title('加权输入a*x1+b*x2的输出');subplot(2,1,2)stem(n,yt);ylabel('振幅')'title('加权输出a*y1+b*y2');0510152025303540-40-2020振幅加权输入a*x1+b*x2的输出0510152025303540-40-2020振幅加权输出a*y1+b*y2２. 用MATLAB 程序仿真步骤1给出的系统，对两个不同的输入序列x(n)和x(n-10)，计算并显示相应的输出序列y3(n)和y4(n)。
clf;n=0:40;n1=0:50a=2;b=-3;x1= cos(2*pi*0.1*n);x2=cos(2*pi*0.4*n);x=2*x1-3*x2;x3=[zeros(1,10) x];num=[2.2403, 2.4908,2.2403];den=[1 -0.4 0.75];ic=[0 0]; %设置零初始条件y1=filter(num,den,x,ic); %计算输入为x1(n)时的输出y1(n)y2=filter(num,den,x3,ic); %计算输入为x2(n)时的输出y2(n)subplot(2,1,1)stem(n,y1);ylabel('振幅');title('时移前');subplot(2,1,2)stem(n1,y2);ylabel('振幅');title('时移后');0510152025303540-30-20-1001020振幅时移前05101520253035404550-40-20020振幅时移后3．用MA TLAB 程序仿真计算下列两个有限长序列的卷积和并显示图形。
目录实验一数字信号的产生和频谱分析实验 (2)一.实验目的 (2)二.实验要求 (2)三.实验原理 (2)四.实验步骤 (3)五.流程图 (3)七.实验结果分析 (6)实验二FIR数字滤波器设计 (7)一.试验目的 (7)二.实验要求 (7)三.实验原理和步骤 (7)四.实验流程图 (7)五.实验波形 (8)六.实验结果分析 (11)实验三IIR数字滤波器设计 (11)一.试验目的 (11)二.实验要求 (11)三.实验原理 (12)四.实验流程图 (12)五．实验波形 (13)六.实验结果分析 (14)实验四模拟调制解调 (14)一.试验目的 (14)二.实验要求 (14)三．实验原理 (15)四.实验流程图 (15)五.实验波形 (16)六.实验结果分析 (17)实验五数字调制解调 (17)一.试验目的 (17)二.实验要求 (17)三.实验原理 (18)四.实验流程图 (18)五.实验波形 (18)六.实验结果分析 (20)实验一数字信号的产生和频谱分析实验一.实验目的1.通过仿真掌握采样定理2.掌握利用FFT进行信号谱分析的原理二.实验要求1. 按照采样定理生成CW信号和LFM信号；2. 画出信号时域波形图和频谱图；3. 生成高斯分布的白噪声；4. 生成一定信噪比的带噪信号，并对其进行谱分析。
三.实验原理1.采样定理:在模拟信号数字化时,需要对模拟信号进行采样,当采样频率fs.max大于信号中最高频率fmax的2倍时(fs.max>=2fmax)，采样之后的数字信号完整地保留了原始信号中的信息，否则会发生频谱混叠,造成最后解调出来的信号失真.一般实际应用中保证采样频率为信号最高频率的5～10倍；采样定理又称奈奎斯特定理。
2.快速傅立叶变换(FFT)算法长度为N的序列的离散傅立叶变换为：N点的DFT可以分解为两个N/2点的DFT，每个N/2点的DFT又可以分解为两个N/4点的DFT。
Matlab数字信号处理实验报告本次实验使用MATLAB进行数字信号处理操作，目的是熟悉MATLAB中数字信号处理的相关工具箱，并进一步理解数字信号处理的基本概念和算法。
一、实验内容1.信号的生成与显示2.时域分析和频域分析3.滤波器设计4.数字滤波器性能分析二、实验步骤在MATLAB中，使用sawtooth函数生成一个锯齿波信号，并使用plot函数进行时域波形的显示。
代码如下：f = 1000;fs = 40000;t = 0:1/fs:0.01;y = sawtooth(2*pi*f*t);plot(t,y);xlabel('Time (s)');ylabel('Amplitude');title('Sawtooth Wave');时域分析包括波形的观察和参数分析，如幅值、均值、方差等。
频域分析则是对信号进行傅里叶变换，得到其频谱图，包括频率分布和强度分布。
%时域分析amp = max(y)-min(y);mean_y = mean(y);var_y = var(y);设计一个低通滤波器，将高于1kHz的频率成分滤掉。
对滤波后的信号进行时域分析和频域分析，比较滤波前后信号的特征参数和频谱特征，并绘制原始信号、滤波后信号及其频谱图。
subplot(2,2,1);plot(t,y);xlabel('Time (s)');ylabel('Amplitude');title('Sawtooth Wave');subplot(2,2,2);plot(t,y_filt);xlabel('Time (s)');ylabel('Amplitude');title('Sawtooth Wave After Filter');subplot(2,2,3:4);plot(f2,fft_y_filt,'r',f,fft_y,'g');xlabel('Frequency (Hz)');ylabel('Amplitude');title('Sawtooth Wave Spectrum Comparison');legend('After Filter','Before Filter');三、实验结果与分析通过生成并显示一段锯齿波信号，并对其进行时域和频域分析，可以得到该信号的关键信息，如幅值、均值、方差和频率分布特性。
MATLAB基础一、MATLAB的工作环境1．命令窗口：是用户和MATLAB系统交互的主要窗口。
在该窗口中，用户可以运行函数，执行MATLAB 的基本操作命令以及对MATLAB系统的参数设置等操作。
在命令提示符>>后输入命令如：t=[1,2,3;4,5,6;7,8,9];完成对t的赋值2．帮助的使用：直接在命令行输入help在菜单栏选择help—matlab help直接按F13．图形窗（Figure）：用于显示绘出的图形。
通常只要执行了任意一种绘图命令，图形窗都会自动产生。
绘图都是在这个图形窗中进行。
如果要在建一个图形窗，则可在命令窗中输入figure命令，MATALB 就回新建一个图形窗口，并自动给它排出序号。
4．文本编辑窗：其作用是用来创建、编辑和调试MATLAB的相关文件（或称程序，即.M文件），它与一般的编辑调试器有相似的功能。
通常，MA TLAB命令编辑有行命令方式和文件两种。
行命令方式，即在命令窗口中一行一行的输入命令，计算机对每一行命令做出反应。
这种方式，只能用于编辑简单的程序，当程序比较复杂的时候，把程序写成一个由多行语句组成的文件（.M文件），让MATLAB来执行这个程序中的全部语句，MA TLAB文本编辑器的功能就是完成编写，修改和调试这种程序。
其进入方式：1）file—new/open—M-file2）直接在命令窗口输入：edit/edit 文件名二、MATLAB的基本语法1．变量及其赋值1）赋值要求在MATLAB中，变量和常量的标识符最长允许19个字符。
MA TLAB内部只有一种数据格式，就是双精度类型，对应于64位二进制。
赋值就是把数赋予代表常量或变量的标识符。
基本格式为：变量=表达式（或数）在MA TLAB中，变量都代表矩阵。
列矢量可被当做只有一行的矩阵，行矢量也可被当做只有一个行的矩阵，标量（或常数）应该做是1×1的矩阵。
输入矩阵时，应遵循一下规则：整个矩阵的值都应放在方括号中；同一行的个元素以逗号或空格分开；不同行的元素以分号隔开。
第一次实验MAT1AB的数字信号处理基础一、实验目的1.掌握在MAT1AB中创建和编写脚本程序的步骤，熟悉在MAT1AB中进行基本的矩阵运算2.了解数字信号在计算机系统中的表示和数字信号处理的基本过程3.掌握用MAT1AB产生数字信号的方法二、实验内容1.熟悉在MAT1AB下创建脚本文件编制程序的方法在MAT1AB中，脚本文件时由一系列的命令构成并储存为.m格式的文件。
通常使用m文件来编写一个完整的仿真程序。
脚本文件的创建，可采用以下两种方式：(1)在菜单栏中选择Fi1e下拉框中New选项，可以新建多种MAT1AB文件，我们编辑仿真程序，通常选择第一项BIankM-Fi1e,即新建一个空的MAT1AB文件。
具体参见下图。
图1(2)采用菜单栏中新建按钮即可新建一个空的MAT1AB文件。
Jk MAT1AB7.9.0(R2009b)Fi1eEditDebugParaJ一一J2.MAT1AB中的数据表示MAT1AB中的基本数据单元为数组矩阵，MAT1AB中的数学运算都是基于矩阵的。
掌握了矩阵运算，就掌握了MAT1AB编程的关键。
MAT1AB中使用到的变量无需事先声明其数据类型，大小等，MAT1AB会自动根据赋值情况进行解析。
比如，可用通过以下命令产生一个矩阵：3.常用序列的MAT1AB实现(1)单位抽样序列。
在MAT1AB中可以用以下函数来实现单位抽样序列function[x,n]=impseq(nθ,n1z n2)%产生x(n)=de1ta(n-n0);n1<=n0<=n2if((nθ<n1)∣(nθ>∩2)∣(n1>n2))error('参数必须满足n1<=nθ<=n2,)endn=[n1:n2];x=[zeros(1,(nO-n1))4∕Zθros(1z(n2-nO))];stem(x);图3(2)单位阶跃序列。
在MAT1AB中可用〃>=0来实现〃(〃一%)。
实验1 离散时间信号的时域分析一、实验目的（1）了解MATLAB 语言的主要特点及作用；（2）熟悉MATLAB 主界面，初步掌握MATLAB 命令窗和编辑窗的操作方法；（3）学习简单的数组赋值、数组运算、绘图的程序编写；（4）了解常用时域离散信号及其特点；（5）掌握MATLAB 产生常用时域离散信号的方法。
二、知识点提示本章节的主要知识点是利用MATLAB 产生数字信号处理的几种常用典型序列、数字序列的基本运算；重点是单位脉冲、单位阶跃、正（余）弦信号的产生；难点是MATLAB 关系运算符“==、>=”的使用。
三、实验内容1. 在MATLAB 中利用逻辑关系式0==n 来实现()0n n -δ序列，显示范围21n n n ≤≤。
（函数命名为impseq(n0,n1,n2)）并利用该函数实现序列：()()()632-+-=n n n y δδ；103≤≤-nn 0212. 在MATLAB 中利用逻辑关系式0>=n 来实现()0n n u -序列，显示范围21n n n ≤≤。
（函数命名为stepseq(n0,n1,n2)）并利用该函数实现序列：()()()20522≤≤--++=n n u n u n y3. 在MATLAB 中利用数组运算符“.^”来实现一个实指数序列。
如： ()()5003.0≤≤=n n x n4. 在MATLAB 中用函数sin 或cos 产生正余弦序列，如：()()2003.0cos 553.0sin 11≤≤+⎪⎭⎫ ⎝⎛+=n n n n x πππ5. 已知()n n x 102cos 3π=，试显示()()()3,3,+-n x n x n x 在200≤≤n 区间的波形。
6. 参加运算的两个序列维数不同，已知()()6421≤≤-+=n n u n x ，()()8542≤≤--=n n u n x ，求()()()n x n x n x 21+=。
Matlab实验报告实验一：1.实验Matlab代码：N=25;Q=0.9+0.3*j;WN=exp(-2*j*pi/N);x=zeros(25,1);format long; %长整型科学计数for k0=1:25x(k0,1)=Q^(k0-1);end;for k1=1:25;X1(k1,1)=(1-Q^N)/(1-Q*WN^(k1-1));end;X1;X2=fft(x,32);subplot(3,1,1);stem(abs(X1),'b.');axis([0,35,0,15]);title('N=25,formular');xlabel('n'); subplot(3,1,2);stem(abs(X2),'g.');axis([0,35,0,15]);title('N=32, FFT');xlabel('n');for(a=1:25)X3(a)=X1(a)-X2(a)end;subplot(3,1,3);stem(abs(X3),'r.');title('difference');xlabel('n');实验结果如图：实验结论：可以看出基2时间抽选的FFT算法与利用公式法所得到的DFT结果稍有偏差，但不大，在工程上可以使用计算机利用FFT处理数据。
2.实验Matlab代码：N = 1000; % Length of DFTn = [0:1:N-1];xn = 0.001*cos(0.45*n*pi)+sin(0.3*n*pi)-cos(0.302*n*pi-pi/4);Xk = fft(xn,N);k=[0:1:N-1];subplot(5,1,1);stem(k,abs(Xk(1:1:N)));title('DFT x(n)');xlabel('k');axis([140,240,0,6])subplot(5,1,2);stem(k, abs(Xk(1:1:N)),'r');%画出sin(0.3npi)-cos(0.302npi-pi/4) axis([140,160,0,6]);title('sin(0.3*pi*n)-cos(0.302*pi*n) ');xlabel('k');subplot(5,1,3);stem(k, 1000*abs(Xk(1:1:N)),'g');%画出0.001*cos(0.45npi)axis([220,230,0,6]);title('cos(0.45*pi*n) ');xlabel('k');subplot(5,1,4);stem(k,0.01*abs(Xk(1:1:N)),'k');%画%sin(0.3npi)-cos(0.302npi-pi/4)axis([140,160,0,6]);title('sin(0.3*pi*n)-cos(0.302*pi*n) ');xlabel('k');subplot(5,1,5);stem(k, 10*abs(Xk(1:1:N)),'m');%画出0.001*cos(0.45npi)axis([220,230,0,6]);title('cos(0.45*pi*n) ');xlabel('k');实验结果如图：实验结论：由上图及过程可知，当DFT变换长度为1000时所得到的谱线非常理想。
目录实验一学习使用MATLAB.................................................... 错误！未定义书签。
一. 实验目的 ............................................................ 错误！未定义书签。
二. 实验内容和要求 (1)三. 实验主要仪器设备和材料 (1)四. 实验方法、步骤及结果测试 (1)五. 实验报告要求 (18)六. 思考题 (18)实验二用FFT做谱分析 (19)一. 实验目的 (19)二. 实验内容和要求 (19)三. 实验主要仪器设备和材料 (20)四. 实验方法、步骤及结果测试 (20)五. 实验报告要求 (40)六. 思考题 (40)实验三用双线性变换法设计IIR数字滤波器 (41)一. 实验目的 (41)二. 实验内容和要求 (41)三. 实验主要仪器设备和材料 (41)四. 实验方法、步骤及结果测试 (41)五. 实验报告要求 (46)六. 思考题 (46)实验四用窗函数法设计FIR数字滤波器 (48)一. 实验目的 (48)二. 实验内容和要求 (48)三. 实验主要仪器设备和材料 (49)四. 实验方法、步骤及结果测试 (49)五. 实验报告要求 (59)六. 思考题 (59)附录一MA TLAB信号处理工具箱函数 (60)一. 表附1-1 波形产生 (60)二. 表附1-2 滤波器分析和实现 (60)三. 表附1-3 线性系统变换 (61)四. 表附1-4 IIR滤波器设计 (62)五. 表附1-5 IIR滤波器阶的选择 (62)六. 表附1-6 FIR滤波器设计 (63)七. 表附1-7 变换 (63)八. 表附1-8 统计信号处理 (64)九. 表附1-9 窗函数 (64)十. 表附1-10 参数化建模 (65)十一. 表附1-11 特殊操作 (65)十二. 表附1-12 模拟原形滤波器设计 (66)十三. 表附1-13 频率变换 (66)十四. 表附1-14 滤波器离散化 (66)十五. 表附1-15 其他 (67)附录二实验中用到的一些子程序 (68)实验内容和要求(1)了解MATLAB 的基本程序设计原则，常量和变量的用法(2)掌握MATLAB中对矩阵进行输入、运算和比较的方法(3)了解循环语句的类型，并掌握循环语句的用法(4)熟悉M文件的作用，并掌握二维图形的绘制一.实验主要仪器设备和材料计算机，MATLAB6.5或以上版本二.实验方法、步骤及结果测试关于MATLAB它是由美国的Math Works 公司推出的一个科技应用软件，它的名字是由MATRIX(矩阵）和LABORATORY（实验室）的前三个字母组合而成MATLAB是一种高性能的、用于工程计算的编程软件，它把科学计算、结果的可视化和编程都集中在一个使用方便的环境中优势在于能很容易求解复数数值问题，速度快且容易扩展创建新的命令和函数主要组成部分：（1）编程语言：以矩阵和数组为基本单位的编程语言（2）工作环境：包括一系列的应用工具，提供编程和调试程序的环境（3）图形处理：包括绘制二维、三维图形和创建图形用户界面（GUI）等（4）数学库函数：包含大量数学函数，也包括复杂功能（5）应用程序接口：提供接口程序，可使MATLAB与其他语言程序进行交互典型特点：（1）语言简洁紧凑，运算符十分丰富，使用方便灵活（2）既具有结构化的控制语言，又能面向对象编程（3）语法限制不严格，程序设计自由度大，可移植性好（4）具有强大的图形功能（5）包含功能强劲的工具箱（6）最重要、最受欢迎的特点是它的开放性（7）典型应用有两方面：数值计算和符号计算建模和动态仿真下面介绍MATLAB的界面、常用命令和使用方法菜单栏和工具栏：位于窗口顶部，用户可以通过它们来执行某些命令命令窗口：位于右边空白部分，用户的数据输入和结果运算，都在此窗口进行，是Matlab 极为重要的部分，也是用户使用最频繁的部分工作台和工具箱：位于主窗口左上部分，双击工具箱或前面的‖+‖号，就能看到工具箱的各项功能工作空间：主窗口的中上部分，可看到Matlab 的各个工作变量，新打开Matlab 时，只能看到系统提供的默认输出变量ans历史命令：主窗口的左下部分，主要保存工作过的变量、表达式等，需要时，用户可以直接提取历史命令在命令窗口中使用当前工作目录：主窗口的中下部分，主要保存在当前工作路径下的图形文件和命令文件二、常量和变量MATLAB 中使用的数据有常量和变量作用标量的实数常量，类似于 C 语言中的整形常量和实形常量，常量的表示方法：十进制数字，可以带正负号和小数点，也可采用科学计数法如：1, 2.5 , 0.0033 , 2 e-7 ，pi , 2+3 i 等变量以其名称在操作语句中第一次合法出现而定义，无需事先定义。
《数字信号处理》Matlab 实验一．离散信号的 FFT 分析1.用Matlab 编程上机练习。
已知：N=2^5。
这里Q=0.9+j0.3。
可以推导出 ，首先根据这个式子计算X(k)的理论值，然后计算输入序列x(n)的32个值，再利用基2时间抽选的FFT 算法，计算x(n)的DFT X(k)，与X(k)的理论值比较（要求计算结果最少6位有效数字）。
解：函数代码：>> function xn()>> format long>> q=0.9+0.3*i;>> wn=exp(-2*pi*i/32);>> xk=(1-q^32)./(1-q*wn.^[0:31])>>xn=q.^[0:31]>> xk1=fft(xn,32)>>diff=xk-xk1具体执行情况：>> function xn()format longq=0.9+0.3*i;wn=exp(-2*pi*i/32);xk=(1-q^32)./(1-q*wn.^[0:31])xk =Columns 1 through 20.5698 + 3.3840i 2.8369 + 8.8597iColumns 3 through 49.3189 - 9.8673i 1.2052 - 3.5439iColumns 5 through 61.8846 -2.0941i 0.8299 - 1.2413i11,011)()()(k k 10nk 10-=--===∑∑-=-=N k QW Q QW W n x k X N NnN N n N N n ，0.9214 - 1.0753i 0.3150 - 0.0644i Columns 9 through 100.9240 - 0.8060i 0.4202 - 0.2158i Columns 11 through 120.8513 - 0.6357i 0.5040 - 0.1701i Columns 13 through 140.6217 - 0.6931i 0.2441 - 0.8978i Columns 15 through 160.9454 - 0.2800i 0.7139 - 0.3158i Columns 17 through 180.6723 - 0.6496i 0.0263 + 0.5093i Columns 19 through 200.5671 + 0.6914i 0.3173 + 0.9841i Columns 21 through 220.8929 + 0.7792i 0.4066 + 0.8452i Columns 23 through 240.5847 + 0.9017i 0.9129 + 0.9283i Columns 25 through 260.0573 + 0.5531i 0.4219 + 0.9562i Columns 27 through 280.3298 + 0.3143i 0.4513 + 0.2638i0.7214 + 0.1879i 0.0933 + 1.7793iColumns 31 through 320.9483 + 1.9802i 0.4932 + 2.6347i>> xn=q.^[0:31]xn =Columns 1 through 21.0000 0.0000 + 0.0000i Columns 3 through 40.0000 + 0.0000i 0.0000 + 0.0000iColumns 5 through 60.0000 + 0.0000i -0.0000 + 0.0000iColumns 7 through 8-0.0000 + 0.0000i -0.0000 + 0.0000iColumns 9 through 10-0.0000 + 0.0000i -0.0000 + 0.0000iColumns 11 through 12-0.0000 - 0.0000i -0.0000 - 0.0000iColumns 13 through 14-0.2000 - 0.4000i -0.3600 - 0.5200iColumns 15 through 16-0.9680 - 0.1760i 0.4816 - 0.3488i0.2381 - 0.2695i 0.2951 - 0.1711i Columns 19 through 200.1169 - 0.4655i 0.4449 - 0.9838i Columns 21 through 220.6955 + 0.2480i 0.5516 + 0.4319i Columns 23 through 240.3669 + 0.5542i 0.9639 + 0.2088i Columns 25 through 260.3049 + 0.9771i -0.5187 + 0.4709i Columns 27 through 28-0.6081 + 0.2682i -0.1278 + 0.5589i Columns 29 through 30-0.4827 + 0.0647i -0.6538 + 0.5134i Columns 31 through 32-0.8425 - 0.4341i -0.1280 - 0.1434i >> xk1=fft(xn,32)xk1 =Columns 1 through 20.5698 + 3.3839i 2.8366 + 8.8599i Columns 3 through 49.3182 - 9.8692i 1.2051 - 3.5439i1.8845 -2.0942i 0.8298 - 1.2413i Columns 7 through 80.9213 - 1.0754i 0.3150 - 0.0645i Columns 9 through 100.9240 - 0.8060i 0.4202 - 0.2158i Columns 11 through 120.8514 - 0.6356i 0.5040 - 0.1701i Columns 13 through 140.6217 - 0.6931i 0.2441 - 0.8977i Columns 15 through 160.9454 - 0.2800i 0.7139 - 0.3159i Columns 17 through 180.6723 - 0.6496i 0.0263 + 0.5093i Columns 19 through 200.5671 + 0.6913i 0.3172 + 0.9840i Columns 21 through 220.8929 + 0.7792i 0.4065 + 0.8452i Columns 23 through 240.5846 + 0.9016i 0.9129 + 0.9283i Columns 25 through 260.0572 + 0.5531i 0.4219 + 0.9563i0.3297 + 0.3144i 0.4512 + 0.2638iColumns 29 through 300.7213 + 0.1879i 0.0932 + 1.7793iColumns 31 through 320.9480 + 1.9802i 0.4928 + 2.6347i>> diff=xk-xk1diff =1.0e-013 *Columns 1 through 20.4625 + 0.8501i 0.9504 - 0.4003iColumns 3 through 40.6010 + 0.4028i 0.4752 + 0.7001iColumns 5 through 60.5502 + 0.8501i 0.4625 + 0.8501iColumns 7 through 80.7751 + 0.9250i 0 + 0.3875i Columns 9 through 100.7751 - 0.4625i 0.3126 - 0.4625iColumns 11 through 12-0.4625 - 0.3126i 0.4625 + 0.3875iColumns 13 through 14-0.9250 + 0.6938i 0.3875 - 0.0781iColumns 15 through 160.3875 - 0.6156i 0 + 0.9641iColumns 17 through 180.9250 - 0.7598i -0.4625 - 0.0422iColumns 19 through 200.4625 + 0.1172i 0.4625 + 0.3094iColumns 21 through 220.9250 + 0.4625i 0.9250 + 0.2313iColumns 23 through 240.3875 + 0.1563i 0.3875 - 0.2313iColumns 25 through 260.8501 -0.9250 - 0.4625iColumns 27 through 280.0127 - 0.7751i 0.7001 - 0.9250iColumns 29 through 300.1626 0.7814 - 0.9250iColumns 31 through 320.4816 + 0.9250i 0.7255 - 0.8501i由以上结果可知，由基2时间抽选的FFT算法所得到的DFT结果与利用公式法所得的理论值稍有偏差，但误差较小，从结果可以看出大概在小数点第15位才开始出现误差，故而用计算机FFT处理数据在精度上是可以接受的。
北京邮电大学DSP之matlab实验1实验题目2. 程序分析实验一(源代码及分析)A=zeros(1,25)c=complex(0.9,0.3)for i=1:1:25A(i)=c^iendn=0:24k=0:24WN=exp(-j*2*pi/25)nk=n'*kWNnk=WN.^nkX=A*WNnk;%X(k)为用定义算法计算的DFT值%下面为32点运算for i=26:1:32A(i)=0endn=0:31k=0:31WN=exp(-j*2*pi/32)nk=n'*kWNnk=WN.^nkX=A*WNnk%下面为32点蝶形运算k1=bin2dec(fliplr(dec2bin([1:32]-1,5)))+1%得到从1到32的倒序顺序B=A(k1)%得到倒序后的A，将其赋给Bfor mm=1:5 %将DFT做m次基2分解，从左到右，对每次分解作DFT运算Nmr=2^mm;u=1; %旋转因子u初始化WN=exp(-j*2*pi/Nmr); %本次分解的基本DFT因子WN＝exp(-i*2*pi/Nmr)for n=1:Nmr/2 %本次跨越间隔内的各次碟形运算for k=n:Nmr:32 %本次碟形运算的跨越间隔为Nmr=2^mmkp=k+Nmr/2; %确定碟形运算的对应单元下标(对称性)t=B(kp)*u; %碟形运算的乘积项B(kp)=B(k)-t; %碟形运算的加法项B(k)=B(k)+t;endu=u*WN; %修改旋转因子，多乘一个基本DFT因子WNendend实验二(源代码及分析)for i=1:1:1000x(i)=0.001*cos(0.45*i*pi)+sin(0.3*i*pi)-cos(0.302*i*pi-4/pi)endX=fft(x);k=0:500;w=2*pi/1000*k;stem(w/pi,abs(X(1:501)))axis([0.25 0.35 490 510])axis([0.3 0.5 0 0.8]) 3. 程序运行结果实验一程序运行结果太长,只附一截图如下:实验二结果的有窗与无窗图如下:4. 总结本次实验是第一次使用matlab,实验过程中遇到了不少的问题,经过和同学们的交流以及通过网络查询,最后得到了结果。
实验一离散时间信号与系统一、实验目的：1、熟悉常见离散时间信号的产生方法；2、熟悉离散时间系统的单位脉冲响应和单位阶跃响应的求解方法；3、熟悉离散时间信号经过离散时间系统的响应的求解方法。
二、实验内容：已知离散时间系统差分方程为y(n)-0.5y(n-1)+0.06y(n-2)=x(n)+x(n-1)，求1、该系统的单位脉冲响应并绘图；2、该系统的单位阶跃响应并绘图；3、已知x(n)=可自己指定用filter函数经过系统的响应并绘图；4、用conv_m函数求系统响应并绘图。
三、实验平台：MA TLAB集成系统四、设计流程：此处写个人自己的设计流程五、程序清单：此处写程序内容六、调试和测试结果：此处写程序的执行结果和实验过程中的调试经过、出现的错误和对应的解决方法七、教师评语与成绩评定此处由老师填写上机操作：实验一离散时间信号与系统实验内容：1．脉冲响应>> b =[1,1]; a = [1,-0.5,0.06];n = [-10:25];>> impz(b,a,n);>> title('Impulse Response'); xlabel('n'); ylabel('h(n)')2.单位阶跃响应>> x = stepseq(0,-10,25); s = filter(b,a,x);Warning: Function call stepseq invokes inexact match d:\MATLAB7\work\STEPSEQ.M.>> stem(n,s)>> title('Step Response'); xlabel('n');ylabel('s(n)')3.>> a=[1,-0.5,0.06];b=[1,1];>> n=-20:120;>> x1=exp(-0.05*n).*sin(0.1*pi*n+pi/3);>> s1=filter(b,a,x1);>> stem(n,s1);;xlabel('n');ylabel('s1(n)');4．>> a=[1,-0.5,0.06];b=[1,1];>> n=-20:120;>> h=impz(b,a,n);>> x1=exp(-0.05*n).*sin(0.1*pi*n+pi/3);>> [y,m]=conv_m(x1,n,h,n);Warning: Function call conv_m invokes inexact match d:\MATLAB7\work\CONV_M.M. >> stem(m,y);title('系统响应');xlabel('m');ylabel('y(m)');实验二离散信号与系统的连续频域分析一、实验目的：1、掌握离散时间信号的DTFT的MATLAB实现；2、掌握离散时间系统的DTFT分析；3、掌握系统函数和频率相应之间的关系。
