文档介绍：function[Fx,Fy]=FDM_circular(X,Y,X_dot,Y_dot)Fx=0;Fy=0;epsilon=;%偏心率==;d=;%轴颈直径dlambda=L/d;m=40;n=30;phi1=0;phi2=2*pi;delta_phi=(phi2-phi1)/m;delta_lambda=2/n;k=1;P0=zeros(m+1,n+1);%zeros就是用来表示线性代数的0矩阵A=zeros(m+1,n+1);B=zeros(m+1,n+1);C=zeros(m+1,n+1);D=zeros(m+1,n+1);E=zeros(m+1,n+1);F=zeros(m+1,n+1);H=zeros(m+1,n+1);fori=1:1:m+1%theta(i)=(i-1)*delta_phi;endforj=1:1:n+1H(i,j)=1+epsilon*cos((i-1)*delta_phi);endendS=0;T=0;ERR=1e-3;GAP=1;whileGAP>ERRk=k+1fori=1:1:m+1forj=1:1:n+1if(i==1)|(j==1)|(i==m+1)|(j==n+1)P(i,j)=0;elseA(i,j)=(1+epsilon*cos((i+1/2-1)*delta_phi))^3;B(i,j)=(1+epsilon*cos((i-1/2-1)*delta_phi))^3;C(i,j)=(d/L)^2*(delta_phi/delta_lambda)^2*(1+epsilon*cos((i-1)*delta_phi))^3;D(i,j)=(d/L)^2*(delta_phi/delta_lambda)^2*(1+epsilon*cos((i-1)*delta_phi))^3;E(i,j)=A(i,j)+B(i,j)+C(i,j)+D(i,j);F(i,j)=6*delta_phi*((1+epsilon*cos((i+1/2-1)*delta_phi))-(1+epsilon*cos((i-1/2-1)*delta_phi)))+12*(delta_phi)^2*(X_dot*cos(phi1+(i-1)*d