MATLAB实验四循环结构程序设计

实验四 循环结构程序设计 1.根据n

3

212

2

2

2

2

1

1

1

1

6

+?+++=

π，求π的近似值。当n 分别取100、1000、10000时，记过是多少？

n=100

程序设计：

n=100;

i=1:n;

f=1./i.^2;

x=sum(f);

y=sqrt(6*x)

运行结果：

x =

3.1321 n=1000

程序设计：

n=1000;

i=1:n;

f=1./i.^2;

x=sum(f);

y=sqrt(6*x)

运行结果：

x =

3.1406 n=10000

程序设计：

n=1000;

i=1:n;

f=1./i.^2;

x=sum(f);

y=sqrt(6*x)

运行结果：

x =

3.1415

2.根据y=1+1

-n 21

5

13

1

+

?++,求： （1）y<3时的最大n 值。

（2）与（1）的n 值对应的y 值。

程序设计：

y=1;n=1;

while(y<3)

n=n+2;

y=y+1/n;

end

y=y-1/n n=(n+1)/2 运行结果： y = 2.9944 n =

57 3.

考虑以下迭代公式：

x x n

n b a

+=

+1 其中a 、b 为正的常数。 （1）编写程序求迭代的结果，迭代的终止条件为105

1-+≤-x x n n ，迭代初值x 0=1.0，迭代次数不超过500次。 程序设计： x=1.0; a=2; b=3; x=a/(b+x); n=1;

while(abs(x-(a/x-b))>

10e-5&n<=500) x=a/(b+x); n=n+1; end x n

运行结果： x = 0.5616 n =

6

（2）如果迭代过程收敛于r ，那

么r 的准确值是

2

42

a

b b

+±

-，当

（a,b ）的值取（1，1）、（8，3）、（10，0.1）时，分别对迭代结果的准确值进行比较。 （a,b ）取（1，1） 程序设计： x=1.0; a=1; b=1; x=a/(b+x);

n=1;

while(abs(x-(a/x-b))>10e-5&n<=5 00)

x=a/(b+x);

n=n+1;

end

x

r1=(-b+sqrt(b^2+4*a))/2

r2=(-b-sqrt(b^2+4*a))/2

运行结果：

x =

0.6181

r1 =

0.6180

r2 =

-1.6180

（a,b）取（8，3）

程序设计：

x=1.0;

a=8;

b=3;

x=a/(b+x); n=1;

while(abs(x-(a/x-b))>10 e-5&n<=500)

x=a/(b+x);

n=n+1;

end

x

r1=(-b+sqrt(b^2+4*a))/2 r2=(-b-sqrt(b^2+4*a))/2 运行结果：

x =

1.7016

r1 =

1.7016

r2 =

-4.7016

（a,b）取（10，0.1）

程序设计：

x=1.0;

a=10;

b=0.1;

x=a/(b+x);

n=1;

while (abs(x-(a/x-b))>10e-5&n<=500) x=a/(b+x); n=n+1; end x

r1=(-b+sqrt(b^2+4*a))/2 r2=(-b-sqrt(b^2+4*a))/2 运行结果： x = 3.1127 r1 = 3.1127 r2 = -3.2127 4.已知

??????

?>+

-=======---3

,2

3

,12,01,13

2

1

321

n n n n f

f

f

f f f f n n n n

求f f

100

1

~

中：

（1）最大值、最小值、各数之和。

程序设计： f(1)=1; f(2)=0; f(3)=1; n=100; max=f(3); min=f(2); sum=2; for i=4:n

f(i)=f(i-1)-2*f(i-2)+f(i-3);

if f(i)>max max=f(i); elseif f(i)

sum=sum+f(i); end sum max min

运行结果：

sum =

-7.4275e+011

max =

4.3776e+011

min =

-8.9941e+011（2）正数、零、负数的个数。

程序设计：

f(1)=1;

f(2)=0;

f(3)=1;

n=100;

pn=2;

nn=0;

zn=1;

for i=4:n

f(i)=f(i-1)-2*f(i-2)+f( i-3);

if f(i)>0

pn=pn+1;

elseif f(i)==0

zn=zn+1;

else

nn=nn+1;

end

end

disp('正数的个数:'),pn disp('负数的个数:'),nn disp('零的个数:'),zn

运行结果：

正数的个数:

pn =

49

负数的个数:

nn =

49

零的个数:

zn =

25.若两个连续自然数的乘积减1是素数，则称这两个连续自然数是亲密数对，该素数是亲密素数。例如，2 3-1=5，由于5

是素数，所以2和3是亲密数对，5是亲密素数。求[2，50]区间内：（1）亲密数对的对数。

（2）与上述亲密数对对应的所有亲密素数之和。

程序设计：

x=1;

i=0;

y=0;

for n=2:50

m=n*(n-1);

l=m-1;

if x==isprime(l);

y=y+x*l;

i=i+1;

end

end

disp('请输入亲密对象的个数'); i

disp('亲密素数之和'); y

请输入亲密对象的个数i =

28

亲密素数之和y =

21066

实验二 Matlab程序设计基本方法1

实验二Matlab程序设计基本方法 覃照乘自092 电气工程学院 一、实验目的： 1、熟悉MATLAB 程序编辑与设计环境 2、掌握各种编程语句语法规则及程序设计方法 3、函数文件的编写和设计 4、了解和熟悉跨空间变量传递和赋值 二、实验基本知识： ◆for循环结构 语法：for i＝初值：增量：终值 语句1 …… 语句n end 说明：1．i＝初值：终值，则增量为1。 2．初值、增量、终值可正可负，可以是整数，也可以是小数，只须符合数学逻辑。 ◆while 循环结构 语法：while 逻辑表达式 循环体语句 end 说明：1、whiIe结构依据逻辑表达式的值判断是否执行循环体语勾。若表达式的值为真，执行循环体语句一次、在反复执行时，每次都要进行判断。若表达 式的值为假，则程序执行end之后的语句。 2、为了避免因逻辑上的失误，而陷入死循环，建议在循环体语句的适当位 置加break语句、以便程序能正常执行。（执行循环体的次数不确定； 每一次执行循环体后，一定会改变while后面所跟关系式的值。） 3、while循环也可以嵌套、其结构如下：

while逻辑表达式1 循环体语句1 while逻辑表达式2 循环体语句2 end 循环体语句3 end ◆if-else-end分支结构 if 表达式1 语句1 else if 表达式2(可选) 语句2 else(可选) 语句3 end end 说明：1.if结构是一个条件分支语句，若满足表达式的条件，则往下执行；若不满足，则跳出if结构。 2．else if表达式2与else为可选项，这两条语句可依据具体情况取舍。 3．注意：每一个if都对应一个end，即有几个if，记就应有几个end。 ◆switch-case结构 语法：switch表达式 case常量表达式1 语句组1 case常量表达式2 语句组2 …… otherwise 语句组n end

c语言实验报告4 循环结构程序设计

《C程序设计》课程实验报告 学院：班级：姓名：学号： 实验设备：计算机1台实验日期：2011年3月1日 实验项目名称循环结构程序设计 实验目的 掌握使用三种循环语句实现循环结构的方法。 实验要求：能用while~ do~while for 三种循环语句实现循环结构，编写简单的程序，掌握这三种循环语句。 实验内容(包括步骤)：1．猴子吃桃问题。猴子第一天摘下若干个桃子，当即吃了一半另一个，以后每天早晨都吃剩下的一半另一个，到第十天早晨再想吃时，就剩一个桃子。问第一天共摘了多少桃子。 要求：分别用三种语句编写程序。 2．编写程序验证下列结论：任何一个自然数n的立方都等于n个连续奇数之和。例如：13=1；23=3+5；33=7+9+11;43=13+15+17+19。 要求：程序对每个输入的自然数计算并输出相应的连续奇数，直到输入的自然数为0时止。 3．编写程序，求1-3+5-7+…-99+101的值。 4．编写程序，求e的值。e ≈ 1+1/1!+1/2!+1/3!+1/4!+…+1/n! （1）用for循环，计算前50项。 （2）用while循环，要求直至最后一项的值小于10-6。 5．编写程序，输出从公元1600年至2000年所有闰年的年号。每输出5个年号换一行。判断公元年是否为闰年的条件是：

（1）公元年数如能被4整除，而不能被100整除，则是闰年。 （2）公元年数如能被400整除也是闰年。 6.编写程序，打印以下图形： * *** ***** ******* ***** *** * 7.输入一行字符，统计输入的字符中字母‘a’出现的次数（包括大小写）。 8.我国1991年有11.6亿人口，要求根据人口平均年增长率，计算从1991年算起经过多少年后我国的人口增加到15亿。 9.有一个四位数F=BBAA，前二位数相同，后二位数相同，此数又是某个自然数T的平方，求此数。 10.分解一个正整数的所有质因数 调试与结果测试：部分程序出现错误经过调整后调试正常，达到题目要求的结果，运行正常。

实验四分支和循环程序设计实验(精)

实验四分支和循环程序设计实验 一、实验要求和目的 1．熟悉汇编语言程序设计结构； 2．熟悉汇编语言分支程序基本指令的使用方法； 3．掌握利用汇编语言实现单分支、双分支、多分支的程序设计方法； 4．了解汇编语言循环程序设计的基本流程； 5．熟悉汇编语言循环基本指令的使用方法； 6．掌握利用汇编语言的循环指令完成循环程序设计方法。 二、软硬件环境 1、硬件环境：计算机系统 windows； 2、软件环境：装有 MASM、DEBUG、LINK、等应用程序。 三、实验涉及的主要知识 在实际应用中，经常根据一些条件来选择一条分支执行。汇编语言的条件判断主要是通过状态寄存器中的状态位、无符号数相减或有符号相减而导致的结果来进行。 1．无条件转移指令 JMP 无条件转移指令 JMP是使程序无条件转移至目标处，又分为段内转移、段间转移。 2．条件转移指令 JXX 条件转移指令可分为三大类： 1）．简单条件转移指令指令。根据单个标志位的状态判断转移条件。标志位指令转移条件意义 CF JC CF=1 有进位/借位 JNC CF=0 无进位/借位 ZF JE/JZ ZF=1 相等/等于 0 JNE/JNZ ZF=0 不相等/不等于 0 SF JS SF=1 是负数 JNS SF=0 是正数 OF JO OF=1 有溢出 JNO OF=0 无溢出 PF

JP/JPE PF=1 有偶数个 1 JNP/JPO PF=0 有奇数个 1 2）．无符号数条件转移指令。 假设在条件转移指令前使用比较指令，比较两个无符号数A，B，指令进行的的操作是 A-B，其转移指令如下： 指令转移条件意义 JA/JNBE CF=0 AND ZF=0 A>B JAE/JNB CF=0 OR ZF=1 A>=B JB/JNAE CF=1 AND ZF=0 AB JGE/JNL SF=OF OR ZF=1 A>=B JL/JNGE SF OF AND ZF=0 A

matlab循环语句

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 顺序结构 顺序结构是由两个程序模块串接构成。一个程序模块是完成一项独立功能的逻辑单元，它可以是一段程序、一个函数，或者是一条语句。 看图可知，在顺序结构中，这两个程序模块是顺序执行的，即先执行<程序

实验3MATLAB程序设计

1，编写M 函数实现求一个数是否为素数，再编写一主程序（脚本文件），要求通过键盘输入一个整数，然后调用判断素数函数，从而确定它是否素数。 x=input('请输入一个整数x:'); if myprime(x) disp('您输入的整数x是一个素数。') else disp('您输入的数x不是一个素数。') end function y=myprime(x) y=1; for i=2:fix(sqrt(x)) if mod(x,i)==0 y=0; end end 2，编写M 函数统计一数值中零的个数，然后编写脚本文件，实现统计从1—2007 中零的总个数。 function num=number0(a) %统计十进制数值中0的个数 sa=num2str(a);%将数值装化为字符串 num=length(find(sa=='0'));% ));%求取字符串中'0’的个数 y=0;

for a=1:2006 num=number0(a); y=num+y; end disp(y) 504 3，编写程序计算x∈[-3,3],字长0.01：并画出曲线x = -3:0.01:3; y=zeros(size(x)); for i = 1:length(x) if -3<= x(i)& x(i)<=-1 y(i)=(-x(i).^2-4*x(i)-3)/ 2; elseif -1<= x(i) & x(i)<=1 y(i)=-x(i).^2+1; elseif 1<=x(:,i)<=3 y(i)=(-x(i).^2+4*x(i)-3)/2; end end plot(x,y) -3-2-10123

循环结构程序设计实验报告

嘉应学院计算机学院 实验报告 课程名称程序设计基础实验名称实验地点 指导老师实验时间提交时间 班级姓名座号 一、实验目的和要求 （1）熟悉掌握用while语句、do…while语句和for语句实现循环的方法。 （2）掌握在程序设计中用循环的方法实现一些常用算法（如穷举、迭代、递推等）。 （3）进一步学习调试程序。 二、实验环境和方法 实验方法： （一）综合运用课本所学的知识，用不同的算法实现在不同的程序功能。 （二）结合指导老师的指导，解决程序中的问题，正确解决实际中存在的异常情况，逐步改善功能。 （三）根据实验内容，编译程序。 实验环境：Windows xp Visual C++6.0 三、实验内容及过程描述 实验步骤： ①进入Visual C++ 6.0集成环境。 ②输入自己编好的程序。 ③检查一遍已输入的程序是否有错（包括输入时输错的和编程中的错误），如发现有错， 及时改正。 ④进行编译和连接。如果在编译和连接过程中发现错误，频幕上会出现“报错信息”， 根据提示找到出错位置和原因，加以改正。再进行编译，如此反复直到不出错为止。 ⑤运行程序并分析运行结果是否合理。在运行是要注意当输入不同的数据时所得结果 是否正确，应运行多次，分别检查在不同情况下结果是否正确。 实验内容：编译以下题目的程序并调试运行。 （1）输入一行字符，分别统计出其中的英文字母、空格、数字和其它字符的个数。 编写程序如下： #include int main() { char c; int letters=0,space=0,digit=0,other=0; printf("请输入一行字符：\n"); while((c=getchar())!='\n')

matlab中循环语句用法

循环结构：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循环

刘卫国版MATLAB程序设计与应用课后实验六八九

实验六 高层绘图操作 %第一题： 程序代码如下： x=linspace(0,2*pi,101); y=(0.5+3*sin(x)./(1+x.^2)).*cos(x); plot(x,y) 01234567 -1 -0.5 0.5 1 1.5 %第二题： %（1） 程序代码如下： x=linspace(-2*pi,2*pi,100); y1=x.^2; y2=cos(2*x); y3=y1.*y2; plot(x,y1,'b-',x,y2,'r:',x,y3,'y--'); text(4,16,'\leftarrow y1=x^2'); text(6*pi/4,-1,'\downarrow y2=cos(2*x)'); text(-1.5*pi,-2.25*pi*pi,'\uparrow y3=y1*y2');

-8 -6 -4 -2 2 4 6 8 -30-20 -10 10 20 30 40 %（2） 程序代码如下： x=linspace(-2*pi,2*pi,100); y1=x.^2; y2=cos(2*x); y3=y1.*y2; subplot(1,3,1);%分区 plot(x,y1); title('y1=x^2');%设置标题 subplot(1,3,2); plot(x,y2); title('y2=cos(2*x)'); subplot(1,3,3); plot(x,y3); title('y3=x^2*cos(2*x)');

-10 10 0510 15202530 35 40y1=x 2 -10 10 -1-0.8 -0.6 -0.4-0.200.20.4 0.6 0.8 1y2=cos(2*x) -10 10 -30-20 -10 10 20 30 40 y3=x 2*cos(2*x) %（3） 程序代码如下： x=linspace(-2*pi,2*pi,20); y1=x.^2; subplot(2,2,1);%分区 bar(x,y1); title('y1=x^2的条形图');%设置标题 subplot(2,2,2); stairs(x,y1); title('y1=x^2的阶梯图'); subplot(2,2,3); stem(x,y1); title('y1=x^2的杆图'); subplot(2,2,4); fill(x,y1,'r');%如果少了'r'则会出错 title('y1=x^2的填充图'); %其他的函数照样做。

C语言循环结构程序设计实验报告

甘肃政法学院 本科生实验报告 ( 五） 姓名:赵明翔 学院:公安技术学院 专业:安全防范工程 班级: 2015级安全防范工程班 实验课程名称:程序设计 实验日期:2016年月日 开课时间：2015学年第二学期

（2）输出所有的水仙花数，所谓水仙花数是指一个3位数，其各位数字李方和等于该数本身。例如，153是水仙花数，因为153=12+53+33程序如下： #include int main() { int i,j,k,n; printf("paicissus numbers are"); for (n=100;n<1000;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); } printf("\n"); return 0; } 运行结果： （3）猴子吃桃问题。猴子第1天摘下若干个桃子，当即吃了一半，还不过瘾，又多吃了一个。第2天早上又将剩下的桃子吃掉一半，又

多吃了一个。以后每天早上都吃了前一天剩下的一半零一个。到第10天早上想再吃时，见只剩下了1个桃子了。求第一天共摘了多少桃子。 程序如下： #include int main() { int day,x1,x2; day=9; x2=1; while(day>0) {x1=(x2+1)*2; x2=x1; day--; } printf("total=%d\n",x1); return 0; } 运行结果： （4）①用牛顿迭代法求方程0634223=-+-x x x 在1.5附近的根. #include #include int main() {double x1,x0,f,f1; x1=1.5; do {x0=x1; f=((2*x0-4)*x0+3)*x0-6; f1=(6*x0-8)*x0+3; x1=x0-f/f1;

Matlab实验六 循环结构程序设计答案

．实验五循环结构程序设计实验内容 （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的编程方法和循环结构程序设计方法，但都不是特别的熟练，都有待加强。

MATLAB程序设计实验报告

MATLAB 程序设计实验报告 一、实验目的 1. 通过实验熟悉MATLAB 仿真软件的使用方法； 2. 掌握用MATLAB 对连续信号时域分析、频域分析和s 域分析的方法，利用绘图命令绘制出典型信号的波形，了解这些信号的基本特征； 3. 掌握用MATLAB 对离散信号时域分析、频域分析和z 域分析的方法，利用绘图命令绘制出典型信号的波形，了解这些信号的基本特征； 4. 通过绘制信号运算结果的波形，了解这些信号运算对信号所起的作用。 二、实验设备 1. 计算机 ： 2. MATLAB R2007a 仿真软件 三、实验原理 对系统的时域分析 信号的时域运算包括信号的相加、相乘，信号的时域变换包括信号的平移、反折、倒相及信号的尺度变换。 （1）信号的相加和相乘：已知信号)(1t f 和)(2t f ，信号相加和相乘记为 )()(1t f t f =)(2t f +；)()(1 t f t f =)(2t f *。 （2）信号的微分和积分：对于连续时间信号，其微分运算是用diff 函数来完成的，其语句格式为：diff(function,’variable’,n)，其中function 表示需要进行求导运算的信号，或者被赋值的符号表达式；variable 为求导运算的独立变量；n 为求导的阶数，默认值为求一阶导数。连续信号的积分运算用int 函数来完成，语句格式为：diff(function,’variable’,a,b)，其中function 表示需要进行被积信号，或者被赋值的符号表达式；variable 为求导运算的独立变量；a,b 为积分上、下限，a 和b 省略时为求不定积分。 （3）信号的平移、翻转和尺度变换 信号的平移包含信号的左移与右移，信号的翻转包含信号的倒相与折叠，平移和翻转信号不会改变信号)(t f 的面积和能量。信号的尺度变换是对信号)(t f 在时间轴上的变化，可使信号压缩或扩展。)(at f 将原波形压缩a 倍，)/(a t f 将原波形扩大a 倍。 ￥ 对系统频率特性的分析

循环程序设计实验

循环程序设计实验 实验目的 1、掌握循环程序的设计方法。 2、掌握比较指令、转移指令和循环指令的使用方法。 3、进一步掌握调试工具的使用方法。 实验预习要求 复习比较指令、条件转移指令和循环指令。 复习循环程序的结构、循环控制方法等知识。 读懂“实验内容”中给出的将十进制数转换为二进制数以及将二进制数转换为十进制数的程序。 根据“实验内容”中给出的流程图和程序框架编写源程序，以便上机调试。 从“实验习题”中任选一道题目，编写源程序，以便上机调试。 实验内容 计算1+2+……n=，其中n通过键盘输入。要求在屏幕上提供如下信息： Please input a number(1627): ；出现此信息后通过键盘输入一个小于628的无符号整数 1+2+…..n=sum；其中n为用户输入的数，sum为所求的累加和 程序运行情况如下图所示（说明：图中所运行程序允许累加和不大于一个32位二进制数所能表示的范围）。编程指导 键盘输入的十进制数如 368在计算机中是以 33H，36H，38H形式存 放的，如何将它们转换 为一个二进制数0B，以 便对累加循环的循环次 数进行控制是本程序首 先要解决的问题。将键 盘输入的十进制数转换为二进制数的程序清单如下： DATA SEGMENT INF1 DB "Please input a number (0-65535):$" IBUF DB 7,0,6 DUP() DATA ENDS CODE SEGMENT ASSUME CS: CODE, DS:DATA START: MOV A X, DATA MOV DS, AX MOV DX, OFFSET INF1 MOV AH, 09H INT 21H MOV DX, OFFSET IBUF ；键入一个十进制数（<65535） MOV AH, 0AH INT 21H MOV CL, IBUF+1 ；十进制数的位数送CX MOV CH, 0 MOV SI, OFFSET IBUF+2 ；指向输入的第一个字符（最高位）

matlab循环结构程序设计

实验四循环结构程序设计 一、实验目的 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 =

Matlab程序设计实验报告

实验七Matlab 程序设计 实验目的： 1、掌握建立和执行M 文件的方法； 2、掌握实现选择结构的方法； 3、掌握实现循环结构的方法。 实验内容： 1. 编写用 5 次多项式拟合函数y=sin(x), x [0, 2 ]的脚本M 文件，要求绘图观察拟合的效果。 function shiyan1 x=0:0.5:2*pi y=sin(x) p=polyfit(x,y,5) x1=0:0.2:2*pi y1=polyval(p,x1) plot(x,y, 'b' ,x1,y1, '*r' x =

Columns 1 through 9 0 0.5000 1.0000 1.5000 2.0000 2.5000 3.0000 3.5000 4.0000 Columns 10 through 13 4.5000 5.0000 5.5000 6.0000 y = Columns 1 through 9 0 0.4794 0.8415 0.9975 0.9093 0.5985 0.1411 -0.3508 -0.7568 Columns 10 through 13 -0.9775 -0.9589 -0.7055 -0.2794 p = -0.0056 0.0881 -0.3967 0.2671 0.8902 0.0029 x1 = Columns 1 through 10 0 0.2000 0.4000 0.6000 0.8000 1.0000 1.2000 1.4000 1.6000 1.8000 Columns 11 through 20

2. 2.2000 2.4000 2.6000 2.8000 3.0000 3.2000 3.4000 3.6000 1.8001 Columns 21 through 30 4.0 4.2000 4.4000 4.6000 4.8000 5.0000 5.2000 5.4000 5.6000 5.8000 Columns 31 through 32 6.0 6.2000 y1 = Columns 1 through 10 0.29 0.1886 0.3786 0.5585 0.7172 0.8461 0.9391 0.9926 1.0048 0.9761 Columns 11 through 20 0.9083 0.8048 0.6701 0.5098 0.3301 0.1381 -0.0590 -0.2538 -0.4389 -0.6073 Columns 21 through 30 -0.7524 -0.8685 -0.9505 -0.9949 -0.9991 -0.9626 -0.8863 -0.7732 -0.6288 -0.4606 Columns 31 through 32

循环程序设计实验

2.4 循环程序设计实验 2.4.1 实验目的 1、掌握循环程序的设计方法。 2、掌握比较指令、转移指令和循环指令的使用方法。 3、进一步掌握调试工具的使用方法。 2.4.2 实验预习要求 1、复习比较指令、条件转移指令和循环指令。 2、复习循环程序的结构、循环控制方法等知识。 3、读懂“2.4.3 实验内容”中给出的将十进制数转换为二进制数以及将二进制数转换为十进制数的程序。 4、根据“2.4.3 实验内容”中给出的流程图和程序框架编写源程序，以便上机调试。 5、从“2.4.4 实验习题”中任选一道题目，编写源程序，以便上机调试。 2.4.3 实验内容 计算1+2+……n=?，其中n通过键盘输入。要求在屏幕上提供如下信息： Please input a number(1 627): ；出现此信息后通过键盘输入一个小于628的无符号整数 1+2+…..n=sum；其中n为用户输入的数，sum为所求的累加和 程序运行情况如下图所示（说明：图中所运行程序允许累加和不大于一个32位二进制数所能表示 的范围）。 1、编程指导 （1）键盘输入的十进 制数如368在计 算机中是以 33H，36H，38H 形式存放的，如 何将它们转换为 一个二进制数 101110000B，以便对累加循环的循环次数进行控制是本程序首先要解决的问题。将键盘输入的十进制数转换为二进制数的程序清单如下： DA TA SEGMENT INF1 DB "Please input a number (0-65535):$" IBUF DB 7,0,6 DUP(?) DA TA ENDS CODE SEGMENT ASSUME CS: CODE, DS:DATA START: MOV AX, DATA MOV DS, AX MOV DX, OFFSET INF1 MOV AH, 09H INT 21H MOV DX, OFFSET IBUF ；键入一个十进制数（<65535） MOV AH, 0AH INT 21H MOV CL, IBUF+1 ；十进制数的位数送CX MOV CH, 0

挺好的——matlab循环语句

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 顺序结构 顺序结构是由两个程序模块串接构成。一个程序模块是完成一项独立功能的逻辑单元，它可以是一段程序、一个函数，或者是一条语句。 看图可知，在顺序结构中，这两个程序模块是顺序执行的，即先执行<程序

MATLAB程序设计及经典例题解析3

MATLAB程序设计 用MATLAB语言编写的程序，称为M文件。M文件可以根据调用方式的不同分为两类：命令文件(Script File)和函数文件(Function File)。 例3-1 分别建立命令文件和函数文件，将华氏温度f转换为摄氏温度c。 程序1：首先建立命令文件并以文件名f2c.m存盘。 clear; %清除工作空间中的变量 f=input('Input Fahrenheit temperature：'); c=5*(f-32)/9 然后在MATLAB的命令窗口中输入f2c，将会执行该命令文件，执行情况为： Input Fahrenheit temperature：73 c =22.7778 程序2：首先建立函数文件f2c.m。 function c=f2c(f) c=5*(f-32)/9 然后在MATLAB的命令窗口调用该函数文件。 clear; y=input('Input Fahrenheit temperature：'); x=f2c(y) 输出情况为： Input Fahrenheit temperature：70 c =21.1111 x =21.1111 3.1.2 M文件的建立与打开 M文件是一个文本文件，它可以用任何编辑程序来建立和编辑，而一般常用且最为方便的是使用MATLAB提供的文本编辑器。

1．建立新的M文件 为建立新的M文件，启动MATLAB文本编辑器有3种方法： (1) 菜单操作。从MATLAB主窗口的File菜单中选择New菜单项，再选择M-file命令，屏幕上将出现MATLAB 文本编辑器窗口。 (2) 命令操作。在MATLAB命令窗口输入命令edit，启动MATLAB文本编辑器后，输入M文件的内容并存盘。 (3) 命令按钮操作。单击MATLAB主窗口工具栏上的New M-File命令按钮，启动MATLAB文本编辑器后，输入M文件的内容并存盘。 2．打开已有的M文件 打开已有的M文件，也有3种方法： (1) 菜单操作。从MATLAB主窗口的File菜单中选择Open命令，则屏幕出现Open对话框，在Open对话框中选中所需打开的M文件。在文档窗口可以对打开的M文件进行编辑修改，编辑完成后，将M文件存盘。 (2) 命令操作。在MATLAB命令窗口输入命令：edit 文件名，则打开指定的M文件。 (3) 命令按钮操作。单击MATLAB主窗口工具栏上的Open File命令按钮，再从弹出的对话框中选择所需打开的M文件。 3.2 程序控制结构 3.2.1 顺序结构 1．数据的输入 从键盘输入数据，则可以使用input函数来进行，该函数的调用格式为： A=input(提示信息，选项)； 其中提示信息为一个字符串，用于提示用户输入什么样的数据。 如果在input函数调用时采用's'选项，则允许用户输入一个字符串。例如，想输入一个人的姓名，可采用命令： xm=input('What''s your name?','s'); 2．数据的输出 MATLAB提供的命令窗口输出函数主要有disp函数，其调用格式为

D实验五 M文件和MATLAB程序设计

实验五M文件和MATLAB程序设计 一、实验目的 matlab作为一种高级计算机语言，不仅可以命令行方式完成操作，也具有数据结构、控制流、输入输出等能力，本次实验通过熟悉和掌握m文件的建立与使用方法，以及函数与控制程序流程语句的使用，使学生具备一定的编程和程序调试能力。 1．掌握M文件的使用方法。 2．掌握if语句和switch语句的使用 3. 掌握循环语句的使用 4. 通过练习理解MATLAB编程方法。 二、实验原理 1．m文件 用matlab语言编写的程序，称为m文件。M文件根据调用方式的不同分为两类，命令文件（Script file）和函数文件（Function file）。区别? 2．程序控制结构 1）顺序结构 2）选择结构 （1）if语句a) 单分支if语句b) 双分支if语句c) 多分支if语句 （2）switch 语句 （3）try语句 3）循环结构 （1）for 语句 （2）while语句 （3）break语句、continue语句、return使用,区别? 3．函数文件 function 输出形参表＝函数名（输入形参表） 注释说明部分 函数体语句 注意事项？ 三、实验要求 1．首先上机练习PPT中各种流程控制语句的有关实例。 2．然后上机练习下面的实验习题。 四、实验习题

1．数论中一个有趣的题目：任意一个正整数，若为偶数，则用2除之，若为奇数，则与3相乘再加上1。重复此过程，最终得到的结果为1。如： 2→1 3→10→5→16→8→4→2→1 6→3→10→5→16→8→4→2→1 运行下面的程序，按程序提示输入n=1,2,3,5,7，8,9等数来验证这一结论。 %classic "3n+1" problem from number theory. while 1 n=input('Enter n,negative quits:'); if n<=0 break end a=n; while n>1 if rem(n,2)==0 n=n/2; else n=3*n+1; end a=[a,n]; end a end Enter n,negative quits: 2． 编程求满足∑=>m i i 11000020的最小m 值。 a=0; i=1; while (a<100000) a=a+pow2(i); i=i+1; end m=i-1 3． 编写一个函数，计算下面函数的值，给出x 的值，调用该函数后，返回y 的值。 function [y]=myfun1(x)

微机原理与汇编语言_实验3_循环程序设计实验汇总

微机原理与汇编语言 实验报告 姓名x x x学号xxxxxx 专业班级计科x班 课程名称微机原理与汇编语言实验日期2014.10.22 实验名称循环程序设计实验成绩 一、实验目的 1、掌握循环程序的设计方法。 2、掌握比较指令、转移指令和循环指令的使用方法。 3、进一步掌握调试工具的使用方法。 二、实验内容 1、实验原理 (1)键盘输入的十进制数如368在计算机中是以33H，36H，38H形式存放的，如何将它们转换为一个二进制数101110000B，以便对累加循环的循环次数进行控制是本程序首先要解决的问题 (2)累加结果为一个16位的二进制数，为了显示结果，必需把它们转换为十进制数 2、实验步骤 （1）流程图 （2）实验源码 .model small .stack ;定义堆栈段 .data ;定义数据段 inf1 db "Please input a number(1-627):$" inf2 db 0ah, 0dh, "1+2+......+$" ibuf db 7, 0, 6 dup(0) obuf db 6 dup(0) .code start: mov ax, @data mov ds, ax mov dx, offset inf1 ;将屏幕输入的十进制数(存放于ibuf)转化为二进制数存储于ax 中 mov ah, 09h int 21h mov dx, offset ibuf mov ah, 0Ah int 21h mov cl, ibuf+1 mov ch, 0 mov si, offset ibuf+2 mov ax, 0 loop1: mov dx, 10 mul dx and byte ptr [si], 0Fh add al, [si] 从键盘输入一个十进制数，并将其转换为二进制数，存放在AX中 开始 结束 素、、 (CX)-1=0 累加循环次数送CX（MOV CX,AX）AX清0（AX中存放累加和） BX送1（BX存放每次循环累加的数，每循环一次，BX值加1） ADD AX，BX INC BX 累加和（在AX中）转换为十进制数并 显示 Y N

矩量法matlab程序设计实例

矩量法m ａtla ｂ程序设计实例： Ha ｌlen 方程求对称振子天线 一、条件与计算目标 已知: 对称振子天线长为Ｌ,半径为a ，且天线长度与波长得关系为,，设,半径a=0、０００0001，因此波数为。 目标： 用H ａll ｅn 方程算出半波振子、全波振子以及不同值得对应参数值。 求:（1）电流分布 （2）E 面方向图 (二维）,H 面方向图（二维）,半波振子空间方向性图(三维） 二、对称振子放置图 图1 半波振子得电流 分布 半波振子天线平行于z 轴放置，在ｘ轴与ｙ轴上得分量都为零，坐标选取方式有两种形式，一般选取图1得空间放置方 式。图１给出了天线得电流分布情况，由图可知，当天线很细时,电流分布近似正弦分布。 三、Ha ｌｌen 方程 得解题思路 ()()()()2 1 ' ' ' ' 12,cos sin sin 'z z i z z z z i z k z G z z dz c kz c kz E k z z dz j ωμ'++=-?? 对于中心馈电得偶极子，Hallen 方程为 ()22'1222 ('),'cos sin sin ,2L L i L L V i z G z z dz c kz c kz k z z j η + -- ++= <<+? 脉冲函数展开与点选配，得到 ()1121 ,''cos sin sin ,1,2,,2n n N z i n m m m m z n V I G z z dz c kz c kz k z m N j η +''=++= =???∑? 上式可以写成 矩阵形式为 四、结果与分析 （1）电流分布