实验5综合应用程序设计
C++语言程序设计形考二实验五

C++语言程序设计形考二实验五上机实验5:数组应用程序设计一、实验内容有如下程序:#include<iostream>using namespacestd;intmain(){__________________ //定义具有6个元素的整形数组dint data,i=0.// data用于暂存每一个输入数据,i记录已插入数据的个数。
cout<<"请输入6个整数:";do{ //输入并处理6个数据cin。
data;int j=i。
//j是候选的插入位置while(j>0){ //若j==0,它就是第一个数据的插入位置,无需进入循环。
if(data>=d[j-1]) __________。
//若找到插入位置,退出循环。
d[j]=d[j-1]。
//否则数据后移一个位置j--;}d[j]=data。
//数据插入到位置j处i++。
//i记录已插入数据的个数}while(________)。
//若缺乏6个数据,连续输入。
for(_______________________) cout<<d[i]<<' ';//表现数组中的6个数据cout<<endl;return 0;}此程序输入6个整数保存到一维数组d中,使得输入的数据在数组中始终按从小到大的顺序排列,最后顺序输出数组中的数据。
例如,若输入是5 7 4 29 3,程序的输出是:2 3 4 5 79.二、实验要求读懂步伐,弥补步伐中下划线处缺失的部分,上机测试,直至产生精确的输出成效。
三、实验结果截图。
python程序设计实验报告五:循环结构程序设计(综合性实验)

python程序设计实验报告五:循环结构程序设计(综合性实验)安徽⼯程⼤学Python程序设计实验报告班级物流192 姓名张⽻学号 3190505221 成绩⽇期 2020.4.23 指导⽼师修宇实验五综合运⽤三种基本结构进⾏程序设计(综合性实验)(⼆学时)【实验⽬的】(1)进⼀步掌握选择与循环结构程序设计的⽅法和技巧;(2)综合运⽤三种基本结构进⾏程序设计;【实验条件】PC机或者远程编程环境【实验内容】1、完成三个编程题。
(1)计算并输出指定范围内的素数素数⼜称质数。
指整数在⼀个⼤于1的⾃然数中,除了1和此整数⾃⾝外,没法被其他⾃然数整除的数。
换句话说,只有两个正因数(1和⾃⼰)的⾃然数即为素数。
⽐1⼤但不是素数的数称为合数。
1和0既⾮素数也⾮合数。
素数在数论中有着很重要的作⽤。
2300多年前,古希腊数学家欧⼏⾥得在《⼏何原本》⼀书中证明了素数有⽆穷多个,如2、3、5、7、11等等。
问题描述:编程实现输出指定范围内的素数,输⼊以逗号间隔,输出以空格间隔。
指定范围由⽤户输⼊指定。
(2)九九乘法表问题描述:请按照如下格式设计⼀个“九九乘法⼝诀表”,并打印出来。
编程要点及提⽰:1. 使⽤for循环或者while循环实现以上功能2. 代码量不超过10⾏3. 在⼀⾏中,每个等式之间使⽤制表符"\t"进⾏分隔;对于每⼀⾏等式之间使⽤print()打印空字符实现,即print("")(3)统计不同字符个数问题描述:⽤户从键盘输⼊⼀⾏字符。
请编写⼀个程序,统计并输出其中英⽂字符、数字符号、空格和其他字符的个数。
输⼊格式:输⼊⼀⾏字符输出:按照顺序输出字母,数字,空格,及其他字符的个数,个数以空格分隔如:【实验记录与分析】(请在填写实验记录与分析结果)1.计算并输出指定范围内的素数输⼊代码,输出结果及分析如下:易错点:1.程序⽆法同时处理两个输⼊数据,在输⼊未处理时。
2.本题中的else语句不能和第⼆个的if对齐。
实验5

实验五第3章MATLAB程序设计第二节程序控制结构程序控制结构有三种:顺序结构、旋择结构、循环结构.任何复杂的程序都可以由这三种基本结构构成.Matlab提供了实现控制结构的语句,利用这些语句可以编写解决实际问题的程序.一、顺序结构顺序结构是指按照程序中语句的排列顺序从上到下依次执行,直到程序的最后一个语句如例1的qiu文件和例2的dd1文件.这是最简单的一种程序结构.一般涉及数据的输入、计算或处理、数据的输出等内容.1. 数据的输入通过input命令来接收从终端输入的内容,它也可以显示文本和提示.命令集14 输入命令input(out ,in)在屏幕上显示出字符串out的文本并等待终端的输入.如果变量in是's',则输入的内容以字符串的形式进行保存,通常MATL AB在保存前要尽可能地求出表达式的值.如果使用格式控制符号如'\n',字符串out可以是若干行.2. 数据的输出可以通过简单地输入变量的名字来显示数字矩阵或者字符串向量的内容,结果将显示出变量的名字和内容.另一种显示变量的值就是使用命令disp.使用它只显示出变量的内容,这是有用的,特别是在字符串的应用中.命令集15显示命令disp(A)显示矩阵A的内容,如果A是字符串,则显示出它的文本.二、循环结构循环结构是按照给定的条件,重复执行指定的语句.Matlab用于实现循环结构的语句有for — end语句和while — end语句.1.for —end语句for-end语句的格式为:for 循环变量=表达式1 : 表达式2 : 表达式3循环体语句end其中表达式1的值为循环变量的初值,表达式2的值为步长,表达式3的值为循环变量的终值.步长为1时,表达式2可以省略.循环语句应用时应注意:● 循环语句可以嵌套使用;● 不能在 for — end 循环体内改变循环变量的值; ● 为了提高代码的运行效率,应尽可能提高代码的向量化程度,避免 for — end 循环的使用,如例8;● 如果预先就知道循环的次数,则可以采用 for — end 循环;否则,如果预先无法确定循环的次数,则可以使用 while — end 循环.2. while — end 循环while-end 语句的一般格式为:while (条件) 循环体语句 end 其执行过程为:若条件成立,则执行循环体语句,执行后再判断条件是否成立,如果不成立则跳出循环.■例1 使用while-end 结构计算1+2+3+…+100 . 解:■例2 计算 Matlab 中 eps 值 (exp1.m) 解:s=0;for i=1:100 s=s+i; end ss=0; i=1;while i<=100 s=s+i; i=i+1; end snum=0; EPS=1; while (1+EPS)>1 EPS=EPS/2; num=num+1; endEPS,num三、选择结构选择结构是根据给定的条件成立或不成立,分别执行不同的语句.Matlab 用于实现选择结构的语句有条件语句:if –end 和 开关语句:switch-end .1.单分支结构 if —end 语句 2.双分支结构 if —else —end 语句 if-end 语句的格式为: if-else-end 语句的格式为:3.多分支结构 if — elseif — else — end 语句 if — elseif — else — end 语句的格式为:if 条件表达式 语句组 endif 条件表达式 语句组1 else语句组2 endif 条件表达式1 语句组1 elseif 条件表达式 2 语句组2 ... ...elseif 条件表达式 m 语句组m else 语句组end■例3求函数1;00;1;x y x x -<⎧⎪= =0⎨⎪ >⎩ ,输入x 值,输出y 值.解:■例4一个三位整数各位数字的立方和等于该数本身则称该数为水仙花数.输出全部水仙花数.解:程序如下运行结果:153 370 371 407function y=fuhao(x) if x<0 y=-1; elseif x==0 y=0; else y=1; endA=[];for m=100:999m1=fix(m/100); %求m 的百位数字 m2=rem(fix(m/10),10); %求m 的十位数字 m3=rem(m,10); %求m 的个位数字 if m==m1*m1*m1+m2*m2*m2+m3*m3*m3 A=[A,m];endend disp(A)4. switch-end 语句根据表达式的不同取值,分别执行不同的语句.其格式如下:■关于 switch-end 语句的几点注解:● Matlab 首先计算表达式的值,然后将它依次与各个case 指令后的检测值进行比较,当比较结果为真时,就执行相应的语句组,然后跳出 switch-end 结构● 如果所有的比较结果都为假,则执行 otherwise 后面的语句组,然后跳出 switch-end 结构 ● otherwise 指令可以不出现● switch 后面的表达式的值可以是一个标量或字符串 ● case 指令后的检测值超过一个时,应用{}括起来switch 表达式 case 值1 语句组1 case 值2 语句组2 ... ... case 值m 语句组m otherwise 语句组n end■例5 阅读程序method=input('请输入方法名:')switch methodcase {'linear','bilinear'}disp('Method is linear')case 'cubic'disp('Method is cubic')case 'nearest'disp('Method is nearest')otherwisedisp('Unknown method.')end■例6用switch-end语句编一程序,对于给定的一个百分制成绩,输出相应的五分制成绩,设:90分以上为…A‟, 80~89分为…B‟,70~79分为…C‟,60~69分为…D‟,60分以下为…E‟.function sctole(score)grade=fix(score/10);switch gradecase {10,9}grade='A'case 8grade='B'case 7grade='C'case 6grade='D'otherwisegrade='E'end四、其他流程控制语句●break语句●continue语句●return语句●try语句●echo命令●keyboard命令●pause命令1.break和continue语句与循环结构相关的语句还有break语句和continue语句.它们一般与if语句配合使用.break语句用于终止循环的执行.当在循环体内执行到该语句时,程序将跳出循环,继续执行循环语句后的下一语句.continue语句控制跳过循环体中的某些语句.当在循环体内执行到该语句时,程序将跳过循环体中所有剩下的语句,继续下一次循环. ■例7求[100,200]之间第一个能被21整除的整数.for n=100:200if rem(n,21)~=0continueendbreakendn for n=100:200if rem(n,21)==0breakendendn2.try语句语句格式为:try语句组1catch语句组2endtry语句先试探性执行语句组1,如果语句组1在执行过程中出现错误,则将错误信息赋给保留的lasterr变量,并转去执行语句组2.这种试探性执行语句是其他高级语言所没有的.■例8矩阵乘法运算要求两矩阵的维数相容,否则会出错.先求两矩阵的乘积,若出错,则自动转去求两矩阵的点乘.function C=rc(A,B)%A,B为矩阵tryC=A*B;catchC=A.*B;endlasterr %显示出错原因3.请求键盘输入命令keyboard暂停运行程序并处于等待键盘输入状态(便于调试及程序运行中修改变量),键入“return”程序继续运行.如运行如下:4. echo 命令(在命令窗口,程序调用前使用)一般情况下,M 文件执行时,文件中的命令不会显示在命令窗口.echo 命令可使文件命令在执行时可见.这对程序的调试及演示很有用.对命令文件和函数文件,echo 的作用稍微有些不同.●对命令文件,其格式如下:echo on 打开命令文件的回应命令 echo off 关闭回应命令●对函数文件,其格式如下:echo file on 使指定的file 文件的命令在执行中被显示出来 echo file off 关闭指定文件的命令在执行中的回应function y=f(x) keyboardy=(5*x-6)/((x-3)*(x-3)+2);>> f K>> x=4; K>> return ans = 4.6667>> y=f(4) K>> return y = 4.6667。
北京理工大学汇编实验五实验报告概要

本科实验报告实验名称:子程序设计实验实验五子程序设计实验(设计性实验)一、实验要求和目的1.熟悉汇编语言程序设计结构;2.熟悉汇编语言子程序设计方法;3.熟悉利用汇编语言子程序参数传递方法;4.熟悉汇编语言字符串处理基本指令的使用方法;5.掌握利用汇编语言实现字符串的输入输出程序设计方法;6.掌握数制转换程序实现方法。
二、软硬件环境1、硬件环境:计算机系统windows;2、软件环境:装有MASM、DEBUG、LINK、等应用程序。
三、实验涉及的主要知识A)子程序知识要点:1、掌握子程序的定义语句;过程名PROC [near/far] 过程体RET过程名ENDP2.子程序结构形式一个完整的子程序一般应包含下列内容:1. )子程序的说明部分在设计了程序时,要建立子程序的文档说明,使用户能清楚此子程序的功能和调用方法.说明时,应含如下内容:.子程序名:命名时要名中见意..子程序的功能:说明子程序完成的任务;.子程序入口参数:说明子程序运行所需参数及存放位置;.子程序出口参数:说明子程序运行结果的参数及存放位置;.子程序所占用的寄存器和工作单元;.子程序调用示例;2、)掌握子程序的调用与返回在汇编语言中,子程序的调用用CALL,返回用RET指令来完成。
.段内调用与返回:调用子程序指令与子程序同在一个段内。
因此只修改IP;.段间调用与返回:调用子程序与子程序分别在不同的段,因此在返回时,需同时修改CS:IP。
3.)子程序的现场保护与恢复保护现场:在子程序设计时,CPU内部寄存器内容的保护和恢复。
一般利用堆栈实现现场保护和恢复的格式:过程名PROC[NEAR/FAR]PUSH AXPUSH BX..PUSH DX...POP DX...POP AXRET过程名ENDP4.子程序的参数传递方法1.寄存器传递参数这种方式是最基本的参数传递方式。
2.存储器单元传(变量)递参数这种方法是在主程序调用子程序前,将入口参数存放到约定的存储单元中;子程序运行时到约定存储位置读取参数;子程序执行结束后将结果也放在约定存储单元中。
C语言程序设计实验实验报告

C语⾔程序设计实验实验报告《C语⾔程序设计实验》实验报告实验九⽂件编程(⼀)班级:_07某某专业1班_学号:______姓名:___实验样例_____⽇期:___2008.6.10_____台州学院数信学院计算机系制作⼀、实验⽬的和要求综合使⽤基本结构、基本算法、函数、指针、结构体、链表等知识设计程序。
⼆、实验指导1、注意:本次实验为综合性实验,需要综合本课程以前所学内容,完成编制⼀个⼩型应⽤程序。
需要复习课本以前所学内容,浏览上课讲解的相关内容。
2、实验时间安排:准备实验⽅案,上课时上机编程,课余完成整个设计并完成实验报告。
3、实验具体要求:以⼀张假想的学⽣期中成绩表为实验对象,编制⼀个⼩型应⽤程序,可以考虑实现以下功能:(1)成绩修改;(2)成绩查询;(3)成绩统计;(4)成绩输⼊。
也可以实现其他⾃定义的功能。
界⾯采⽤简易菜单形式,可以分级显⽰。
实验需要进⾏以下步骤:需求分析、程序设计、程序调试、程序运⾏,并对程序运⾏效果进⾏分析,提出改进⽅案。
初步锻炼学⽣的程序设计能⼒、程序调试能⼒、对实验结果进⾏综合分析的能⼒。
4、实验成绩评价标准:参见成绩评定表。
四、实验报告内容(参见样本,请另附页)实验报告必须包含以下内容:1、对你设计的程序功能进⾏的较详细描述。
2、画出程序测试⽤的学⽣期中成绩表(包含:表头、若⼲⾏包含学⽣成绩等信息的数据⾏)。
3、给出程序清单(包含结构体的声明与定义、所有的函数声明、所有的全局变量、所有的函数定义、必要的注释。
注:如果程序过长,可以不⽤给出所有的函数定义,只给出主要的函数定义)。
4、运⾏结果(写出运⾏时的主界⾯、各模块的实现界⾯等,不作具体要求,但要求能清晰地体现出程序的功能)。
5、对你的设计作⼀个评价(可以包含你所花的时间、精⼒,你遇到的困难及如何解决的,你的设计过程所遇到的其他事情;必须包含你对设计的客观评价,以及提出的改进⽅案)。
1、程序功能描述:(略)2、成绩表样表:3、程序清单:(红⾊部分是和上⼀个程序的区别之处)#includestruct student /*期中成绩表的结构声明*/{int num;char name[20];char sex;float YuWen;float ShuXue;float YingYu;float totle;};#define MAX_STU_NUMBER 100 /*宏定义最⼤学⽣数为100⼈*/struct student Stu[MAX_STU_NUMBER]; /*⼀次性定义100个学⽣的存储空间*/int STU_NUMBER; /*存放学⽣数*/char *FileName = "chengji.txt"; /*成绩存放的⽂件*//*所有的函数声明:*/int SelectMainMenu(); /*显⽰并选择主菜单。
《Java程序设计》上机实验报告 实验五 流、文件及基于文本的应用

信息科学与工程学院《Java程序设计》上机实验报告专业班级姓名学号主要实验记录及个人小结(包括部分实验源程序、调试结果及实验结果分析等)2.编写一个程序,能将一个源程序中的空行及注释去掉,并将处理后的内容写到目标文件中。
(源文件名称和目标文件名称均要求以命令行参数形式输入)import java.io.BufferedReader;import java.io.File;import java.io.FileNotFoundException;import java.io.FileReader;import java.io.FileWriter;import java.io.IOException;import java.io.PrintWriter;public class CopyFile {static String delComm(String s){if(s==null)return s;int pos=s.indexOf("//"); //去掉以//开始的注释if(pos<0)return s;return s.substring(0,pos);}public static void main(String args[]){String infname="CopyFile.java";String outfname="CopyFile.txt";if(args.length>=1)infname=args[0];if(args.length>=2)outfname=args[1];try{File fin=new File(infname);File fout=new File(outfname);BufferedReader in = new BufferedReader(new FileReader(fin));PrintWriter out = new PrintWriter(new FileWriter(fout));int cnt=0;String s = in.readLine(); //写字符串while(s != null){主要实验记录及个人小结(包括部分实验源程序、调试结果及实验结果分析等)cnt ++;s=delComm(s); //去掉以//开始的注释out.println(cnt+":\t"+s);//写出s=in.readLine(); //读入}in.close(); //关闭缓冲读入流及文件读入流的连接out.close();}catch (FileNotFoundException e1){System.out.println("File not found");}catch(IOException e2){e2.printStackTrace();}}}主要问题的出错现象:1)运行程序后,无对应的输出结果出现2)源程序编写过程中,出现方法和类的调用无法完成.出错原因:1)函数调用时候,相关函数输入变量错误2)方法名的书写不符合规范,使得方法变为变量.解决方法及效果:1) 编写和调用类的方法时,注意参数的类型和数量的对应.2)规范命名.常量用大写字母,变量用小写字母开始,类以大写字母开始.实验小结:1) 为每个类分配正确的职责.尽量保持类功能单一. 不要在一个类中做太多事.2) BufferedReader的方法:readLine():String ,当他的返回值是null时,就表示读取完毕了。
C++语言程序设计形考二实验五

上机实验5:数组应用程序设计一、实验内容有如下程序:#include <iostream>using namespace std;int main(){__________________ //定义具有6个元素的整形数组dint data,i=0; // data用于暂存每一个输入数据,i记录已插入数据的个数。
cout<<"请输入6个整数:";do{ //输入并处理6个数据cin>>data;int j=i; //j是候选的插入位置while(j>0){ //若j==0,它就是第一个数据的插入位置,无需进入循环。
if(data>=d[j-1]) __________; //若找到插入位置,退出循环。
d[j]=d[j-1]; //否则数据后移一个位置j--;}d[j]=data; //数据插入到位置j处i++; //i记录已插入数据的个数}while(________); //若不足6个数据,继续输入。
for(_______________________) cout<<d[i]<<' ';//显示数组中的6个数据cout<<endl;return 0;}此程序输入6个整数保存到一维数组d 中,使得输入的数据在数组中始终按从小到大的顺序排列,最后顺序输出数组中的数据。
例如,若输入是5 7 4 29 3,程序的输出是:2 3 4 5 7 9。
二、实验要求读懂程序,补充程序中下划线处缺失的部分,上机测试,直至产生正确的输出结果。
三、实验结果截图(请将实际上机运行的结果截图并上传)。
Python语言程序设计实验

Python语言程序设计实验一、Python语言简介1.Python语言的发展历程2.Python语言的特点–解释型语言–拥有丰富的标准库3.Python语言的应用领域–Web开发–自动化运维二、Python程序设计基础1.变量和数据类型–浮点数类型–字符串类型2.控制结构–if-elif-else语句–for循环–while循环–函数的定义与调用–匿名函数(lambda)–函数的递归调用3.模块和包–模块的概念–模块的搜索路径4.文件操作–文件的打开与关闭–文件读写操作–文件常用方法三、Python高级特性1.面向对象编程–构造函数(__init__)–析构函数(__del__)2.异常处理–异常的概念–异常的捕获与处理–异常的传递–自定义异常3.列表推导式–列表推导式的概念–常用列表推导式–列表推导式的嵌套–生成器的概念–生成器的使用–生成器的迭代器–装饰器的概念–装饰器的使用–带参数的装饰器四、Python常用库•math库•datetime库•json库•urllib库•shutil库2.第三方库–requests库–numpy库–pandas库–matplotlib库–scrapy库–Django库五、Python实验教程1.Python编程环境搭建–安装Python–使用IDE(如PyCharm、VSCode)–安装与使用文本编辑器(如Sublime Text、Atom)2.Python编程规范–代码格式化–注释与文档3.编写与运行Python程序–编写Python源代码–保存源代码–运行Python程序4.Python编程实践–计算器程序–学生信息管理系统–数据分析与可视化–网页自动化测试5.Python项目实战–选择合适的项目–项目需求分析–项目设计与开发–项目测试与部署六、实验要求与评价1.实验报告2.代码质量–代码可读性3.实验评价通过以上知识点的学习与实践,学生可以掌握Python语言的基本语法、编程技巧和项目开发流程,为以后的学习和工作打下坚实的基础。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
实验5. 综合应用程序设计
一、 实验目的
1.理解命令.m 文件与函数.m 文件的异同;
2.掌握定义和调用matlab 函数的方法。
二、 实验环境
1. 计算机
2. matlab7.1
三、 实验说明
3. 熟悉matlab7.1运行环境,正确操作,认真观察;
4. 实验学时:2学时;
5. 学会使用help ;
6. 编写程序前先写出程序的流程图框架;
7. 自主编写调试程序直到正确为止,必要时参考相关资料,; 8. 保存整理实验结果,提交实验报告。
四、 实验内容
1. 设两个系统各自有5个冗余件可供选择,可靠性与成本如下表(数据来自上海大众4S 店)所示。
则对成本与可靠性的集成研究可化为以下两个系统中四个冗余件的优化组合。
(注:尽量不用循环,用矩阵结构,保存文件名为s5_学号_01.m )
图 调角器系统与电机系统可靠性的并串联结构
表 调角器系统与电机系统可靠度及成本数据
R 1i
R 1i ’ R 2j ’
R 2j
目标:
Min C = C1 + C2 约束条件:
C1 = 65.470×
()
2
1i
R
i
∑
(i取1,2,3,4,5中任两个数)
C2 = 119.281×
()
2
2j
R
j
∑
(j 取1,2,3,4,5中任两个数)
R1 =1-
2
1i
(1R)
i
∏-
(i取1,2,3,4,5中任两个数)
R2 =1-
2
2j
(1R)
j
∏-
(j取1,2,3,4,5中任两个数)
R1 ≥ 0.9825 , R2 ≥ 0.9819
五、实验程序及结果
(1)程序:
R1i=[0.961 0.920 0.876 0.792 0.723];
R2j=[0.951 0.820 0.788 0.712 0.657];
A=[R1i(1)+R1i;R1i(2)+R1i;R1i(3)+R1i;R1i(4)+R1i;R1i(5)+R1i];%两个R1i相加
B=[R2j(1)+R2j;R2j(2)+R2j;R2j(3)+R2j;R2j(4)+R2j;R2j(5)+R2j];%两个R2j相加
C=(1-R1i)'*(1-R1i);%两个(1-R1i)相乘
D=(1-R2j)'*(1-R2j); %两个(1-R2j)相乘
R1=1-C;%R1的值
R2=1-D;%R2的值
C1=65.470*(A.*(R1>=0.9825).*(R2>=0.9819)); %在满足(R1>=0.9825)&(R2>=0.9819)的条件下求C1,C2
C2=119.281*(B.*(R1>=0.9825).*(R2>=0.9819));
Cm=C1+C2;
C=Cm-diag(diag(Cm)) %令Cm中的对角线上元素值为0,即排除选择的是同一个冗余件的情况。
%此处Cm就是C=C1+C2的所有可能情况
C=C(find(C~=0));%把Cm中的非零元素提取出来,组成新的矩阵Cm
Min=min(C)
(2)结果:
C =
0 334.3957 327.6980 313.1332 302.0553
334.3957 0 0 0 0
327.6980 0 0 0 0
313.1332 0 0 0 0
302.0553 0 0 0 0
Min =
302.0553
2.定义一个函数文件,功能是采用重叠相加法完成无限长序列与有限长序列的线性卷积,然后用该函数文件计算序列x[k]=2k+1,0≤k≤18与序列h[k]={1,2,3,4;k=0,1,2,3}的线性卷积结果。
(注:可能用到的函数conv,保存文件名为s5_学号_02.m)
(1)函数文件
%% 重叠相加法,实现长短序列的线性卷积
function y=overlap_add(x,h,M)
% x为较长的输入序列,h为较短的系统函数,M为分段大小
N=length(h); % 序列h(n)的长度
Lx=length(x); % 序列x(n)的长度
if N>M % 算法要求 N<=M
M=N+1;
end
L=M+N-1; % 用每段循环卷积计算线性卷积所需点数
T=ceil(Lx/M); % 分段数,ceil向上取整
t=zeros(1,N-1); % 缓存序列初始化
x=[x,zeros(1,(T+1)*M-Lx)]; % 最后一个不足M的分段补零
y=zeros(1,(T+1)*M); % 生成输出序列y(n)
for i=0:T
xi=i*M+1;
x_seg=x(xi:xi+M-1); % 低点数计算时的分段x(n)
y_seg=lin_conv(x_seg,h,L); % 计算 i 分段和 h 的循环线性卷积
y_seg(1:N-1)=y_seg(1:N-1)+t(1:N-1); % 与前一段卷积的后N-1位重叠相加
t(1:N-1)=y_seg(M+1:L); % 缓存序列更新
y(xi:xi+M-1)=y_seg(1:M); % 每卷积一段输出M个点
end
y=y(1:Lx+N-1) ; % 最终输出序列
end
%% 利用FFT循环卷积计算线性卷积的函数,可放在另一个m文件中被调用
function y=lin_conv(x1,x2,L)
% 利用FFT循环卷积计算线性卷积
y1=fft(x1,L);
y2=fft(x2,L);
Yk=y1.*y2;
y=ifft(Yk);
end
(2)程序
k=0:18;
x=2*k+1
h=[1,2,3,4]
M=input('请输入将x分段后各段的长度M:')
foverlap=overlap_add(x,h,M)
(3)结果
x =
Columns 1 through 18
1 3 5 7 9 11 13 15 17 19 21 23 25 27 29 31 33 35
Column 19
37
h =
1 2 3 4
请输入将x分段后各段的长度M:7
M =
7
xf =
1 3 5 7 9 11 13
15 17 19 21 23 25 27
29 31 33 35 37 0 0
y =
1 5 14 30 50 70 90 95 83 52
15 47 98 170 190 210 230 221 181 108
29 89 182 310 330 311 251 148 0 0
f =
Columns 1 through 18
1 5 14 30 50 70 90 95 83 5
2 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 15 47 98 170 190 210 230 221 181 108 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 29 89 182 310
Columns 19 through 22
0 0 0 0
0 0 0 0
330 311 251 148
foverlap =
Columns 1 through 18
1 5 14 30 50 70 90 110 130 150 170 190 210 230 250 270 290 310
Columns 19 through 22
330 311 251 148。