MATLAB应用第四章-程序设计初步知识研讨.pptVIP

第四章 MATLAB程序设计初步;4.1 M文件简介;启动M文件;4.2 脚本文件和函数文件;edit;脚本文件示例 % 注释行 % M文件示例 % “flower petal” % 以下为代码行 % 计算 theta=-pi:0.01:pi; rho(1,:)=2*sin(5*theta).^2; rho(2,:)=cos(10*theta).^3; rho(3,:)=sin(theta).^2; rho(4,:)=5*cos(3.5*theta).^3; for k=1:4 % 图形输出 subplot(2,2,k) polar(theta,rho(k,:)) end disp(程序运行结束！);Eg 4-2续;4.2.2 M函数文件 ;第一步：在editor内输入下列命令 行并存为“average.m” 001 function y=average(x) 002 % AVERAGE 求向量元素的均值 003 % 语法： 004 % Y=average(X) 005 % 其中，X是向量，Y为计算得到向量元素的均值 006 % 若输入参数为非向量则出错 007 008 % 代码行 009 [m,n]=size(x) 010 % 判断输入参数是否为向量 011 if(~((m==1)|(n==1))|(m==1n==1)) 012 % 若输入参数不是向量，则出错 013 error(Input must be a vector) 014 end 015 % 计算向量元素的均值 016 y=sum(x)/length(x);函数定义行 001 function y=average(x) 包括 关键字function 函数输出参数y 函数名称average 函数输入参数x 函数名称定义要求 必须以字符开头，后面可以用字符、数字和下划线的组合构成函数名称 MATLAB对函数名称的长度有限定 函数的M文件名称最好和函数名称保持一致，若不一致，则调用函数时需要使用文件名称而非函数名称;在线帮助 M函数文件的在线帮助为紧随在函数定义行的注释行 注意，这些%符号须顶格输入，否则matlab不认 002 % AVERAGE 求向量元素的均值 003 % 语法： 004 % Y=average(X) 005 % 其中，X是向量，Y为计算得到向量元素的均值 006 % 若输入参数为非向量则出错;若在MATLAB命令行窗口中键入指令lookfor average 则显示 ;M语言代码 008 % 代码行 009 [m,n]=size(x) 010 % 判断输入参数是否为向量 011 if(~((m==1)|(n==1))|(m==1n==1)) 012 % 若输入参数不是向量，则出错 013 error(Input must be a vector) 014 end 015 % 计算向量元素的均值 016 y=sum(x)/length(x) 代码行需要完成具体的算法，实现用户的具体功能;2 M函数的类型; date ans = 20-Sep-2008 date=@()datestr(now); date() ％必须有空括号 ans = 20-Sep-2008 15:38:56 date date = @()datestr(now);function[avg,med]=newstats(u) % 主函数 % NEWSTATS计算均值和中间值 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);％将v按升序进行排列 if rem(n,2)= =1 m=w((n+1)/2); else m=(w(n/2)+w(n/2+1))/2; end;局部变量：在M函数内部声明并使用的变量 只在本函数范围内有效 仅能在函数调用执行期间被使用，一旦函数结束运行，则变量占用的内存空间自动释放，变量的数值也不存在了 例4-22 局部变量的例子 function local % LOCAL察看局部变量的例子 x=rand(2,2); y=zeros(2,2); z=函数中的变量; u={x,y,z}; dis