MATLAB程序设计方法.pptVIP

处理数据程序如下： score=0.1*datas(:,1)+0.3*datas(:,2)+0.6*datas(:,3); N=length(score); %统计总人数 II=find(score60);F0=length(II); %统计不及格人数 II=find(score=60score70);F1=length(II); II=find(score=70score80);F2=length(II); II=find(score=80score90);F3=length(II); II=find(score=90);F4=length(II); fansu=[F0,F1,F2,F3,F4] %显示各分数段人数 format bank R=fansu/N %显示各百分比 MATLAB程序设计方法全文共34页，当前为第31页。 牟合方盖的几何图形 绘制柱面x2 + y2 = R2与柱面x2 + z2 =R2所围立体在x-y平面上半部分曲面。由第二个方程解出z，得 利用极坐标变换： h=2*pi/100; t=0:h:2*pi;r=0:0.05:1; x=r*cos(t); y=r*sin(t); z=sqrt(1-x.^2); meshz(x,y,z) MATLAB程序设计方法全文共34页，当前为第32页。 练习与思考题 1.叙述古希腊三大难题之一的化圆为方问题,并用实验数据显示误差。 2. 用一句话说明循环控制语句的特点。 3.说明for循环while循环语句的相同和不同之处。 4.修改“试商法判断素数”的程序，产生100以内的全部素数。 5.用help命令获取primes的使用方法，用求素数函数primes()产生100以内的全部素数。 MATLAB程序设计方法全文共34页，当前为第33页。 练习与思考题 1.修改摩托车问题程序文件为函数文件，只输出摩托车跑了多少趟数据。 2.hanoi问题中,n个盘需要移动多少次才可以完成任务？ 3.是否可用mesh()方法绘正六面体 4. 分析正六面体绘图方法 B=[0 0 0;1 0 0;1 1 0;0 1 0;0 0 1;1 0 1;1 1 1;0 1 1] fac=[1 2 3 4;1 2 6 5;1 4 8 5;7 8 5 6;7 3 2 6;7 3 4 8] patch(faces,fac,vertices,B,faceColor,c) MATLAB程序设计方法全文共34页，当前为第34页。 */32 数学科学学院：汪小平 wxiaoping325@163.com 第二章 MATLAB程序设计方法 表达式和运算符 命令文件(程序文件) 程序设计中的流程控制 函数文件的编写格式 ? ? ? ? * MATLAB程序设计方法全文共34页，当前为第1页。 表达式——由变量、运算符、函数、数字组成 球坐标变换 ————表达式和运算符———— 用赋值语句创建变量 变量名=数据 变量名=表达式 x=R*cos(theta)*cos(fai); y=R*cos(theta)*sin(fai); z=R*sin(theta); ——纬度 ——经度 MATLAB程序设计方法全文共34页，当前为第2页。 例2.1 格林威治天文台建于1675年，其经度为零,纬度为北纬510。计算格林威治天文台在地心直角坐标系中的坐标并在球面上标记 [X,Y,Z]=sphere(24); mesh(X,Y,Z),hold on fai=0; theta=51*pi/180; x0=cos(theta)*cos(fai); y0=cos(theta)*sin(fai); z0=sin(theta); plot3(x0,y0,z0,ro,linewidth,2) view(124,30) MATLAB程序设计方法全文共34页，当前为第3页。 例2.2 四边形柱面绘图原理和方法。 设四边形顶点为(-1, -1), (1, -1), (1, 1), (-1, 1). 设柱面高为1,创建Z坐标矩阵 创建柱面的X坐标和Y矩阵 MATLAB程序设计方法全文共34页，当前为第4页。 cylinder(R，N)创建单位高度绕Z轴旋转曲面，R为母线,N表示圆圈上点数,N缺省默认值为 20。 Cylinder %单位高度二十边形柱面 [x,y,z]=cylinder([1,1],4); mesh(x,y,z)