MATLAB 编程技巧

MATLAB 编程技巧

ID:41037291

大小:138.50 KB

页数:5页

时间:2019-08-14

MATLAB 编程技巧_第1页
MATLAB 编程技巧_第2页
MATLAB 编程技巧_第3页
MATLAB 编程技巧_第4页
MATLAB 编程技巧_第5页
资源描述:

《MATLAB 编程技巧》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库

1、MATLAB编程技巧一、尽量避免使用循环:循环语句及循环体经常被认为是~MATLAB编程的瓶颈问题。改进这样的状况有两种方法:(1)尽量用向量化的运算来代替循环操作。我们将通过如下的例子来演示如何将一般的循环结构转换成向量化的语句。〖例3.19〗考虑下面无穷级数求和问题:如果我们只求出其中前有限项,比如100,000项之和(要精确地求出级数的和,无需求100000项,几十项往往就能得出满意的精度。这里主要是为了演示循环运算向量化的优越性。),则可以采用下面的常规语句进行计算>>tic,s=0;fori=1:100000,s=s+(1/2^i+1/3^i);en

2、d,s,tocs=1.5000elapsed_time=1.9700如果采用向量化的方法,则可以得出下面结果。可以看出,采取向量化的方法比常规循环运算效率要高得多。>>tic,i=1:100000;s=sum(1./2.^i+1./3.^i),tocs=1.5000elapsed_time=0.3800(2)在必须使用多重循环的情况下,如果两个循环执行的次数不同,则建议在循环的外环执行循环次数少的,内环执行循环次数多的。这样也可以显著提高速度。〖例3.20〗考虑生成一个5x10000的Hilbert长方矩阵,该矩阵的定义是其第i行第j列元素为h_{i,j}=1

3、/(i+j-1)。我们可以由下面语句比较先进行i=1:5的循环和后进行该循环的耗时区别,其效果和前面分析的是一致的。>>ticfori=1:5forj=1:10000H(i,j)=1/(i+j-1);endendtocelapsed_time=8.6800>>tic,forj=1:10000fori=1:5J(i,j)=1/(i+j-1);endendtocelapsed_time=25.7000一、大型矩阵的预先定维给大型矩阵动态地定维是个很费时间的事。建议在定义大矩阵时,首先用MATLAB的内在函数,如zeros()或ones()对之先进行定维,然后再进行

4、赋值处理,这样会显著减少所需的时间的。再考虑例3-20中的问题,如果输入下面的命令>>ticH=zeros(5,10000);fori=1:5forj=1:10000H(i,j)=1/(i+j-1);endendtocelapsed_time=1.0400则采用预先定维的方法,再结合向量化的方法,我们可以给出下面的MATLAB语句。>>ticH=zeros(5,10000);fori=1:5H(i,:)=1./[i:i+9999];endtocelapsed_time=0.060可见,预先定维后,所需要的时间显著地减少了。可以看出,同样一个问题,由于采用了有效

5、的措施,所需的时间就可以从25.7秒减少到0.06秒,亦即效率提高了428倍。二、对二重循环这样的特殊问题,我们还可以使用meshgrid()函数构造两个5x10000矩阵i和j,从而直接得出H矩阵,更进一步地加快速度。>>tic,[i,j]=meshgrid(1:5,1:10000);H=1./(i+j-1);tocelapsed_time=0一、优先考虑内在函数矩阵运算应该尽量采用MATLAB的内在函数,因为内在函数是由更底层的编程语言C构造的,其执行速度显然快于使用循环的矩阵运算。二、采用有效的算法在实际应用中,解决同样的数学问题经常有各种各样的算法。例

6、如求解定积分的数值解法在MATLAB中就提供了两个函数quad()和quad8(),其中后一个算法在精度、速度上都明显高于前一种方法。所以说,在科学计算领域是存在“多快好省”的途径的。如果一个方法不能满足要求,可以尝试其他的方法。三、应用Mex技术虽然采用了很多措施,但执行速度仍然很慢,比如说耗时的循环是不可避免的,这样就应该考虑用其他语言,如C或Fortran语言。按照Mex技术要求的格式编写相应部分的程序,然后通过编译联接,形成在MATLAB可以直接调用的动态连接库(DLL)文件,这样可以显著地加快运算速度。MATLAB图形上的文字修饰字符对象及其属性 文

7、字标注是图形修饰中的重要因素,它可以是用户在窗口上随意添加的字符说明,还可以是坐标轴对象中所用到的刻度标志等。字符对象的常用属性如下:·Color属性:字符的颜色。该属性的属性值是一个1x3颜色向量。·FontAngle属性:字体倾斜形式。如正常'normal'和斜体'italic'等。·FontName属性:字体的名称。如'TimesNewRoman'与'Courier'等。·FontSize属性:字号大小。默认以pt为单位,属性值应该为实数。·FontWeight属性:字体是否加黑。可以选择'light'、'normal'(默认值)、'demi'和'bol

8、d'4个选项,其颜色逐渐变黑。·Hor

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。