《MATLAB编程与作图》由会员分享,可在线阅读,更多相关《MATLAB编程与作图(47页珍藏版)》请在金锄头文库上搜索。
1、数学实验MATLAB编程与作图数据的输出格式 MATLAB用十进制数表示一个常数,具体可 采用日常记数法和科学记数法两种表示方法 。 数据输出时用户可以用format命令设置或改 变数据输出格式。format命令的格式为: format 格式符 注意,format命令只影响在屏幕上的显示数 据输出格式,而不影响数据的计算和存储 ,MATLAB的数据存储和运算总是以双精度 进行的。可用菜单File/Preferences中的Command Window 里的Numarical Format来设定数据显示方式,有short( 小数点后4位), long(15位数字), bank(小数点后2位),
2、short e(5位科学计数法), long e(15位科学计数法), rat(最接近的有理数).也可在Command Window中用命令设置,例如:format longformat short :5位定点格式,如 1.4142 format long :15位定点格式,如 1.41421356237310 format short e :5位浮点,如1.4142e+000 format long e :15位浮点,如 1.414213562373095e+000 format short g :5位定点或浮点取优,如1.4142 format long g:15位定点或浮点取优,如 1.
3、4142135623731 format hex :十六进制数,如3ff6a09e667f3bcd format bank :货币格式,如1.41 format rat :分数之比,如1393/985 format +:若为正数则显示+,若为负数则显示-,若为0 则显示空如 format short %设置定点格式,显 示5位 pians =3.1416 format long %设置定点格式,显示 15位 pians =3.14159265358979关系与逻辑运算1. 关系操作符2. 逻辑运算符MATLAB流程控制 循环语句 for, while 条件转移 if end, if elsei
4、f else end 开关语句 switch case 注释语句 % 中断语句 break 暂停语句 pause 回显语句 echo on/off选择结构1. if语句在MATLAB中,if语句有3种格式。1.单分支if语句语句格式为:if 条件语句组end当条件成立时,则执行语句组,执行完之 后继续执行if语句的后继语句,若条件不成立 ,则直接执行if语句的后继语句。2.双分支if语句语句格式为:if 条件语句组1 else语句组2end当条件成立时,执行语句组1,否则执行语句组 2,语句组1或语句组2执行后,再执行if语句的 后继语句。例1. 计算分段函数值程序如下: x=input(请输
5、入x的值:);if x=10y=cos(x+1)+sqrt(x*x+1);elsey=x*sqrt(x+sqrt(x);endy 请输入x的值:y =47.1893121.m3.多分支if语句 语句格式为: if 条件1语句组1else if 条件2语句组2else if 条件m语句组melse语句组m+1 end例2. 输入一个字符,若为大写字母,则输出其 后继字符,若为小写字母,则输出其前导字符 ,若为数字字符则输出其对应的数值,若为其 他字符则原样输出。c=input(请输入一个字符,s);if c=A end end a=1.0000 0.5000 0.3333 0.2500 0.20
6、000.5000 0.3333 0.2500 0.2000 0.16670.3333 0.2500 0.2000 0.1667 0.14290.2500 0.2000 0.1667 0.1429 0.12500.2000 0.1667 0.1429 0.1250 0.1111例5. 已知求y的表达式,当n=100时,求y的值。 在实际MATLAB编程中,为提高程序的执行速度,常 用向量运算来代替循环操作:n=100; i=1:n;f=1./i.2;y=sum(f)程序如下:y=0;n=100;for i=1:ny=y+1/i/i;endyy =2047/1252 for语句更一般的格式 :fo
7、r 循环变量=矩阵表达式循环体语句end执行过程是依次将矩阵的各列元素赋给 循环变量,然后执行循环体语句,直至 各列元素处理完毕。实际上,“表达式1: 表达式2:表达式3”是一个仅为一行的矩阵 (行向量),因而列向量是单个数据。例6. 已知5个学生4门功课的成绩,求每名学生 的总成绩。 程序如下:try.ms=0; a=65,76,56,78;98,83,74,85;76,67,78,79;98,58,42,73;67,89,76,87 for k=as=s+k; end disp(s);a =65 76 56 78 98 83 74 85 76 67 78 79 98 58 42 73 67
8、 89 76 87 275 340 300 271 319 2. while语句while语句一般用于不能事先确定循环次数 的情况 while语句的一般格式为:while (条件)循环体语句end其执行过程为:若条件成立,则执行循环 体语句,执行后再判断条件是否成立,如 果不成立则跳出循环(如图3.5所示)。 While循环和for循环的区别在于,while循环结构的 循环体被执行的次数不是确定的,而for结构中循环体的执行次数是确定的。例7. 用循环求解例8. 用循环求解最小的 m例9. 用while循环语句编写一个计算1000以 内的Fibonacci数:5.mf=1 1;i=1;whil
9、e f(i)+f(i+1)1000f(i+2)=f(i)+f(i+1);i=i+1;endf=1 1 2 3 5 8 13 21 34 55 89 144 233 377 610 987例10. 产生一个长度为的调和数列。程序如下:X=zeros(1,6);i=1;while i=6X(i)=1/i;i=i+1;endformat rat %使用分数表示数值disp(X)Columns 1 through 5 1 1/2 1/3 1/4 1/5 Column 6 1/6 与循环结构相关的还有一个break语句,当在循 环体内执行到该语句时,程序将跳出循环。该 语句一般与if语句配合使用。例11
10、. 用while语句实现例5。6.m程序如下:y=0; i=1;while 1f=1/i/i;y=y+f;if i=100break;endi=i+1;endy3. 循环的嵌套 如果一个循环结构的循环体又包括一个循环 结构,就称为循环的嵌套,或称为多重循环 结构。可以按照嵌套层数,分别叫做二重循 环、三重循环等。处于内部的循环叫作内循 环,处于外部的循环叫作外循环。 在设计多重循环时,要特别注意内、外循环 之间的关系,以及各语句放置的位置,不要 搞错。例12. 用筛选法求某自然数范围内的全部 素数。程序如下:m=input(m=); p=2:m;for i=2:sqrt(m)n=find(re
11、m(p,i)=0p(n)=;endpm=100 p =Columns 1 through 5 2 3 5 7 11 Columns 6 through 10 13 17 19 23 29 Columns 11 through 15 31 37 41 43 47 Columns 16 through 20 53 59 61 67 71 Columns 21 through 25 73 79 83 89 97 7.m小结:MATLAB程序的基本组成结构 说明 清除命令:清除workspace中的变量和图形( clear,close) 定义变量:包括全局变量的声明及参数值的设定 逐行执行命令:指MA
12、TLAB提供的运算指令或工具箱 提供的专用命令 控制循环 : 包含for,if then,switch,while等语句 逐行执行命令 end 绘图命令:将运算结果绘制出来 当然更复杂程序还需要调用子程序,或与simulink以 及其他应用程序结合起来。二维图形三维图形图形处理实例特殊二、三维图形MATLAB提供多种图形功能,使你的数据或函数可 视化,使数据不再是枯燥乏味的。使用MATLAB的图形 函数,可以绘制二维或三维的数据图形和函数图形, 如数据的散点图、直方图、茎干图、饼图、阶梯图和 面积图等。使数据可视化的基本步骤是:1.准备好数据;2.选择适用的绘制图形函数;3.选择窗口和位置;4
13、.编辑图形标注和说明;5.输出或保存图形。MATLAB提供了丰富的绘图功能help graph2d可得到所有画二维图形的命令 help graph3d可得到所有画三维图形的命令MATLAB作图是通过描点、连线来实现的,故在 画一个曲线图形之前,必须先取得该图形上的 一系列的点的坐标(即横坐标和纵坐标),然 后将该点集的坐标传给MATLAB函数画图. 命令为: plot(X,Y,S)plot(X,Y) 画实线 plot(X,Y1,S1,X,Y2,S2,X,Yn,Sn) 将多条线画在一起X,Y是向量,分别表示点集的横坐标和纵坐标线型y 黄色 . 点 - 连线 m 洋红 o 圈 : 短虚线 c 蓝绿色 x x-符号 -. 长短线 r 红色 + 加号 - 长虚线1.曲线图例 在0,2 用红线画sin x,用绿圈画cos x. x=linspace(0,2*pi,30); y=sin(x); z=cos(x); plot(x,y,r,x,z, g0)解 MATLAB liti12.符号函数(显函数、隐函数和参数方程)画图 (1) ezplotezplot(x(t),y(t),tmin,tmax) 表示在区间tminttmax绘制参数方程 x=x(t),y=y(t)的函数图.ezplot(f(x),a,b) 表示在axb绘制显函数f=f(x)的函数图.ezplot(f(x,y),xm