MALAB简单程序设计

MATLAB程序设计方法

表达式和运算符命令文件(程序文件)程序设计中的流程控制练习与思考题表达式——由变量、运算符、函数、数字组成球坐标变换————表达式和运算符————用赋值语句创建变量变量名=数据变量名=表达式

x=R*cos(theta)*cos(fai);y=R*cos(theta)*sin(fai);z=R*sin(theta);——纬度——经度例2.1格林威治天文台建于1675年，其经度为零,纬度为北纬510。计算格林威治天文台在地心直角坐标系中的坐标并在球面上标记[X,Y,Z]=sphere(24);mesh(X,Y,Z),holdonfai=0;theta=51*pi/180;x0=cos(theta)*cos(fai);y0=cos(theta)*sin(fai);z0=sin(theta);th=linspace(-pi/2,pi/2,100);x=cos(th);y=zeros(size(x));z=sin(th);plot3(x,y,z,'b',x0,y0,z0,'ro','linewidth',2)view(124,30)+

加-

减

*

乘

.*点乘/

右除

./点右除\

左除

.\点左除^

矩阵幂

.^点幂1.算术运算符2.关系运算符

<

小于

<=

小于等于

>

大于

>=

大于等于

==

等于

~=

不等于3.逻辑运算符&与(and)|或(or)~非(not)————表达式和运算符————运算优先级①函数运算②算术运算③关系运算④逻辑运算exp()、log()、sin()、abs()、fix()、···.^、^、.*、./、*、/、+、–、；<、<=、>、>=、==、~=&、|表达式(1+fix(pi))*mod(2,4)+2*3^2结果为ans=26；三角形任意两边之和大于第三边的逻辑表达式

(A)a+b>=c|a+c>=b|b+c>=a(B)a+b<=c|a+c<=b|b+c<=a(C)a+b>c|a+c>b|b+c>a(D)a+b>c&a+c>b&b+c>a

命令文件——MATLAB命令的有序集合。文件执行——对文件中命令进行批处理，即从第一条命令开始按顺序执行，直到最后一条命令。如果中间某条命令出错，则中断并输出错误信息M文件分为命令文件和函数文件两种①在编辑窗口中编写;②保存并对文件命名；③命令窗口键入文件名运行；④观察运行结果；·······························例2.2成都地理位置:北纬30o,东经104o;北京地理位置:北纬40o,东经116o。计算两城市距离近似值。由两城市经纬度可计算出地心直角坐标系的坐标P1(x1，y1，z1)，P2(x2，y2，z2)向径之间的角度

球面短程线计算公式

①输入经纬度数据和地球半径；②转换两城市的经纬度为地心直角坐标数据；③提取两个点的向径坐标；④计算向径间的夹角和短程线长度并输出计算结果。datas=[30,104;40,116];%经纬度数据R=6400;theta=datas(:,1)*pi/180;%弧度制fai=datas(:,2)*pi/180;x=R*cos(theta).*cos(fai);%直角坐标转换y=R*cos(theta).*sin(fai);z=R*sin(theta);Op1=[x(1),y(1),z(1)];%第一城市坐标Op2=[x(2),y(2),z(2)];%第二城市坐标d=R*acos(Op1*Op2'/(R*R))%计算并显示距离

按顺序录入程序文件(文件名:distance.m)在命令窗口运行程序distance，屏幕将显示：d=1563.42条件控制;循环控制;错误控制;终止运行控制条件控制——有选择地运行程序块1.if/elseif语句例2.8判润年程序————程序设计中的流程控制————year=input('inputyear:=');n1=year/4;n2=year/100;n3=year/400;ifn1==fix(n1)&n2~=fix(n2)

disp('是润年')elseifn1==fix(n1)&n3==fix(n3)

disp('是润年')else

disp('不是润年')end润年条件有二①能被4整除，但不能被100整除；②能被4整除，又能被400整除。1.衬f迈or循环语法：fo随r变量=初值：步长见：终孤值可执阔行语锈句en兔d——侄——程序债设计芳中的宾循环巡寿控制——冈——例2.揪3裴波朽拉奇昼数列n=陆in叙pu吼t(北'i扇np公ut助n络:=霸')尊;f(预1)厦=1卧;f葱(2雄)=扰1;fo作rk=像3:承nf(胳k)优=f悄(k纷-1丝式)+践f(腔k-抓2)已;en凳df(痛n)f(分44疑)=谦70贪14伯08另73磁3例2.予3程序争功能为如下：把边庄长为2以原叨点为接中心干的正拉方形兆旋转pi献/2叉4，将其尊压缩(r=幼0.苏89),重复24次并康绘图xy=[斯-1扮-士1;裹1个-1楼;1镇1繁;-建1挠1;迎-1镜-招1]吸;A=闪[c知os组(p满i/椅24昼)尸-冲si正n(抹pi描/2德4)赤;抛%创建徐正交缎矩阵si权n(赴pi腊/2级4)撑co兄s(膜pi觉/2傻4)转];x=太xy闪(:撑,1在);高y=嫩xy缸(:论,2油);难%提取咱坐标嘉数据li粪ne丹(x塑,y昨),撑pa聚us蚊e(钥1)杆%画线医并暂屠停一脑秒fo掉r羡k=烫1:威24xy=.各89惧*xy*A互’;盘%旋转刊并压挣缩x=枯xy嫩(:庙,1饥);让y=休xy窑(:购,2纯);li轰ne碗(x耳,y邮),厘pa切us摄e(抹1)en东d例2.脑4空间肤四面培体底粗面是跟一三忘角形，将三妄角形降三个生顶点械与四裙面体加顶部构的顶搞点相怒连接以就形凶成三逐角形躁锥面。利用MA牛TL先AB的帧汗动画吐方法窄演示孤一个逢四面哗体旋裹转过途程。M谢=诵mo既vi森ei执n(洽16色);[X正,Y报,Z套]=踪蝶cy语li锈nd确er使([检1,猪0]彻,3厨);敲%创建泪四面登体坐墙标A=技[c旁os厌(p折i/猎8)闭-踏si替n(器pi块/8辅);惊si柳n(盟pi速/8屋)过co凯s(肿pi蒜/8熔)]天;fo室r弄k=宅1:淋16me政sh闪(X颈,Y扭,Z泻),塌ax叨is疾o降ff月%绘图co百lo虏rm头ap写([届0厕0杜1]念),荷vi蛮ew醒(-奸30睡,4泻0)M(脱:,愧k)评=ge扔tf阵ra件me;畏%帧设腹计xy=A葱*[众X(棉1,萄:)葛;Y夕(1偏,:贩)]洪;酱%旋转x=版xy膛(1领,:默);效X=痕[1衬;0嫁]*世x;子%坐标怀提取熊扩充药矩阵y=演xy挎(2此,:作);槽Y=世[1镇;0厘]*击y;en厨dmo俊vi舞e(拐M,绣5)宾%帧播洁放2.售w震hi红le循环将循捐环体耕中语江句循吸环执郊行直咳到条着件不斧成立消为止。语法：wh哗il科e条件被表达约式可执顿行语租句en页d条件督表达添式一与般由着变量需、数置字、茧逻辑穴运算缸、关察系运算符驾和一茄般运卖算符框组成爱，以平判断虫循环葵的进世行和鲜停止；只要予表达特式的停值(逻辑蚁值)结果匠为正丛确(非0)，循环爆继续；凳直到企表达精式值勇为0，循笔环停斜止。——眨——程序盘设计信中的边循环腔控制——恒——例2.浸5述3n+状1问题.n=她in写pu踢t(批‘i祸np鲁ut锁n=馅’)蒜;%输入婚数据wh袍il水en~膏=1r=亡re奏m(宽n,薄2)缺;%求n/宾2的余跨数if孝r批=雄=规0n=恰n/拼2%第一阔种操撕作el饭sen=美3*兔n+耕1%第二完种操羞作en灵den捎d对任赏一自吧然数n，按如魂下法仅则进江行运顷算:若n为偶捞数,则将n除2;若n为奇远数，则将n乘3加1。将运盯算结珍果按上面晌法则浑继续匹运算,重复量若干扩次后究结果卡最终谎是1.n=防51葱6,京8选,海4,叛2岛,否1——巩——程序惧设计疑中的稻流程钻控制——照——3.碍c逢on卵ti季nu理e命令通常驶用于fo脑r或wh途il龟e循环激语句陪中,与if语句挂一起寨使用,跳过库本次别循环,去执暴行下规一轮绩循环4.卸b坡re砖ak命令通常丸用于fo收r或wh哨il贷e循环的语句季中,与if语句焦一起盗使用,中止躬本次宽循环,跳出苍最内阴层循贫环If条件型表达光式，br雹ea京k,钟e驱nd5.er聋ro脾r(民‘m慢es叫sa舍ge联’)显示及文本me架ss贫ag岔e，并中欺断程挥序执游行If条件上表达管式,er档ro浴r(痰‘m加es中sa虹ge哭’),铁en汁d——敲——程序才设计皱中的剪流程具控制——肉——如果境一个连正整拌数n只能督被1和它符自身姻整除鸣，则鲁称这瓶个数饲为素城数(或质窃数)。例哪如：2、3、5、7、11、··谎··垦··例2.捆6用试悄商法替判别剧素数n=铅in稍pu殖t(豪'i利np钳ut温n殃:=侦')泡;fo巩r钉k=救2:惜n-棍1if谋m环od辽(n符,k绞)=克=0巾,b赠re续ak灿,e沈nd饱%中止崭循环en董dif输k俱<前n-顾1di玻sp('不是炼素数')el季sedi帝sp('是素搞数')en仿d——产——程序耻设计阁中的峰流程烛控制——肝——例2.踩7海伦公式其中，p牢=(a捕+寺b闹+骂c)/拉2,其中，a,缺b,线c分别仪为三室条边届长。>>音h瓶er贝oin仿pu棵t成a:伶=3in胃pu糖t中b:垫=4in袭pu垂t日c:言=8??如?笋Er摘ro倘r浪us众in疲g初==立>油he饲ro输入先错误,三角渠形两荒边之犁和大巷于第勾三边if色a仆+b凭<c旱|丸a舅+c吊<b若|压b巨+c吵<aer源ro历r(券'输入些错误,三角签形两洋边之银和大乓于第确三边')en烛dfu免nc彻ti坝on返回变量穗列表＝函数单名(输入变量胞列表)%注解钥说明输入治变量沉检测,输出廊变量缘瑞检测函数诊体——隔——函数时文件概的编述写——浙——fu朋nc那ti插on披Y屡=y关ou盯ng贸(n换)ifna弓rg哄in==协0,歇n=基3;爪en跟dY=池ey父e(挠n)幕;Y邪(:旅,1犹)=锅on魔es尝(n渐,1隶);fo美r丸k=已3:勒nY(狐k,悉2:歇k-吐1)混=Y绝(k棵-1完,1亡:k淘-2坟)+咽Y(宰k-尝1,挡2:乓k-帮1)辛;en朴d例2.8杨辉套三角鬼形红、绿两运队从销相距10朽0公里责的地糊点同扭时出罗发相抗向行沿军红队届速度搅为10糠(公里/小时)绿队事速度鞠为8(公里/小时)开始株时,通讯小员骑载摩托毫从红昌队出鞋发为健行进铁中的屡两队传递维消息铲。摩他托车圣的速尽度为60西(公里/小时)往返费于两市队之间.每遇捉一队,立即况折回下驶向训另一考队.当两对队距登离小于0.德2公里插时,摩托晌车停兰止。谱计算迟通讯董员驾慈驶摩轻托车姜跑了多扮少趟(从一待队驶逃向另位一队令为一榆趟)。例2.遍8摩托路车问感题将红队、疗绿队定和摩处托车财假设干为A、B、C三个犬点.A点初加始位饱置A=榜0,速度va=1络0漏(运动搅向右);B点初辈始位能置B=械10朴0,速度vb=8仇(运动嫂向左);C点初始遇位置C=庭0,速度vc=6提0课(f董=±收1表示运动方向)当C向右树运动答时,名C、B相遇蛋时间:tk=(B秆–享A炒)/棵(8料+6低0)当C向左扛运动冠时,A、C相遇鸦时间:tk=(B材–役A团)/拣(6惧0+