如果您无法下载资料,请参考说明:
1、部分资料下载需要金币,请确保您的账户上有足够的金币
2、已购买过的文档,再次下载不重复扣费
3、资料包下载后请先用软件解压,在使用对应软件打开
Matlab编程入门用Matlab语言编写的程序称为M文件。M文件是一个文本文件,可以用任何文本编辑器来建立和编辑,通常使用Matlab自带的M文件编辑器。在命令窗口中输入f2cs,即可执行该M文件。function输出形参列表=函数名(输入形参列表)%注释说明部分(可选)函数体语句(必须)例1编写函数文件求半径为r的圆的面积和周长。函数文件如下:function[s,p]=fcircle(r)%calculatetheareaandperimeterofacircleofradiir%r,s,p分别为圆的半径、面积以及周长s=pi*r*r;p=2*pi*r;functiony=myfun1(x)%returnsXwithrowpreservedandcolumns%flippedintheleft/rightdirection%%X=123becomes321%456654%num_col=size(x,2);%取x的列数y=x(:,num_col:-1:1);functiony=myfun2(x)%计算分段函数n=length(x);fork=1:nifx(k)<-1y(k)=x(k);elseifx(k)>=1y(k)=exp(1-x(k));elsey(k)=x(k)^3;endend函数调用的一般格式程序控制结构有三种:顺序结构、选择结构和循环结构。任何复杂的程序都由这三种基本结构组成.例:disp(X)例求一元二次方程ax2+bx+c=0的根。程序如下:functionx=myfun3(a,b,c)a=input(‘请输入系数a:');b=input('请输入系数b:');c=input('请输入系数c:');d=b*b-4*a*c;x=[(-b+sqrt(d))/(2*a),(-b-sqrt(d))/(2*a)];disp(['x1=',num2str(x(1)),',x2=',num2str(x(2))]);数据输出fprintf数据输出fprintf选择结构是根据给定的条件成立或不成立,分别执行不同的语句。Matlab用于实现选择结构的语句有if语句和switch语句。单分支结构ifexpression1(条件1)statements1(语句组1)elseifexpression2(条件2)statements2(语句组2)......elseifexpressionm(条件m)statementsm(语句组m)elsestatements(语句组)end例:编写函数,计算functiony=myfun4(x)n=length(x);fori=1:nifx(i)>1y(i)=x(i)^2;elseifx(i)>-1y(i)=1;elsey(i)=3+2*x(i);endend根据表达式的不同取值,分别执行不同的语句。Matlab首先计算expression的值,然后将它依次与各个case指令后的检测值进行比较,当比较结果为真时,就执行相应的语句组,然后跳出switch结构。例某商场对顾客所购买的商品实行打折销售,标准如下(商品价格用price来表示):price<200没有折扣200≤price<5003%折扣500≤price<10005%折扣1000≤price<25008%折扣2500≤price<500010%折扣5000≤price14%折扣输入所售商品的价格,求其实际销售价格。price=input('请输入商品价格');switchfix(price/100)case{0,1}%价格小于200rate=0;case{2,3,4}%价格大于等于200但小于500rate=3/100;casenum2cell(5:9)%价格大于等于500但小于1000rate=5/100;casenum2cell(10:24)%价格大于等于1000但小于2500rate=8/100;casenum2cell(25:49)%价格大于等于2500但小于5000rate=10/100;otherwise%价格大于等于5000rate=14/100;endprice=price*(1-rate)%输出商品实际销售价格循环结构是按照给定的条件,重复执行指定的语句。Matlab用于实现循环结构的语句有for语句和while语句forvariable=expressionstatement(循环体)end例一个三位整数各位数字的立方和等于该数本身则称该数为水仙花数。输出全部水仙花数。程序如下:form=100:999m1=fix(m/100);%求m的百位数字