实验四循环结构设计实现 一、实验目的 1. 掌握while,do-while循环语句的使用与区别; 2. 掌握循环条件、循环体、循环终止等循环要素, 3. 理解循环的执行过程。 二、实验内容 1、上机验证课堂实例：求1+2+3+…+100的和。要求分别用while和do-while和for语句实现。 2、打印出所有的“水仙花数”。 算法提示： 1）水仙花数是一个3位数，被判断的范围是[100，999]。 2）如何分离出百、十、个位数。可参考教材P115例5.16。 3、验证课堂例题：打印“九九乘法表”。 4、公元钱五世纪，我国古代数学家张丘建在《算经》一书中提出了“百鸡问题”：鸡翁一值钱五，鸡母一值钱三，鸡雏三值钱一。百钱买百鸡，问鸡翁、鸡母、鸡雏各几何？ 数学模型：设i代表母鸡数，j代表公鸡数，k代表小鸡数。用凑数法解决问题。 i+j+k=100 (1) 5*i+3*j+k\3=100 (2) 算法提示： S1． i从1到20依次一一取值； S2．对每一个固定的i值，j都要从1到33依次一一取一遍值； S3．对每一个固定的i值及每一个固定的j值，按公式k=100-i-j取得k值； S4．做下列验证工作：S3中所取的一组i，j，k是否满足条件(2) 5*i+3*j+1/3*k=100？ 若满足，则输出这组解i，j，k，然后转2；否则直接转2。 当i已取到20，j也取到33时整个任务就结束。 三、仪器、设备、材料 微机 四、实验准备 1．理论知识预习及要求 ①表达式的书写规范； ②while语句和do-while语句的格式、功能、执行过程； ③while语句和do-while语句间的区别； ④循环结构程序设计的简单算法。 2．实验指导书预习及要求 对程序执行流程不清楚时,请先画出程序的流程图. 上机前先思考编出程序,并分析结果，上机时调试,并写出实验结果。 3．其他准备
实验6 循环结构程序设计 一．实验目的 1．掌握在设计条件型循环结构时，如何正确地设定循环条件。 2．掌握如何正确地控制计数型循环结构的循环次数。 3．练习并掌握选择结构与循环结构的嵌套、多重循环的应用。 4．掌握在程序设计中用循环的方法实现一些常用算法，加强调试程序的能力。二．实验要求 1．复习while、do-while、for语句和continue、break语句。 2．能够用流程图表示实验题目的算法。 3．能够独立调试运行实验题目。 4．本实验要求4学时完成。 三．实验内容和步骤 题目1：分别用while循环和for循环计算：1+2+3+…+100=？ 要求：写出程序并上机查看运行结果。 题目2：输入一行字符,分别统计出其中大小写英文字母、空格、数字和其他字母的个数。要求：程序填空并写出2组运行结果。 01#include
matlab 基本语句 1.循环语句for for i=s1:s3：s2 循环语句组 end 解释：首先给i赋值s1；然后，判断i是否介于s1与s2之间；如果是，则执行循环语句组,i=i+s3（否则，退出循环.）；执行完毕后，继续下一次循环。 例：求1到100的和，可以编程如下： sum=0 for i=1:1:100 sum=sum+i end 这个程序也可以用while语句编程。 注：for循环可以通过break语句结束整个for循环. 2.循环语句while 例：sum=0;i=1; while(i<=100) sum=sum+i;i=i+1; end 3.if语句 if(条件) 语句 end if(条件) 语句 else 语句 end if(条件) 语句 elseif 语句 end 4.关系表达式：
=,>,<,>=,<=,==(精确等于) 5.逻辑表达式：|(或),&(且) 6.[n,m]=size(A)(A为矩阵) 这样可以得到矩阵A的行和列数 n=length(A)，可以得到向量A的分量个数；如果是矩阵，则得到矩阵A的行与列数这两个数字中的最大值。 7.！后面接Dos命令可以调用运行一个dos程序。 8.常见函数： poly():为求矩阵的特征多项式的函数，得到的为特征多项式的各个系数。如 a=[1,0,0;0,2,0;0,0,3]，则poly(a)＝1 -6 11 -6。相当于poly(a)＝1入^3+(-6)入^2+11入+(-6)。 compan():可以求矩阵的伴随矩阵. sin()等三角函数。 MATLAB在数学建模中的应用(3) 一、程序设计概述 MATLAB所提供的程序设计语言是一种被称为第四代编程语言的高级程序设计语言，其程序简洁，可读性很强，容易调试。同时，MATLAB的编程效率比C/C++语言要高得多。 MATLAB编程环境有很多。常用的有： 1.命令窗口 2.word窗口 3.M-文件编辑器，这是最好的编程环境。 M-文件的扩展名为“.m”。M-文件的格式分为两种： ①λ M-脚本文件，也可称为“命令文件”。 ② M-函数文件。这是matlab程序设计的主流。λ 保存后的文件可以随时调用。 二、MATLAB程序结构 按照现代程序设计的观点，任何算法功能都可以通过三种基本程序结构来实现，这三种结构是：顺序结构、选择结构和循环结构。其中顺序结构是最基本的结构，它依照语句的自然顺序逐条地执行程序的各条语句。如果要根据输入数据的实际情况进行逻辑判断，对不同的结果进行不同的处理，可以使用选择结构。如果需要反复执行某些程序段落，可以使用循环结构。 1 顺序结构 顺序结构是由两个程序模块串接构成。一个程序模块是完成一项独立功能的逻辑单元，它可以是一段程序、一个函数，或者是一条语句。 看图可知，在顺序结构中，这两个程序模块是顺序执行的，即先执行<程序
一、实验目的 1.利用MATLAB 实现循环卷积。 2.比较循环卷积与线性卷积的区别。 二、实验条件 PC 机，MATLAB7.0 三、实验内容 1）循环卷积的定义：两个序列的N 点循环卷积定义为： )0()()()]()([1 0N n m n x m h n x n h N k N N <≤-=?∑-= 利用MATLAB 实现两个序列的循环卷积可以分三个步骤完成： (1)初始化：确定循环点数N ，测量输入2个序列的长度。 (2)循环右移函数：将序列x(n)循环右移，一共移N 次（N 为循环卷积的循环次数），最后将每次循环成的新序列组成一个矩阵V 。 (3)相乘：将x(n)移位后组成的矩阵V 与第二个序列h(n)对应相乘，即得循环卷积结果。程序如下： 程序一： clear;close all ; N=10; x1=[6 15 -6 3 5 7 0 1]; x2=[7 1 2 9 4 3 20 6]; xn1=length(x1); xxn1=0:xn1-1; xn2=length(x2); xxn2=0:xn2-1; subplot(3,1,1); stem(xxn1,x1); subplot(3,1,2); stem(xxn2,x2); x11=fft(x1,N);
x12=fft(x2,N); y11=x11.*x12; y1=ifft(y11,N); subplot(3,1,3); n=0:length(y1)-1; stem(n,y1,'.'); title('循环卷积的结果'); xlabel('n');ylabel('y1(n)'); 运行后所得图形如下： 观察所得的循环卷积结果发现并没有呈现周期性的序列，因此将程序做下列改变。程序二： clear;close all; N=40; x1=[6 15 -6 3 5 7 0 1]; x2=[7 1 2 9 4 3 20 6]; x2=[x2,x2,x2,x2]; xn1=length(x1); xxn1=0:xn1-1; xn2=length(x2); xxn2=0:xn2-1; subplot(3,1,1);
嘉应学院计算机学院 实验报告 课程名称程序设计基础实验名称实验地点 指导老师实验时间提交时间 班级姓名座号 一、实验目的和要求 （1）熟悉掌握用while语句、do…while语句和for语句实现循环的方法。 （2）掌握在程序设计中用循环的方法实现一些常用算法（如穷举、迭代、递推等）。 （3）进一步学习调试程序。 二、实验环境和方法 实验方法： （一）综合运用课本所学的知识，用不同的算法实现在不同的程序功能。 （二）结合指导老师的指导，解决程序中的问题，正确解决实际中存在的异常情况，逐步改善功能。 （三）根据实验内容，编译程序。 实验环境：Windows xp Visual C++6.0 三、实验内容及过程描述 实验步骤： ①进入Visual C++ 6.0集成环境。 ②输入自己编好的程序。 ③检查一遍已输入的程序是否有错（包括输入时输错的和编程中的错误），如发现有错， 及时改正。 ④进行编译和连接。如果在编译和连接过程中发现错误，频幕上会出现“报错信息”， 根据提示找到出错位置和原因，加以改正。再进行编译，如此反复直到不出错为止。 ⑤运行程序并分析运行结果是否合理。在运行是要注意当输入不同的数据时所得结果 是否正确，应运行多次，分别检查在不同情况下结果是否正确。 实验内容：编译以下题目的程序并调试运行。 （1）输入一行字符，分别统计出其中的英文字母、空格、数字和其它字符的个数。 编写程序如下： #include
循环结构：for语句 格式： for 循环变量=表达式1：表达式2：表达式3 循环体 end 【注】：表达式1：循环变量初值， 表达式2：步长，为1时，可省略； 表达式3：循环变量终值。 或： for循环变量=矩阵表达式 循环体 end 【注】：执行过程是依次将矩阵的各列元素赋给循环变量，然后执行循环体语句，直至各列元素处理完毕。 2 while语句： 格式： while（条件） 循环体 end 【注】：条件成立时，执行循环体 3
break语句&& continue语句： break：破坏，破坏循环，终止循环的进行，跳出循环，程序将执行循环语句的下一语句。 continue：继续，循环继续，程序将跳过循环体中剩下的语句，继续下一次循环。 4 循环的嵌套—多重循环结构 5 选择结构：if-else语句 格式： if 表达式 程序模块 end 或 if 表达式 程序模块1 else 程序模块2 end 6 switch语句： 格式：
switch 表达式 case 数值1 程序模块1 case 数值2 程序模块2 case 数值3 程序模块3 ...... otherwise 程序模块n end 执行过程：首先计算表达式的值， 然后将其结果与每一个case后面的数值依次进行比较， 如果相等，则执行该case的程序模块； 如果都不相等，则执行otherwise模块中的语句。 switch语句可以替代多分支的if语句，而且switch语句简洁明了，可读性更好。 7 matlab中一些基本知识： END 注意事项 for循环可以通过break语句结束整个for循环
实验报告 班级：电信13-1班学号：130******** 姓名：谢朗星成绩： 实验四循环结构 一、实验目的 1．掌握循环结构C程序的编写和调试方法。 2．掌握循环结构C程序中使用while语句和do while语句。 二、实验内容和步骤 1.调试下面三个程序，分析执行的结果。 1. #include
} 在该程序中，先判断i是否<=100，如果是，在循环体先执行了sum=sum+i，接着再执行i++，然后输出当次运行的结果，然后再返回判断i是否<=100，接着再重复上步骤，直到i>100，结束该程序。 3.#include
．实验五循环结构程序设计实验内容 （1）编写程序实现求两矩阵A*B,A.*B. A*B >> A=input('输入矩阵A:'); B=input('输入矩阵B:'); C=zeros(size(A,1),size(B,2)); for ii=1:size(A,1) for jj=1:size(B,2) for k=1:size(A,2) C(ii,jj)=C(ii,jj)+A(ii,k)*B(k,jj); end end end C 输入矩阵A:[1 2;3 4] 输入矩阵B:[5 6;7 8] C = 19 22 43 50 A.*B >> A=input('请输入矩阵A:'); B=input('请输入矩阵B:'); if size(A)==size(B) [m,n]=size(A); for i=1:m for j=1:n C(i,j)=A(i,j)*B(i,j); end end disp('输出A.*B:') C
else disp('无效') end 请输入矩阵A:[1 2 3] 请输入矩阵B:[4 5 6] 输出A.*B: C = 4 10 18 （2）编写程序实现求矩阵的转置。 >> A=input('请输入矩阵A:'); [m,n]=size(A); B=zeros(n,m); for i=1:m; for j=1:n; B(i,j)=A(j,i); end end disp(B) 请输入矩阵A:[1 2 3;4 5 6;7 8 9] 1 4 7 2 5 8 3 6 9 小结实验心得体会 通过本次实验，我慢慢理解了循环控制结构的执行过程，掌握了一些循环结构程序设计方法，但还不是很熟练。我了解了一些MATLAB的编程方法和循环结构程序设计方法，但都不是特别的熟练，都有待加强。
实验三循环结构 一．实验目的 1．熟悉用while语句，do-while语句和for语句实现循环的方法。 2．掌握循环嵌套的用法。 3．掌握在程序设计中用循环的方法实现各种算法（如迭代，递推，穷举等）。4．掌握break和continue语句的用法。 二．实验要点 在C语言中主要用以下方法实现循环： 1.while语句格式：while （表达式）语句 2.do- while语句格式：do 语句 while (表达式)； 3.for语句 for语句的一般格式： for （表达式1；表达式2；表达式3）语句 for (初值表达式；条件表达式；增值表达式) 语句 三．实验要求 1．上机前编写好以下程序。 2．上机输入和调试自己所编的程序并存在自己的U盘上。 3．检查实验结果是否正确。 4．上机结束后，整理实验报告，以“学号—姓名-第3次实验”命名的文件上传到作业服务器。 四．实验内容 上机调试运行以下程序。 1．下面的程序段是从键盘输入的字符中统计数字字符的个数,用换行符结束循环.请填空. #include”stdio.h” void main(){ int n=0; char c; c=getchar(); while (【】) { if ( 【】) n++; c=getchar(); } printf(“n=%d”,n);}
键盘输入：abcdef12345k
实验四循环结构程序设计 一、实验目的 1．掌握利用for语句实现循环结构的方法。 2．掌握利用while语句实现循环结构的方法。 3．熟悉利用向量运算来代替循环操作的方法。 二、实验内容 1. 根据pi*pi/6=1/1^2+1/2^2+1/3^2+……+1/n^2,求pi的近似值。当n分别 取100，1000，10000时，结果是多少？ 2. 根据y=1+1/3+1/5+……+1/（2n-1），求 （1）y<3时的最大n值 （2）与（1）的n值对应的y值 4. 已知{f1=1，n=1；f2=0，n=2；f3=1，n=3；fn=fn-1-2fn-2+fn-3，n>3 求f1-f100中 （1）最大值，最小值，各数之和。 （2）正数，零、负数的个数。 三、实验步骤(程序和结果) 1. 方法1：y=0;n=100; >> for i=1:n y=y+1/i/i; end >> y y = 1.6350 >> pi=sqrt(6*y) pi = 3.1321 y=0;n=1000; for i=1:n y=y+1/i/i; end >> y
y = 1.6439 pi=sqrt(6*y) pi = 3.1406 y=0;n=10000; >> for i=1:n y=y+1/i/i; end >> y y = 1.6448 pi=sqrt(6*y) pi = 3.1415 方法2：n=100;i=1:n; >> f=1./i.^2; >> y=sum(f) y = 1.6350 >> pi=sqrt(6*y) pi = 3.1321 n=1000;i=1:n; >> f=1./i.^2; >> y=sum(f) y =
实验五循环结构程序设计 【目的与要求】 熟练掌握while、do—while和for三种循环语句的应用。 【上机内容】 【示例】只要求运行前4个例题，并在/* */中说明该语句的作用或使用的方法等；典型例题仅供参考，不作要求。 【例1】求5！ main() { int n,t; n=1; t=1; while(t<=5) { n=n*t; t=t+1;/**/ } printf(“%d”,n); } 【例2】求和s=1!+2!+3! main() { int n,s=0,t=1; for(n=1;n<=3;n++) { t=t*n;/*计算阶乘*/ s=s+t;/* */ } printf(“%d”,s); } 【例3】求和s= 3+33+333 main() { int s=0,t=3,i; for(i=1;i<=3;i++) { s=s+t; t=10*t+3;/*生成3,33,333等数据*/ } printf(“%d”,s); } 【例4】打印500以内的“水仙花数”。“水仙花数”是一个三位数，其各位数立方和等于该数本身。main()
int i,j,k,n; for(n=100;n<=500;n++) { i=n/100;/* */ j=n/10-i*10; /* */ k=n%10; /* */ if(n= =i*i*i+j*j*j+k*k*k) /* */ printf(“ %d”,n); } } 【典型示例】 【例1】计算两个数的最大公约数；分析：首先，随机输入两个数m，n(默认m>n)；其次，算法：使k为m除以n的余数，如果m能被n整除，则k值为0，n为这两个数的最大公约数，否则，使k代替n，n代替m，重复以上过程，直到k值为0。 #include"stdio.h" main() { int m,n,k,result; printf("Enter two numbers:"); scanf("%d,%d",&m,&n); if(m>0&&n>0) /*限定两个正整数*/ { do { k= m % n; if(k= =0) result=n; else { m = n; n=k; } }while(k>0); /*循环取余求出最大公因子*/ printf("The greatest common divistor is:%d\n",result); } else printf("Nonpositive values not allowed\n"); } 【例2】打印下列图案 * *** ***** ******* main()
Matlab 基本语句 1.循环语句for for i=s1:s3：s2 循环语句组 end 解释：首先给i赋值s1；然后，判断i是否介于s1与s2之间；如果是，则执行循环语句组,i=i+s3（否则，退出循环.）；执行完毕后，继续下一次循环。 例：求1到100的和，可以编程如下： sum=0 for i=1:1:100 sum=sum+i end 这个程序也可以用while语句编程。 注：for循环可以通过break语句结束整个for循环. 2.循环语句while 例：sum=0;i=1; while(i<=100) sum=sum+i;i=i+1; end 3.if语句 if(条件) 语句 end if(条件) 语句 else 语句 end if(条件) 语句 elseif 语句 end 4.关系表达式：
=,>,<,>=,<=,==(精确等于) 5.逻辑表达式：|(或),&(且) 6.[n,m]=size(A)(A为矩阵) 这样可以得到矩阵A的行和列数 n=length(A)，可以得到向量A的分量个数；如果是矩阵，则得到矩阵A的行与列数这两个数字中的最大值。 7.！后面接Dos命令可以调用运行一个dos程序。 8.常见函数： poly():为求矩阵的特征多项式的函数，得到的为特征多项式的各个系数。如 a=[1,0,0;0,2,0;0,0,3]，则poly(a)＝1 -6 11 -6。相当于poly(a)＝1入^3+(-6)入^2+11入+(-6)。 compan():可以求矩阵的伴随矩阵. sin()等三角函数。 MATLAB在数学建模中的应用(3) 一、程序设计概述 MATLAB所提供的程序设计语言是一种被称为第四代编程语言的高级程序设计语言，其程序简洁，可读性很强，容易调试。同时，MATLAB的编程效率比C/C++语言要高得多。 MATLAB编程环境有很多。常用的有： 1.命令窗口 2.word窗口 3.M-文件编辑器，这是最好的编程环境。 M-文件的扩展名为“.m”。M-文件的格式分为两种： ①λ M-脚本文件，也可称为“命令文件”。 ② M-函数文件。这是matlab程序设计的主流。λ 保存后的文件可以随时调用。 二、MATLAB程序结构 按照现代程序设计的观点，任何算法功能都可以通过三种基本程序结构来实现，这三种结构是：顺序结构、选择结构和循环结构。其中顺序结构是最基本的结构，它依照语句的自然顺序逐条地执行程序的各条语句。如果要根据输入数据的实际情况进行逻辑判断，对不同的结果进行不同的处理，可以使用选择结构。如果需要反复执行某些程序段落，可以使用循环结构。 1 顺序结构 顺序结构是由两个程序模块串接构成。一个程序模块是完成一项独立功能的逻辑单元，它可以是一段程序、一个函数，或者是一条语句。 看图可知，在顺序结构中，这两个程序模块是顺序执行的，即先执行<程序
while语法： while expression statements end 说明：while expression, statements, end 计算一个表达式，并在该表达式为true 时在一个循环中重复执行一组语句。表达式的结果非空并且仅包含非零元素（逻辑值或实数值）时，该表达式为true。否则，表达式为false。 示例代码如下： function [sum] = summation(ratio, head, top) sum = 0; while (head <= top) sum = sum + ratio ^ head; head = head + 1; end end
假设ratio = 2，head = 0，top = 63 matlab循环语句for怎么用？ matlab中for语句使用方法和应用实例 for循环语句 1、一般格式为： for x（循环变量）= array（数组） commands（执行的循环代码） end 2、array可以是一个数字，也可以是数组，例如输入：for a=5 for a=1:5 for a=1:1:5（以1为步长到5）
只不过在a=1:5和a=1:1:5时，会显示之间的结果，a=5时只显示最后结果。a的变动就是第一次循环a=1，第二次循环a=2，第三次循环a=3，第四次循环a=4，第五次循环a=5。 3、(commands)就是命令，其中的命令行可以很多很多，最常见的就是调用上面说的a变动比如： for a=10 %循环10次 s=a+1 %循环语句 end %结束 上面的a不需要再指定，a的变化就是上面说的先是1，然后是2，3，…那么s就是先s=1+1=2，然后再s=2+1=3，s=3+1=4，…，s=10+1=11，循环结束，就是一个连续加s的指令，最后=11。 4、for语句可以嵌套的，和C一样 for a=5 %第一个for循环
C语言程序设计实验教学（5） 【实验目的】通过程序设计实现，掌握while结构、do-while结构和for结构各种类型的循环结构，完成各种循环程序的设计和实现。 【实验要求】同一个题目尝试使用多种循环控制结构来实现，并分析其设计和实现的差别和难易程度。 【实验课时】6.0 【实验内容】 1、完成课堂实例的实现。 （1）计算s＝1+2+…+100。（累加型） （2）请输入数n，计算n！。（连乘型） （3）输出所有的水仙花数。（范围型） （4）判断n是否为素数。（反向思维型） （5）输出100-200之间能同时被3和5整除的数。(用continue语句) （6）输入任意n个数，输出其中最大数和最小数，并输出它们在序列中的位置。 （7）输出以下图形（1）： （8）输出以下图形（2）： （9）输出以下图形（3）： 图形（1）图形（2）图形（3） （10）编写输出如下字母塔的程序： A A B A A B C B A A B C D C B A
……………… …………………… ………………………… A B C D E F G H I J K L M N O P Q R S T U V W X W V U T S R Q P O N M L K J I H G F E D C B A A B C D E F G H I J K L M N O P Q R S T U V W X Y X W V U T S R Q P O N M L K J I H G F E D C B A A B C D E F G H I J K L M N O P Q R S T U V W X Y Z Y X W V U T S R Q P O N M L K J I H G F E D C B A 2、计算多项式7+72+73+······+7n 的值，n 由键盘输入。 3、求 ∑=201! n n （即求1！+2！+3！+……+20！） 4、输入2个正整数m,n ，求其最大公约数和最小公倍数。 5、请编程序打印出以下数列：1、1、2、3、5、8、13、…的前40项。每行输出4个数。（斐波那契数列） 6、求100~200间的全部素数。每行输出6个数。 7、“鸡兔同笼”是一个中国古算题，最早出现在《孙子算经》中。笼中有若干只鸡和兔子，它们共有88个头，244只脚，问鸡和兔各有多少只？ 8、求下列分数序列的前20项之和。 12，23，35，58，813，1321，…… 9、输入x 值，按照公式计算cos(x)的值，直到最后一项小于10-6为止。 ,! 6x !4x 2!x 1)x (cos 6 42 +-+-= 10、找出1000之内的所有完数，一个数如果恰好等于它的因子之和，这个数就称为“完数”。例如，6的因子为1、2、3，而6=1+2+3，因此6是“完数”。输出时按下面的格式输出其因子： 6 its factors are 1,2,3 补充作业：（以下题目任选2题） 1、我国古代有一道有名的数学难题，称为“百鸡问题”：鸡翁一，值钱五；鸡母一，值钱三；鸡雏三，值钱一。百钱买了百鸡，问鸡翁、鸡母、鸡雏各买了多少只？ 2、猴子吃桃问题。猴子第一天摘下若干桃子，当即吃了一半，还不过瘾，又多吃了一个。第二天早上又将剩下的桃子吃掉一半，又多吃了一个。以后每天早上都吃了前一天剩下的一半零一个。到第10天早上再想吃时，就
MATLAB编程如何不用for循环--以DFT变换矩阵为例 缘起：大家都知道MATLAB中用for循环编写的代码执行起来效率不高，如何用矩阵和向量的运算提升效率对每一个用matlab的人来说都是很有必要的，但是此项功夫高手一般不愿意给初学者讲，此功夫是高手和低手的分水岭，高手们更是拿此功夫在初学者面前炫耀。本人当初怀着很恭敬的心向高手请教，高手笑笑说这要我自己编。出于让后来人受益，帮助和我一样无助的求知者。本人今天话了一天时间将此问题研究下，并且将代码毫无保留的公布出来。希望大家能够受益，阿弥陀佛！ 上述W矩阵的第一列代表直流成分，第二列到最后一列是信号的交流成分，可以看出倍频关系！我以前不知道DFT 可以通过矩阵表示。注意matlab中dftmtx实现上述W矩阵的时候没有用1/sqrt(N) 进行归一化！可以通过dftmtx(2)验证，没有1/sqrt(2)。 例1 DFT matrix 是Hermitian的 Nfft = 8; xn=rand(1,Nfft); y=dftmtx(Nfft)*xn.' %结果是个列向量y=dftmtx(Nfft)*x.' 和y=fft(x,Nfft)是等价的 y = xn*dftmtx(Nfft) %结果是个行向量 y=fft(xn,Nfft) %结果是个行向量 %dft变换公式，n代表时域采样点，k代表频域采样点 Y(k)=sum(x.*exp(-j*2*pi*n*k)) 相应的，dftmtx(Nfft)产生的矩阵中，第k行，n列元素=exp(-j*2*pi*k*n/Nfft)，与x.'相乘正好对应fft变换后的每个频点值。
例1：双重循环求DFT N = 8; x=rand(1,N); for k=0:N-1 sum=0; % 注意每个X(k)的值不应该受上次计算的影响 for m=0:N-1 w=exp(-j*2*pi/N.*m.*k); %DFT matrix 的每一行的元素是不同的 sum=sum+x(m+1).*w; % 这个循环实际上是计算DFT matrix的每一行与信号x的内积 end X(k+1)=sum; %Matlab下标从1开始 end % 注意X 是个行向量，是个数组，我以前不知道X(k)的循环赋值的结果是个行向量 y=fft(x,N) %验证 例2：单循环求DFT %% x 是行向量的Version，结果X也是行向量 N = 8; x=rand(1,N); m=0:N-1; sum=0; for k=0:N-1 %核心思想是内积运算：x躺着，后面的必须站着 X(k+1)= x*exp(-j*2*pi/N.*m'.*k) ; % 以前这样的表达物理含义是不明确的--> X(k+1)= x*exp(-j*2*pi/N.*m.*k).' ; end y=fft(x,N) %验证 %% x 是列向量的Version，结果X是行向量，fft(x)的结果是列向量 N = 8; x=rand(1,N)'; m=0:N-1; sum=0; for k=0:N-1 % 核心思想是内积运算：x躺着，后面的必须站着 X(k+1)= exp(-j*2*pi/N.*m.*k)*x ; % 以前这样的表达物理含义是不明确的--> X(k+1)= x*exp(-j*2*pi/N.*m.*k).' ; end y=fft(x,N) %验证结果y是个列向量，X是y的转置。因为X按下标赋值的结果是个行向量！ 例3：不用循环求DFT %% x 是行向量的Version，结果X也是行向量 N = 8; xn=rand(1,N); %一次去掉2个循环，不要试图一次去掉一个 n=[0:1:N-1]; k=[0:1:N-1]; WN=exp(-j*2*pi/N); nk=n'*k; %看看如何不用for循环来实现，向量的外积N-by-1 multiply 1-by-N WNnk=WN.^nk; %WNnk就是dftmtx(8)，dftmtx原来是这样构造的 Xk=xn*WNnk; % 究竟nk这个向量外径是n'*k还是k'*n要看内积在哪一维执行
for 循环语句 for语句作用是按照预先设定的循环步骤重复执行某语句段, 其语法: for循环控制变量=存储着该变量依次所取值的一个向量 循环语句块, 本块重复执行的次数由上方向量的长度决定, 每次执行, 循环控制变量依次取该向量中的值. end 注意, 如果你matlab的for 循环语句有其独有的特点, 它使用一个向量来控制循环, 循环次数由向量的长度来决定, 而每次循环都依次从向量中取值. 这使得Matlab循环更灵活多样, 其循环变量取值可以不按照特定的规律; 但是另一方面, Matlab 的for循环也有独特的限制, 当次循环中改变循环变量赋值, 不会代入下次循环, 所以除非在其中用break提前退出, Matlab循环的次数是预先定好的. 举个例子: for a=1:2:7 , 将循环4次, a的取值依次是1, 3, 5, 7. 另一个例子: for a=[1, 5, 3, 4, 6], 这个循环将被执行5次, 循环控制变量a的取值依次为: 1, 5, 3, 4, 6. 如果还对for循环这个特性不太熟悉, 试试下面这段Matlab代码: for a=[1 5 3 4 6] disp(['第', num2str(find(a==[1 5 3 4 6])) , '次循环, a的取值为: ',num2str(a)] ); end While循环 与for循环不同, while循环不预先指定好循环次数, 只要符合条件循环就一直执行下去, while语句的语法: while判断条件 循环语句块 end 此处的'判断条件'和if语句中的那个一样, 要求这个变量或者表达式最终得到一个逻辑型标量, 每次循环之前, while语句会判断这个条件是否满足, 如果满足则开始循环模块, 否则跳过整个循环语句. 在循环语句块中控制循环退出有两种办法, 其一是直接或间接地改变'判断条件'的值, 使之为'false', 其二是在循环块中执行'break' 语句直接退出循环. 一个例子, 每次循环变量a都将增加1, 我们欲控制使a大于10时终止循环, 下面分别采用上述提到的两种方法控制循环. approach 1 a=0; while a<=10 a=a+1; disp(a); end approach 2: a=0; while 1 %由于判断条件是'1', 永远为'true', 所以如果不在循环块中设置跳出条件,循环将永久性进行下去! a=a+1; disp(a); if a>=10 break; end
读书破万卷下笔如有神 上机实验五循环结构程序设计 姓名：学号: 一．目的要求 1．掌握用for、while、do-while语句实现循环的基本技巧； 2．掌握循环次数不确定的解题技巧； 3．掌握break语句、continue语句在循环语句中的作用以及不同之处； 4．掌握多重循环的编程技巧。 二．实验内容 【实验题1】. 程序填空：分别用for、while、do-while语句编程，输出菲波那契数列的前20项，要求每行输出10项。 源程序1：用for语句实现 编译、连接并运行该程序，则显示 思考：line 5和line 8行中printf()函数的格式控制符均为”m”, 假如均改为”%d”, 输出结果有何不同？
输出的数据无间隔 源程序2：用while语句实现 #include
读书破万卷下笔如有神 编译、连接并运行程序，输入0.5, 则显示：；输入 1 ，则显示：；输入1.7，则显示：。 略 line 21中的表达式“f*=i”，等同于表达式“f=f*i ”;同理，line 11中的表达式“s= s+item”也可替换为“s+=item ”. 【实验题3】程序填空：验证哥德巴赫猜想——任何一个大于6的偶数均可表示为两个素数之和。例如：6=3+3， 8=3+5，10=5+5，…，18=7+11。要求将6-100之间的偶数均表示为素数之和，一行输出5组。【分析】 1. 本程序采用双层循环：外层循环负责对6～100之间的偶数逐个输出加法式子；内层循环负责为一个特定的偶数k, 找出相应的素数a和b. 2. 对于偶数k（k>=6）,符合条件的两个素数a和b必然位于区间[3, k-1].(为什么不从2开始？) 3. 由于k=a+b，通过内层循环试探性地找出第一个加数a后，第二个加数b可用推算而得，再调用函数prime()来判断a、b是否是素数，如果两者均为素数，则停止寻找素数，并输出k=a+b. 源程序如下：
实验四 循环结构程序设计 教材P364 一、实验目的 1. 掌握利用for 语句实现循环结构的方法。 2. 掌握利用while 语句实现循环结构的方法。 3. 熟悉利用向量运算来代替循环操作的方法。 二、实验内容 1. 根据2222π11116123n =++++"2，求的近似值。当n 分别取100、1000、10000时，结果是多少？ π要求：分别用循环结构和向量运算（使用sum 函数）来实现。 答：程序设计： (1) 循环结构 clear all; close all; clc; for n=[100,1000,10000] sum=0; for x=1:n sum=sum+1/(x.^2); end value=sqrt(6*sum) end (2) 向量运算 clear all; close all; clc; for n=[100,1000,10000] x=1:n; value=sqrt(6*sum(1./x./x)) end 运行结果： value = 3.1321 value = 3.1406 value = 3.1415 2. 根据11113521 y n =++++?"，求： (1) 时的最大n 值。 3y <(2) 与(1)的n 值对应的y 值。 答：程序设计：
clear all; close all; clc; y=0; n=1; while y<3 y=y+1/(2*n-1); n=n+1; end n=n-1; y=y-1/(2*n-1) n=n-1 运行结果： y = 2.9944 n = 56 验证： clear all; close all; clc; n=56; i=1:n; f=1./(2*i-1); y=sum(f) 运行结果： y = 2.9944 3. 考虑以下迭代公式： 1n n a x b x +=+ 其中、b 为正的常数。 a (1) 编写程序求迭代的结果，迭代的终止条件为5110n n x x ?+?≤，迭代初值，迭代次数不超过500次。 0 1.0x =(2) 如果迭代过程收敛于r ，那么r 的准确值是2 b ?±，当的值取、(8、时，分别对迭代结果和准确值进行比较。 (,)a b (1,1),3)(10,0.1)答： (1) 程序设计： clear all; close all; clc; a=2; b=3; x0=1.0; x1=0; n=1; y=abs(x1-x0); while n<=500 & y>10^(-5)