数值分析及其MATLAB实现第二篇MATLAB快速入门word第一章13.docVIP

1.3 MATLAB的命令和窗口环境 MATLAB是一个标准的Windows界面，可以利用菜单中的命令完成对命令窗口的操作。它的使用方法与Windows的一般应用程序相同，参看图1-8。下面将对菜单进行介绍。 图1-8 1.3.1 M文件 M文件有两种类型：文本M文件和函数M文件。 （一）文本M文件 一个比较复杂的程序常常要作反复的调试，这时你不妨建立一个文本文件并把它储存起来，可以随时调用进行计算。建立文本文件可以在File菜单中选择New，再选择M-file，这时MATLAB将打开一个文本编辑窗口，在这里输入命令和数据。储存时文件名遵循MATLAB变量命名的原则，但必须以 m 为扩展名，其一般形式为 M文件名.m 如 juzhen.m, pp.m 等。 值得注意的是，文本M文件中的变量都是全局变量，在执行过程中, 文本M文件中的命令可以使用目前工作区中的变量，它所产生的变量也将成为工作区的一部分。例如我们把下面的程序保存为名为budijifen1.m的文本M文件： function budijifen1 x=-4*pi:0.014*pi;a=1; for C=-3:3 y =1/2*(-sin(a*x)+a*x)/a+C plot(x,y) hold on end grid hold off xlabel(自变量 X) ylabel(因变量 Y) title(不定积分的积分曲线族) legend(不定积分的积分曲线族) 然后在MATLAB工作窗口中输入文件名： budijifen1 运行后屏幕显示计算结果(略)和图形(见图1-9)。 再如，将下面的一个生成矩阵的程序保存为名为juzhen.m的文本M文件： for i=1:k for j=1:n b(i,j)=7/(2*i+3*j-6); end end a=rats(b) 然后在MATLAB工作窗口中输入： k=3;n=4; juzhen,b 图1-9 名为budijifen1.m的文本M文件的图形 运行后屏幕显示运行结果如下： a = -7 7/2 7/5 7/8 7 7/4 1 7/10 7/3 7/6 7/9 7/12 b = -7.0000e+000 3.5000e+000 1.4000e+000 8.7500e-001 7.0000e+000 1.7500e+000 1.0000e+000 7.0000e-001 2.3333e+000 1.1667e+000 7.7778e-001 5.8333e-001 （二） 函数M文件 函数M文件是另一类M文件，我们可以根据需要建立自己的函数文件，它们能够像库函数一样方便地调用，从而极大地扩展MATLAB的能力。如果对于一类特殊的问题，建立起许多函数M文件，就能最终形成独立的工具箱。 函数M文件的第一行有特殊的要求，其形式必须为 function 因变量= 函数名（自变量） 其他各行为从自变量计算因变量的语句，并最终将结果赋予因变量。而这个M文件的文件名必须是函数名.m。下面给出函数文件的一个简单例子。 如果我们要在某个程序中要调用函数e，就必须建立并保存下面的名为fun.m的函数M文件。 function y = fun(x) y = exp((-x.^2)./2)./(sqrt(2*pi)); 然后在MATLAB工作窗口中输入程序： [Q1,FCNT14] = quad(@fun,0,pi/2,1.e-14,3) % quad是一个计算数值积分的程序名 运行后屏幕显示计算的部分结果如下： ………………………………………………………………………… 551 1.5441340226 1-002 0.0015978928 553 1.5574651747 1-002 0.0015650585 Q1 = 0.44188501721659 FCNT14 = 553 符号 % 后面是注释部分，给出了这个程序的说明，MATLAB执行时忽略这些内容。实际上MATLAB的所有函数都有这样一段注释。 函数M文件有多个因变量时，要用[]将它们括起来，请看下面数值积分程序quad.m的部分内容: function [Q,fcnt] = quad(funfcn,a,b,tol,trace,varargin) %QUAD Numer