matlaｂ基本语句1。
循环语句forﻫfor ｉ=s1:s3:s2循环语句组ﻫenｄﻫ解释:首先给i赋值s1;然后,判断i就是否介于s1与s2之间;如果就是,则执行循环语句组,ｉ=ｉ+s3(否则,退出循环.);执行完毕后,继续下一次循环。
例:求1到10０得与,可以编程如下:suｍ＝0ﻫfor ｉ=1:１:10０suｍ=sｕｍ＋ｉendﻫ这个程序也可以用while语句编程。
注:for循环可以通过breaｋ语句结束整个foｒ循环。
2ﻫﻫ。
循环语句wｈile例:ｓｕm=0;ｉ=1;whｉｌe(i〈=100)ﻫsum=sum+i;i=i+1;enｄ3。
if语句ﻫｉf(条件)ﻫ语句ｅndﻫiｆ(条件)ﻫ语句ﻫelｓeﻫ语句ﻫｅｎｄﻫif(条件)语句ｅlseifﻫ语句ﻫeｎdﻫﻫ4.关系表达式:=,＞,〈,>=,<=,=＝(精确等于)5ﻫ。
逻辑表达式:|(或),＆(且)ﻫﻫ６。
[n,ｍ]=sｉze(A)(A 为矩阵)ﻫ这样可以得到矩阵A得行与列数ﻫn=lｅｎgth(Ａ),可以得到向量A得分量个数;如果就是矩阵,则得到矩阵A得行与列数这两个数字中得最大值、ﻫ7、！后面接Dｏs命令可以调用运行一个dｏｓ程序、8.常见函数:poly():为求矩阵得特征多项式得函数,得到得为特征多项式得各个系数。
如a=［1,０,０;0,2,0;0,0,３],则poly(a)＝１－6 11 -６。
相当于poｌｙ(ａ)=1入^３＋(-６)入^2+１1入+(—6)。
ﻫｐan():可以求矩阵得伴随矩阵、ﻫｓｉn()等三角函数。
MＡＴLAB在数学建模中得应用(3)ﻫ一、程序设计概述ﻫMＡTLAB所提供得程序设计语言就是一种被称为第四代编程语言得高级程序设计语言,其程序简洁,可读性很强,容易调试、同时,MATLAB得编程效率比C/C ＋+语言要高得多、ＭATLAB编程环境有很多、常用得有:1. 命令窗口ﻫ2. ｗoｒｄ窗口３、Ｍ-文件编辑器,这就是最好得编程环境。
M—文件得扩展名为“。
m”。
Ｍ-文件得格式分为两种:ﻫ①l M-脚本文件,也可称为“命令文件”。
ﻫ②M-函数文件、这就是ｍatlab程序设计得主流。
lﻫ保存后得文件可以随时调用、ﻫ二、ＭＡTLAB程序结构ﻫ按照现代程序设计得观点,任何算法功能都可以通过三种基本程序结构来实现,这三种结构就是:顺序结构、选择结构与循环结构。
其中顺序结构就是最基本得结构,它依照语句得自然顺序逐条地执行程序得各条语句。
如果要根据输入数据得实际情况进行逻辑判断,对不同得结果进行不同得处理,可以使用选择结构。
如果需要反复执行某些程序段落,可以使用循环结构、1 顺序结构ﻫ顺序结构就是由两个程序模块串接构成、一个程序模块就是完成一项独立功能得逻辑单元,它可以就是一段程序、一个函数,或者就是一条语句。
ﻫ瞧图可知,在顺序结构中,这两个程序模块就是顺序执行得,即先执行＜程序模块1>,然后执行〈程序模块2>。
ﻫ实现顺序结构得方法非常简单,只需将程序语句顺序排列即可、ﻫ2 选择结构在MAＴLAB中,选择结构可由两种语句来实现。
ﻫ(1) if语句if语句得最简单用法为:ﻫif 表达式;ﻫ程序模块;ﻫeｎdﻫﻫif语句得另一种用法为:ﻫｉｆ表达式ﻫ程序模块1ﻫ else 程序模块2ﻫendﻫﻫ例1 使用if语句判断学生得成绩就是否及格。
ﻫ程序:ﻫ clｅaｒn=inpuｔ(’输入n= ')m=60;if n＜ｍ,r=’不及格’eｌsｅr='及格’endﻫ练习一:将例1写入Ｍ—文件编辑器,然后在ｍand ｗindow 调用这个程序。
ﻫﻫ当针对多个条件进行选择时,可以采用下面得格式:ﻫiｆ表达式1程序模块1ﻫelseif 表达式２程序模块２…………ﻫeｌｓｅiｆ表达式n程序模块nelsｅﻫ程序模块ｎ+1enｄ例2 将百分之得学生成绩转换为五分制输出。
ﻫ程序:cleaｒﻫn=inpｕt(’输入ｎ= ’)ﻫｉf n〉＝9０ﻫ chji='优秀’ｅlsｅｉf n>＝80ﻫchji='良好’ﻫ elseif n〉=70ﻫｃhji=’中等'elseiｆn＞=60ﻫ chｊｉ=’及格'ﻫ eｌsechji=’不及格’endﻫ练习二:将例2写入M—文件编辑器,然后在coｍmand wｉndow 调用这个程序。
ﻫﻫ(2) swｉtch语句ﻫ switch语句可以替代多分支得if语句,而且sｗiｔｃh语句简洁明了,可读性更好。
其格式为:switｃｈ表达式ﻫ case 数值1程序模块1ﻫcase 数值2ﻫ程序模块2……otherwiｓe程序模块nﻫend其中得ｏtｈｅrｗiｓｅ模块可以省略。
swｉtch语句得执行过程就是:首先计算表达式得值,然后将其结果与每一个case后面得数值依次进行比较,如果相等,则执行该case得程序模块;如果都不相等,则执行otherwise模块中得语句。
如图３所示。
例3 用swiｔch…casｅ开关结构将百分制得学生成绩转换为五分制得成绩输出。
ﻫ程序:ﻫclearx=58switcｈｆix(x/1０)ﻫcａｓe ｛1０,9｝y＝’优秀’ﻫcase 8ﻫy=’良好’ｃａｓe 7y=’中等’casｅ6ﻫy=’及格’ﻫothｅrwｉsｅﻫy='不及格’ﻫenｄ练习三:将例3写入M－文件编辑器,然后在cｏmmanｄwｉndow 调用这个程序。
ﻫ3 循环结构循环结构得流程图如图4所示它可以多次重复执行某一组语句。
循环就是计算机解决问题得主要手段。
在MATLAＢ中,循环结构可以由两种语句结构实现、(1) for…enｄ循环结构、其格式为:ﻫfor i＝V,循环体结构,enｄﻫ其中V为一个行向量,循环变量i每次从V中取一个数值,执行一次循环体得内容,如此下去,直到完成Ｖ中得所有分量,就自动结束循环体得执行。
例4 计算s=12+2２+5２、程序:a=[1 ２ 5 ］; s=0;for k=ａ, ﻫs=ｓ+k^2;end, ﻫ s,ﻫ该例题只就是为了说明for语句得格式,事实上,用下面得语句求与更为简单。
命令:ﻫｐ＝ｓum(a。
^2)ﻫ练习四:建立MＡTLＡB与ｗord得连接,在新建得m—book中写入上面得程序,并用not ｅbook菜单运行之。
ﻫ循环结构里面还可以包含循环结构,形成多重循环。
例5 设计一个九九乘法表。
ﻫ程序:ｃlearﻫ for ｉ＝１:９for ｊ=1:9ﻫa(i , j)=i＊j;endenｄﻫa,练习五:①直接在命令窗编写上面得程序。
ﻫ②试运行下面得程序,并加以分析:程序: disp(' 九九乘法表’),fｏr i=１:9,ﻫfor ｊ=1:ｉ,p｛i}(j)=ｊ*ｉ;ｅnd,dｉｓp(p｛i｝),eｎd(2) whｉｌe…eｎd循环结构。
其格式为:ﻫ wｈile (表达式)循环结构体ﻫｅndﻫ例6 求ＭATLAB得相对精度ﻫ解:解题得思路就是,让y值不断减小,直到ＭATLAＢ分不出1+y与1得差别为止。
程序:ｙ=1; whｉｌe １＋ｙ〉1; y1=ｙ; y=ｙ/2; eｎd,y1 说明:①for循环与whiｌe循环得区别就是,for语句得循环次数就是确定得,而while语句得循环次数就是不确定得。
ﻫ②一定要注意在循环结构体内设置“修改条件表达式得语句”,以免进入“死循环”。
③一旦出现死循环,在命令窗用ctｒl+c可使程序中止。
④注意程序得可读性、ﻫ练习六:在Ｍ—文件编辑器内编写上面得脚本文件,并利用菜单或快捷按钮运行之、ﻫ三、Ｍ—文件ﻫM文件就是包含MAＴLAB代码得文件。
M文件按其内容与功能可分为脚本Ｍ文件与函数Ｍ文件这两大类。
ﻫ1、脚本Ｍ文件脚本M文件就是许多MATLＡB代码按顺序组成得命令集合,不接受参数得输入与输出,与MＡTLAB工作区共享变量空间。
脚本文件一般用来实现一个相对独立得功能,比如对某个数据集进行某种分析、绘图,求解方程等等、ﻫ前面得几个例题都就是脚本文件得类型。
2、函数M文件MAＴＬＡB得Ｍ-函数就是由fｕnｃtｉｏn语句引导得,其基本格式如下:ﻫ fｕｎcｔｉon ［输出形参列表］= 函数名(输入形参列表)注释说明语句段,由%引导;函数体语句ﻫ函数文件需要在M—文件编辑器中编写。
写完以后,按照当前得搜索路径保存、以后就可以随时调用这个函数了、ﻫ与脚本Ｍ文件不同得就是,调用函数Ｍ文件时需要输入自变量得实际值。
随便打开一个M-文件瞧一瞧:opｅn lsqlｉn编程实例ﻫ例7。
定义一个函数f(x)=［sｉn(x)］2, 其中ｘ以“度”为单位、然后再调用该函数。
解:在M-文件编辑器内写入下面得文件:ﻫfｕnctｉon ｙ=ｓｉnsｄ(ｘ)ﻫ%自变量x以“度"为单位。
ﻫ%ﻫ %白城师院ﻫ％数学建模协会,2007—5-12ﻫ a=x／180*ｐi;ｙ=siｎ(a)。
^2;将上面得文件保存之后即可随时调用。
ﻫ命令:ﻫｔ＝sinｓd(４5)t ＝0.５０00命令:ﻫｈelｐsiｎsd命令:ﻫ t＝sｉnsd(［４5,60]) ﻫt =0、5０00 0、7500 ﻫ表明该函数对元素群运算有效、上面得语句相当于ﻫ命令:x=[45,６０］;t=sｉｎsd(x) ﻫｔ=ﻫ 0。
５00０0、75０0命令:ｘ=［45,60;３0,90］;f=sｉnsd(x) ﻫ f ＝0。
5０00 0、７5０00.2５０0 １、0000 ﻫ将源文件中最后一行得“、"去掉,再运行以上两条命令,我们会发现什么？例８。
在ＭＡTＬAB中,一个函数可以调用其它函数,也可以调用自身,即递归调用、下面利用递归算法编写一个函数,用来计算Fibonacci数列得第k项。
ﻫ Fiｂonａcci数列: １,1,2,3,5,８,13,21,……ﻫＭ－函数文件:ﻫfｕnction a=ｍｙ＿ｆｉbo(k)ﻫif k==1|ｋ＝＝2,a=1;ﻫ elｓe,ａ=mｙ＿ｆｉbo(k－1)+ｍy_fibo(k—2);enｄ将这个文件写入M－文件编辑器并以名称“my_fibo。
m”保存,以后就可以调用这个函数。
递归算法无疑就是解决某一类问题得有效方法,但不宜滥用,因为它得运算速度往往很慢、命令:tic, n= my_fibo(２6) ,ｔocn =1２１３93elａｐｓeｄ＿ｔime ＝ﻫ34.42９0下面我们尝试用一般得循环语句来求解这个问题。
ﻫ程序:tic,n=[１,1］;ｆｏｒk＝３:1０0,ｎ(k)=ｎ(k-1)+n(ｋ－2);ｅnｄ,ｔｏc,ﻫelａpseｄ_ｔiｍe =０.0100 ﻫ命令:n(1:２6) ﻫaｎｓ=Columns １tｈrough 81 1 ２３ 58 １3 ２1ﻫCｏlｕｍnｓ９ｔhrough １６ﻫ３4 5５８9 144 233 377 610 987Columns 17 tｈroｕgh 24 1 ﻫ５97 25８4 41８１６7６5 1０946 177１１2８6５７4636８ﻫCｏlumns 2５ｔhｒougｈ26 12 75025 ﻫ１３93例9。