腐蚀电化学极化曲线拟合程序.docVIP

曲线拟合及其MATLAB程序 实验任务：用最小二乘法将一组离散数据进行曲线拟合，并用MATLAB做出曲线。 实验要求： 第三次作业的最后一题求函数关系并作出曲线，此组x和y的实验数据如下表所示。 xi 1 2 3 4 5 6 7 8 9 10 yi 6.00 3.50 2.67 2.75 2.00 1.80 1.70 1.63 1.56 1.50 实验内容： 在在MATLAB工作窗口输入程序 x=[1 2 3 4 5 6 7 8 9 10]; y=[6.00 3.50 2.67 2.75 2.00 1.80 1.70 1.63 1.56 1.50]; plot(x,y,r*), legend(实验数据(xi,yi)) xlabel(x), ylabel(y), title(实验数据点(xi,yi)的散点图) 运行后屏幕显示数据的散点图（略）。 （3）编写下列MATLAB程序计算在处的函数值，即输入程序 syms a1 a2 a3 a4; x=[1 2 3 4 5 6 7 8 9 10]; fi=a1.*x.^3+ a2.*x.^2+ a3.*x+ a4 运行后屏幕显示关于a1, a2, a3和a4的线性方程组 fi = [a1+a2+a3+a4, 8*a1+4*a2+2*a3+a4, 27*a1+9*a2+3*a3+a4, 64*a1+16*a2+4*a3+a4, 125*a1+25*a2+5*a3+a4, 216*a1+36*a2+6*a3+a4, 343*a1+49*a2+7*a3+a4, 512*a1+64*a2+8*a3+a4, 729*a1+81*a2+9*a3+a4, 1000*a1+100*a2+10*a3+a4] 编写构造误差平方和的MATLAB程序 y=[6.00 3.50 2.67 2.75 2.00 1.80 1.70 1.63 1.56 1.50]; fi =[a1+a2+a3+a4,8*a1+4*a2+2*a3+a4,27*a1+9*a2+3*a3+a4,64*a1+16*a2+4*a3+a4,125*a1+25*a2+5*a3+a4,216*a1+36*a2+6*a3+a4,343*a1+49*a2+7*a3+a4,512*a1+64*a2+8*a3+a4,729*a1+81*a2+9*a3+a4, 1000*a1+100*a2+10*a3+a4]; fy=fi-y; fy2=fy.^2;J=sum(fy.^2) 运行后屏幕显示误差平方和如下 J = (a1+a2+a3+a4-6)^2+(8*a1+4*a2+2*a3+a4-7/2)^2+(27*a1+9*a2+3*a3+a4-267/100)^2+(64*a1+16*a2+4*a3+a4-11/4)^2+(125*a1+25*a2+5*a3+a4-2)^2+(216*a1+36*a2+6*a3+a4-9/5)^2+(343*a1+49*a2+7*a3+a4-17/10)^2+(512*a1+64*a2+8*a3+a4-163/100)^2+(729*a1+81*a2+9*a3+a4-39/25)^2+(1000*a1+100*a2+10*a3+a4-3/2)^2 为求使达到最小，只需利用极值的必要条件。，得到关于的线性方程组，这可以由下面的MATLAB程序完成，即输入程序 syms a1 a2 a3 a4 J=(a1+a2+a3+a4-6)^2+(8*a1+4*a2+2*a3+a4-7/2)^2+(27*a1+9*a2+3*a3+a4-267/100)^2+(64*a1+16*a2+4*a3+a4-11/4)^2+(125*a1+25*a2+5*a3+a4-2)^2+(216*a1+36*a2+6*a3+a4-9/5)^2+(343*a1+49*a2+7*a3+a4-17/10)^2+(512*a1+64*a2+8*a3+a4-163/100)^2+(729*a1+81*a2+9*a3+a4-39/25)^2+(1000*a1+100*a2+10*a3+a4-3/2)^2; Ja1=diff(J,a1); Ja2=diff(J,a2); Ja3=diff(J,a3); Ja4=diff(J,a4); Ja11=simple(Ja1), Ja21=simple(Ja2), Ja31=simple(Ja3), Ja41=simple(Ja4), 运行后屏幕显示J分别对a1, a2 ,a3 ,a4的偏导数如下 Ja11 =6050*a4+3956810*a1+441650*a2+50666*a3-497579/50 Ja21 =770*a4+441650*a1+50666*a2+6050*a3-66