Matlab第5次课 2、三维曲线和曲面 3、M文件及程序调试
Matlab第5次课 • 2、三维曲线和曲面 • 3、M文件及程序调试
255三维曲线和曲面 1、空间曲线绘制plot3 格式:plot(x2y;z’r) 例子 z=0:0.14*pi; X-COS(Z y=sin(z) gure plot(x, y, z)
2.5.5三维曲线和曲面 • 1、空间曲线绘制—plot3 • 格式:plot(x,y,z,’r’) • 例子 • clc • z=0:0.1:4*pi; • x=cos(z); • y=sin(z); • figure • plot3(x,y,z)
15 10 0.5 0.5 -0.5 -0.5
2、空间曲面的绘制 函数mesh和surf用来绘制三维曲面。 X=8:0.5:8; y-X Ⅹ1=ones(SiZe(y)* Yl-y*ones(size(x R=sqrt(Xl. * Xl+Yl. Y1) z=sin(r). R mesh(z
2、空间曲面的绘制 • 函数mesh和surf用来绘制三维曲面。 • clc • x=-8:0.5:8; • y=x'; • X1=ones(size(y))*x; • Y1=y*ones(size(x)); • R=sqrt(X1.*X1+Y1.*Y1); • z=sin(R)./R; • mesh(z)
26M文件及程序调试 26.1主程序文件 注释,以“%”对本程序的特性以及关键的部 份进行说明。 1)一般以clc、 clear、 close all等语句开始,清 除工作空间中原有的变量和图形,以避免其它 已执行的程序残留数据对本程序的影响 2)程序的主体 如果程序中要用到全局变量,要在这里说明 格式: global变量名1变量名2变量名3
2.6 M文件及程序调试 • 2.6.1 主程序文件 • 注释,以“%”对本程序的特性以及关键的部 份进行说明。 • 1)一般以clc、clear、close all等语句开始,清 除工作空间中原有的变量和图形,以避免其它 已执行的程序残留数据对本程序的影响。 • 2)程序的主体 • 如果程序中要用到全局变量,要在这里说明。 • 格式:global 变量名1 变量名2 变量名3…
般地,若程序中用到流程控制语句, 都要缩进,并且要注意end的对应 3)文件命名要应按前面介绍的标志符来 命名。文件名中不能有汉字,也不能以 数字开头
• 一般地,若程序中用到流程控制语句, 都要缩进,并且要注意end的对应。 • 3)文件命名要应按前面介绍的标志符来 命名。文件名中不能有汉字,也不能以 数字开头
函数文件及其应用 函数文件是M文件重要的组成部分,M语言函数 文件能够接受用户的输入参数,进行计算,并将 计算结果作为函数的返回值返回给调用者 格式 function[x,y= myfun(a,b,c)--数定义行 %函数功能说明行 %帮助文件,说明如何使用该函数 函数体
函数文件及其应用 • 函数文件是M文件重要的组成部分,M语言函数 文件能够接受用户的输入参数,进行计算,并将 计算结果作为函数的返回值返回给调用者。 • 格式 • function [x,y]=myfun(a,b,c)-----函数定义行 • %函数功能说明行 • %帮助文件,说明如何使用该函数 • 函数体
例子 计算一数组的平均值及找出其中间值。 主程序: cIc clear close all n=input(n=) X=1:N mean, mid]=newstats(x) omean %omid X1=[5.8,1,3.5,2263,3.8] meanl, mid1]-newstats(x1
例子 • 计算一数组的平均值及找出其中间值。 • 主程序: • clc • clear,close all • N=input('n=') • x=1:N; • [mean,mid]=newstats(x) • %mean • %mid • x1=[5.8,1,3.5,2.2,6.3,3.8] • [mean1,mid1]=newstats(x1)
函数 function avg, med=newstats(u) °n= ength(u avg-mean(u, n) med=median(u, n) function a=mean(v, n) a=sum(v)/n; function m=median(v, n) w=sort(v) if rem(n, 2)==1 mw(n+1)2) else m(w(n/2)+wn/2+1)/2 ent
• 函数 • function [avg,med]=newstats(u) • n=length(u); • avg=mean(u,n); • med=median(u,n); • function a=mean(v,n) • a=sum(v)/n; • function m=median(v,n) • w=sort(v); • if rem(n,2)==1 • m=w((n+1)/2); • else • m=(w(n/2)+w(n/2+1))/2; • end
程序调试 M文件编辑器不仅是个文件编辑器,而 且是个可视化的调试开发环境。 在这里可以对主程序以及函数文件进行 调试,以排查程序的错误。 般地,程序错误有两类: 语法错误 ·2、运行错误
程序调试 • M文件编辑器不仅是个文件编辑器,而 且是个可视化的调试开发环境。 • 在这里可以对主程序以及函数文件进行 调试,以排查程序的错误。 • 一般地,程序错误有两类: • 1、语法错误; • 2、运行错误