当变量之间为非线性相关时,可用非线性最小二乘数据拟合高斯.pptVIP

例题1: 某筛分试验测得筛孔尺寸与累计产率间的关系如表所示，试建立模型。 用MATLAB编程 首先建立一个′ep′函数文件。 function f=ep(b,x)； f=(b(1) * exp(－b(2) * x) 输入 x=[50 25 13 6 3 0.5]’; y=[1.5 11 30 55 70 88]; b0=[1 1]; %初始值 [b,r,j]=nlinfit(x,y, ′ep′,b0);b [yp,d]=nlperdci(′ep′,x,b,r,j); ci=nlparci(b,r,j) %输出预测值、残差及置信区间 nlintool(x,y, ′ep′,b0); %绘制非线性拟合曲线图3.14 plot(x,r, ′――′) %绘残差曲 结果 得 b= 91.4158 0.0855 ci= 90.2085 92.6231 0.0828 0.0883 移动垂直于x轴的虚线可以显示y随x的变化。数学模型为 y=91.4158e-0.0855x 例2 已知浮煤累计灰分及产率为x和y，试进行非线性拟合。 x=[3.46 7.33 9.47 10.48 11.19 11.79 12.78 20.50]； y=[10.69 56.84 76.89 82.15 84.70 86.32 88.45 100]； 现采用含有双曲正切的函数拟合。 输入 x=x’； b0=[0.2 0 -0.4 -0.1 5]；%初始值 [b,r,j]=nlinfit(x,y, ′mf′,b0)； [yp,d]=nlpredci (′mf′,x,b,r,j); nlintool(x,y, ′mf′,b0) 结果 b= 0.3414 0.0136 -0.3797 -0.3365 6.2948 即 y=100*(0.3414+0.0136x-0.3797tanh(-0.3365(x-6.2948))) 由计算出的yy可知，拟合效果很好。 第三讲 线性规划 线性规划问题是目标函数和约束条件均为线性函数的问题 补充知识 :解方程或方程组 用solve 符号解代数方程 solve (′eq1′，′eq2′, ′eqn′, ′var1,var2,…,varn′) eqn为符与方程，varn为符号变量。 求解x2+xy+y=3，x为参数。 solve(′x ^ 2+x * y+y=3′, ′x′) ans= [-1/2 * y+1/2 * (y ^ 2-4 * y+12) ^ (1/2)] [-1/2 * y-1/2 * (y ^ 2-4 * y+12) ^ (1/2)] 解au2+υ2=0, u-υ=1。 [u,v]=solve (′a * u ^ 2+v ^ 2=0′, ′u-v=1′) u= [1/2/(a+1) * (-2 * a+2 * (-a) ^ (1/2))+1] [1/2/(a+1) * (-2 * a-2 * (-a) ^ (1/2))+1] v= [1/2/(a+1) 8 (-2 * a +2 * (-a) ^ (1/2))] [1/2/(a+1) 8 (-2 * a -2 * (-a) ^ (1/2))] 2 解线性方程组 MATLAB用函数linsolve求解线性方程组Ax=b，要求A的列数等于b的行数；也可用矩阵除等方法求解。linsolve的语法格式为 x=linsolve(A,b) 例题 解线性方程组 a=[1 2 1；2 2 3；-1 -3 0]； % rank(a)=3，有惟一解 b=[0 3 2]’; x=linsolve(a,b) x= 1 -1 1 3 非线性方程 非线性方程的标准形式为f(x)=0 函数 fzero 格式 x = fzero (fun,x0) %用fun定义表达式f(x)，x0为初始解。 x = fzero (fun,x0,options) [x,fval] = fzero(…) %fval=f(x) [x,fval,exitflag] = fzero(…) [x,fval,exitflag,output] = fzero(…) 说明 该函数采用数值解求方程f(x)=0的根。 非线性方程解例题 求x^3-2*x-5=0 的根 解： fun=x^3-2*x-5; z=fzero(fun,2) %初始估计值为2 结果为 z = 2.0946 非线性方程组求解 非线性方程组的标准形式为：F(x) = 0 其中：x为向量，F(x)为函数向量。 函数 fsolve 格式 x = fsolve(fun,x0) %用fun定义向