您所在位置： &
&nbsp&&nbsp&nbsp&&nbsp
matlab程序语言设计第4章.ppt 159页
本文档一共被下载：
次 ,您可全文免费在线阅读后下载本文档。
下载提示
1.本站不保证该用户上传的文档完整性，不预览、不比对内容而直接下载产生的反悔问题本站不予受理。
2.该文档所得收入(下载+内容+预览三)归上传者、原创者。
3.登录后可充值，立即自动返金币，充值渠道很便利
需要金币：200 &&
matlab程序语言设计第4章
你可能关注的文档：
··········
··········
再次产生更大的向量进行检验 && v=2.1+6*rand(1,100); var2=std(v).^2 var2 =
2.9614 && myvar2=variance(v) myvar2 =
2.9614 这说明调试后的程序正确。
语言结构与调试函数
　　这一节将给出MATLAB的语言结构与调试函数。先简要列出这些函数，如表4.3所示，然后分类列出各个函数的使用说明。
语言结构和调试
MATLAB编程语言 　　1．function，script 　　功能：MATLAB的函数文件和脚本文件。
　　说明： 　　以function开头的M文件(文件名后缀为M)是MATLAB的函数文件。MATLAB的M文件分为两类：函数文件和脚本文件(也称程序文件)，其区别在于函数文件应该通过输入/输出变量传递参数，函数文件中使用的变量不会进入基本工作空间；而脚本文件中的变量会完全带入到基本工作空间，基本工作空间中的任意变量也可为脚本文件中的语句所用。
　　函数文件必须以function开头，而脚本文件可直接写出各种MATLAB语句。例如有一函数文件stat.m，其内容为 function [mean,stdev]=stat(x) n=length(x) mean=sum(x)/n; stdev=sqrt(sum(x-mean).^2/n); 　　这时函数名为stat，它应与文件名同名。函数中使用的变量都是局部变量，它只有借助于输入/输出变量才能读取或返回到基本工作空间。
　　我们还可以在同一个函数文件中定义子函数，但子函数只能在这一个函数文件中才能调用，它不能为其它函数文件或脚本文件所调用。例如，上述stat函数可采用自带子函数方法编写，这里return为结束引用命令。 function [mean,stdev] = stat(x) n = length(x); mean = avg(x,n); stdev = sqrt(sum1(x,avg(x,n))/n); return
function s=sum1(x,y) s=sum((x-y).^2); return ? function mean = avg(x,n) mean = sum(x)/n; return
2．global 功能：定义全局变量。 格式： global
　　说明： 　　global X Y Z可将变量X、Y、Z定义成全局变量。一般而言，每个MATLAB函数都有自己的局部变量，它们与其它函数中的局部变量无关，也与基本工作空间中的变量无关，因此它们可与基本工作空间和其它函数文件采用同名的变量，其内容之间也没有关系。函数与基本工作空间之间的参数主要依靠输入/输出变量传递。然而，如果将某一变量声明为全局变量，则只要在某函数中也将它声明为全局变量，那么在该函数中就可以存取这一变量。例如，在基本工作空间中，已声明一矩阵a为全局变量，则可编写下列函数文件：
function y=abc(x) global a m=mean(a)； y=x*m;
　　注意，利用clear global variable可从工作空间中清除指定的全局变量variable；利用clear variable可从当前工作空间中清除变量variable的全局连接。当然，这不会影响到全局变量的值。
　　3．persistent 　　功能：定义永久变量。 　　格式： 　　persistent
X Y Z 　　说明： 　　persistent X Y Z可将变量X、Y、Z定义成永久变量。persistent函数只用于函数文件中，使在每次调用时保持变量的值不变。 　　永久变量只有在从内存中清除M文件或已改变M文件时才能清除。为此，要想一直保持某个函数中定义的永久变量，应采用mlock锁定相应的M文件。
　　4．nargchk 　　功能：检查输入变量数。 　　格式： 　　　　msg=nargchk(low，high，number)
　　说明： 　　nargchk通常在M文件内使用，用来检查输入变量数的有效性。 　 　　msg=nargchk(low，high，number)用于检查输入变量数number(可由nargin函数获得)是否有效，当number小于low或大于high时，则得到一个包含出错信息的字符串msg；否则，msg为空矩阵。例如，有一函数foo：function f=foo(x,y,z) 　　error(nargchk(2,3,nargin)) 则在MATLAB下输入foo(1)时会产生出错信息： && foo(1) ??? Error using ==&gt
正在加载中，请稍后...一些MATLAB的编程规范总结1.0版
已有 12037 次阅读
|个人分类:|系统分类:
好的编程习惯，可以提高编程效率，不仅可以使代码容易修改，也容易给别人看懂，便于交流。我们不仅要写出“给机器读懂的代码”，也写出“给人看得懂的代码”。本文根据一些目前搜索到的文献和自己的一些使用心得，整理出这个文档，大家可以根据经验提出自己的心得，相互促进，共同提高。&1，命名规则（1）变量a，使用英文。变量名最重要的原则就是，一看就知道这个变量是什么意思。由于程序大部分是字母，所以最好还是有英文字母表示意思比较靠谱。如果英文稍微差点，一下子不知道变量的英文怎么写，可以查字典，现在在线词典很多，这样还可以使自己的英文水平提高，二则可以提高程序可读性。比如，我们要表达“风速”，“windSpeed”比“fengsu”要好，时间久了，你再看“fengsu”的的时候，还必须得拼一下，搞不好就弄成“风俗”。b，变量名应该以小写字母开头的大小写混合形式。例如，maturityDay,anthesisDay。注意：有些人喜欢用“下划线”把词语隔开，maturity_day,athesis_day,这样也一目了然，但是在MATLAB中，不推荐这么用，因为下划线会在Tex解释程序中为下标的转换符，比如：xlabel(“maturity_day”)，在图中会显示成“maturityday”。c，用大范围意义的变量名应该带有意义的名称，但下范围的变量名可以用短变量名。例如：stepSize：大范围意义 & & dt（δt）：小范围d，前缀n用在作为数值对象申明例如：nRowse，遵循的一个有关复数变量的惯例所有变量名要么为单数形式，要么为复数形式。两个变量只是最后相差一个字母s加以区别的情况应该避免，apple和apples区别2个变量，看着头大。可以用一个代表复数的标示，如：appleArrayf，代表单个实体数据的变量可以加以后缀No或者是前缀i例如，tableNo,employeeNo,g，循环变量应该以i,j,k为前缀例如，iFiles,jPositions。不使用i,j,的原因是：因为这两个在matlab里面是虚数。h，嵌套循环，循环变量要以字母顺序for iFile = 1：nFilesfor jPosition = 1：nPositions…endendi，避免否定式的布尔变量命名。例如：命名，在使用判断的时候，~isNotFound，搞半天才知道啥意思。所以不适用否定式布尔变量命名。j，缩写形式，即使是通常的大写缩写，也应该于小写字母混合使用。采用：html，isUsaSpecific，checkTiffFromat（）避免使用：hTML，isUSASpecific，checkTIFFFormat（）&（2）常数a，命名常数（包括全局变量）应该采用大写字母，用下划线分割单词。例如：MAX_ITERATIONS，COLOR_REDb，参数可以以某些通用类型名作为前缀COLOR_RED，COLOR_GREEN，COLOR_BLUE&（3）结构体a，结构体的命名应该以一个大写字母开头这是区别一般变量，例如：ParameterSetb，结构体的命名应该是隐性的，并且不需要包括字段名用Segment.length.避免Segment.segmentLength&（4）函数a，函数名应该采用小写字母，且将函数名与它的m 文件名保存为相同b，函数名应该有具体的意义避免短的函数名，这经常使得其名字含糊不清。例如：采用： computetotalwidth（）避免： compwid（）c，单输出变量的函数可以根据输出参数命名例如，mean,stdd，没有输出变量或者返回值为句柄的函数应该根据其实现的功能命名例如：plote，前缀get/set 应该作为访问对象或者属性的保留前缀例如：getobj（.）f，前缀 compute应该用在计算某些量的函数的地方例如：ComputWeightedAverage()g，前缀find可以用在那些具有查询功能的函数的地方例如：FindOldestRecord(.);FindHeaviestElement(.);h，前缀initialize可以用在对象或者是概念（concept）建立的地方例如：InitializeProblemState(.);i，前缀is应该用在布尔函数的命名的地方例如：IsOverPriced(.)j，避免无意识地覆盖有时候我们取的名字，可能在MATLAB中含有这个函数名了，可以用exist检查是否含有了。k，避免缩写采用：ComputeArrivalTime(.)避免：CompArr(.).l，考虑使得名字可以拼读&&&2，文件与程序结构（1）M文件a，模块化编写一个大程序的最好的方法是将它以好的设计分化为小块（通常采用函数的方式）。这种方式通过减少为了理解代码的作用而必须阅读的代码数量使得程序的可读性、易于理解性和可测试性得到了增强。超过编辑器两屏幕的代码都应该考虑进行分割。并且设计规划很好的函数也使得它在其他的应用中可用性增强了。b，确保交互过程清晰函数通过输入输出参数以及全局变量与其他代码交互通信。使用参数几乎总是比使用全局变量清楚明了。采用结构可以避免那种一长串儿的输入输出参数的形式。c，分割所有的子函数和所有的函数都应该只把一件事情做好。每个函数应该隐藏（hide）一些东西。d，利用现有的函数开发一个有正确功能的、可读的、合理灵活性的函数在一项有重大意义的任务。或许寻找一个现成的提供了要求的部分、甚至全部功能的函数应该更快也更具有正确性。e，任何在多个m文件中出现的代码块都应该考虑用函数的形式封装起来f，子函数只被另外一个函数调用的函数应该作为一个子函数写在同一个文件中。这使得代码更加利于理解与维护。g，测试脚本为每一个函数写一个测试脚本。这样可以提高初期版本的质量和改进版本的可靠性。&（2）输入与输出a，编写输入/输出模块输出要求可以无需特别注意就可以根据变化而改变，输入的格式与内容根据变化的时候经常很混乱。找到处理输出的地方进行改善，提高其可维护性。避免将输入/输出部分的代码与计算功能的代码混淆在一起，单个函数的预处理的时候除外。各种功能混合的函数的可再用性一遍很小。b，格式化输出使得其易于利用如果输出很大可能是人工阅读，那么就让输出采用易于越多的描述性的方式。如果输出更多的可能是通过其他软件调用而不是人，那么应该使得输出易于解析。如果这以上两种情况都很重要，将输出表达成易于解析的格式，并编写一个格式化输出的函数用来产生一个人工可读的输出版本。&&3，基本语句（1）变量与常数a，变量不应该重复使用通过确保所有的概念都只有唯一的意义可以加强代码的可读性，以及通过消除误解的定义可以减少错误的可能。b，同种类型的相近的变量可以在同一个语句中定义c，不相近的变量应该不要在同一个语句中定义d，注意在文件开始部分的注释中为重要变量编写文档MATLAB不需要变量申明，这种信息就可以在注释中提供。示例：％pointArray Points are in rows with coordinates in columns.e，注意在语句行注释的最后为常数编写文档示例： THRESHOLD = 10；％Maximumnoise level found by experiment。&（2）全局变量a，应该尽量少地使用全局变量&（3）循环语句a，循环变量应该在循环开始前立即被赋值。这可以提高循环的速度，有助于防止循环没有执行所有的可能索引而产生的虚假值。示例：refult= zeros( nEntries，1)；forindex = 1：nEntriesresult（index）= foo（index）；endb，在循环中应该尽量少用break与 continue这些结构可以与goto 相比较，只有当他们可以证明用这些结构可以比他们相应的结构化部分有更好的可读性的时候，才可以使用。c，在嵌套式循环的时候应该在end行加上注释在长的嵌套循环的end命令行添加注释可以有助于弄明白哪些语句在那个循环体内、在此处之前已经完成了哪些功能。（4）条件语句a，应该避免复杂的条件表示式，而采用临时逻辑变量进行替代通过对表达式指定逻辑变量，使得程序更能够自为文档，使得程序结构更易于阅读与调试。示例：避免使用：if（value&=lowerLimit）& (values&=upperLimit)&~ismember(value,valueArray)…………end而应该用如下的方式代替：isValid = （value=lowerLimit）& (values&=upperLimit)；isNew =~ismember(value,valueArray)if ( isValid &isNew)…………endb，在if else 结构的时候，发生较频繁的事件应该放在if 部分，例外情况放在else部分这样通过将例外情况排除在常规执行路径之外可以提高程序的可读性。示例：fid = fopen（fileName）；if（fid~ = -1）…else…endc，条件表达式if 0 是应该避免的，除非在对临时程序块进行注释的时候如果确信表达式在程序正常执行的时候不会发生，首选的方法是采用编辑器的块注释。d，一个 switch语句应该包含otherwise条件将otherwise情况遗漏在外是一种通常错误，这或许会导致不可预测结果。正确示例：switch（condition）case ABC处理语句；case DEF处理语句；otherwise处理语句；ende，switch 变量应该通常是字符串字符串在这种情况下能够很有效，通常他们比采用列举值的形式意义更丰富。（5）其它a，采用附加说明MATLAB对于操作运算有个优先级的文档，但是谁愿意记住它们的具体内容呢？如果在某些地方有任何疑问，采用附加说明使得表达清楚。特别是在扩展的逻辑表达式的时候尤其有用。b，尽量在表达式中少用数字。可能会改变的数字应该用常数代替如果一个数字它的本身没有明确的意义，采用将它命名为常数可以加强程序的可读性。并且，改变参数的定义比改变文件中所有的相应出现地方的数字要容易得多。c，浮点常数应该在小数点前面写上一个阿拉伯数据这是坚持数学习惯的语法要求，而且，0.5 比.5 更具有可读性，因为.5 很有可能被误认为是整数5。d，浮点数的比较应该要小心二进制表达可能导致麻烦，如下面的例子所示：shortSide = 3；longSide = 5；otherSide = 4；longSide^2 ＝＝（shortSide^2+otherSide^2）ans =1scaleFactor = 0.01；（scaleFactor*longSide）^2 = =((scaleFactor*shortSide)^2+(scaleFactor*otherSide)^2ans =0e，close all在一般的M文件开头，写这三个东西，可以初始化MATLAB，把打开的表，工作空间清理。&4，排版、注释与文档（1）排版a，应该将代码内容控制在前80列之内对于一个编辑器、终端仿真器、打印机、调试器以及文件的通常列数是80列，因此通常几个人的程序共享的时候，大家通常将内容控制在前80列之内。在程序员之间传递文件的时候，避免无意识的分行可以增强程序代码的可读性。MATLAB的M文件编排，有这个线提示。b，在恰当的地方应该将行进行切分当语句长度超过80 列的限制的时候应该切分行。通常：在一个逗号或者空格之后进行断开；&在一个操作符之后断开；在表达式开始前的地方重新开始新的一行；示例：totalSum = a +b +c +…d+e；function (param1，param2，…param3)setText（['Long line split'…'into two parts.']）；c，基本缩排应该是3或者4个空格好的缩排或许是唯一的一个展现程序结构的好方法。1个空格是缩排太小而不能够强调出代码的逻辑分层；2个空格的缩排在为了减少因为嵌套循环超过80列而切分行的断裂的时候被建议采用，而MATLAB通常没有太多太深的循环嵌套。大于4个空格的缩排使得因为行切分的机会增大而使得代码的可读性变差。4个空格的缩排是MATLAB编辑器的缺省设置，在以前的一些版本，缺省缩排是3个空格。d，应该与MATLAB编辑器的缩排一致MATLAB编辑器提供了使得代码结构清晰的缩排，并且与C++与Java推荐使用的缩排方式相一致。e，通常情况下，一行代码应该只包含一个可执行语句f，短的单个if，for 或者while语句可以写在一行。if（condition），statement；endwhile（condition），statement；endfor iTest = 1：nTest，statement；endg，空白空格在 =，&，与| 前后加上空格在指定的字符前后加上空格可以增强其可视化的分割提示，明显地将语句左右两部分分开。在二值逻辑操作符前后加上空格可以使得复杂的表达式清晰。示例：simpleSum = firstTerm+secondTerm；h，常规的操作符两边可以加上空格这种方式是有争议的。部分人认为它可以增强其可读性。示例：simpleAverage= (firstTerm + secondTerm) / two；1 : nIterationsi，逗号后面可以加上空格这些空格可以增强可读性。有些程序员为了避免切分行而不采用这种方式。示例： foo(alpha, beta, gamma)，也可以foo(alpha,beta,gama)j，分号或者同一行多条指令的逗号之后应该加上一个空格字符空格加强可读性。示例：if （pi&1），disp( 'Yes' ), endk，关键字后面应该加上空格这种方式有助于区分关键字与函数。l，一个块（block）内部的一个逻辑组语句应该通过一个空白行将其分隔开在块的逻辑单元之间加入空白行可以增强代码的可读性。m，块（blocks）之间应该用多行空白行分隔一种方式是采用3 隔空白行。采用大的间隔来与块内分隔相区别，使得在文件中，块看起来非常明显。另外一种方式是采用注释符号后面跟多个诸如*或者－。n，通过排列成行列整齐的方式来加强可读性代码排列成行列整齐的形式可以使得切分表达式容易阅读与理解。这种排版也有助于揭示错误。示例： weithedPopulation = (doctorWeight *nDoctors) + …(layerWeight* nLawyers) + …(chiefWeight* nChiefs)；（2）注释注释的目的是为代码增加信息。注释的典型应用是解释用法、提供参考信息、证明结果、阐述需要的改进等。经验表明，在写代码的同时就加上注释比后来再补充注释要好。a，注释不能够改变写得很糟糕的代码效果注释不能够弥补因为代码命名不当、没有清晰的逻辑结构等造成的缺陷。存在这样缺陷的代码应该重写。b，注释文字应该简洁易读一个糟糕的或者是无用的注释反而会影响读者的正常理解。N.Schryer提到：“如果代码与注释不一致，那么或许两者都是错误的。”一个通常更重要的是注释应该讲的是“为什么”（Why）和“怎么做”（how），而不是“是什么”（what）。c，函数的注释写法（尽量使用英文）如下例子：%================================================================% 功能： &求圆孔的夫琅禾费衍射光强分布% 参数： &CircleHoleFD为圆孔结构体，包含圆孔衍射相关信息；% & & & & &theta为衍射场的次波方向，可以为向量，求取各方向的光强% 返回值： I为衍射光强分布% 主要思路：使用夫琅禾费单缝衍射公式计算% 备注： & 入射角只考虑一个维度的% 调用方法：见CalcCircleHoleFD_Test文件% 日期： &
20:37%================================================================function I =CalcCircleHoleFD(CircleHoleFD, theta)。。。（具体编程开始）&（3）文档a，文档规范化作为有用的文档应该包含一个对如下内容的可读性的描述：代码打算干什么（要求），它是如何工作的（设计），它依赖于什其他什么函数以及怎么被其他代码调用（接口），以及它是如何测试的等。对于额外的考虑，文档可以包含解决方案的选择性的讨论以及扩展与维护的建议。b，首先考虑书写文档一些程序员相信的方法是：“代码第一，回答问题是以后的事情。”而通过经验，我们绝大多数人知道先开发设计然后再实现可以导致更加满意的结果。如果将测试与文档留在最后，那么开发项目几乎不能够按期完成的。首先书写文档可以确保其按时完成甚至可能减少开发时间。c，修改一个专业的对代码修改进行管理和写文档的方法是采用源程序控制工具。对于很简单的工程，在函数文件的注释中加入修改历史比什么都不做要好。&参考文献： & &Anonymous，Good Matlab Programming Practices for the Non-Programmer & &Richard Johnson，MATLAB Programming Style Guidelines & &Fayssal El Moufatich，MATLAB Programming Tips & Tricks:Part I Functional Programming & &Genial 译，MATLAB 编程风格指南了凡春秋：&
转载本文请联系原作者获取授权，同时请注明本文来自何亮科学网博客。链接地址：
上一篇：下一篇：
当前推荐数：5
评论 ( 个评论)
扫一扫，分享此博文
作者的精选博文
作者的其他最新博文
热门博文导读
Powered by
Copyright &用matlab画函数频域图,求程序 以下是一个函数的程序,求教将其频域图画出的程序,用matlab编程.clear d1;clear d2;clear d3;clear l1;clear l2;%输入输出对比图fs=480000;C=2.37;A=1;B=1;G=1;H=1;j1=besselj(1,2.37);j2=besselj(2,2.37);f0=60000;t0=0.003;t=0:1/fs:t0;D=1;fc=5000;signal=D*cos(2*pi*fc*t);% signal=signal+D*cos(2*pi*fc*t+rand(1)*2*pi);F0=cos(signal);I=A+B*cos(C*cos(2*pi*f0*t)+signal);画出函数I的频域图 一楼的朋友请先帮我进行一下傅里叶变换,就是这个我不知道该怎么弄 ,在matlab里运行有错误：Error using ==> plotVectors must be the same lengths.
请问你是要画出I的频域图吗,如果是试试下面这段程序是不是你要的I_fft=fftshift(abs(fft(I)));f=linspace(-0.5,0.5,length(t));plot(f,I_fft)
sinα-2cosα=0移向 sinα=2cosα tana=2 sina=2/杠号5 cosa=1/杠号5 2sinαcosα=4/5(sinα)^2=4/5结果是8/5
cos(2(x-pai/4))=2 (cos(x-pai/4))^2-1=2×(√2/10)^2-1=-24/25cos(2(x-?/4)) = cos(2x-pai/2)=sin2xsin2x=-24/25x(pai/2,3pai/4)则2x（pai,3pai/2）cos2x=-√(1-sin^2(2x))=-7/25sin(2x+pai/3)=sin2xcos pai/3+cos2x sin（pai/3）=-(48+7√3)/50
设f(x)是定义在R上的奇函数,且y=f(x)的图像关于直线x=1/2对称若f(x)是定义在R上的奇函数,且y= f(x)的图像关于x=1/2 对称,则f(1)+f(2)+f(3)+f(4)+f(5)=________
f(x)是定义在R上的奇函数f(0)=0 y= f(x)的图像关于x=1/2 对称f(1)=0 f(-1)=-f(1)=0 y= f(x)的图像关于x=1/2 对称f(2)=0 f(-2)=-f(2)=0 ……f(1)+f(2)+f(3)+f(4)+f(5)=____0____
函数f（x）的定义域是（1,9】.函数g（x）=f（x）+f（x的平方）的定义域是多少?我算出来（1,3】对么.
g（x）=f（x）+f（x^2）因为函数f（x）的定义域是（1,9】.又因为 x^2符合函数关系f所以 1
其他相关问题