腐蝕電化學極化曲線擬合程序

22曲線擬合及其MATLAB程序實驗任務：用最小二乘法將一組離散數據進行曲線擬合，并用MATLAB做出曲線。實驗要求：第三次作業(yè)的最后一題求函數關系并作出曲線，此組x和y的實驗數據如下表所示。xi12345678910yi6.003.502.672.752.001.801.701.631.561.50實驗內容:在在MATLAB工作窗口輸入程序>>x=[12345678910];y=[6.003.502.672.752.001.801.701.631.561.50];plot(x,y,'r*'),legend('實驗數據(xi,yi)')xlabel('x'),ylabel('y').title('實驗數據點(xi,yi)的散點圖')運行后屏幕顯示數據的散點圖(略)。6-實驗數據點網汕的散點團6-實驗數據點網汕的散點團5.5*實驗數iSfxi.yi)4.54A3.52.5*-TOC\o"1-5"\h\z**卡IIIIIIT*：?2345678910編寫下列MATLAB程序計算f(x)在(x,)處的函數值，即輸入程序>>symsa1a2a3a4;x=[12345678910];fi=al.*x43+a2.*x.A2+a3.*x+a4運行后屏幕顯示關于al,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.003.502.672.752.001.801.701.631.561.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.A2;J=sum(fy.A2)運行后屏幕顯示誤差平方和如下J=(a1+a2+a3+a4-6)A2+(8*a1+4*a2+2*a3+a4-7/2)A2+(27*a1+9*a2+3*a3+a4-267/100)A2+(64*a1+16*a2+4*a3+a4-11/4)A2+(125*a1+25*a2+5*a3+a4-2)A2+(216*a1+36*a2+6*a3+a4-9/5)A2+(343*a1+49*a2+7*a3+a4-17/10)A2+(512*a1+64*a2+8*a3+a4-163/100)A2+(729*a1+81*a2+9*a3+a4-39/25)A2+(1000*a1+100*a2+10*a3+a4-3/2)A2J=0為求《，役，行3使J達到最小，只需利用極值的必要條件叫。(k二1,2,3,4)，得到關于a1,a2,S'a4的線性方程組，這可以由下面的MATLAB程序完成，即輸入程序>>symsa1a2a3a4J=(a1+a2+a3+a4-6)A2+(8*a1+4*a2+2*a3+a4-7/2)A2+(27*a1+9*a2+3*a3+a4-267/100)人2+(64*al+16*a2+4*a3+a4-ll/4)人2+(125*al+25*a2+5*a3+a4-2)人2+(216*al+36*a2+6*a3+a4-9/5)人2+(343*al+49*a2+7*a3+a4-l7/10)A2+(512*a1+64*a2+8*a3+a4-163/100)A2+(729*a1+81*a2+9*a3+a4-39/25)A2+(1000*a1+100*a2+10*a3+a4-3/2)A2;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分別對al,a2,a3,a4的偏導數如下Jall=6050*a4+3956810*a1+441650*a2+50666*a3-497579/50Ja21=770*a4+441650*a1+50666*a2+6050*a3-66681/50Ja31=110*a4+50666*a1+6050*a2+770*a3-10679/50Ja41=20*a4+6050*a1+770*a2+110*a3-2511/50解線性方程組Ja11=0,Ja21=0,Ja31=0,Ja41=0,輸入下列程序>>A=[3956810,441650,50666,6050;441650,50666,6050,770;50666,6050,770,110;6050,770,110,20];B=[497579/50,66681/50,10679/50,2511/50];C=B/A,f=poly2sym(C)運行后屏幕顯示擬合函數f及其系數C如下C=-0.01760.3744-2.64527.9627f=-5066427105276917/288230376151711744*xA3+6744376644764349/18014398509481984*xA2-2978180218619021/1125899906842624*x+1120645707277323/140737488355328故所求的擬合曲線為f(x)=-0.0176x3+0.3744x2-2.6452x+7.9627.編寫下面的MATLAB程序估計其誤差，并作出擬合曲線和數據的圖形，輸入程序>>xi=[12345678910];y=[6.003.502.672.752.001.801.701.631.561.50];n=length(xi);f=-0.0176.*xi.A3+0.3744.*xi.A2-2.6452.*xi+7.9627;x=1:1:10;F=-0.0176.*xi43+0.3744.*xi42-2.6452.*xi+7.9627;fy=abs(f-y);fy2=fy.A2;Ew=max(fy),E1=sum(fy)/n,E2=sqrt((sum(fy2))/n)plot(xi,y,'r*'),holdon,plot(x,F,'b-'),holdofflegend('數據點(xi,yi)','擬合曲線y=f(x)'),xlabel('x'),ylabel('y'),tit