清华数学实验第2章MATLAB程序设计方法I.pptVIP

MATLAB 程序设计方法 I;表达式——由变量、运算符、函数、数字组成 球坐标变换;例2.2 格林威治天文台建于1675年，其经度为零,纬度为北纬510。计算格林威治天文台在地心直角坐标系中的坐标并在球面上标记 ;+ 加 - 减 * 乘 .* 点乘 / 右除 ./ 点右除 \ 左除 .\ 点左除 ^ 矩阵幂 .^ 点幂;运算优先级;————命令文件————;命令文件——MATLAB命令的有序集合。 文件执行——对文件中命令进行批处理，即从第一 条命令开始按顺序执行，直到最后一条命令。如果 中间某条命令出错，则中断并输出错误信息;例2.3直线平行于Z轴沿x-y平面上的四边形移动，形成四边形柱面。利用矩阵方法绘制四边形柱面. ;x=[-1,1,1,-1,-1]; %确定四边形X坐标 y=[-1,-1,1,1,-1]; %确定四边形Y坐标 X=[1;1]*x ;Y=[1;1]*y; %创建X和Y坐标矩阵 Z=[0;1]*ones(1,5); %创建柱面Z坐标矩阵 mesh(X,Y,Z) %绘图 colormap([0 0 0]) %设定黑色 axis off %去掉坐标框架;例2.7 成都地理位置:北纬30o,东经104o;北京地理位置:北纬40o,东经116o。计算两城市距离近似值。 ;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)) %计算并显示距离 ;条件控制; 循环控制; 错误控制;终止运行控制;1. for 循环 语法： for 变量 =初值：步长：终值 可执行语句 end;例2.11 程序功能如下：把边长为2以原点为中心的正方形旋转pi/24，将其压缩(r=0.89),重复24次并绘图 ;例2.12空间四面体底面是一三角形，将三角形三个顶点与四面体顶部的顶点相连接就形成三角形锥面。利用MATLAB的帧动画方法演示一个四面体旋转过程。 ;2. while循环 将循环体中语句循环执行直到条件不成立为止。 语法： while 条件表达式 可执行语句 end; 例2.14 3n + 1 问题. ;3. continue 命令;如果一个正整数 n 只能被1和它自身整除，则称这个数为素数(或质数)??例如：2、3、5、7、11、······;例2.17海伦公式;立方倍积问题实验。作一立方体,使其体积等于单位正方体的两倍.对边长取不同位数,显示计算结果误差。 ;练习与思考题