数模程序资料.docVIP

数学建模程序资料 [例1] 输入命令： v=eye(3,4) %3×4对角线为1的矩阵。 运行后输出结果如下： v= 1 0 0 0 0 1 0 0 0 0 1 0 [例5] 作下面函数的图形. . 解 用以下程序实现： x=-7.5:0.5:7.5; y=x; [X,Y]=meshgrid(x,y); R=sqrt(X.^2+Y.^2)+eps; Z=sin(R)./R; mesh(X,Y,Z) 运行后输出图1-6 图1-6 例5的图 [例6] 作螺旋线 . 解 用以下程序实现 t=0:pi/50:10*pi; plot3(sin(t),cos(t),t) 运行后输出图1-7。 图1-7 例6的图 [例8] 作曲面在时的图形。 解 作曲面的图形在MATLAB工作区输入以下程序实现： x=-500:20:500;y=x; [X,Y]=meshgrid(x,y); Z=X.^2-Y.^2; mesh(X,Y,Z) 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)。 图1-9 名为budijifen1.m的文本M文件的图形 函数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 函数的插值方法 例6.2.2 求函数e在上线性插值多项式，并估计其误差. 解 输入程序 X=[0,1]; Y =exp(-X) , l01= poly(X(2))/( X(1)- X(2)), l11= poly(X(1))/( X(2)- X(1)), l0=poly2sym (l01), l1=poly2sym (l11), P = l01* Y(1)+ l11* Y(2), L=poly2sym (P), 运行后输出基函数l0和l1及其插值多项式的系数向量P和插值多项式L为 l0 = l1 = P = -x+1 x -0.6321 1.0000 L = -1423408956596