第13讲 最佳一致逼近法.ppt

《第13讲 最佳一致逼近法.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、3.2.4最佳(一致)逼近法3.2.4最佳(一致)逼近法3.2.4最佳(一致)逼近法3.2.4最佳(一致)逼近法最佳(一致)逼近法应用实例Matlab程序：>>fork=1:1:4phi1(k)=50*(1-cos(k*pi/5))end>>phi1=[0,phi1,100]>>fork=1:1:4kphi(k)=-25*(1-cos(k*pi/5))end>>kphi=[0,kphi,-50]>>phi2=100/9*(10.^(-1*kphi/50)-1)phi2=10/9*acos(kph

2、i/50+1)>>phi=(phi1+phi2)/2>>kphi=-50*log10(1+9/100*phi)kphi=50cos(9/10*phi-1)>>f0=cos(phi*pi/180)>>f1=sin(phi*pi/180)>>f2=cos(kphi*pi/180)>>f3=sin(kphi*pi/180)>>f5=cos((kphi-phi)*pi/180)>>f6=sin((kphi-phi)*pi/180)>>x0=[111111]'>>fun=inline('[x(1)+x(3

3、)+x(5)+(-x(1)*x(3)-x(2)*x(4))-1*x(6);0.9943*x(1)+0.107*x(2)+0.986*x(3)-0.166*x(4)+x(5)+0.963*(-x(1)*x(3)-x(2)*x(4))-0.271*(x(2)*x(3)-x(1)*x(4))+1*x(6);….]','x')>>[xn,fval,flag]=fsolve(fun,x0)>>x0=[2;3];>>fun=inline('[2*x(1)-x(2)^2-exp(-x(1));-x(1)^3+

4、x(1)*x(2)-exp(-x(2))]','x')>>[xn,fval,flag]=fsolve(fun,x0)作业（任选一题）