1/14数据显示格式字符串变量操作图形文件的输入/输出宝石切割实验MATLAB 使用入门 III2/14计算机浮点数表示法 6.378137 103 6.378137 e +003 2.99792458 105 2.99792458e +005 例1.24近似计算地球表面积,并以不同格式输出数据R=6378.137; S=4*pi*R2format long,Sformat bank,Sformat short 短格式 5.1121e+008 format long 长格式 5.112078933958109e+008 format bank 银行格式 511207893.40数据显示格式3/14MATLABMATLAB内部常数内部常数常 数返 回 值 ans默认变量名，保存最近的结果 eps浮点数相对精度(2.2204e-016) realmax最大浮点数(1.7977e+308) realmin最小浮点数(2.2251e-308) pi圆周率（3.1416） i,j虚数单位 inf无限大 NaN不合法的数值，非数值数据显示格式4/14例1.25立方倍积问题:求作一立方体,使其体积为已知 立方体的二倍.如果已知立方体体积为V0=1,要作的立 方体体积为V=2,则所求立方体高度为 21/3在命令窗口中直接使用命令h=vpa(2(1/3),40) 数据结果为h = 1.259921049894873190666544360283296555281MATLAB具有一种称为VPA类型的数据显示格式 .VPA的全称是可变精度浮点算法（Variable precision arithmetic）.显示S的D位数方法如下R=vpa(S，D)5/14字符串变量：将字符串用单引号括起来赋值给变量拼接字符串命令: strcat(T1,T2,T3) 将整数转换成字符命令: int2str( )例1.26 三个名人Euler，Elizabeth，Plato职业分别 是mathematician，movie star，philosopher，编写 程序正确联接他们的名字和职业并输出 n1=Euler;n2=Elizabeth;n3=Plato； p1=mathematician;p2=movie star; p3=philosopher； s1=strcat(n1,-,p1), s2=strcat(n2,-,p2), s3=strcat(n3,-,p3)6/14例1.27 生肖问题08年是鼠年还是狗年 ？n=input(input n:=); S=鼠牛虎兔龙蛇马羊猴鸡狗猪; k=rem(n-4,12)+1; s=S(k); s=strcat(int2str(n), 年是, s,年)input n:=2006s= 2006年是狗年 Input n:=2008s=2008年是鼠年7/14复变函数图形设有复平面上单位圆域内变化的变量 以 u = z 的实部函数绘图并输出图形文件 r=linspace(0,1,20); theta=linspace(-pi,pi,25); z=r*exp(i*theta); x=real(z); y=imag(z); mesh(x,y,x),hold on colormap(0 0 1) mesh(x,y,-ones(size(x) axis off8/14输出图形文件 1.在图形窗口用鼠标点 击窗口左上方的菜单 栏“file”，选择下拉菜 单中的“Export”. (在高 版本中选择“Save As”)2.在对话框中选择文 件类型“bmp”格式， 将图形文件命名为 riman. 9/14宝石加工厂需要将石料内已知位置 上的一块精品切割出来.石料尺寸为 191412(cm3) ,需要切割的精品 尺寸为542(cm3),且二者的左侧 面、前面、底面相互平行,距离分别 为6cm,7cm,9cm。宝石加工问题：切割加工费用:垂直切割3元/cm2 ,水平切割5元/cm2 。试确定切割的工作流程,使切割的总费用最少。 (为了减少旋转刀具的次数，要求同向切割连续两次 后再旋转刀具)10/1419 14 12 5 4 2S1=2(1412)5 14 12 5 4 25 4 12 5 4 25 4 2 5 4 2S2=2(512)S3=2(54)11/14全部方案比较石料尺寸(cm) 石料:长宽高=a1a2a3； 精品:长宽高= b1b2b3；六种方案的切割面积计算a2a3 b1a3 b1b2a2a3 b1a2 b1b3 a1a3 b2a3 b1b2 a1a3 a1b2 b2b3 a1a2 a2b3 b1b3 a1a2 a1b3 b2b3(1,2,3) (1,3,2) (2,1,3) (2,3,1) (3,1,2) (3,2,1)1 左、右 2 前、后 3 上、下12/14a1=19;a2=14;a3=12;b1=5;b2=4;b3=2; p=1 2 3;1 3 2;2 1 3;2 3 1;3 1 2;3 2 1; f1=3;3;5;f2=3;5;3;f3=5;3;3; A=a2*a3,b1*a3,b1*b2;a2*a3,b1*a2,b1*b3;a1*a3,b2*a3,b1*b2;a1*a3,a1*b2,b1*b3;a1*a2,a2*b3,b1*b3;a1*a2,a2*b3,b2*b3; q(1)=A(1,:)*f1;q(2)=A(2,:)*f2;q(3)=A(3,:)*f1; q(4)=A(4,:)*f2;q(5)=A(5,:)*f3;q(6)=A(6,:)*f3; qmin,index=min(q); operate=p(index,:) expense=2*qminMATLAB程序(mlab0.m)13/14练习与思考题1.设地球半径为6400(km),试分析“坐地日行八万里 ”这句话是否有道理2.如何显示出圆周率的100位数字？ 3.分析复变函数计算中 z=r*exp(i*theta);x=real(z);y=imag(z); 的功能作用4. 2007年农历是丁亥年,排列天干(甲乙丙丁戊己庚辛 壬癸)/地支(子丑寅卯辰巳午未申酉戍亥)，用字符串 数组编写程序,要求输入年份,输出农历纪年14/145. 下面程序是绘抛射曲线族,分析每行命令功能alpha=pi*(1:20)/40; v0=100;g=9.8; Taim=2*v0*sin(alpha)/g; T=(0:16)*Taim/16; X=v0*T*diag(cos(alpha); Y=v0*T*diag(sin(alpha)-g*T.2/2; plot(X,Y,r)