西南交大MATLAB 语言程序设计.ppt

4.2 inline 函数和匿名函数 inline 函数，可以免去文件 f=inline(‘sin(x.^2+y.^2)’,’x’,’y’) MATLAB 7.0 5 二维图形绘制 二维图形绘制基本语句 构造向量: 例：选项为红色点划线且每个转折点用五角星表示 ‘r-.pentagram’ 例： x=[-pi : 0.05: pi]; % 以 0.05 为步距构造自变量向量 y=sin(tan(x))-tan(sin(x)); % 求出各个点上的函数值 plot(x,y) plot(x,y,r-.pentagram) x=[-pi:0.05:-1.8,-1.801:.001:-1.2, -1.2:0.05:1.2,... 1.201:0.001:1.8, 1.81:0.05:pi]; % 以变步距方式构造自变量向量 y=sin(tan(x))-tan(sin(x)); % 求出各个点上的函数值 plot(x,y) % 绘制曲线 6 三维图形绘制 6.1三维曲线绘制 stem3, fill3, bar3 等 例：参数方程 t=0:.1:2*pi; % 构造 t 向量，注意下面的点运算 x=t.^3.*sin(3*t).*exp(-t); y=t.^3.*cos(3*t).*exp(-t); z=t.^2; plot3(x,y,z), grid % 三维曲线绘制 stem3(x,y,z); hold on; plot3(x,y,z), grid 6.2 三维曲面绘制 一般曲面绘制 mesh( )绘制网格图,surf( )绘制表面图。 其他函数，光照下 surfl( ), surfc( )，瀑布型waterfall( ) 等高线绘制 contour( ), contour3( ) 例:Butterworth 滤波器 [x,y]=meshgrid(0:31); n=2; D0=200; D=sqrt((x-16).^2+(y-16).^2); % 求距离 z=1./(1+D.^(2*n)/D0); mesh(x,y,z), % 计算并绘制滤波器 axis([0,31,0,31,0,1]) % 重新设置坐标系，增大可读性 surf(x,y,z) % 绘制三维表面图 contour3(x,y,z,30) ％三维等高线图，30等高线条数 7. 多项式、插值与数据拟合 多项式MATLAB命令 插值 Lagrange插值 Hermite插值 Runge现象和分段插值 分段插值 样条插值的MATLAB表示 数据拟合 多项式拟合 函数线性组合的曲线拟合方法 最小二乘曲线拟合 B样条函数及其MATLAB表示 关于多项式MATLAB命令 一个多项式的幂级数形式可表示为： 幂系数：在MATLAB里，多项式用行向量表示，其元素为多项式的系数，并从左至右按降幂排列。 例： 被表示为 p=[2 1 4 5] poly2sym(p) ans = 2*x^3+x^2+4*x+5 Roots: 多项式的零点可用命令roots求的。 例： r=roots(p) 得到 r = 0.2500 + 1.5612i 0.2500 - 1.5612i -1.0000 所有零点由一个列向量给出。 polyval: 可用命令polyval计算多项式的值。 例： 计算y(2.5) c=[3,-7,2,1,1]； xi=2.5; yi=polyval(c,xi) yi = 23.8125 如果xi是含有多个横坐标值的数组，则yi也为与xi长度相同的向量。 c=[3,-7,2,1,1]; xi=[2.5,3]; yi=polyval(c,xi) yi = 23.8125 76.0000 polyfit:给定n+1个点将可以唯一确定一个n阶多项式。利用命令polyfit可容易确定多项式的系数。 例： x=[1.1,2.3,3.9,5.1]; y=[3.887,4.276,4.651,2.