评阅标准和注意事项.ppt

function y=f(x) y=x^3*exp(x)+sqrt(x)*tan(x); 在Command窗口中用以下命令来调用该函数： y=f(3) 这种调用更加符合数学上的习惯，但函数表达式需要单独编写成一个文件。 （2）用function格式的M函数文件来定义 在Editor窗口中编写函数文件： f =inline(x^3*exp(x)+sqrt(x)*tan(x), x) f = Inline function: f(x) = x^3*exp(x)+sqrt(x)*tan(x) 执行结果： 要求点3上的函数值，可以使用与function格式下相同的命令： y=f(3) （3）用inline命令来定义 sym、function 和 inline 所定义的函数，其适用范围最主要的区别为： 前者用于解析运算（对绘图比较有用），后两者则一般用于数值运算。 解析运算 数值运算 解析方法求方程根： 因式分解法 求根公式法 数值方法求方程根（牛顿切线法）： x0 x1 x2 x x0 为初始搜索点，迭代序列收敛于根 x，求解成功 x0 x1 x2 x 初始搜索点选择不当时，迭代序列可能不收敛，导致求解失败 x y=f (x) x01 x02 x03 x04 x05 插值 线性拟合 多项式拟合 指数拟合 观 测 点 观 测 值 插值与拟合 例：设有6个数据点，其横、纵坐标分别用向量a、b表示： a=[1 2 5 7 11 12] b=[33.78 32.65 26.25 25.55 19.24 16.65] 1、直接以这6点坐标来绘制曲线； 2、通过插值，加密作图点，再重新绘制曲线。 解： a=[1 2 5 7 11 12]; b=[33.78 32.65 26.25 25.55 19.24 16.65]; plot(a,b, -o) %直接绘图 aa=1:0.5:12; % 以0.5的步长来得到插值点 bb=interp1(a,b,aa, spline); % 用“样条插值”方法求得插值 hold on, plot(aa,bb, ‘r.-’) %在保留上幅图的情况下绘制新图 一维插值问题举例（更多内容见课件《第六讲》） 例：将上题中的6个数据点，分别作 1、线性拟合； 2、3次多项式拟合。 解： a=[1 2 5 7 11 12]; b=[33.78 32.65 26.25 25.55 19.24 16.65]; plot(a,b, o) % 绘散点图 p= polyfit(a,b,1); y1=p (1)*a+p (2); % 线性拟合 hold on; plot(a,y1,r) % 绘制直线图形 q= polyfit(a,b,3); y2= q(1)*a.^3+q(2)*a.^2+q(3)*a+ q(4); % 3次多项式拟合 plot(a,y2,k) 多项式拟合举例（更多内容见课件《第六讲》） 例：已知 x = 0:0.1:3，y =[0.1200 0.1826 0.2421 0.2977 0.3485 0.3940 0.4337 0.4670 0.4937 0.5135 0.5264 0.5323 0.5313 0.5237 0.5098 0.4900 0.4647 0.4345 0.4000 0.3618 0.3205 0.2770 0.2319 0.1858 0.1396 0.0939 0.0494 0.0066 -0.0338 -0.0714 -0.1057]，并已知该组数据满足 ，其中 为待定系数。用最小二乘拟合求出这些待定系数。 一般的最小二乘拟合举例 问题一般化：设有数据组 (xi, yi)（i =1, 2, …, n），且已分析出这组数据满足 某一函数原型： y =f (a, x) （a=(a1, a2,…, am)为待定系数向量） 求函数中的待定系数，以确定整个函数。 即函数的结构已确定，只是其中的系数尚未确定（多项式拟合是此类问题的特殊情况）。 解：输入MATLAB命令 x=0:0.1:3; y=[0.1200 0.1826 0.2421 0.2977 0.3485 0.3940 … 0.4337 0.4670 0.4937 0.5135 0.5264 0.5323 0.5313 … 0.5237 0.5098 0.4900 0.4