matlab函数主要是计算与坐标差相关的矩阵 `xx` 和 `yy` 及其衍生矩阵

function [C,W,x,y]=GDQ(N,M)
% function [C,W,y]=GDQ(N,n1,m1,n2,m2)

Deta=1e-4; % 0; % 
Deta2=0; % 1e-4; % 
Deta3=0; % 1e-4; % 
% pi=3.1415926;
% if (n1==3 && m1==3)||(m1==3 && n2==3)||(n2==3 && m2==3)||(m2==3 && n1==3)
%     Alfa=0.5;     % FF--含自由角点
% else
%     Alfa=1;       % 不含自由角点
% end
%% ---------------------------------------------------------------
% x=zeros(N);
% y=zeros(N);

for i=1:N
    x(i)=( 1-cos(pi*(i-1)/(N-1)) )/2;
end

% for i=2: N-1
% %     xx=( 1-cos(pi*(i-1)/(N-1)) )/2;               % 思考此处 i、xx后面+Deta与最后+Deta区别
% %     x(i)=(1-Alfa)*( 3*xx^2 - 2*xx^3 ) + Alfa*xx;  % 参考 C Shu书籍：Differential Quadrature and Its Application in Engineering
% %     x(i)=(1-Alfa)*( (1-cos(pi*xx))/2 ) + Alfa*xx; % 自己构造
%     x(i)=x(i) +Deta;
% end

%% ---------------------------------------------------------------
y=x;
for i=1:N
    for j=1:N
        xx(i,j)=[x(i)-x(j)]; 
        yy(i,j)=[y(i)-y(j)];
    end
end

Xx=xx';
Yy=yy